Detecting and Mapping Salt-Affected Soil with Arid Integrated Indices in Feature Space Using Multi-Temporal Landsat Imagery

: Salinity systems are well known as extreme environmental systems that occur either naturally or by certain human activities, in arid and semiarid regions, which may harm crop production. Soil salinity identification is essential for soil management and reclamation projects. Information derived from space data acquisition systems (e.g., Landsat, ASTER) is considered as one of the most rapid techniques in mapping Salt-Affected Soil (SAfSoil). The current study tested the previously proposed salinity indices on the northern Nile Delta region, Egypt. The results indicated that most of the indices were not suitable to detect the SAfSoil in the area, due to the interaction between the bare soils, salts, and urbanization. To resolve this issue, the current work suggested a new index for detecting and monitoring the SAfSoil in the Nile Delta region. The newly proposed index takes into consideration plant health, the salt crust at the surface of the soils, as well as urbanization. It facilitates the mapping processes of SAfSoil in the area compared to any other previously proposed index. In this respect, multi-temporal Landsat-7 and 8 satellite data, acquired in 2002, 2016, and 2021, were used. The new index was prepared using the 2002 data and verified using the 2016 and 2021 data. Field measurements and data collected during 2002, 2016, and 2021 were utilized as ground truth data to assess the accuracy of the results obtained from the proposed index. The evaluation of the results indicated that the accuracy assessment for 2002, 2016, and 2021 images was 94.58, 96.08, and 95.68%, respectively. Finally, the effectiveness of using remote sensing in detecting and mapping SAfSoil is outlined.


Introduction
Soil salinity is a primary land degradation process that reduces biological production in arid and semi-arid climates, requiring comprehensive monitoring and management.Recent advances in remote sensing technology allowed soil salinity to be recognized and monitored reliably [1][2][3].Soil salinity is known as a dynamic process, which arises as a result of natural processes and human activities, particularly in arid and semi-arid environments as a result of insufficient rainfall, an excessive evaporation of shallow groundwater, and a highly soluble concentration of salt [3,4].Salinity alters the chemical and physical properties of soils, which results in a considerable reduction in crop productivity [5][6][7].
Saline soils are defined by the USA Salinity Laboratory Staff [8] as a soil with an electrical conductivity (EC) value > 4 (dS m −1 ), pH < 8.5, and exchangeable sodium percentage (ESP) < 15.An increase in soil salinity affects both soil and water quality [9][10][11][12], which, in turn, affects plant growth, production, and may eventually leads to soil deterioration [13][14][15][16][17][18].Soil salinity is usually evaluated by collecting soil samples and measuring electrical conductivity in the laboratory.This method is well-known to be time-consuming and expensive due to the number of samples required to understand the spatial variability in the soil salinity for a given area [19][20][21][22][23][24].The high geographic variability of soil salinity over short distances necessitates the collection of multiple soil samples, which results in such an expensive and time-consuming process.Remote sensing technologies, when compared to traditional approaches, provide considerable cost-effective advantages in detecting soil salinity for vast regions with high precision [1,21,25].Satellite images, along with advanced processing techniques, provide a better tool in detecting and mapping soil salinity due to the spectral demeanor of the salt lineaments on the Earth's surface, which is a reasonable indicator for assigning salinity [24,[26][27][28][29][30].It is important to keep in mind that many factors influence soil spectral reflectance, e.g., the soil physicochemical properties (i.e., wetness or moisture-content, organic matter (OM), texture, clay content, and/or types and degree of surface roughness) [30][31][32][33][34], the salt crust being mixed with other soil constituents, or being invisible beneath them [35,36].Thus, it is sometimes challenging to predict soil salinity from reflectance [37].Consequently, several studies used vegetation cover (plant health) as an indicator for soil salinity to avoid soil reflectance complications [22].Scattered (sparse) vegetation on the soil surface, which can function as a sign of high salinity, have also been used as an essential indicator to identify pieces of land affected by salinity using flower reflectance [38].Additionally, the photosynthetic characteristics of unhealthy vegetation can boost the reflectance of visible light and reduce reflectance from the Near Infra-Red (NIR) range [39].This behavior has been recognized in various plants, which undergo salinity fatigue [29].Several Vegetation Indices (VIs), such as the Normalized Difference Vegetation Index (NDVI) and Soil-Adjusted Vegetation Index (SAVI), were also exploited as main indicators for estimating soil salinity [40][41][42].Moreover, numerous scientists have developed and used other sophisticated salinity indices along with the Normalized Difference Salinity Index (NDSI) and Salinity Index (SI) [36].
Several authors used remote sensing data analysis to investigate the soil salinity at different locations by investigating the spectral properties of the saline soils in visible (VIS) and infrared (IR) wavelengths [21,[37][38][39].Their primary method was mainly to create regression models between ground sample locations and remote sensing datasets.Three approaches have mainly been used to map the soil salinity: (1) model development between reflection values of spectral bands in satellite images and measured soil salinity values [37,40], (2) directly estimating soil salinity by calculating salinity indices [21,39,40], and (3) indirectly estimating soil salinity by calculating vegetation indices [21,41,42].
Due to the unique reflection of electromagnetic energy from salt crusts on the soil surface, soil salinity is detected and monitored at numerous scales using remote sensing technologies [4,39].However, some factors such as soil texture, moisture, organic matter, and iron oxide content [43], as well as vegetation type and development period [44], interfere with the energy reflected or adsorbed from saline soils under natural conditions, reducing the success of remote sensing in determining surface salinity.Furthermore, the resolution of remote sensing images has a substantial impact on the efficacy of remote sensing-based soil salinity determination [45][46][47].Nonetheless, various authors claimed that employing remote sensing techniques to determine soil salinity has several advantages [1,4,39].This viewpoint is widely acknowledged as a low-cost soil salinity monitoring tool, particularly in basins with insufficient drainage [48].However, developments in remote sensing technology have yet to produce an optimum data type and interrogation method combination for assessing soil salinity under a variety of ambient situations [49,50].
Multi and hyperspectral images have been previously used in soil salinity mapping, but the results obtained from these studies were mixed [36].Consequently, researchers interested in the prediction and delineating of soil salinity usually accompanied their research with the indices of flora, i.e., NDVI, SAVI, Radar Vegetation Index (RVI), Salt Balance Index (SBI), Green Vegetation Index (GVI), and Humidity Index (WI) [41,45,46,57,[60][61][62]. In the Landsat image, a significant relationship was revealed between the reflectance in the blue, green, red, and near-infrared bands, and the EC values, as well as brightness (BI) and Humidity (WI) indices [63][64][65].
Kafr El-Sheikh Governorate is located in the far north of Egypt.The entire area is located in the north of the Nile Delta and overlooks the Mediterranean Sea with a coastal extension of 100 km in length and is bordered to the west by the Rashid branch with a length of 85 km until its mouth in the Mediterranean Sea.The land cover/land use of the region varies between plants, sabkhas, urban areas, and water bodies [66].This diversity and its intensity affect the ability of remote sensing data to identify saline soils.With reference to the area landscape and its geology structure, the soils consist of alluvial plain (flood plain), coastal plain (marine plain), and lacustrine deposits (mixed of flood plain and marine plain) [17].The saline soils refer to lands with inherited salts, while the salt-affected soils refer to acquired salts formed by nature, climate, and human activity.The current research proposes a Soil Salinity Index for arid and semi-arid conditions (SIA), which is tested along the northern Nile Delta region, Egypt, where different levels of soil salinities have been recorded.Landsat Thematic Mapper (TM) imagery is used in the development process of the SSI, which relies on: (1) NDVI, to extract flora fraction, WI, to extract the waterlogging, and the alkalinity index to extract the alkalinity; and (2) examination of the relationships between the direct and indirect indicators in relation to their variation in the saline soil.

Materials and Methods
Study area: Kafr El-Sheikh Governorate is one of the north Nile Delta governorates of Egypt.It is located between longitudes 30°20′5″ and 31°20′10″ East, and latitudes 30˚50′35″ and 31°55′32″ North, with a total area of 377,798 ha (Figure 1).It is characterized by imperfect to well-drained soil as it was formed within the alluvium of the Nile delta, which has been deposited in semiarid climates [66].The hydrology of the area is affected by both the River Nile irrigation canals, and the seepage from the Mediterranean Seawater, but overall, the groundwater quality of the area is mostly good [66][67][68].
The area slopes gently to very gently toward the north.It has a semi-arid climate and receives a total annual rainfall of about 72.44 mm, which spreads mostly along 6 months/year.The annual mean and total evaporations are 89.76 and 1077.08 mm, respectively.The maximum and minimum mean annual temperatures are about 32.3 °C and 18.7 °C, respectively.The annual mean relative humidity is about 56.78%.Additionally, the annual mean wind speed is about 3.9 km h −1 [17].
The main economic activity is agriculture and fishing.The governorate is considered an agricultural area.It is famous for the cultivation of rice, beets, cotton, and wheat.It is distinguished by its industrial activity, as it manufactures dairy, oils, soap, fodder, rice, gin, and beet sugar, in addition to the spread of many fish farms [66].The natural environments in the region can be classified into three main types: agricultural and urban environments, coastal environments, and wetlands.The center and south of the province cover the sediments of the modern geological age (the Holocene era), which are dark brown formations composed of deposits of clay and sandy clay.These sediments are deposited over the ancient marine sediments (under-delta formations) that date back to the Pleistocene era.They are yellow in color and consist of coarse and fine sand and pebbles consisting of quartz or igneous and metamorphic rocks.As for the northern coastal zone, it is a low sandy coast consisting of soft brittle sediments belonging to the Pleistocene and Holocene periods.The 2016 and 2021 soil samples were air-dried, softly crushed, sieved through a 2 mm sieve, and stored for analysis.Soil parameters were determined in triplicate in accordance with SSSA-Methods of Soil Analysis: Part 3 Chemical Methods, 5.3 (1996) https://acsess.onlinelibrary.wiley.com/doi/book/10.2136/sssabookser5.3(accessed on 25 May 2022).A Perkin Elmer PinAAcle 500 Flame Atomic Absorption Spectrometer was used.EC was determined using soil saturated water extract according to the methods of USDA [70].
In order to quantitatively measure land surface salinity, a set of multi-temporal Landsat images, digital topographic (scale of 1: 25,000), geologic, and hydrologic maps were exploited to obtain an accurate soil classification map.The multi-temporal Landsat-7 and 8 images were acquired in 2002, 2016, and 2021, respectively, with both atmospheric and radiometric corrections applied.Periodic noise was automatically removed using the Fast Fourier transform (FFT) method.The satellite overpasses as well as the path/row for the Landsat images were, respectively, as follows: (LE07_L1TP_177038_20020617_20200916_02_T1, LC08_L1TP_177038_20160818_20200906_02_T1, and LC08_L1TP_177038_20210731_20210804_01_T1) (Table 2).
Table 2. List of the satellites used, the types of sensors captured, and their bandwidths.The imagery was taken from https://www.usgs.gov(accessed on 30 June 2021).The imagery was processed and analyzed using ENVI 5.1.In the initial analysis of the current study, the indices of Khan et al. [71], Farifteh [72], Setia et al. [73], and Dehni and Lounis [74] (Table 3) were selected, respectively, to check their ability in mapping the SAfSoil in the area.A mask was applied to the vicinity of Kafr El-Sheikh Governorate (Figure 1), as well as the Borollus Lake, to maximize the delineation accuracy of SSI, as well as enhance the extraction process of urban areas.Initial results were acquired using the integration between the Normalized difference built-up index (NDBI), supervised classification, and visual interpretation where urban areas were subtracted from SAfSoil.The soil salinity layer was then extracted using the previously published indices (Table 3).Finally, SAfSoil distribution was verified against the map produced from the geostatistical techniques of the collected soil samples.All the previously extracted layers were then combined to produce the final land-cover map of the area (Figure 3).

Salinity Indices Salinity Indices References
Normalized Differential Salinity Index

R NIR NDSI R NIR
Dehni and Lounis, 2012 [74] Thematic maps for the soil salinity were developed using ArcGIS 10.5 software.Thematic maps of the study area were generated using Inverse Distance Weighted (IDW) interpolation as recommended by AbdelRahman et al. [75] and Yao [76].IDW interpola-tion determines cell values using a linearly weighted combination of a set of sample points.
The following equation was used for IDW interpolation methods as well-discussed by Yao et al. [58]  * ( ) =  ( ) where the Z(xi) data value of locations is used to generate the variable Z value of x0 at the unsampled location, the Z(xi) value is assigned by the weight wi, and n is the number of closest neighboring data points used for estimation.
where di is the distance between the estimated point and the observed point.
where xi and xi + h are sampling locations separated by a distance h, and Z(xi) and Z(xi + h) are the observed values of variable Z at the corresponding locations.
The least squares method that was used to estimate the linear regression is the following equation: After AbdelRahman [66], the research area was divided into three categories (flood plain, lacustrine deposits, and coastal plain).Geomorphology and geology are depicted in Figure 3: 1. Nile silt covers the most land, consisting of fine-grained sediment (silt and clay) deposited on a flood plain by floodwaters that cannot be confined inside the stream channel.2. Quaternary marine deposits (such as carbon deposits and gypsum metal) are found on lower sections of the existing land surface, and terrestrial deposits are found on higher parts of the present sea floor.3. Sabkha deposits, which are flat and extremely saline patches of sand or silt lying slightly above the water table and typically having soft nodules and enterolithic veins of gypsum or anhydrite, are salt-flat soils.4. Sand dunes are wind-created sand ridges that can be seen in deserts or near lakes and oceans.5. Wind-stabilized dunes migrate inward from the lake.Movement is halted if the sand dune is stabilized by vegetation, such as this dune slope.
Results obtained from these indices (Table 3) were carefully examined as they were not efficient in detecting the SAfSoil in the area, as they were either inconclusive or there was a major overlap between the SAfSoil and other classes in the area (e.g., urban areas, waterlogging).
According to Richard's criteria, the major predominant soil types in the area are Nonsaline-nonalkaline soil with area is about 65%, whereas saline soil is about 7%, while alkaline soil is about 4% whilst saline-alkaline soil is about 5% of the total area.Soils are alluvial due to silt deposits (clayey) resulting from the flow of the Nile River.In this region, diverse soil types are observed, which also vary in their production potential.In general, these soils have been described as a complex of salinity ranges in the topographic units and overlapped with the results of the approached index.On the basis of texture, they are categorized into two ranges: light (coarse textured) soils in the north of Borollus lake along the coastal shoreline and heavy soils (clayey black types) in the rest of the area.Therefore, the laboratories soil analyses help in the delineation of saline, sodic, and Saline-Alkaline (sodic) Soils.Accordingly, it is possible to identify nonsaline or alkaline areas.These are the boundaries that were compared with the areas obtained from the new proposal, and the accuracy of the results was measured, which are shown in the tables (5, 6 and 7) discussed later.
During field work, stressed vegetation was used as an indirect sign for the presence of salt in the soils, where salt-affected soils are usually characterized by poorly vegetated areas.The selected remote sensing indices were used to discriminate and map salt-affected soils.Salinity index was the best combined ratio of selected bands to compute the targeted areas.
Detailed examination of Landsat data shows that the maximum and minimum observed wavelength values of SAfSoil ranged between 10.5 and 12.5, and 0.76 and 0.90 μm, respectively, while the urban areas maximum and minimum wavelength were well represented in the ranges between 2.08 and 2.35, and 10.5 and 12.5 μm, respectively.These observations were used to produce much-suitable SAfSoil indices for the study area.Additionally, the study area is mostly covered with vegetation, water lakes, and waterlogging locations, which means involving other indices for vegetation and water is essential for the good separation between the different classes in the area.The current proposed new indices take into account the enhancement for both vegetation and water following the indices of Deering and Rouse [77], and Mcfeeters [78], respectively (Table 4).NDWI and NDVI were used to distinguish vegetation from urban areas.After that, both features were used as a comparison with the new index, taking into account the improvements achieved compared to the old indices, thus ensuring that the results obtained are not affected by unfavorable categories.

Data Processing
Before developing the proposed soil salinity index SIA, the soil samples collected were analyzed to make an interpolation as representative as possible of the salinization process.The analysis of the variability and spatial structure of these data was performed under Arc GIS (Spatial Analyst).This module, which integrates the IDW (inverse distance weighting) interpolation method, offers the possibility to analyze the spatial variability in the data.This algorithm was used to interpolate the 2002, 2016, and 2021 soil data, including sodium adsorption ratio (SAR) and soil electrical conductivity (EC).Several studies have recommended this interpolation procedure as it allows reliable results according to AbdelRahman [68].Landsat data were used to develop the land cover map following the approach mentioned in AbdelRahman [75].Finally, the factorial thematic maps produced were overlaid in a GIS environment to evaluate the SIA map.

Atmospheric Correction
The initial image processing investigation was to perform atmospheric correction for the Landsat data.The atmospheric correction algorithm was effectively and successfully used in soil salinity studies, particularly in arid and semi-arid environments.Once the atmospheric correction process was finalized, the Landsat digital number values (DN) were converted into their corresponding ground surface reflectance values.Using the field-measured GPS, the location of each soil sampling location was identified in the satellite images and their corresponding spectral reflectance values were extracted.The field GPS receiver had a locational accuracy of + or − a meter.It is worth noting here that only the spectral reflectance values of the selected Landsat bands were used in this study.

Index Development
To identify and observe how the saline soil's spectral reflectance differs from nonsaline soils, it was necessary to plot the spectral reflectance of various levels of soil salinities.This is believed to aid in the process of developing a soil salinity index particularly suited for the study area.
All soils, with their varying levels of salinities, peaked in their reflectance at the short-wave infrared reflectance band.Saline soils had higher reflectance values within the visible and blue*red/green parts of the spectrum in comparison to nonsaline soils.This was likely due to the high intensity of reflectance from such surfaces, particularly due to the whitish color of salt covering a saline pixel within the study area.The saline soil's reflectance at this band was extremely low in comparison to nonsaline soils, suggesting strong absorption.The close-up of saline and nonsaline soils' reflectance in the Red*NIR/Green regions and demonstrated this further.
While both saline and nonsaline soils all had varying degrees of response in the Blue − Red/Blue + Red regions, it appeared that a key difference between the two was within the Blue − Red/Blue + Red regions.One key noticeable difference in the spectral reflectance between saline and nonsaline soils was within the three parts' combinations of all bands.For this, the study developed a SIA, which took into account this difference in spectral response.
It should be expected that this index would be significantly higher than unity in its magnitude for saline soils and quite low in magnitude for nonsaline soils.High-salinity soils will portray higher values for the term equation part1-part3 and lower values for part 2, leading to index values higher than unity.To assess the accuracy of this index in estimating soil salinity, it is imperative to compare it with some of the widely used salinity indices such as those mentioned in Table 3.
The proposed index relies mostly on the visible and NIR portion of the electromagnetic spectrum, as illustrated below: The proposed index used three parts of the equation, where the first part used blue, red, and green bands, as illustrated in Figures 4-6.This part of the equation illustrates the highly saline soils (>4 dS m −1 ) in an accurate way, but, unfortunately, there was an overlap between highly saline soils in red (at the sand beach of the shoreline) and dense urban areas, while it succeeded in extracting the saline sodic soils, presented in light green along with the waterlogged areas.Furthermore, the model succeeded in recognizing the saline soils with high accuracy, but there was an overlap with waterlogged, dry sabkhas, and low-density urban areas (in yellow).The model observed the nonsaline and nonsodic areas in high accuracy (dark green).
For the second part of the equation, as shown in Figure 5, the model succeeded in recognizing the waterlogged, fish farms, and dry sabkhas in an accurate way in red.In addition, for the highly saline and highly sodic areas that are presented in orange, the model succeeded in configuring slightly to moderately saline, moderately sodic soils (yellow color), while the light green color was nonsaline-nonsodic soils, and the dark green color indicated the saline highly sodic soils, but it overlapped with dense urban areas.For the third part of the model, as seen in Figure 6, the model succeeded in recognizing the highly saline highly sodic soils but with an overlap with the dense urban areas (red color).While the orange color illustrated slightly saline highly sodic soils, the yellow color illustrated the saline sodic soils with sabkhas.The light green color illustrated the saline sodic soils mixed with fishponds, and the dark green color demonstrated nonsaline-nonsodic soils.

Empirical Models Development
Statistical and mathematical models were developed to model vegetation canopies, soil moisture, and soil salinity.Research in this area has resulted in considerable improvements in estimating soil salinity, its spatial distribution, and the interaction between electromagnetic radiation and saline soils.In an attempt to develop a new empirical model, the SIA values for the field soil sampling locations were regressed with the lab-estimated salinity EC measurements.Various regression models were developed, including linear and logarithmic; however, the exponential relation portrayed the highest R-square values.The exponential relation is ideal with soil salinity indices due to their low saturation levels.The developed semi-empirical models were later applied for the Landsat for generating soil salinity and soil salinity class maps.

Results and Discussion
Results obtained from the indices [71][72][73][74] (Table 3) were not successful in mapping the SAfSoil in the study area (Figures 5 and 7-9).The indices of Khan et al. [71] and Farifteh [72] (Table 3) failed to separate between the different types of soil degradation (e.g., waterlogging, saline soil, and SAfSoil).It also failed to separate between the SAfSoil and the soil sealing.Its ability for water bodies extraction was also limited (Figures 8 and  9).There was also a considerable overlap between the SAfSoil and the urban areas.On the other hand, the indices of Setia et al. [73] and Dehni and Lounis [74] (Table 3) were inconclusive where the different classes were mostly indistinguishable (Figures 10 and  11).Results from our proposed indices to an image acquired in 2002 (Figure 12) clearly separated the SAfSoil in the study area.It was able to distinguish between the saline and nonsaline soil, as well as urban areas and waterlogging locations.These results were verified using data collected during fieldwork conducted during the same year of 2002.While applying the proposed index, five classes were extracted, namely: waterlogging, saline-alkaline soil, nonsaline nonalkaline soil, urban areas, and saline soils.The classification accuracy assessment was 97.8% with a kappa coefficient of 0.88 (Table 5).To further assess the accuracy of the results, the analysis was performed on satellite data acquired during 2016 and 2021 with a time span of 10 years between the two datasets, so we can better see the changes.Field visits to some key sites were carried out to collect sufficient field data for the validation of the results.The results were also compared with those obtained from 2002 to see if the changes were consistent.The classified images with the same five classes are listed in Figures 13 and 14 for the years 2016 and 2021, respectively, with the tables showing the results from their accuracy assessment listed in Tables 6 and 7, respectively.8).The waterlogging class represented 7.5%, 8.4%, and 9.9% of the total area for the years 2002, 2016, and 2021, respectively, showing a clear increase over time.This can be attributed to the intensive agricultural practices, as well as the misuse of irrigation water in the area.The saline soils classes represented 4.1%, 5.2%, and 7% of the total area for the years 2002, 2016, and 2021, respectively, showing a clear increase over time as well.This can be clearly observed by the large patch of salt-affected area located along the shoreline of the Mediterranean Sea and at the southern boundary of Borollus Lake.The saline-alkaline soils represented 3.8%, 2.8%, and 0.7% of the total area for the years 2002, 2016, and 2021, respectively, while the nonsaline nonalkaline class represented 67.3%, 65.6%, and 64.5% of the total area for the years 2002, 2016, and 2021, respectively, showing a clear decrease over time.Lastly, the urban areas represented 4.1%, 5.3%, and 5.9% of the total area for the years 2002, 2016, and 2021, respectively.According to Madani [80], examining the Landsat images, it was found that saline soils had a high spectral reflectance in the red and NIR bands and seemed to have a bright image signature.Salinity prediction models were created using these factors to infer soil salinity across a vast region.In contrast, Mehrjardi [81] discovered that band 3 (red band) had the strongest association with EC among the Landsat ETM+ bands 1-5 and 7, and that a regression model built to tie EC to band 3 and the exponential relation was determined to be the best kind of model based on that finding.Salinity and sodicity may arise for a variety of causes, including geological formation or human and climate change effect.For the observation of soil salinity and sodicity model, ground water level, landform, land management type, soil texture, and land cover were employed.The utilization of Landsat ETM+ data of the research region indicated bare patches with a high reflectance owing to high salinity and/or salt-efflorescence on the soil surface, according to our findings, also in agreement with the finding of Abbas and Khan [82].Elnaggar and Noller [65] observed that salt-affected soils with salt encrustation at the surface were smoother than nonsaline surfaces and had a high reflectivity in the visible and near-infrared bands.They also discovered that in Landsat bands 1, 2, 3, and 4, there was a substantial link between EC values and reflectance.
In soil salinity research, the thermal band has shown to be an effective technique.It plays an important role in distinguishing saline soils, particularly in locations where there is bare soil or scant vegetation.Verma et al. [64], for example, proved that adding the Landsat TM thermal band to the visible-NIR bands helps address spectral similarity difficulties with saline soils.The thermal band was also employed to distinguish between salt-and sodium-affected soils.Metternicht and Zinck [83] discovered that including the thermal band improved salt and sodium detection.Furthermore, AlaviPanah and Goossens [84] discovered that adding the thermal band to the best Landsat TM visible-NIR band combination has excellent promise for distinguishing saline soil from gypsiferous soil in an Iranian case study.This research corroborated the findings of Goossens et al. [85] who found that the TM thermal band plays a critical role in distinguishing gypsiferous from salty soils.Furthermore, Furby et al. [86] claimed that multispectral satellite sensors confounded reflectance because nonsaline soils were mistaken with bare, exces-sively saline places.It was found that Landsat's spectral resolution was inadequate as the difference between the spectra of salty land and wet land was not sufficient to enable spectral separation [87], as well as the fluctuating spectral response of saline soil [88].The proposed index illustrated a high correlation with the measured EC and ESP in the lab with an accuracy of 96%, which is not in agreement with Bannari et al. [89], who presented three alternative salinity indices, SI-1, SI-2, and SI-3, based on the EO-1 ALI spectral bands, to distinguish between modest and moderate soil salinity and sodicity in Morocco.Despite the fact that SI-3 had the greatest correlation (46.9%), the data from this index were insufficient to offer exact information.As a result, they created two more soil salinity and sodicity indices.Their findings showed that these SSSI indices were likely to improve the identification accuracy in low-and medium-salinity regions as they the highest correlation (52.9%) with ground EC measurements.To measure soil salinity in Pakistan, Abbas and Khan [82] proposed an integrated technique based on geographical analysis of both ground and satellite data.To identify soil salinity, remote sensing data-based salinity indices and a Principal Component Analysis (PCA) were created.Their findings revealed that, when compared to ground data, S3 offered the most promising outcome of the six salinity indices.
The SIA model outperformed all other salinity indices models and showed the highest R-square (0.90) relationship between the EC salinity measures and the SIA in 2021.While the model showed and R-square value of 0.88 for the data of year 2016 followed by value of 0.83 for 2002.

Conclusions
The identification of the spatial extent of SAfSoil is a very important task for crop manufacturing.This drove many authors to develop efficient remote sensing algorithms (indices) to map this type of soil (e.g., [71][72][73][74]).Results of mapping of the SAfSoil using the previously proposed indices were not conclusive when applied to the study area in the northern Nile Delta region, Egypt.The current study proposed a new SAfSoil index, which took into account different parameters (e.g., vegetation and water saturation).The proposed index, on the other hand, was very successful in delineating the SAfSoil as well as nonsaline soil and waterlogging in the study area.The newly proposed SAfSoil index was tested using Landsat-7 and Landsat-8 data acquired in 2002, 2016, and 2021, providing near accurate results.These results were verified using field measurements performed during the same time frames of the images.
The newly developed SAfSoil index (SIA) achieved a much greater result when compared to the previously published indices.It can also be used worldwide in locations of similar environments.

Figure 1 .
Figure 1.Location of Kafr El-Sheikh Governorate area and land use/land cover units.

Figure 2 .
Figure 2. Historic climate data from 2005 to 2015 of Kafr El-Sheikh Governorate area.

Table 3 .
Salinity indices as proposed by different authors (where R = Red, G = Green, B = Blue, NIR = Near-Infrared, and SI = Salinity Index).

Figure 3 .
Figure 3. Soil profile location and land use land cover, geomorphology, geology, and landscape of the area.

Table 4 .
Vegetation and water indices as proposed by different authors (where R = Red, NIR = Near-Infrared, SWIR = Shortwave-Infrared, NDVI = Normalized difference vegetation index, and NDWI = Normalized difference water index).

Figure 6 .
Figure 6.3rd part of the equation (Blue − Red/Blue + Red) applied on 2021.The image resulting from the proposed indices was tested using the satellite image and field data obtained during 2002.The validation of the indices was carried out using the field data collected during 2016 and 2021 from different key locations in the area.The field data of 2016 and 2021 comprised 20 to 60 samples and GPS locations from the five different classes we used.The suggested formula of Congalton and Green [79] was used to calculate the Kappa coefficient and estimate the overall accuracy of our results.A flowchart showed the main steps of the study work, illustrated in Figure 7.

Figure 7 .
Figure 7. Flowchart showing the main steps of the study.

Figure 8 .
Figure 8.A 2002 map showing the SAfSoil results from the indices of Khan et al. [71].

Figure 9 .
Figure 9.A 2002 map showing the SAfSoil results from the indices of Farifteh [72].

Figure 10 .
Figure 10.A 2002 map showing the SAfSoil results from the indices of Setia et al. [73].

Figure 11 .
Figure 11.A 2002 map showing the SAfSoil results from the indices of Dehni and Lounis [74].

Figure 12 .
Figure 12.A 2002 map showing the salt-affected soils in the study area.

Figure 13 .
Figure 13.A 2016 map showing the salt-affected soils in the study area.

Figure 14 .
Figure 14.A 2021 map showing the salt-affected soils in the study area.

Table 1 .
Analysis of soil samples used to verify the results of the proposed Index.

Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) for the Years 2016 and 2021
* ETM+ Band 6 is acquired at 60 m resolution, but products are resampled to 30 m pixels.** TIRS bands are acquired at 100 m resolution, but are resampled to 30 m in delivered data product.

Table 5 .
Accuracy assessment using error matrix in 2002 based on soil survey of 20 soil profiles.

Table 6 .
Accuracy assessment using error matrix in 2016 based on soil survey of 20 soil profiles.

Table 7 .
Accuracy assessment using error matrix in 2021 based on soil survey of 20 soil profiles.Looking at the classification results for 2002, 2016, and 2021, we can easily see how the different classes changed with time (Table

Table 8 .
Output classes from the proposed SI.