Modeling and Assessment of Land Degradation Vulnerability in Arid Ecosystem of Rajasthan Using Analytical Hierarchy Process and Geospatial Techniques

: Wind erosion is a major natural disaster worldwide, and it is a key problem in western Rajasthan in India. The Analytical Hierarchy Process (AHP), the Geographic Information System (GIS), and remote sensing satellite images are effective tools for modeling and risk assessment of land degradation. The present study aimed to assess and model the land degradation vulnerable (LDV) zones based on the AHP and geospatial techniques in the Luni River basin in Rajasthan, India. This study was carried out by examining important thematic layers, such as vegetation parameters (normalized difference vegetation index and land use/land cover), a terrain parameter (slope), climatic parameters (mean annual rainfall and land surface temperature), and soil parameters (soil organic carbon, soil erosion, soil texture, and soil depth), using the Analytical Hierarchical Process (AHP) and geospatial techniques in the Luni River basin in Rajasthan, India. The weights derived for the thematic layers using AHP were as follows: NDVI (0.27) > MAR (0.22) > LST (0.15) > soil erosion (0.12) > slope (0.08) > LULC (0.06) > SOC (0.04) > soil texture (0.03) > soil depth (0.02). The result indicates that nearly 21.4 % of the total area is prone to very high degradation risks; 12.3% is prone to high risks; and 16%, 24.3%, and 26% are prone to moderate, low, and very low risks, respectively. The validation of LDV was carried out using high-resolution Google Earth images and ﬁeld photographs. Additionally, the Receiver Operating Characteristic (ROC) curve found an area under the curve (AUC) value of 82%, approving the prediction accuracy of the AHP technique in the study area. This study contributes by providing a better understanding of land degradation neutrality and sustainable soil and water management practices in the river basin.


Introduction
Land is an essential part of the supportive ecosystem of the Earth. Humans and animals are dependent on the ecosystem services, such as food, fiber, and shelter, provided by this natural resource [1]. Resource degradation is one of the hottest topics of research at present due to its alarming threat to biodiversity [2]. Land degradation is nothing but the unsustainable over-exploitation of resources resulting in decreasing productivity gion to evaluate and develop a land degradation vulnerability map. Therefore, an integrated approach of AHP and GIS has immense potential to accurately estimate the wind-erosion-prone areas and identification of most priority areas for suitable erosion control measures. Hence, the objectives of the present study were as follows: (i) to characterize the vegetative, soil, climatic, and terrain parameters of the basin and (ii) to identify the area's most vulnerable to land degradation using AHP along with RS and GIS techniques. The results of the present study help to provide important information for soil and water management practices, land use plans, and environmental sustainability for the study area.

Study Area
Luni, the largest river basin in Rajasthan (India), is located between latitude 24°30′ N to 27°10′ N and longitude 70°49′ E to 75°04′ E ( Figure 1). Geographically, it is located adjacent to the western part of the Aravalli Range and occupies an area of 68,939 km 2 . The basin covers the districts of Ajmer, Jodhpur, Nagaur, Barmer, Jalore, and Sirohi in part and the district of Pali in full. The elevation ranges from −5 to 1613 m from the mean sea level (MSL), with a decreasing gradient from the northern to the southern part of the basin. In this basin, rainfall occurs mostly from June to September (>90%), and it ranges from as high as ~1000 mm in the southern and eastern parts to as low as ~234 mm in the western part of the basin [39]. The annual maximum potential evapotranspiration in the river basin is 1850 mm, and it is the highest (260 mm) in May and the lowest (77 mm) in December [39]. The study area, with the characteristics of low rainfall, high temperatures, and low relative humidity, has an arid to hyper-arid climate [40].

Data Acquisition
In the present study, nine input layers were used to identify the area most vulnerable to land degradation, namely, the normalized difference vegetation index (NDVI) from a moderate-resolution imaging spectroradiometer (MODIS), land surface temperature (LST) from MODIS, mean annual rainfall from the Climate Hazards Group Infrared Precipitation with Station data (CHIRPS), land use/land cover (LULC), slope, soil texture, soil erosion, soil depth, and soil organic carbon. The United States Geological Survey

Data Acquisition
In the present study, nine input layers were used to identify the area most vulnerable to land degradation, namely, the normalized difference vegetation index (NDVI) from a moderate-resolution imaging spectroradiometer (MODIS), land surface temperature (LST) from MODIS, mean annual rainfall from the Climate Hazards Group Infrared Precipitation with Station data (CHIRPS), land use/land cover (LULC), slope, soil texture, soil erosion, soil depth, and soil organic carbon. The United States Geological Survey (USGS) website was used to download the shuttle radar topography mission (SRTM) DEM data (https://earthexplorer.usgs.gov (accessed on 29 June 2022)). A land use/land cover (LULC) map was downloaded from ESRI [41]. MODIS NDVI, MODIS LST, and CHIRPS rainfall were downloaded for the 21-year period of 2001-2021 using Google Earth Engine (GEE). Land degradation, particularly wind erosion and desertification, cannot be assessed during a short period of time, such as a few years or less; it should be evaluated over a long period of time [42]. So, we used long-term mean data (21 years) to evaluate land degradation vulnerability. Soil organic carbon data were downloaded from SoilGrids Land 2023, 12, 106 4 of 20 (https://soilgrids.org/ (accessed on 25 June 2022)). Additionally, data on soil texture, soil erosion, and soil depth were collected at a scale of 1:250,000 from the ICAR-National Bureau of Soil Survey and Land Use Planning (NBSS&LUP), Nagpur. All thematic layers were converted to the WGS 1984 datum and Universal Transverse Mercator (UTM) 43N coordinate system. Table 1 provides an overview of various datasets and their specifications. The downloaded SRTM DEM was reprojected to the UTM coordinate system (43N) and filled in ArcGIS 10.8. A slope map of the basin was then generated using the filled DEM.

Climatic Parameters
Land degradation is caused by interactions between soil, climate, and land-use patterns. In the present study, CHRIPS-based rainfall products (https://developers.google. com/earth-engine/datasets/catalog/UCSB-CHG_CHIRPS_DAILY (accessed on 22 April 2022) at a spatial resolution of 5 km were downloaded for the period of 2001 to 2021 using GEE. The total annual rainfall was computed by calculating the sum of all images for the respective year. These twenty years of annual rainfall data were used for the computation of the long-term mean annual rainfall. The data were reprojected to the UTM 43N coordinate system in ArcGIS 10.8. The downloaded rainfall data were resampled to a 30 m cell size by using the bilinear interpolation method. The bilinear interpolation technique is the most widely used interpolation technique for resampling raster images [35,38,43]. The thematic layer of rainfall was categorized into five subclasses, namely, 234-401, 401-524, 524-685, 685-941, and 941-1654 mm, using the natural break method in ArcGIS. We used the natural break method based on previous research that used it for various hazard mapping applications, such as [44][45][46][47].
In this study, MODIS global LST data (MOD11A2) (https://developers.google.com/ earth-engine/datasets/catalog/MODIS_006_MOD11A2 accessed on 23 April 2022) were downloaded for 2001 to 2021 using GEE. The mean annual LST was computed by calculating the mean of all images of the respective year. These mean annual data of all years were used for the computation of long-term (21 years) mean LST. Using Equation (1), data were converted to degrees Celsius ( • C): The datasets were subsequently reprojected and resampled in the ArcGIS environment. The thematic layer was divided into five subclasses, namely, 27-34, 34-36, 36-37, 37-39, and 39-43 • C, using the natural break method in ArcGIS.

Vegetation Parameters
In the present study, NDVI and LULC were considered important parameters for evaluating the land degradation of the river basin. An LULC map was downloaded for the year 2018 from ESRI and resampled to 30 m in the ArcGIS environment. NDVI data were Land 2023, 12, 106 5 of 20 obtained from NASA Land Processes Distributed Active Archive Center's MODIS products (MOD13Q1) (https://developers.google.com/earthengine/datasets/catalog/MODIS_00 6_MOD13Q1 accessed on 23 April 2022) at spatial and temporal resolutions of 250 m and 16 days, respectively, for a period of 21 years (2001-2021) using GEE. The mean annual NDVI was computed from all images of a particular year. These mean annual data of all years were used for the computation of long-term (21 years) mean NDVI. The vegetative greenness of an area can be assessed using the NDVI, a dimensionless index that varies from −1 to +1. A high NDVI value indicates healthy vegetation, while a low value signifies stressed vegetation. The data were reprojected to UTM 43N as previously described and resampled to a 30 m cell size by using the bilinear interpolation method in ArcGIS. Finally, the layer was categorized into five subclasses, namely, −0.039-0.19, 0.19-0.24, 0.24-0.29, 0.29-0.36, and 0.36-0.60, using the natural break method in ArcGIS.

Soil Parameters
Soil texture, soil erosion class, and soil depth data were taken from ICAR-NBSS&LUP, Nagpur, and resampled to 30 m in the ArcGIS environment. Seven soil texture classes were identified, namely, clay loam, fine loam, loam, loamy skeletal, sandy, sandy skeletal, and rock. For soil erosion, four classes were identified, namely, slight, moderate, severe, and very severe erosion. Six classes were reported for soil depth and named, viz., rock, <25, 25-50, 50-75, 75-100, and >100 cm. Soil organic carbon data were downloaded from soil grids and subsequently reprojected to the UTM 43N coordinate system and resampled to a 30 m cell size. The layer was classified into five subclasses, namely, 0-52, 52-67, 67-93, 93-140, and 140-315 decigrams/kg. Prior to the AHP weightage assignment, all input layers should have a uniform resolution. Therefore, all layers were resampled to be comparable to DEM (30 m). A detailed flowchart of the methodology is displayed in Figure 2.

Analytical Hierarchical Process and Weightage Assignment
The multi-criteria decision analysis (MCDA) using the AHP technique is the most popular and well-known method for identifying land degradation vulnerability [48]. AHP is a pairwise comparison assessment theory, where Saaty's scale of relative im-

Analytical Hierarchical Process and Weightage Assignment
The multi-criteria decision analysis (MCDA) using the AHP technique is the most popular and well-known method for identifying land degradation vulnerability [48]. AHP is a pairwise comparison assessment theory, where Saaty's scale of relative importance is used to compare parameters to one another ( Table 2) [49]. We need comparisons and a scale of numbers that indicate how important one parameter is in comparison to another in terms of the criterion being compared in order to draw organized conclusions about priorities. Intermediate between scale 1 and 3 3 Moderate importance 4 Intermediate between scale 3 and 5 5 Strong importance 6 Intermediate between scale 5 and 7 7 Very strong importance 8 Intermediate between scale 7 and 9 9 Extreme importance The consistency ratio (CR) was used to validate the decision about the pairwise comparison of the various thematic layers and their subclasses [50]. The following equation was used to calculate the CR: where RCI indicates the random consistency index, and Saaty's standard is used to calculate its values (Table 3); CI stands for the consistency index, which was calculated using the following equation: where λ max indicates the principal eigenvalue, and n indicates the total number of input layers used in the LD assessment. In the weighted overlay analysis using AHP, a CR value of ≤0.10 is acceptable. If the CR is greater than 0.10, the decision must be revised to identify the source of the inconsistency, and it must be resolved so ensure CR value below 0.10.

Generating Land Degradation Vulnerability Map
Thematic layers and their subclasses were given AHP-based weightages in ArcGIS to delineate the land degradation vulnerability (LDV) zones. The LDV map in the current study was generated using the following equation: where NDVI, MAR, LST, SE, S, LULC, SOC, ST, and SD indicate the normalized difference vegetation index, mean annual rainfall, land surface temperature, soil erosion, slope, land use/land cover, soil organic carbon, soil texture, and soil depth, respectively; Cwi and Scwi represent the class weight and subclass weight, respectively. Finally, the LDV map was categorized into five classes, viz., very low, low, moderate, high, and very high, using the quantile breaks in ArcGIS. Using the ultra-high definition (UHD) Google Earth imageries of 2022, the LDV map was validated at ten randomly selected sites. Moreover, an ROC curve was generated from the fifty randomly selected sites using Google Earth to validate the results. The ROC curve was used to estimate the area under the curve (AUC), which has values between 0.5 and 1. AUC values closer to 1 indicate excellent model performance whereas, values closer to 0.5 indicate poor prediction accuracy. Validation was also performed using field photographs to determine the agreement with the AHPgenerated LDV map.

Input Thematic Layers and Their Variabilities
The mean NDVI values of the study area range from −0.03 to 0.60, and they were divided into five subclasses, viz., −0.03-0.19, 0.19-0.24, 0.24-0.29, 0.29-0.36, and 0.36-0.60 ( Figure 3A). The highest area (27.9%) falls under subclass −0.03-0.19, followed by 0.24-0.29 (27.4%), and the lowest area falls under subclass 0.36-0.60 (4.9%). The low NDVI class was assigned higher weights, while the high-value class was assigned lower weights. The study of LULC established the baseline information for activities such as the thematic mapping and change detection analysis of earth over time. The study area's LULC map represents cropland (46.9%), current fallow land (21.2%), land with shrub/scrub (20.3%), forest (shrub/scrub/degraded) (1.4%), other wastelands (8.6%), built-up land (urban/rural) (0.5%), deciduous forest (0.8%), water bodies (0.5%), and gullied/ravines (0.01%), and it is depicted in Figure 3B. Most of the area (68.7%) is used for agriculture, which spread acrossthe study area, mostly covering the northeastern and southeastern parts rather than the western part. Around 2.15% area of the river basin is covered by various types of forests, including deciduous and degraded forests, which lie in the southeastern part, and the central part of the study area is covered by land with shrub/scrub (20.3%).
The mean annual rainfall was categorized into five subclasses, namely, very low (234-401 mm), low (401-524 mm), moderate (524-685 mm), high (685-941 mm), and very high (941-1654 mm) rainfall, occupying approximately 42, 30, 20, 5, and 2% of the total geographical area, respectively ( Figure 3C). The study area receives very high precipitation in a small portion of the south and southeastern parts, whereas the western part receives very low rainfall. In the present study, the low rainfall subclass was assigned high weightage and vice versa.
The land surface temperature of the river basin was divided into five subclasses, namely, 27-34, 34-36, 36-37, 37-39, and 39-43 • C ( Figure 3D). About 65 % of the area of the basin falls under the fourth and fifth subclasses of LST (37-39, 39-43 • C), lying in the western and central parts of the study area. The third class (36-37 • C) of LST is distributed mostly in the northern and southeastern parts of the river basin and occupies around 23% of the total study area. Only around 12% of the area falls under the first subclass (27-34 • C), lying mostly in the eastern part of the river basin. A high LST was assigned high weights and vice versa.
The study area was classified into five subclasses of SOC content, namely, 0-52, 52-67, 67-93, 93-140, and 141-315 decigrams/kg ( Figure 4A). The maximum area (42.1%) of the study area falls under subclass 52-67 decigram/kg, followed by subclass 67-93 (26.9%), whereas the lowest area (1.6%) falls under the 141-315 decigram/kg subclass. The northwestern part of the study area has a low SOC as compared to the southeastern part, which is an important factor for the stability of the soil structure. The soil erosion map of the study area is depicted in Figure 4B, and four erosion classes were identified. The study area is mainly composed of moderate soil erosion (56.3%), followed by severe erosion (35.2%), and very severe erosion occurs the least (4.2%). Severe and very severe erosion are mainly concentrated in the western part, and moderate erosion is found in the eastern part of the study area. The soil depth map is presented in Figure 4C, and six subclasses were identified according to the ICAR-NBSS&LUP, Nagpur, India. The highest area falls under subclass >100 cm (59%), and the lowest area falls under subclass <25 cm, with an area of 1.9%. The soil in the study area was classified into different texture classes, namely, clay/fine loam, loamy skeletal, sandy skeletal, and sandy soils ( Figure 4D). The study area is mainly composed of sandy soil (49.5%), followed by clay loam (24.4%), with loam representing the least amount (1.1%). The western part of the study area is dominated by sandy soil, whereas the eastern part is dominated by fine loam and clay loam texture.
high (941-1654 mm) rainfall, occupying approximately 42, 30, 20, 5, and 2% of the total geographical area, respectively ( Figure 3C). The study area receives very high precipitation in a small portion of the south and southeastern parts, whereas the western part receives very low rainfall. In the present study, the low rainfall subclass was assigned high weightage and vice versa.
The land surface temperature of the river basin was divided into five subclasses, namely, 27-34, 34-36, 36-37, 37-39, and 39-43 °C ( Figure 3D). About 65 % of the area of the basin falls under the fourth and fifth subclasses of LST (37-39, 39-43 °C), lying in the western and central parts of the study area. The third class (36-37 °C) of LST is distributed mostly in the northern and southeastern parts of the river basin and occupies around 23% of the total study area. Only around 12% of the area falls under the first subclass (27-34 °C), lying mostly in the eastern part of the river basin. A high LST was assigned high weights and vice versa.  The study area's slope was divided into six subclasses, namely, <3, 3-8, 8-15, 15-30, 30-50, and >50% ( Figure 5). Almost 62% of the area of the basin falls under the slope subclass of <3%, which is nearly flat and mainly concentrated in the northern and eastern parts of the basin. The slope of the remaining 38% of the area falls under subclass >3%. The eastern and southeastern parts are mainly composed of high slopes.

Land Degradation Vulnerability
Before the thematic layer integration, consistency ratios for each thematic layer (Table 4), normalized matrix (Table 5), and thematic layer subclass (Table 6) were computed. The findings show that the judgment matrices used in the analyses were reasonably consistent and accurate (CR ≤ 0.10). The NDVI, LST, rainfall, topography, and pedological parameters were computed and the final weightage assigned to each parameter was Land 2023, 12, 106 9 of 20 computed using the AHP model. The weighted overlay approach was used to integrate the thematic layers according to their associated weights. AHP found that NDVI had the highest priority in land degradation vulnerability zone (LDVZ) identification, with a weight value of 0.27, followed by rainfall (0.22), LST (0.15), erosion (0.12), slope (0.08), LULC (0.06), SOC (0.22), texture (0.22), and depth (0.22). In this study, with the help of the AHP-and the GIS-based modeling approach, five LDVZs, namely, very low, low, moderate, high, and very high, were identified ( Figure 6). The findings show that 50.3% of the river basin area (34.72 lakh ha) is highly vulnerable to land degradation, and it was classified in the very low to low class, accounting for nearly half of the river basin (Figure 7). These very low to low classes of LDV cover the eastern, southern, and northeastern parts of the study area. About 11.02 lakh ha (16.0%) of the river basin area falls under the moderate class of land degradation, mostly covering the central to southeastern parts of the study area. The high and very high classes of land degradation vulnerability cover about 23.19 lakh ha (33.7%) of the river basin area. The western and somewhat central portions of the study region mostly fall under these classes. These two classes particularly indicate wind erosion, which is one of the most severe forms of land degradation.
Land 2023, 12, x FOR PEER REVIEW 9 of 21 The study area was classified into five subclasses of SOC content, namely, 0-52, 52-67, 67-93, 93-140, and 141-315 decigrams/kg ( Figure 4A). The maximum area (42.1%) of the study area falls under subclass 52-67 decigram/kg, followed by subclass 67-93 (26.9%), whereas the lowest area (1.6%) falls under the 141-315 decigram/kg subclass. The northwestern part of the study area has a low SOC as compared to the southeastern part, which is an important factor for the stability of the soil structure. The soil erosion map of the study area is depicted in Figure 4B, and four erosion classes were identified. The study area is mainly composed of moderate soil erosion (56.3%), followed by severe erosion (35.2%), and very severe erosion occurs the least (4.2%). Severe and very severe erosion are mainly concentrated in the western part, and moderate erosion is found in the eastern part of the study area. The soil depth map is presented in Figure 4C, and six subclasses were identified according to the ICAR-NBSS&LUP, Nagpur, India. The highest area falls under subclass >100 cm (59%), and the lowest area falls under subclass <25 cm, with an area of 1.9%. The soil in the study area was classified into different texture classes, namely, clay/fine loam, loamy skeletal, sandy skeletal, and sandy soils ( Figure  4D). The study area is mainly composed of sandy soil (49.5%), followed by clay loam (24.4%), with loam representing the least amount (1.1%). The western part of the study area is dominated by sandy soil, whereas the eastern part is dominated by fine loam and clay loam texture.  The study area's slope was divided into six subclasses, namely, <3, 3-8, 8-15, 15-30, 30-50, and >50% ( Figure 5). Almost 62% of the area of the basin falls under the slope subclass of <3%, which is nearly flat and mainly concentrated in the northern and eastern parts of the basin. The slope of the remaining 38% of the area falls under subclass >3%. The eastern and southeastern parts are mainly composed of high slopes.

Land Degradation Vulnerability
Before the thematic layer integration, consistency ratios for each thematic layer (Table 4), normalized matrix (Table 5), and thematic layer subclass (Table 6) were computed. The findings show that the judgment matrices used in the analyses were reasonably consistent and accurate (CR ≤ 0.10). The NDVI, LST, rainfall, topography, and pedological parameters were computed and the final weightage assigned to each parameter was computed using the AHP model. The weighted overlay approach was used to integrate the thematic layers according to their associated weights. AHP found that NDVI had the highest priority in land degradation vulnerability zone (LDVZ) identification, with a weight value of 0.27, followed by rainfall (0.22), LST (0.15), erosion (0.12), slope (0.08), LULC (0.06), SOC (0.22), texture (0.22), and depth (0.22). In this study, with the help of the AHP-and the GIS-based modeling approach, five LDVZs, namely, very low, low, moderate, high, and very high, were identified ( Figure 6). The findings show that 50.3% of the river basin area (34.72 lakh ha) is highly vulnerable to land degradation, and it was classified in the very low to low class, accounting for nearly half of the river basin ( Figure  7). These very low to low classes of LDV cover the eastern, southern, and northeastern parts of the study area. About 11.02 lakh ha (16.0%) of the river basin area falls under the moderate class of land degradation, mostly covering the central to southeastern parts of the study area. The high and very high classes of land degradation vulnerability cover about 23.19 lakh ha (33.7%) of the river basin area. The western and somewhat central

Validation of Land Degradation Vulnerability Zones
The validation of the LDVZs of the study area was conducted using two approaches, namely, field photographs and Google Earth images. Figure 8 indicates good agreement between the LDV map of the study area with the field photographs. In the second approach, validation was carried out with ten randomly selected high-resolution Google Earth images (Figure 9). In addition to this, fifty points were randomly selected in Google Earth to plot the ROC curve, and the AUC value was found to be 82% (Figure 10). The AHP model predicts land degradation vulnerability zones in the research area with reasonable accuracy, as shown by the ROC curve.

Validation of Land Degradation Vulnerability Zones
The validation of the LDVZs of the study area was conducted using two approaches, namely, field photographs and Google Earth images. Figure 8 indicates good agreement between the LDV map of the study area with the field photographs. In the second approach, validation was carried out with ten randomly selected high-resolution Google Earth images (Figure 9). In addition to this, fifty points were randomly selected in Google Earth to plot the ROC curve, and the AUC value was found to be 82% (Figure 10). The AHP model predicts land degradation vulnerability zones in the research area with reasonable accuracy, as shown by the ROC curve.

Validation of Land Degradation Vulnerability Zones
The validation of the LDVZs of the study area was conducted using two approaches, namely, field photographs and Google Earth images. Figure 8 indicates good agreement between the LDV map of the study area with the field photographs. In the second approach, validation was carried out with ten randomly selected high-resolution Google Earth images (Figure 9). In addition to this, fifty points were randomly selected in Google Earth to plot the ROC curve, and the AUC value was found to be 82% (Figure 10). The AHP model predicts land degradation vulnerability zones in the research area with reasonable accuracy, as shown by the ROC curve.

Discussion
Assessments of land degradation status are important for sustainable agricultural planning and development in an ecosystem. Since land degradation has recently accelerated, spatial risk mapping of land degradation is required in order to help authorities make trustworthy and fair decisions on ecosystem rehabilitation or restoration and investment priorities [51]. In the hot, arid western Rajasthan, wind erosion causes serious land degradation and has an impact on the crop production system [19]. For sustainable natural resource management, the identification of areas sensitive to wind erosion is urgently required. Therefore, using an AHP-and GIS-based modeling technique, the current research was carried out to identify land degradation vulnerability zones in the Luni River basin.
In the present investigation, NDVI, MAR, LST, erosion class, slope, LULC, SOC, soil depth, and soil texture were taken into consideration for land degradation vulnerability mapping. The most important layers for assessing vulnerability to land degradation were determined to be NDVI and MAR. It has long been known that NDVI is a valuable tool for assessing the vegetation cover of any study area, and it is widely acknowledged in research that a decline in NDVI is an indication of land degradation and is directly related to climatic factors [52,53]. Rainfall is the most significant factor for wind erosion, as low-rainfall zones are highly vulnerable to sand movement and, consequently, wind erosion [54,55]. LST is a crucial parameter in arid and semi-arid zones due to its direct relationship with soil moisture availability, and it is also a crucial parameter in this study because of its indirect relationship with the vegetation cover of the research area [56,57]. Increased land degradation and decreased vegetation greenness could be the effects of an increasing LST. LULC is a key contributor to land degradation since it indicates both natural and manmade changes in land surfaces [58]. SOC is an important soil parameter that acts as a binding agent and stabilizes the soil structure [59]. Therefore, the depletion of soil organic carbon can reduce soil aggregate stability and, as a consequence, increase soil erodibility [59,60]. Soil texture plays a significant role in the susceptibility of a soil

Discussion
Assessments of land degradation status are important for sustainable agricultural planning and development in an ecosystem. Since land degradation has recently accelerated, spatial risk mapping of land degradation is required in order to help authorities make trustworthy and fair decisions on ecosystem rehabilitation or restoration and investment priorities [51]. In the hot, arid western Rajasthan, wind erosion causes serious land degradation and has an impact on the crop production system [19]. For sustainable natural resource management, the identification of areas sensitive to wind erosion is urgently required. Therefore, using an AHP-and GIS-based modeling technique, the current research was carried out to identify land degradation vulnerability zones in the Luni River basin.
In the present investigation, NDVI, MAR, LST, erosion class, slope, LULC, SOC, soil depth, and soil texture were taken into consideration for land degradation vulnerability mapping. The most important layers for assessing vulnerability to land degradation were determined to be NDVI and MAR. It has long been known that NDVI is a valuable tool for assessing the vegetation cover of any study area, and it is widely acknowledged in research that a decline in NDVI is an indication of land degradation and is directly related to climatic factors [52,53]. Rainfall is the most significant factor for wind erosion, as low-rainfall zones are highly vulnerable to sand movement and, consequently, wind erosion [54,55]. LST is a crucial parameter in arid and semi-arid zones due to its direct relationship with soil moisture availability, and it is also a crucial parameter in this study because of its indirect relationship with the vegetation cover of the research area [56,57]. Increased land degradation and decreased vegetation greenness could be the effects of an increasing LST. LULC is a key contributor to land degradation since it indicates both natural and manmade changes in land surfaces [58]. SOC is an important soil parameter that acts as a binding agent and stabilizes the soil structure [59]. Therefore, the depletion of soil organic carbon can reduce soil aggregate stability and, as a consequence, increase soil erodibility [59,60]. Soil texture plays a significant role in the susceptibility of a soil surface to wind erosion [61,62]. Sandy soils are intrinsically more erodible than fine-textured soil due to the lack of clay and silt, which are necessary for physical crusting and soil aggregation [63,64]. Similarly, soil depth is also an important parameter, as shallow soil is more vulnerable to soil loss than deep soil [65]. We used slope as an input in the present study because topographic relief and surface fragmentation both increase with an increase in slopes, which results in a significant increase in surface roughness [66]. So, high-slope areas are less vulnerable to wind erosion.
Finally, weights were assigned to the thematic layers based on their importance, i.e., NDVI (0.27), MAR (0.22), LST (0.15), erosion class (0.12), slope (0.08), LULC (0.06), SOC (0.04), soil texture (0.03), and soil depth (0.02). Since NDVI is the most significant indicator of vegetal degradation, we identified it as the first most crucial parameter in the hierarchy. The mean annual rainfall was selected as the second most important parameter because low rainfall leads to a dry climate and, consequently, more wind erosion. The land surface temperature was chosen as the third layer in the hierarchy because regions with higher LSTs are assumed to have less vegetative cover and an arid climate compared to areas with lower LSTs. The soil erosion class was chosen as the fourth element in the hierarchy because soil erosion is a direct indicator of land degradation in the study area. Lower ranks in the hierarchy were given to the remaining layers. As a result, a pairwise comparison matrix was used to compare all layers. Before the integration of the thematic layers, consistency ratios for each thematic layer and its subclasses were computed. The computed CR value was ≤0.10, demonstrating the accuracy of all the parameters' predictions regarding how they affect soil erosion.
The land degradation vulnerability (LDV) map generated using AHP-and GIS-based modeling was categorized into five subclasses, viz., very low, low, moderate, high, and very high, occupying 26%, 24.3%, 16%, 12.3%, and 21.4% of the river basin area, respectively. The result shows that the very low to low category of LDV covered half (50.3%) of the study area, predominantly covering the eastern and southern parts of the basin. These categories are mainly associated with high vegetative cover (NDVI >0.24), high to very high rainfall (685-1654 mm), LST <39 • C, deep soils with a clay loam to loam texture and a high organic carbon content, and slight-to-moderate soil erosion. Tolche et al. 2021 [35] also reported that deep soils with acidic-to-neutral pH were very slightly and slightly vulnerable to LD. High NDVI with optimized LULC can improve SOC, which acts as a binding agent and stabilizes soil structure and reduces soil erosion losses [67]. Similarly, many researchers have reported that very low to low categories of LDV have good vegetative coverage with open forests, very low vegetative degradation, adequate rainfall, and well-drained soils [38,43]. About 16% of the study area falls under moderate LDV, mainly concentrated in the central part of the basin. This region is characterized by less vegetative coverage, the scrub/shrub class of LULC, moderate soil erosion, medium rainfall (401-524 mm), medium LST, and sandy soil. Mzuri et al. 2021 [68] and Tolche et al. 2021 [35] also reported that a moderate level of LD was mainly associated with less vegetative cover, insufficient rainfall, a moderate temperature, and sparse shrubland.
The high and very high LDV classes covered about 33.7% of the river basin area mostly comprising the western part of the basin. These two classes showed high to very high severity of land degradation, such as sand dunes and wind erosion, with sparse vegetation conditions. This study area zone had very low vegetation cover (NDVI < 0.19), poor rainfall (<400 mm), a higher LST (39-43 • C), and sandy to sandy skeletal texture soils with low soil organic carbon (0-52 decigram/kg) and severe soil erosion. In the present study, NDVI and rainfall were the most sensitive input parameters. Rukhovich et al. (2021) [69] found that areas with low NDVI values are more prone to land degradation, where the probability of assessing LD using indicators such as NDVI computed from Landsat imagery was 87.5%. Similar to this, it was determined that the assessment of LD was more consistent with spatial and temporal variations in NDVI, especially in arid regions, than with other hydro-climatic indicators [70]. Due to the high temperatures in these locations, the little available precipitation will evaporate, preventing vegetation growth and reducing rain use effectiveness [71,72]. Similar findings were obtained in various studies on LDV for India's semi-arid ecosystem, where low vegetation cover, little rainfall, and high temperatures make the area extremely vulnerable to LD [38,43]. The prioritization of the river basin was carried out to identify degraded land with high erosion activity so that suitable conservation measures can be taken to minimize LD [73]. The study area can be categorized into three priority zones: high priority (the very high and high LDV class), medium priority (the moderate LDV class), and low priority (the low and very low LDV class).
High-resolution Google Earth images and field photographs were used to validate the land degradation vulnerability zones, and the results are in excellent agreement with the model-based AHP-GIS approach. In earlier studies, similar types of validation approaches have been used by many researchers with respect to land degradation [38,43] and soil erosion hazards [74]. Moreover, the ROC curve has been used to assess the accuracy of AHP models [75][76][77]. ROC curves are commonly used to evaluate test accuracy, as they plot a test's true positive rate versus its false positive rate [78,79]. Vulnerable to land degradation and non-vulnerable areas were selected from Google Earth images and used them to test the performance of AHP in river basins [38]. The AUC value of the ROC curve in the current investigation was observed to be 82%. Therefore, it could be determined that the AHP model produces satisfactory results in estimating the zones vulnerable to LD in the river basin.
In this research, the AHP-and GIS-based modeling approach showed its usefulness for the assessment of vulnerability to LD by combining different parameters. As expected, arid and semi-arid terrains are generally sensitive to LD; thus, it is important to consider both scientific and policy approaches [80]. The methodology followed in the present study could be used as a tool to guide decision makers in the prioritization of the river basin. The results of this study suggest that very high and high LDV zones should be considered hotspots for initiatives suitable for soil and water conservation measures and sustainable land resource management.

Conclusions
In the present study, nine thematic layers, namely, NDVI, LULC, rainfall, LST, slope, SOC, soil erosion, soil texture, and depth, were taken into consideration for LDVZ identification and mapping in the Luni River basin. The analysis revealed that 33.7% of the area falls under high to very high vulnerability, followed by 16% of the area falling under moderate vulnerability and 50.3% of the area falling under low to very low vulnerability. The validation of the LDVZ with Google Earth images and field photographs clearly showed that the remote sensing data combined with AHP adequately distinguished the sites prone to land degradation in the study basin in a cost-effective and time-efficient manner. Additionally, the ROC curve analysis, with an area under the ROC curve value of 82%, validated the AHP method's potential to accurately estimate the LD vulnerability zones in the research area. Soil and water conservation structures that are appropriate could be suggested for the regions that are extremely and very highly prone to LD. Multi-criteria decision analyses and geospatial techniques can be used as tools for prioritization management in order to achieve LD neutrality in arid and semi-arid regions.