Quantitative Assessment of Environmental Sensitivity to Desertification Using the Modified MEDALUS Model in a Semiarid Area

Iran is mainly located in the arid and semiarid climate zone and seriously affected by desertification. This is a severe environmental problem, which results in a persistent loss of ecosystem services that are fundamental to sustaining life. Process understanding of this phenomenon through the evaluation of important drivers is, however, a challenging work. The main purpose of this study was to perform a quantitative evaluation of the current desertification status in the Segzi Plain, Isfahan Province, Iran, through the modified Mediterranean Desertification and Land Use (MEDALUS) model and GIS. In this regard, five main indicators including soil, groundwater, vegetation cover, climate, and erosion were selected for estimating the environmental sensitivity to desertification. Each of these qualitative indicators is driven by human interference and climate. After statistical analysis and a normality test for each indicator data, spatial distribution maps were established. Then, the maps were scored in the MEDALUS approach, and the current desertification status in the study area from the geometric mean of all five quality indicators was created. Based on the results of the modified MEDALUS model, about 23.5% of the total area can be classified as high risk to desertification and 76.5% classified as very high risk to desertification. The results indicate that climate, vegetation, and groundwater quality are the most important drivers for desertification in the study area. Erosion (wind and water) and soil indices have minimal importance.


Introduction
Desertification is a global ecological threat, which is defined as the depletion of biological and economic productivity of both rainfed and irrigated agricultural land, rangeland, forest, and shrubland, through a combination of various human and natural factors [1][2][3][4]. The process of desertification occurs in climates with severe humidity deficits, i.e., arid, semiarid, and dry sub-humid areas [5,6]. Studies confirm that between 10 and 20% of arid and semiarid areas in the world [7,8] and approximately 40% of the total Earth's land surface suffer from desertification [9,10]. According to the Forest, Range, and Watershed Management Organization of Iran, more than 80% of Iran's land is directly or indirectly affected by desertification [11].
Considering the complexity of interactions between the driving forces of desertification, assessment of land sensitivity to degradation is one of the most important research efforts in different parts of the world [12,13]. During the last decades, various approaches and models have been developed to assess sensitivity of a region to desertification and Sustainability 2021, 13, x FOR PEER REVIEW 3 of 20

General Applied MEDALUS Framework
Desertification is a complex land degradation process that eventually leads to the depletion of land productivity [24]. To assess the sensitivity of desertification, the MED-ALUS model was used based on five main quality indices specific to the environmentally sensitive areas (ESAs) including soil, groundwater, vegetation, climate, and erosion. These indices were computed based on a comprehensive data set shown in Figure 2. Based on the modified MEDALUS model, the value of each indicator was determined between 100 and 200, where 100 was assigned to areas of the least sensitivity and 200 to areas with the highest sensitivity to desertification [38]. Each index was categorized into five classes based on sensitivity to desertification including very low, low, moderate, high, and very high (Table 1) [32]. After calculating the indices, ArcGIS software (10.3) (Esri, Redlands, CA, USA) was used to depict the spatial distribution of each index and indicator throughout the study area. Inverse distance weighting (IDW) was applied for interpolation. Finally, the total sensitivity of desertification in the Segzi plain was calculated from the geometric mean of indices, according to:

MEDALUS = (SQI × GWGI × VQI × CQI × EQI)
where SQI is soil quality index, GWQI is groundwater quality index, VQI is vegetation quality index, CQI is climate quality index, and EQI is erosion quality index. The calculation of each index is described below.

General Applied MEDALUS Framework
Desertification is a complex land degradation process that eventually leads to the depletion of land productivity [24]. To assess the sensitivity of desertification, the MEDALUS model was used based on five main quality indices specific to the environmentally sensitive areas (ESAs) including soil, groundwater, vegetation, climate, and erosion. These indices were computed based on a comprehensive data set shown in Figure 2. Based on the modified MEDALUS model, the value of each indicator was determined between 100 and 200, where 100 was assigned to areas of the least sensitivity and 200 to areas with the highest sensitivity to desertification [38]. Each index was categorized into five classes based on sensitivity to desertification including very low, low, moderate, high, and very high (Table 1) [32]. After calculating the indices, ArcGIS software (10.3) (Esri, Redlands, CA, USA) was used to depict the spatial distribution of each index and indicator throughout the study area. Inverse distance weighting (IDW) was applied for interpolation. Finally, the total sensitivity of desertification in the Segzi plain was calculated from the geometric mean of indices, according to: where SQI is soil quality index, GWQI is groundwater quality index, VQI is vegetation quality index, CQI is climate quality index, and EQI is erosion quality index. The calculation of each index is described below.   The SQI is a main factor in the desertification process, which affects the soil state due to strength of cohesion between soil particles, water retention capacity, texture, and structure [39,40]. To evaluate SQI, 50 samples in total were collected from a soil depth of 0-20 cm from five different land units ( Table 2). The location of soil samples and land units is shown in Figure 3. Next, the SQI was computed based on nine indicators including pH, soil electric conductivity (EC), soil texture (T), soil gypsum (GYP), soil organic matter (OM), chloride (Cl -), sodium (Na + ), bicarbonate (HCO3 -), and sodium adsorption ratio (SAR).  [38].

Characteristics
Name No. Physiographical features are related to Quaternary Period. The process of water erosion has made them worn and their remains appear as complications.
High plateaus and terraces 1 Flat plains have been created due to accumulation of fine sediments transported from the foot of mountains.
Piedmont alluvial plains 2 The river creates this unit during overflows. River alluvial plains 3 Low or concave slopes are located in the center of the intermountain plains. Lowlands 4 Foothill slopes. Coarse and fine sediments are transferred to these lands by streams from the catchment area.
Gravelly alluvial fans 5  The SQI is a main factor in the desertification process, which affects the soil state due to strength of cohesion between soil particles, water retention capacity, texture, and structure [39,40]. To evaluate SQI, 50 samples in total were collected from a soil depth of 0-20 cm from five different land units ( Table 2). The location of soil samples and land units is shown in Figure 3. Next, the SQI was computed based on nine indicators including pH, soil electric conductivity (EC), soil texture (T), soil gypsum (GYP), soil organic matter (OM), chloride (Cl − ), sodium (Na + ), bicarbonate (HCO 3 − ), and sodium adsorption ratio (SAR). Table 2. Characteristics of the physiographic units in the study area [38].

Characteristics Name
No.
Physiographical features are related to Quaternary Period. The process of water erosion has made them worn and their remains appear as complications.
High plateaus and terraces 1 Flat plains have been created due to accumulation of fine sediments transported from the foot of mountains.
Piedmont alluvial plains 2 The river creates this unit during overflows. River alluvial plains 3 Low or concave slopes are located in the center of the intermountain plains. Lowlands 4 Foothill slopes. Coarse and fine sediments are transferred to these lands by streams from the catchment area.
The SQI was calculated based on the geometric mean of indicators as: The attributed weights for each soil indicator are presented in Table 3.

Groundwater Quality Index (GWQI) and Indicators
Groundwater is the primary source of water supply in arid and semiarid regions [41]. Reduced groundwater quality is one of the major causes of land degradation in these re-  (Table 2) and soil sampling.
The SQI was calculated based on the geometric mean of indicators as: The attributed weights for each soil indicator are presented in Table 3. Table 3. Weighted scores and classes of soil indicators [38].

Groundwater Quality Index (GWQI) and Indicators
Groundwater is the primary source of water supply in arid and semiarid regions [41]. Reduced groundwater quality is one of the major causes of land degradation in these regions [42]. The assessment of groundwater quality index can help policymakers to develop sustainable management of semiarid areas. Therefore, to evaluate GWQI, 17 samples in total were collected from wells in the study area. The WGQI was calculated based on five indicators including: sodium adsorption ratio (SAR), electric conductivity (EC), chloride (Cl − ), pH, and total dissolved solids (TDS). The GWQI was assessed based on the geometric mean of indicators according to: The weights attributed to each water indicator are presented in Table 4.

Vegetation Quality Index (VQI) and Indicators
Vegetation cover is an important biotic land component in the evaluation of desertification in semiarid areas [43]. Vegetation cover reduces runoff by promoting water infiltration. Besides, the vegetation cover enriches the surface soil with organic matter and improves its structure and cohesion, thus providing a high capacity to reduce the intensity of soil erosion [36]. In this study, the VQI was assessed by three different indicators including drought resistance (DR), fire risk (FR), and plant cover (PC). The erosion protection, drought resistance, and fire risk indicators were derived from the land use/land cover map, and the plant cover indicator was extracted from the normalized difference vegetation index (NDVI). The weighted scores and classes of vegetation cover indicators are presented in Table 5. The VQI was calculated by:

Climate Quality Index (CQI) and Indicators
The spatiotemporal distribution of precipitation and evapotranspiration are the most important parameters in hydrology and ecology [44,45]. Precipitation plays a significant role in land degradation and desertification, as it determines runoff and soil water content [43]. The CQI was calculated based on three indicators including: precipitation (P), potential evapotranspiration (ETp), and aridity index (AI). Two climatic stations, located in the study area, were used to assess the precipitation and evapotranspiration indicators. Furthermore, aridity index was computed based on annual precipitation and potential evapo-transpiration. The weighted scores and classes of CQI indicators are presented in Table 6. The CQI was calculated according to:

Erosion Quality Index (EQI) and Indicators
The EQI was assessed using wind and water erosion indicators. The wind erosion and water erosion indicators were computed using Iranian Research Institute of Forest and Rangelands (IRIFR) [46] and the Pacific South-West Inter-Agency Committee (PSIAC) [47] methods, respectively. To calculate the wind erosion, nine parameters that affect the wind erosion process were considered. These parameters included lithology, morphology, wind velocity, soil characteristics, type and plant cover percentage, wind erosion features, soil moisture, type and distribution of sand dunes, land use, and land management. To evaluate the water erosion, the following parameters were used: lithology, soil, climate, runoff, morphology, vegetation cover, and land use. Then, the wind and water erosion indicators were categorized into five classes as shown in Table 7. Ultimately, the EQI was computed as:

Soil Quality Index
The summary statistics of soil parameters such as minimum, maximum, mean, and standard deviation are presented in Table 8. The EC values ranged from 2.5 to 174.7 ds/m in the study area. The Cl − , HCO 3 − , and Na + concentrations varied from 17.5 to 2766 meq/L, 5 to 30 meq/L, and 4 to 10,946 meq/L, respectively. The average pH of the soil was 8, which is a weak alkaline soil with a small spatial variation. Moreover, the average organic matter and gypsum of soil samples was 0.5% and 9.5%, and their maximum was 14.6 and 16.5, respectively. The spatial distribution of soil indicators and soil quality index based on MEDALUS is shown in Figures 4 and 5, respectively. The results show that 18.9% (143.8 km 2 ), 54.6% (415.1 km 2 ), and 26.5% (201.6 km 2 ) of the study area can be classified into low, moderate, and high sensitivity to desertification, respectively. According to the spatial distribution of soil quality index, the central part of the region has high degradation of SQI due to the accumulation of salts in the soil, and very low content of organic matter (Figure 4h). This area also has a high gypsum content (Figure 4f), which can cause soil to lose its structure and become more susceptible to water and wind erosion [38]. Due to the significant influence of SQI on vegetation growth, this index also indirectly indicates desertification. According to Figure 4, the Na + , EC, Cl − , and SAR have a similar spatial distribution with high values in the central part of the region corresponding to a desert area. However, other districts, especially agricultural fields, are in the southwest of the study area. Moreover, the soil quality indicators show that texture and organic matter with 125 and 157 as mean value, respectively, are the least and the most sensitive indicators to desertification.

Soil Quality Index
The summary statistics of soil parameters such as minimum, maximum, mean, and standard deviation are presented in Table 8. The EC values ranged from 2.5 to 174.7 ds/m in the study area. The Cl -, HCO3 -, and Na + concentrations varied from 17.5 to 2766 meq/L, 5 to 30 meq/L, and 4 to 10,946 meq/L, respectively. The average pH of the soil was 8, which is a weak alkaline soil with a small spatial variation. Moreover, the average organic matter and gypsum of soil samples was 0.5% and 9.5%, and their maximum was 14.6 and 16.5, respectively. The spatial distribution of soil indicators and soil quality index based on MEDALUS is shown in Figures 4 and 5, respectively. The results show that 18.9% (143.8 km 2 ), 54.6% (415.1 km 2 ), and 26.5% (201.6 km 2 ) of the study area can be classified into low, moderate, and high sensitivity to desertification, respectively. According to the spatial distribution of soil quality index, the central part of the region has high degradation of SQI due to the accumulation of salts in the soil, and very low content of organic matter (Figure 4h). This area also has a high gypsum content (Figure 4f), which can cause soil to lose its structure and become more susceptible to water and wind erosion [38]. Due to the significant influence of SQI on vegetation growth, this index also indirectly indicates desertification. According to Figure 4, the Na + , EC, Cl − , and SAR have a similar spatial distribution with high values in the central part of the region corresponding to a desert area. However, other districts, especially agricultural fields, are in the southwest of the study area. Moreover, the soil quality indicators show that texture and organic matter with 125 and 157 as mean value, respectively, are the least and the most sensitive indicators to desertification.

Groundwater Quality Index
A summary statistical description of groundwater quality parameter Table 9. The groundwater in this area is slightly alkaline to neutral, as pH to 8, with a mean of 7.5. The electrical conductivity (EC) varies from 1465 t with a mean of 4752 µS/cm. The total dissolved solids ranges from 1025 to average concentration of Clis 32.6 meq/L. The spatial distribution of groundwater quality index and its indicato Figures 6 and 7, respectively. The results of GWQI show that groundwat

Groundwater Quality Index
A summary statistical description of groundwater quality parameters is presented in Table 9. The groundwater in this area is slightly alkaline to neutral, as pH ranges from 7.1 to 8, with a mean of 7.5. The electrical conductivity (EC) varies from 1465 to 12,648 µS/cm, with a mean of 4752 µS/cm. The total dissolved solids ranges from 1025 to 8853 mg/L. The average concentration of Clis 32.6 meq/L. The spatial distribution of groundwater quality index and its indicators are shown in Figures 6 and 7, respectively. The results of GWQI show that groundwater in the area is

Groundwater Quality Index
A summary statistical description of groundwater quality parameters is presented in Table 9. The groundwater in this area is slightly alkaline to neutral, as pH ranges from 7.1 to 8, with a mean of 7.5. The electrical conductivity (EC) varies from 1465 to 12,648 µS/cm, with a mean of 4752 µS/cm. The total dissolved solids ranges from 1025 to 8853 mg/L. The average concentration of Cl − is 32.6 meq/L. The spatial distribution of groundwater quality index and its indicators are shown in Figures 6 and 7, respectively. The results of GWQI show that groundwater in the area is not in a good condition, as 6.5% (49.5 km 2 ) of the area can be classified as moderate and 90.9% (691.8 km 2 ) as highly sensitive to desertification. A small percentage (2.6%) of the study area is in the very high class based on the GWQI. not in a good condition, as 6.5% (49.5 km 2 ) of the area can be classified as moderate and 90.9% (691.8 km 2 ) as highly sensitive to desertification. A small percentage (2.6%) of the study area is in the very high class based on the GWQI.  The results demonstrate that almost the whole study area can be classified with high sensitivity to desertification based on groundwater indicators (Figure 7). The EC and TDS have the greatest effect on desertification. The SAR with a mean of 120 mean value has the lowest effect on land degradation in the study area.

Vegetation Quality Index
Because of the dry climate along with alkalinity and salinity of the soil, the situation for vegetation is not good. The areas in the very high desertification class lack vegetation because of the high amount of saline and alkaline soil. Among the vegetation quality indices, plant cover (mean 194) has the largest impact for the desertification in the area. Fire risk (mean 156) is the least sensitive indicator to desertification (Figure 8). The results demonstrate that almost the whole study area can be classified with high sensitivity to desertification based on groundwater indicators (Figure 7). The EC and TDS have the greatest effect on desertification. The SAR with a mean of 120 mean value has the lowest effect on land degradation in the study area.

Vegetation Quality Index
Because of the dry climate along with alkalinity and salinity of the soil, the situation for vegetation is not good. The areas in the very high desertification class lack vegetation because of the high amount of saline and alkaline soil. Among the vegetation quality indices, plant cover (mean 194) has the largest impact for the desertification in the area. Fire risk (mean 156) is the least sensitive indicator to desertification (Figure 8). The results demonstrate that almost the whole study area can be classified with high sensitivity to desertification based on groundwater indicators (Figure 7). The EC and TDS have the greatest effect on desertification. The SAR with a mean of 120 mean value has the lowest effect on land degradation in the study area.

Vegetation Quality Index
Because of the dry climate along with alkalinity and salinity of the soil, the situation for vegetation is not good. The areas in the very high desertification class lack vegetation because of the high amount of saline and alkaline soil. Among the vegetation quality indices, plant cover (mean 194) has the largest impact for the desertification in the area. Fire risk (mean 156) is the least sensitive indicator to desertification (Figure 8). According to the vegetation map (Figure 9) of this area, 261.2 km 2 high class and 483.7 km 2 (63.6%) in very high risk of desertification. How percentage (2%) of the study area is in the low class for vegetation quali

Climate Quality Index
In terms of climate quality index, the whole area is in the very hig desertification. Rainfall average in the study area is very small (mean 1 relation to potential evapotranspiration (mean 2920 mm per year). The c (precipitation, potential evapotranspiration, and aridity indicators) for t in Figure 10, are not significantly different. Figure 11 demonstrates that classified as very high in terms of precipitation, potential evapotranspir According to the vegetation map ( Figure 9) of this area, 261.2 km 2 (34.3%) are in the high class and 483.7 km 2 (63.6%) in very high risk of desertification. However, a negligible percentage (2%) of the study area is in the low class for vegetation quality index.

Climate Quality Index
In terms of climate quality index, the whole area is in the very high class of risk for desertification. Rainfall average in the study area is very small (mean 125.7 mm/year) in relation to potential evapotranspiration (mean 2920 mm per year). The climate indicators (precipitation, potential evapotranspiration, and aridity indicators) for the region, shown in Figure 10, are not significantly different. Figure 11 demonstrates that the entire area is classified as very high in terms of precipitation, potential evapotranspiration, and aridity

Climate Quality Index
In terms of climate quality index, the whole area is in the very high class of risk for desertification. Rainfall average in the study area is very small (mean 125.7 mm/year) in relation to potential evapotranspiration (mean 2920 mm per year). The climate indicators (precipitation, potential evapotranspiration, and aridity indicators) for the region, shown in Figure 10, are not significantly different. Figure 11 demonstrates that the entire area is classified as very high in terms of precipitation, potential evapotranspiration, and aridity indicators. Moreover, the results of climate quality indicators show that potential evapotranspiration (mean 197) is the most sensitive indicator to desertification (Figure 10).
indicators. Moreover, the results of climate quality indicators show that potential evapotranspiration (mean 197) is the most sensitive indicator to desertification ( Figure 10).   Figures 12 and 13 show maps of erosion indicators and erosion quality index. Based on wind erosion, almost 204 km 2 (27%) of the study area are located in the low, 535 km 2 (70%) in moderate, and 21 km 2 (3%) in the high class of desertification ( Figure 12). More- Figure 11. Spatial distribution of climate quality index. Figures 12 and 13 show maps of erosion indicators and erosion quality index. Based on wind erosion, almost 204 km 2 (27%) of the study area are located in the low, 535 km 2 (70%) in moderate, and 21 km 2 (3%) in the high class of desertification ( Figure 12). Moreover, in terms of water erosion indicator, the area classified in very low, low, and moderate is 16.5 km 2 (2.5%), 590 km 2 (77.5%), and 154 km 2 (20%), respectively ( Figure 12). About 18.6% (1414.5 km 2 ) of the total area is located in the low class, regarding erosion quality index, and the moderate and high classes each occupy approximately 40% of the study area. Figure 11. Spatial distribution of climate quality index. Figures 12 and 13 show maps of erosion indicators and erosion quality index. Based on wind erosion, almost 204 km 2 (27%) of the study area are located in the low, 535 km 2 (70%) in moderate, and 21 km 2 (3%) in the high class of desertification ( Figure 12). Moreover, in terms of water erosion indicator, the area classified in very low, low, and moderate is 16.5 km 2 (2.5%), 590 km 2 (77.5%), and 154 km 2 (20%), respectively ( Figure 12). About 18.6% (1414.5 km 2 ) of the total area is located in the low class, regarding erosion quality index, and the moderate and high classes each occupy approximately 40% of the study area.

Environmentally Sensitive Areas to Desertification
Based on the results, desertification sensitivity can be described as presented in Table  10.

Environmentally Sensitive Areas to Desertification
Based on the results, desertification sensitivity can be described as presented in Table 10.
After averaging the indices, the desertification sensitivity map was produced based on MEDALUS results. Figure 14 shows the distribution of ESAs in the study area, and the highly sensitive areas to desertification in Segzi plain are found in the central part, where soil, groundwater, and erosion quality are low. The southwest part of the study area is classified as high sensitivity to desertification because it has better quality of soil, water, and vegetation. The result illustrates that approximately 23.5% (178 km 2 ) of the area is located in high and 76.5% (582 km 2 ) in very high class of desertification. The weighted mean intensity of desertification for the total area of Segzi is 156, which is in the very high class of desertification risk. In other words, the entire study area is highly sensitive to desertification. The results illustrate that climate, vegetation, and groundwater quality are the most important indicators affecting the desertification process. Soil and erosion quality has the least impact on desertification in the area, with a mean of 138 and 140, respectively.  Climate quality index is most important for desertification in this region because of high rates of evapotranspiration, along with low precipitation. Similar results were reported by researchers working on desertification in different parts of the world [10,48,49]. High evapotranspiration rates cause salinization of soils. At the same time, soils are easily erodible [50]. The salinity causes decreased vegetation density that further strengthens desertification. Lack of vegetation results in direct wind impact on soil surface, causing faster-drying soil, high amount of soil erosion, and desertification.
The central parts of the area are characterized by excessive soil salinity and alkalinity (mean EC = 79.84 ds/m) with little vegetation. Due to the fine-grade soil, wind erosion occurs during the dry season and water erosion during the wet season. In the northwest, west, and southwest regions, desertification is less intense than in other regions. In these areas, there is slight wind erosion in most years because of agricultural crops. The existing crop cover and decreased wind speed, which directly controls wind erosion, prevent the soil from drying out and indirectly cause erosion control in the region.
In the Segzi plain, the effect of wind erosion is much greater than water erosion. One of the main reasons for this is that the flat area inhibits low water erosion. Instead, strong winds (average 8 m/s) increase wind erosion intensity.

Conclusions
Desertification is one of the major environmental problems in arid and semiarid areas. In this study, the modified MEDALUS model was used to improve the understanding Climate quality index is most important for desertification in this region because of high rates of evapotranspiration, along with low precipitation. Similar results were reported by researchers working on desertification in different parts of the world [10,48,49]. High evapotranspiration rates cause salinization of soils. At the same time, soils are easily erodible [50]. The salinity causes decreased vegetation density that further strengthens desertification. Lack of vegetation results in direct wind impact on soil surface, causing faster-drying soil, high amount of soil erosion, and desertification.
The central parts of the area are characterized by excessive soil salinity and alkalinity (mean EC = 79.84 ds/m) with little vegetation. Due to the fine-grade soil, wind erosion occurs during the dry season and water erosion during the wet season. In the northwest, west, and southwest regions, desertification is less intense than in other regions. In these areas, there is slight wind erosion in most years because of agricultural crops. The existing crop cover and decreased wind speed, which directly controls wind erosion, prevent the soil from drying out and indirectly cause erosion control in the region.
In the Segzi plain, the effect of wind erosion is much greater than water erosion. One of the main reasons for this is that the flat area inhibits low water erosion. Instead, strong winds (average 8 m/s) increase wind erosion intensity.

Conclusions
Desertification is one of the major environmental problems in arid and semiarid areas. In this study, the modified MEDALUS model was used to improve the understanding of the degradation/desertification process based on a comprehensive set of indices in the Segzi plain. The model is a valuable method to assess the desertification risk in arid and semiarid areas. Five important quality indices including soil, groundwater, vegetation, climate, and erosion were used in the investigation. Each index was weighted and calculated by a combination of sub-indicators. They were then combined to evaluate final desertification sensitivity index.
The results, based on all five quality indices, demonstrated that the study area is highly exposed to desertification processes, especially the central parts of the plain. The results also showed that climate, vegetation, and groundwater are the most important quality indices to describe the desertification processes in the study area. These indices present a high to a very high sensitivity to desertification in 100, 97.9, and 90.9% of the study area. The results are in line with some previous studies [10,27,32,33,36,39,45,51], while several surveys stress that one of the main quality indices to describe desertification processes is the management quality index [3,12,30], which should be considered in future studies. According to the produced desertification map, almost 178 km 2 (23.5%) of the area are classified in the high and 582 km 2 (76.5%) in the very high classes of desertification risk. Hence, the outcomes demonstrate that the study area is very sensitive to land degradation and desertification. Therefore, the developed model can be used to evaluate the quantity, intensity, and severity of different desertification processes in the study area.
The MEDALUS model is a useful method to evaluate desertification that can help decision makers to support sustainable management of the land that is sensitive to degradation processes like the Segzi plain. Climate variation, vegetation cover, and groundwater quality should be continuously monitored to inform sustainable land management strategies. Therefore, the approach used in this investigation can provide information to help decision makers for the region to combat desertification more effectively. The model can also be used as a useful tool in other regions with similar conditions. The use of the MEDALUS approach to study desertification processes was delimited by available data such as groundwater quality and erosion from mainly wind and water. The results are representative for arid areas such as in the studied catchment. Future research would benefit from including more data types in the process study such as management systems and socioeconomic data.  Acknowledgments: The authors are grateful to the United States Geological Survey (USGS) team for providing the Landsat images via https://earthexplorer.usgs.gov/ (accessed on 13 May 2020). We also would like to thank the Centre for Middle Eastern Studies at Lund University and the Ministry of Science, Research and Technology of the Islamic Republic of Iran. The authors also thank two anonymous reviewers and the journal editor for constructive comments on the manuscript.