Soil Water Erosion Vulnerability and Suitability under Different Irrigation Systems Using Parametric Approach and GIS, Ismailia, Egypt

: Preserving the sustainable agriculture concept requires identifying the plant response to the water regime and rationing the water for irrigation. This research compares different irrigation designs coupled with a parametric evaluation system on soil water erosion and soil suitability to assess the sites vulnerable to soil erosion based on a soil water erosion model (ImpelERO) in an area of 150.0 hectares, Ismailia Governorate, Egypt. Land suitability maps are prepared using the Geographic Information System (GIS), and the soil properties are analyzed and evaluated for the different surface, sprinkler, and drip irrigation methods. The results show that the sprinkler and drip irrigation strategies are more practical irrigation methods and additional environment friendly than surface irrigation for enhancing land productivity. Moreover, the principle acumen for creating use of the surface irrigation on this space is for lowering the soil salinity. Land capability index for surface irrigation ranges from 20.5 to 72.2% (permanently not suitable N2 to moderately suitable S2); and the max capability index (Ci) for drip irrigation was 81.3% (highly suitable-S1), while the mean capability index (Ci) was 42.87% (Currently not suitable-NI). The land suitability of the study area using sprinkler irrigation was ranked as highly suitable (S1), moderately suitable (S2), marginally suitable (S3), and currently not suitable (N1). Thus, the obtained data indicated that applying drip irrigation (trickle irrigation) was the most efﬁcient system compared to the sprinkle and surface irrigation systems. To identify the soil, water erosion vulnerability, and soil optimal management strategies for the agricultural parcel in that region, the ImpelERO model (soil erosion vulnerability/impact/management) was applied. Erosion risk classes ranged from V2 (small) to V3 (moderate), that that region categorized as small-sensitive to water erosion by alfalfa, to moderate-sensitive to water erosion by olive. The results of soil losses varied from 7.1 to 37.9 t ha − 1 yr − 1 with an average of 17.7 t ha − 1 yr − 1 . Thus, guarantee efﬁcient water use and soil suitability for food production in the future will require the use of an efﬁcient irrigation system.


Introduction
Water scarcity is the lack of water resources that sufficiently meet the demands of water usage in the Middle East and North Africa (MENA) regions. Due to the rapidly growing population in Egypt, it is ranked among the top ten countries by 2025 with water

Study Area Description
Ismailia is an Egyptian city located in the middle of the Suez Canal territory-the entire area of Ismailia governorate is 5067 km 2 . The study area is bounded by longitudes 32 •

Study Area Description
Ismailia is an Egyptian city located in the middle of the Suez Canal territory-the entire area of Ismailia governorate is 5067 km 2 . The study area is bounded by longitudes 32°8′29″ and 32°9′46″ E, latitudes 30°29′8″ and 30°28′11″ N, covering a surface area of roughly 150 hectares ( Figure 1). The analysis of meteorological data of the recent time coupled with the back-history reflects the weathering status, as well as the soil genesis. The studied region is characterized by a long dry summer with a high temperature and evaporation rate, and a short, mild winter with rare rainfall in autumn. The utmost amount of rainfall was 15.7 mm/month in February. The highest humidity percent is 75% that recorded in January (Meteorological Station of Port Said/El Gamil, 2017-2018, 31°28′ N and 32°23′ E). Ismailia climate is distinguished into three forms (cold winter for an extended time, intermediate with light rains, moderate summer with humidity of 75%), that generally dominant by moderate climate the whole year. Rainfalls do not exceed 50 mm/year, mainly in winter seasons-therefore, the limited rainfall coupled with the rise of evaporation rate results in the drought phenomena cases in that region. The presence of Ismailia canal water and groundwater was used as the main irrigation sources for the Ismailia lands [30,31].
According to Reference [32], the area under investigation is completely covered by deposits of Quaternary sediments that possibly mask earlier tectonic deformations. The Miocene bedrock exposures southward are composed of sandstone, clay, and limestone. The area around the Ismailia Canal consists of two hydrogeological units: The shallow and the deeper aquifers. The salinity of groundwater in the deeper aquifer is low and rarely exceeds 1500 mg/L, while the groundwater salinity of shallow aquifer ranges from 340 to 7650 mg/L [33].
The morphological properties greatly affect the hydraulic conditions of the shallow aquifer around the Ismailia canal. Sand accumulations as dunes bound these pools and separate from the gravelly and sandy slopes in the southwest of the studied area. The natural hydrogeological conditions of the study area are poor, due to the saline water content of the underlying sediments, which are characterized by poor hydraulic properties and moderate to high salt content [32]. The analysis of meteorological data of the recent time coupled with the back-history reflects the weathering status, as well as the soil genesis. The studied region is characterized by a long dry summer with a high temperature and evaporation rate, and a short, mild winter with rare rainfall in autumn. The utmost amount of rainfall was 15.7 mm/month in February. The highest humidity percent is 75% that recorded in January (Meteorological Station of Port Said/El Gamil, 2017-2018, 31 • 28 N and 32 • 23 E). Ismailia climate is distinguished into three forms (cold winter for an extended time, intermediate with light rains, moderate summer with humidity of 75%), that generally dominant by moderate climate the whole year. Rainfalls do not exceed 50 mm/year, mainly in winter seasonstherefore, the limited rainfall coupled with the rise of evaporation rate results in the drought phenomena cases in that region. The presence of Ismailia canal water and groundwater was used as the main irrigation sources for the Ismailia lands [30,31].
According to Reference [32], the area under investigation is completely covered by deposits of Quaternary sediments that possibly mask earlier tectonic deformations. The Miocene bedrock exposures southward are composed of sandstone, clay, and limestone. The area around the Ismailia Canal consists of two hydrogeological units: The shallow and the deeper aquifers. The salinity of groundwater in the deeper aquifer is low and rarely exceeds 1500 mg/L, while the groundwater salinity of shallow aquifer ranges from 340 to 7650 mg/L [33].
The morphological properties greatly affect the hydraulic conditions of the shallow aquifer around the Ismailia canal. Sand accumulations as dunes bound these pools and separate from the gravelly and sandy slopes in the southwest of the studied area. The natural hydrogeological conditions of the study area are poor, due to the saline water content of the underlying sediments, which are characterized by poor hydraulic properties and moderate to high salt content [32].

In-Situ Field Investigations
The soil survey was accomplished at the detailed soil survey level, and the soil profiles were prepared every 250.0 m. Twenty-one soil profiles were chosen to represent the lands under reclamations based on a detailed soil survey ( Figure 2). The soil profiles were wide open to a depth of 150 cm. Soil profiles were expected to reflect the wide variations in soil texture and soil salinity. A morphological description of the soil profiles was undertaken according to the criteria established by FAO, Guidelines for soil description [34]. Soil erosion was described according to Reference [35], based on category and degree. Soil samples, collected from different profile layers, were air-dried and passed through a 2 mm sieve. The rock fragments (>2 mm diameter) were determined volumetrically to measure their sizes and percentages from a total sample, while the fine soil (<2 mm) was subjected to physical and chemical analyses as indicated by Reference [36]. The Soil Particle size distribution was carried out for the sandy soils by the dry sieving method. CaCO 3 contents were determined using Collin's calcimeter. Soil pH, soil salinity, and soluble cations and anions were determined in the soil extract. Sodium and potassium ions were measured photometrically. Calcium, magnesium, chloride, carbonates, and bicarbonates were determined titrimetrically. Cation exchange capacity (CEC) of soil was determined using the sodium and ammonium acetate method, and exchangeable sodium percentage (ESP) was determined using the ammonium acetate method [36]. Organic matter content was determined using the Walkley and Black titration method [36].

In-Situ Field Investigations
The soil survey was accomplished at the detailed soil survey level, and the soil profiles were prepared every 250.0 m. Twenty-one soil profiles were chosen to represent the lands under reclamations based on a detailed soil survey ( Figure 2). The soil profiles were wide open to a depth of 150 cm. Soil profiles were expected to reflect the wide variations in soil texture and soil salinity. A morphological description of the soil profiles was undertaken according to the criteria established by FAO, Guidelines for soil description [34]. Soil erosion was described according to Reference [35], based on category and degree. Soil samples, collected from different profile layers, were air-dried and passed through a 2 mm sieve. The rock fragments (>2 mm diameter) were determined volumetrically to measure their sizes and percentages from a total sample, while the fine soil (<2 mm) was subjected to physical and chemical analyses as indicated by Reference [36]. The Soil Particle size distribution was carried out for the sandy soils by the dry sieving method. CaCO3 contents were determined using Collin's calcimeter. Soil pH, soil salinity, and soluble cations and anions were determined in the soil extract. Sodium and potassium ions were measured photometrically. Calcium, magnesium, chloride, carbonates, and bicarbonates were determined titrimetrically. Cation exchange capacity (CEC) of soil was determined using the sodium and ammonium acetate method, and exchangeable sodium percentage (ESP) was determined using the ammonium acetate method [36]. Organic matter content was determined using the Walkley and Black titration method [36].

Parametric Approach for Soil Suitability
To assess the land appropriateness for various irrigation systems techniques, the parametric assessment framework portrayed by Reference [37] was applied. These qualities are appraised and used to compute the capacity record for capability index for irrigation (Ci) as indicated by the formula: (Ci) =   100 100 100 100 100 where: Capability index for irrigation = Ci; soil texture rating = A; soil depth rating = B; CaCO3 status = C; electro-conductivity rating = D; drainage rating = E; and slope rating = F. These factors and capability index and also the corresponding suitability class ratings are described by Reference [24], and presented in Table 1.

Parametric Approach for Soil Suitability
To assess the land appropriateness for various irrigation systems techniques, the parametric assessment framework portrayed by Reference [37] was applied. These qualities are appraised and used to compute the capacity record for capability index for irrigation (Ci) as indicated by the formula: where: Capability index for irrigation = Ci; soil texture rating = A; soil depth rating = B; CaCO 3 status = C; electro-conductivity rating = D; drainage rating = E; and slope rating = F. These factors and capability index and also the corresponding suitability class ratings are described by Reference [24], and presented in Table 1. The land assessment was carried out according to the topography and soil characteristics of the region. Topographical features consisted of slope and soil properties, such as texture, depth, salinity, drainage, and calcium carbonate. Under the parametric method, the land evaluation is assessed using a numerical index. In this classification system, first of all, a degree, the rate of which is between 0 and 100, is attributed to any property characteristic by comparison with the tables of land requirements.

Assessment of Soil, Water Erosion Vulnerability
Soil erosion continues to be an essential concern for the improvement of sustainable agricultural management systems. The ImpelERO model (soil erosion vulnerability/impact/management) is used to quantify the soil erosion vulnerability and the result of the soil erosion impact on crop production and on the optimal management strategies (https://evenor-tech.com/microleis/microlei/microlei.aspx) Figure 3.  The land assessment was carried out according to the topography and soil characteristics of the region. Topographical features consisted of slope and soil properties, such as texture, depth, salinity, drainage, and calcium carbonate. Under the parametric method, the land evaluation is assessed using a numerical index. In this classification system, first of all, a degree, the rate of which is between 0 and 100, is attributed to any property characteristic by comparison with the tables of land requirements.

Assessment of Soil, Water Erosion Vulnerability
Soil erosion continues to be an essential concern for the improvement of sustainable agricultural management systems. The ImpelERO model (soil erosion vulnerability/impact/management) is used to quantify the soil erosion vulnerability and the result of the soil erosion impact on crop production and on the optimal management strategies (https://evenor-tech.com/microleis/microlei/microlei.aspx) Figure 3.  Land quality will be described by land characteristics associated with soil parameters. A set of management features, which essentially the crop properties and cropping practices Sustainability 2021, 13, 1057 6 of 20 will be considered. This is a simple attribute of agricultural management that can be estimated and can be used as a means of describing management qualities. ImpelERO Model parameters: -Region: That the evaluated region is a set of field-units. Moreover, it can be named evaluating-scenario. -Field-unit: To identify each spatial unit that is analyzed and will be evaluated in the ImpelERO model. The field-unit is composed of climate, soil, and management components [27]. -Soil: The set of soil characteristics that referred to as the soil types which will define the field-unit [27]. -Climate: The set of land characteristics (LC's) referred to as the climate, which defines a field-unit [27]. -Management: Management characteristics (MC's), that mainly consider the crop properties and cultivation practices [27].

Morphological Description and Analytical Results
The study area is characterized as a sandy plain and alluvial deposits. The dominant fruit trees cultivated in this area are (apple, mango, olive, and citrus) besides fodders (alfalfa) and the vegetable crop (tomato). The detailed morphological descriptions of the soil profiles were recorded based on the basis outlined by References [34,35,38]. The topography of the landscape generally ranges between flat (<0.5%) and almost flat (0.5-2.0%) with a nearly-level sloping surface (0.5-1.0%), according to Reference [34]. Very few scattered desert shrubs are found near the new reclamation areas. The common features of this soil type with a depth of 150 cm [35], ranged from sand to sandy clay loam, and clay loam texture with slightly calcium carbonate content and moderately well-drained. Four profiles were selected as a guide from the collected profiles; in the following Tables 2 and 3, an overview of morphological description and analytical data is given. The dominant dry soil color, shown in Table 2, was yellow (10 YR 7/6) and very pale brown (10 YR 7/4). On moistening, these colors become brownish-yellow (10 YR 6/6) and yellowishbrown (10 YR 6/4). Soil structure is generally massive throughout the successive layers of single-grain surface, and dry consistency is losses, soft to slightly hard to increase soil compaction downward from the soil profiles. Being non-sticky and non-plastic, it agrees well with coarse soil texture. With the exception of fine profile texture layers (13), whose wet consistency has been changed to sticky and plastic. Erosion may be defined as water or wind erosion in accordance with Reference [35], and the degree of erosion varies from slight to moderately severe. Table 3 indicates that these soils are non-calcareous to moderately calcareous, where the content of calcium carbonate varies widely from 0.8 to 9.75%, and the surface layers detect the highest carbonate content. As indicated by pH values that range from 7.7 to 9.7, the soil response ranges from slightly alkaline to strongly alkaline. Electrical conductivity values vary from 0.33 to 9.75 dSm −1 based on soil texture, suggesting non-saline to very high saline values.

Weighted Mean of Soil Characteristics
The physical, chemical, and physical-chemical weighted average of the soil properties are determined and shown in Table 3. Based on the formation of the parent material, the results obtained indicate that soil texture varies from area to area (such as profiles 13 and 16) and subsequent layers within each soil profile. A combination of alluvial soils and sand dune soils is the parent content between soil profiles. Soil texture throughout the entire depth is variable textural composition, from sandy clay to coarse sand at the soil surface.
Soil salinity values as represented in Figure 4 ranged widely between 0.33 and 9.70 dSm −1 (non-saline to very strongly saline), with an average value of 2.9 dSm −1 (moderately saline) according to Reference [35]. The relatively low EC values in these soils with coarse sand textured and low content of fine fractions of silt and clay are mainly affected by the particle size distribution of texture. The lowest values of EC exist in the top surface layer, which is dominated by coarse texture. The electrical conductivity weighted mean of the soil saturation extract varies from 0.90 to 8.74.0 dSm −1 , indicating non-saline to strongly saline condition; the lowest values are mostly detected in the coarse texture profiles. The distribution of the weighted mean of soil salinity (EC) values showed in Table 3. Salinity is one of the main edaphic factors which limits the distribution of plant communities in their natural habitats and which is causing increasingly severe agricultural problems [14,15]. particle size distribution of texture. The lowest values of EC exist in the top surface layer, which is dominated by coarse texture. The electrical conductivity weighted mean of the soil saturation extract varies from 0.90 to 8.74.0 dSm −1 , indicating non-saline to strongly saline condition; the lowest values are mostly detected in the coarse texture profiles. The distribution of the weighted mean of soil salinity (EC) values showed in Table 3. Salinity is one of the main edaphic factors which limits the distribution of plant communities in their natural habitats and which is causing increasingly severe agricultural problems [14,15]. Soil reaction pH is an important parameter, due to its effects on nutrient availability, microbial activity, and plant growth. From the data presented in Table 3, the weighted mean of soil reaction varies considerably from 7.94 to 9.63, indicating moderately alkaline to strongly alkaline soil reaction [38], as shown in Figure 5. Soil profiles under study are slightly calcareous to moderately calcareous according to Reference [35], where the CaCO3 content ranges from 0.9% to 9.75% with an average value of 4.4%. In addition, weighted mean values of CaCO3 content in soil profiles varied from 1.72% to 7.94%, as presented in Figure 6. The highest CaCO3 contents are mostly detected in the surface layers and increased with fine soil texture. Soil reaction pH is an important parameter, due to its effects on nutrient availability, microbial activity, and plant growth. From the data presented in Table 3, the weighted mean of soil reaction varies considerably from 7.94 to 9.63, indicating moderately alkaline to strongly alkaline soil reaction [38], as shown in Figure 5. Soil profiles under study are slightly calcareous to moderately calcareous according to Reference [35], where the CaCO 3 content ranges from 0.9% to 9.75% with an average value of 4.4%. In addition, weighted mean values of CaCO 3 content in soil profiles varied from 1.72% to 7.94%, as presented in Figure 6. The highest CaCO 3 contents are mostly detected in the surface layers and increased with fine soil texture.    Concerning the levels of organic matter values in the upper soil layers, Table 4 shows that organic matter content varied from very low to low, and ranges from 0.27 to 0.82%. The bulk density values vary from 1.52 g/cm 3 to 1.59 g/cm 3 -these contents were coinciding very well with texture classes. The cation exchange capacity is used as a parameter for the fertilizers buffering capacity. Cation exchange capacity (CEC) of soils is closely related to the content of clay content. The exchange characteristics of the soils under study dictate that CEC weighted mean values are low and coincide well with soil texture [39], being in the Concerning the levels of organic matter values in the upper soil layers, Table 4 shows that organic matter content varied from very low to low, and ranges from 0.27 to 0.82%. The bulk density values vary from 1.52 g/cm 3 to 1.59 g/cm 3 -these contents were coinciding very well with texture classes. The cation exchange capacity is used as a parameter for the fertilizers buffering capacity. Cation exchange capacity (CEC) of soils is closely related to the content of clay content. The exchange characteristics of the soils under study dictate that CEC weighted mean values are low and coincide well with soil texture [39], being in the range from 3.49 to 27.16 meq/100 g soil with an average of 15.3 meq/100 g soil and coincide very well with texture classes (Figure 7). In addition, the essential substantial water high-quality popular on crop productiveness is that the water salinity hazard. That the quantity of water transpired via a crop was once immediately associated to yield; therefore, irrigation water with excessive EC reduces the yield potential. Groundwater salinity (ECw) ranges from 2500 to 5500 µmhos (Figure 8). According to permissible limits for classes of irrigation water [40], the study classification for EC indicated that the water quality varies from Class 4 (doubtful) to Class 5 (unsuitable). For groundwater, approximately 40% may be used with no troubles, and 40% may be used with occur moderate problems.  Desertification is also termed soil erosion of dry land. Human-made issues comprise the plurality of ways of soil degradation. While some physical parameters are involved, an important aspect is misuse. Poor land management and irrigation, including used surface irrigation, increase the pace of land degradation [41].    Desertification is also termed soil erosion of dry land. Human-made issues comprise the plurality of ways of soil degradation. While some physical parameters are involved, an important aspect is misuse. Poor land management and irrigation, including used surface irrigation, increase the pace of land degradation [41].

Parametric System and Land Suitability for Different Irrigation Methods
The applied irrigation methods; surface (gravity), sprinkler, and drip (trickle) irrigation, were classified as highly suitable (S1), moderately suitable (S2), marginally (S3), currently not suitable (N1), and permanently not suitable (N2), according to [24], and the results of the land suitability in the study area under these systems presented in Table 5 and Figure 9. Comparing the three irrigation systems revealed that using drip irrigation (trickle irrigation) was more efficient. The land capability index for surface irrigation ranges from 20.5 to 72.2% (N2 to S2), as shown in Figure 9. Soils that are permanently not suitable (N2) for surface irrigation cover 115 hectares (78.2%), whereas the moderately suitable soils (S2) cover 12 hectares (8.1%), as displayed in (Table 6). The results showed that the land suitability of sprinkler irrigations is highly suitable (S1), moderately suitable (S2), marginally suitable (S3), and currently not suitable (N1), as shown in (Figure 10). By analyzing the land suitability maps for sprinkler irrigation, it is evident that a highly suitable area can only be observed in some parts (0.2 ha) Table 6. Furthermore, Sprinkler irrigation revealed the second priority investigated to cultivation by field crops, fodder, and some vegetables. The max capability index (Ci) for drip irrigation was 81.3% (highly suitable-S1), while the mean capability index (Ci) was 42.87% (currently not suitable-NI).  Figure 9. Land suitability classes for surface irrigation. Figure 9. Land suitability classes for surface irrigation. Highly suitable (S1); moderately suitable (S2); marginally (S3); currently not suitable (N1), and permanently not suitable (N2).   Data represented in Figures 9-11 and Table 5, indicated that applying a drip irrigation system was the most efficient option compared to sprinkle and surface irrigation systems. Drip irrigation suitability gave more irrigable areas compared to the surface irrigation practice because of the topographic (slope), soil (depth and texture) consistent with main data [42]. Thus, our results are convenient with [43,44] which investigated the land suitability applying the parametric methods for the surface and drip irrigation in Southern Ankara, Turkey. In addition, the obtained results showed that 51.2% of the studied area were highly suitable for drip irrigation method, whereas 13.1% was highly suitable for surface irrigation methods. Moreover, the results revealed that the main limitation factors were soil texture and soil depth for both irrigation systems. These results are in agreement with Reference [18], which investigated the land suitability using the parametric for surface, sprinkler, and drip irrigation in North West of Egypt. The results obtained showed that the study of various types of irrigation techniques showed that sprinkler and drip irrigation systems were more efficient and effective than surface irrigation for increasing land productivity. Depth of soil, soil texture, CaCO 3 , drainage, and slope were the key limiting factors for the use of surface irrigation techniques in Siwa, Western Northern Egypt [19].
with Reference [18], which investigated the land suitability using the parametric for surface, sprinkler, and drip irrigation in North West of Egypt. The results obtained showed that the study of various types of irrigation techniques showed that sprinkler and drip irrigation systems were more efficient and effective than surface irrigation for increasing land productivity. Depth of soil, soil texture, CaCO3, drainage, and slope were the key limiting factors for the use of surface irrigation techniques in Siwa, Western Northern Egypt [19]. Figure 11. Land suitability classes for drip irrigation. Figure 11. Land suitability classes for drip irrigation.

Assessment of Soil Erosion Vulnerability and Potential Impact on Crop Productivity
Sustainability modeling analysis for soil erosion should include not only vulnerability prediction, but also address impacts and response assessment in an integrated way. In the present study, an integrated neural network-based model, ImpelERO, was applied to identify the vulnerability index, erosion risk class, and soil loss rates. Our results revealed that the soil vulnerability indices ranged from 0.21 to 0.34, and 0.44 for sandy clay, sandy clay loam, and sand texture, respectively. Moreover, the main data showed that the elevation ranges from 20.0 to 64.0 m above sea level ( Figure 12)-these results indicated that the elevation and soil texture are the main limiting factors for soil loss rates. Data showed in Table 7 indicated that the values of erosion risk classes range from V2 (small) to V3 (moderate), as shown in Figure 13, which categorize the region as small-sensitive by alfalfa, to moderate-sensitive to erosion by olive. The values of soil losses varied between 7.1 to 37.9 t ha −1 yr −1 with an average of 17.7 t ha −1 yr −1 .
The overall approach of ImpelERO was applied for the three different soil textures selected to quantify the soil erosion vulnerability with several crops, the impact of soil erosion on crop production, and the optimal management strategies. The potential impact of soil erosion on crop productivity was interesting that the projection in a productivity reduction in 2050 will be 6.8, 2.3, and 0.9% for fine sand, sandy clay loam, and sandy clay, respectively. Thus, farming practices can be widely modified to protect environmental qualities.
clay loam, and sand texture, respectively. Moreover, the main data showed that the elevation ranges from 20.0 to 64.0 m above sea level ( Figure 12)-these results indicated that the elevation and soil texture are the main limiting factors for soil loss rates. Data showed in Table 7 indicated that the values of erosion risk classes range from V2 (small) to V3 (moderate), as shown in Figure 13, which categorize the region as small-sensitive by alfalfa, to moderate-sensitive to erosion by olive. The values of soil losses varied between 7.1 to 37.9 t ha −1 yr −1 with an average of 17.7 t ha −1 yr −1 .    Figure 13. Soil water erosion vulnerability using ImpelERO model at the study area.
The overall approach of ImpelERO was applied for the three different soil textures selected to quantify the soil erosion vulnerability with several crops, the impact of soil erosion on crop production, and the optimal management strategies. The potential impact of soil erosion on crop productivity was interesting that the projection in a productivity reduction in 2050 will be 6.8, 2.3, and 0.9% for fine sand, sandy clay loam, and sandy clay, Figure 13. Soil water erosion vulnerability using ImpelERO model at the study area.

Conclusions
Evaluating the land suitability for different irrigation methods using the parametric evaluation systems is an essential strategy for the decision-makers from the viewpoint of water-saving. We found that the land suitability of 78.2% (109 hectares) will be improved through applying sprinkler and drip irrigations compared to the surface irrigation method. Nevertheless, the main limiting factors in choosing the relevant irrigation methods in this region were the soil salinity. Thus, the areas belong to highly salinity content in the soil surface layers and high CaCO 3 content; it is more relevant to use the surface irrigation method. That could help in leaching salts out of the soil profile. The exercise proved that soil characteristics, land capability, and suitability for crop systems, combined with tools of land evaluation and GIS maps, is a powerful tool that will be used as a decision support system.
Thus, the recommendation of the research revealed that it is more efficient to use the drip irrigation method to irrigate the soil as it ensures the sustainable use of the land for agriculture. Whereas, sprinkler irrigation is considered the second priority in the study area, the field crops, fodder, and some vegetables can be grown in priority lands in sprinkler irrigation.
In addition, for the soil water erosion, the sustainability modeling analysis should not consider only the vulnerability prediction, but also the response assessment. The integrated neural network-based model, ImpelERO, was successfully used to identify the soil water erosion vulnerability index, and soil loss rates. The soil water erosion potential impact on the crop productivity was interesting that the projection in productivity reduction in 2050 will be 6.8, 2.3, and 0.9% reductions for fine sand, sandy clay loam, and sandy clay, respectively. Thus, the relevant farming and irrigation practices can be widely modified to protect future environmental qualities.