Monitoring the Variations of Soil Salinity in a Palm Grove in Southern Algeria

: Soil salinity is considered the most serious socio-economic and environmental problem in arid and semi-arid regions. This study was done to estimate the soil salinity and monitor the changes in an irrigated palm grove (42 ha) that produces dates of a high quality. Topsoil samples (45 points), were taken during two di ﬀ erent periods (May and November), the electrical conductivity (EC) and Sodium Adsorption Ratio (SAR) were determined to assess the salinity of the soil. The results of the soil analysis were interpolated using two geostatistical methods: inverse distance weighting (IDW) and ordinary Kriging (OK). The e ﬃ ciency and best model of these two methods was evaluated by calculating the mean error (ME) and root mean square error (RMSE), showing that the ME of both interpolation methods was satisfactory for EC ( − 0.003, 0.145) and for SAR ( − 0.03, − 0.18), but the RMSE value was lower using the IDW with both data and periods. This can explain the accuracy of the IDW interpolation method. This model showed a dominance of soil salinity distribution in the South and South-East of the study area during the ﬁrst season, and for the second season, the salts were concentrated in the middle of the area. Several factors could interact in this variation such as the topographic direction of the water ﬂow and the aridity of the climate (evaporation). From this study emerges the need to maintain a better management of agricultural water and soils, avoiding salt accumulation, to ensure a good yield and the sustainability of agriculture in arid environments.


Introduction
Soil salinization is the process in which a level of salt accumulates in both the surface and the subsurface of soils [1]. It is a global environmental concern, although it mainly affects arid and semi-arid areas. It also causes soil degradation through sterilization [2] and threatens the sustainable development of oasis ecosystems and the future development of agriculture with scarce availability of resources (water and soil). This can occur naturally or due to poor management practices caused by humans [3].
The spatial variability of soil salinity relies on several factors: soil (characteristics, permeability, depth), water quality (availability, salinity, groundwater), local topography, climatic factors (low precipitation, high temperature and evaporation) and mainly mismanagement (irrigation and drainage system) [4][5][6][7][8][9][10][11], particularly in arid regions [12,13]. The region is part of the Great Desert, which is characterized by short and rigorous winters, while the summers are long and hot. The average annual temperature exceeds 25 °C. July is the hottest month (41.50 °C), while January is the coldest (6.20 °C). Rain is rare and irregular, with average annual rainfall ranging from 100 to 200 mm/year, with evaporation of 2000 mm per year [43].
The study area is characterized by an almost equal topography, with a difference of altitude not exceeding 1.5 m between the highest and the lowest points of the palm grove (highest point in the NW to the lowest point in the SE); the slope is under 0.5%. The soils in the region are mainly classified as Psamment in the Soil Taxonomy [44]; these are soils with a sandy texture, which consist of unconsolidated sand deposits or Arenosols, considering the World Reference Base for Soil Resources (WRB) [45]. After intensive cultivation for long periods, the agricultural soils of these oases can be considered Anthrosols (WRB).
The slope of the palm grove goes from the northwest (NW) to the southeast (SE) of the study area, determining the water flow. There, a non-permanent natural watercourse where water intermittently runs only a few days a year is located. This system acts as natural drainage for the area and the oasis.
The soil is irrigated with water from the Continental Interlayer. Each plot is fed by a tertiary network of open concrete canals (Seguias) that carry the water to the entrances of the plots. Irrigation is organized in rotation with one water turn per week, and farmers have water determined according The region is part of the Great Desert, which is characterized by short and rigorous winters, while the summers are long and hot. The average annual temperature exceeds 25 • C. July is the hottest month (41.50 • C), while January is the coldest (6.20 • C). Rain is rare and irregular, with average annual rainfall ranging from 100 to 200 mm/year, with evaporation of 2000 mm per year [43].
The study area is characterized by an almost equal topography, with a difference of altitude not exceeding 1.5 m between the highest and the lowest points of the palm grove (highest point in the NW to the lowest point in the SE); the slope is under 0.5%. The soils in the region are mainly classified as Psamment in the Soil Taxonomy [44]; these are soils with a sandy texture, which consist of unconsolidated sand deposits or Arenosols, considering the World Reference Base for Soil Resources (WRB) [45]. After intensive cultivation for long periods, the agricultural soils of these oases can be considered Anthrosols (WRB).
The slope of the palm grove goes from the northwest (NW) to the southeast (SE) of the study area, determining the water flow. There, a non-permanent natural watercourse where water intermittently runs only a few days a year is located. This system acts as natural drainage for the area and the oasis.
The soil is irrigated with water from the Continental Interlayer. Each plot is fed by a tertiary network of open concrete canals (Seguias) that carry the water to the entrances of the plots. Irrigation is organized in rotation with one water turn per week, and farmers have water determined according to the number of palm groves. The recommended standard for a palm grove is more than 70 L/s·ha, and the daily cycle does not exceed 7 days [46]. In our case, the amount of irrigation does not exceed 50 L/s·ha, and the daily cycle is in some cases 8 days. In summer, the amount decreases up to 50%. This quantity does not cover the needs of the plants. In addition, ditches and canals are in poor condition to provide good drainage. The drainage canals are about 1 m deep, located around the orchards, forming the drainage system. However, the abandonment of the ditches leads to stagnation of drainage water due to poor management of the system.
In order to improve the organic fertility of the soil, farmers add around 20 to 30 kg of sheep manure and cover it with sand at the bottom of each date palm to prepare for cultivation.

Soil Sampling and Analysis
To explore the spatial distribution of soil salinity, soil samples were taken according to systematic planning covering the entire study area and ensuring at least one sample per plot ( Figure 1). Samples were taken beside the palm tree, approximately 1.5 m from the stipe. 45 soil-sampling points were taken in 2018 in two periods: May and November. These two months coincide with a period of major agricultural work (beginning of the fruiting season and end of fruit maturation) [47]. At each point, we took 0-30 cm depth samples (topsoil).
The samples obtained were transferred to the laboratory and air dried, crushed and sieved at 2 mm. In the fine earth, electrical conductivity was determined in water extraction (EC; ratio 1:2.5 w/v) following the procedure described by the United States Salinity Laboratory (USDA) [28]. This method consists of mixing the sample with a lot of water to obtain a very strongly diluted solution. It allows for maximum extraction of salts as opposed to saturated paste extracts, the soil/water ratio remaining constant whatever the texture of the sample [48]. Soluble cations (Na + . Ca 2+ and Mg 2+ ) were analyzed using standard USDA procedures [49].

Characteristics of Irrigation Water
The sustainability of agricultural yield in oases depends on the quality and quantity of irrigation water. In our study area, soils are irrigated with water from the Continental Interlayer. Table 1 shows the mean physico-chemical characteristics of the water used for irrigation and drainage water, taken in June ( Figure 2). Knowing that all the soils were irrigated with the same water, the salinity may be influenced by other factors like farm management and environmental conditions.

Statistical Analysis
The data were subjected to a standard analysis to obtain descriptive statistics, specifically the mean, minimum and maximum, median, variance, standard deviation (SD), coefficient of variation (CV), kurtosis and skewness of each parameter. To identify the normal distribution of the data, skewness is the most common statistical parameter used with values ranging from −1 to +1. For this preliminary analysis, the normality of the data was assessed before using geostatistics to obtain prediction maps. The normality of each dataset was additionally verified by the (QQ plot) test to ensure normal distribution.

Predictive Mapping
The IDW and OK geostatistical techniques were used to determine the spatial variability of soil salinity in both periods (spring and autumn). Moreover, after that, we compare the two methods by using ME and RMSE.

The IDW Interpolation
The Inverse Distance Weighting (IDW) interpolation is a spatial prediction technique [50] commonly used in Geosciences. This tool attributes each entry point to a local influence that decreases with distance and calculates the values of prediction for an unknown interpolated point by weighting the medium of known data point values. This method can be used if enough sample points have a distribution that occupies the area on a local scale. It weights the points nearest to the prediction point to those more distant. IDW is an exact and convex interpolation method that adapts only to the continuous model of spatial variation.

Statistical Analysis
The data were subjected to a standard analysis to obtain descriptive statistics, specifically the mean, minimum and maximum, median, variance, standard deviation (SD), coefficient of variation (CV), kurtosis and skewness of each parameter. To identify the normal distribution of the data, skewness is the most common statistical parameter used with values ranging from −1 to +1. For this preliminary analysis, the normality of the data was assessed before using geostatistics to obtain prediction maps. The normality of each dataset was additionally verified by the (QQ plot) test to ensure normal distribution.

Predictive Mapping
The IDW and OK geostatistical techniques were used to determine the spatial variability of soil salinity in both periods (spring and autumn). Moreover, after that, we compare the two methods by using ME and RMSE.

The IDW Interpolation
The Inverse Distance Weighting (IDW) interpolation is a spatial prediction technique [50] commonly used in Geosciences. This tool attributes each entry point to a local influence that decreases with distance and calculates the values of prediction for an unknown interpolated point by weighting the medium of known data point values. This method can be used if enough sample points have a distribution that occupies the area on a local scale. It weights the points nearest to the prediction point to those more distant. IDW is an exact and convex interpolation method that adapts only to the continuous model of spatial variation.
A general form of predicting an interpolated value Z at a given point x based on samples Zi = Z(x i ) for i = 1, 2, . . . , N using IDW is an interpolating function: As defined by Shepard [50], the equation presented above designates a simple IDW weighting function, where x stands for the predicated location of an unknown interpolated point, x i is the known data point, d is the distance from x i to x, n refers to number of points used in interpolation, and p resembles an arbitrary real positive number known as the distance-decay or the power parameter (normally α = 2 in the standard IDW). Note that in the standard IDW, parameter α is a specified constant value defined by the user for all unknown interpolated points [51]. The selected interpolation technique is commonly used for this type of study and is the most common scatter-point interpolation method. It is based on the fundamental assumption that the interpolation surface should be influenced most by near points and least by far points. The interpolation surface is a weighted average of scattering points, and the weight assigned to each scattering point decreases as the distance between the interpolation point and the scattering point increases. In fact, values at unknown points are calculated as a weighted average of the values available at known points [52,53].

The Ordinary Kriging Interpolation
Kriging is a group of geostatistical methods for interpolating the values of different regional variables at an unobserved location from observations of its value at nearby locations, consisting of ordinary Kriging, universal Kriging, indicator Kriging, co-Kriging, and others [54][55][56].
Ordinary Kriging (OK) is a commonly used method [57], and it was applied in this study. The OK method plays a special role because it is compatible with a stationary model, involves only the variogram, and it is in fact the form of Kriging that is most often used [58,59].
This method estimates the values of a target variable in unmeasured locations, where Z*(x 0 ) is at an unsampled location, using measured values Z(x i ) (i = 1, 2, . . . , N), [60,61] as follows: λi returns the weight attributed to the i observation. According to [62], weights are assigned to each sample so that the variance of the estimate is minimized and estimates are unbiased.
Semivariograms are considered a basic tool for assessing the spatial distribution of soil properties. Based on the theory of regionalized variables and intrinsic assumptions [63], a semivariogram is expressed as follows: where γ(h) is the semivariance, h is the lag distance, Z is the parameter of the soil property, N(h) is the number of pairs of locations separated by a lag distance h, Z(x i ), and Z(x i + h) are values of Z at positions x i and x i + h [64].

Comparison of Methods
For better selection of the method used to determine the spatial variability of soil salinity, both interpolation methods have been evaluated to verify the compatibility of the relative performance of IDW and OK [65].
In order to evaluate the accuracy of the determined model, two statistical indices, such as the Mean Error (ME) and the Root Mean Square Error (RMSE), were calculated using the following equations: -Mean error: -Root mean square error: where Z(x i ) is the observed value at location i, Z*(x i ) is the predicted value at location i, and n is the sample size. Squaring the difference at any point gives an indication of the magnitude of differences, in such a way that a value of RMSE close to zero illustrates the accuracy of the prediction of the model. It is assumed that if the variogram model is correct, ME should be almost zero [66][67][68].
According to [68][69][70], the indicator of spatial dependence can be classified according to the nugget/sill ratio. The variable considers a strong spatial, a moderate spatial and a weak spatial dependence if the ratio is ≤25%, 25-75% and ≥75%, respectively.
The choice of an appropriate method depends on the particularities of the data and the type of spatial model desired. The accuracy of prediction of ordinary Kriging (OK) and inverse distance weighted (IDW) was evaluated by using the electrical conductivity and SAR data.
All data processing and analysis was carried out using ArcGis 10.2 © software and statistical processing with Microsoft Excel ©.

Soil Salinity
As a common rule, high concentrations of salts in the irrigation water and the nature of the hyper-arid climate, combined with many other factors, significantly influence soil quality, causing an increased presence of salts in the root zone and leading to a reduction in crop growth and agricultural yields [71,72]. However, the quality of the irrigation water in our agrosystem is admissible, and the salinity maybe associated with salts remaining in the soil after a long period of cultivation. This may be a consequence of climatic conditions and the accumulation of salts in the topsoil through capillarity during periods where no irrigation is applied (end of autumn and winter).
Statistical analysis of the results ( Table 2) reveals that in the first season, the EC value varies between 0.15 dS/m and 7.27 dS/m with a coefficient of variation of 104.84, and the SAR value varies between 1.33 and 16.98 with a coefficient of variation of 77.03. In the second season, the EC value varies from 0.21 dS/m to 9.39 dS/m with a coefficient of variation of 129.94, and the SAR value varies from 0.3 to 30.93 with a coefficient of variation of 131.53. The Coefficient of Variation (CV) can reflect the degree of dispersion of random variables: a CV < 10% generally denotes weak variability, while 10% < CV < 100% denotes moderate variability, and CV > 100% denotes strong variability [73], demonstrating that both conditions exhibit strong variation. In our case, this indicates that both periods exhibit strong variability, except for the SAR in the first period, which exhibits moderate variation. Standard deviation and CV were greater in the second period, indicating more differences associated with the location of sampling points, showing the highest variability in soil salinity.  Based on the classes determined by Durand (1983) [74], who considers five salinity classes (Table 3), more than a third of the soil samples would be considered non-saline soils. The analysis of this table (Table 3) allows us to observe that the percentages of saline soils in the first and second periods were 41.49% and 34.8%, respectively (yield can be affected), whereas non-saline soils represent 54.17% and 58.68% (productive soil), and 4.34% and 6.52% represented moderately saline soils (soil that can be used for tolerant crops). Moreover, these results show that, even when using the same irrigation water and the same soil sample points, differences in salinity occurred. Even more, important changes were found in the classes with high salinity, maybe due to the movement of salts. However, we can notice that approximately 55-60% of the soils remained in the salinity classes of non-saline soil and slightly saline soil during both periods.
It is important to understand the spatial distribution of salinity in order to check where the salts have been accumulated and, in this sense, make better decisions for managing the system.

Distribution of Soil Salinity Using IDW and OK Interpolation
In this study, interpolation IDW (inverse distance weighting) and OK (ordinary Kriging) were used to estimate the spatial distribution of soil salinity.
IDW is one of the most common techniques for interpolating scattering points and is based on the fundamental assumption that the interpolation area should be influenced most by near points and least by far points. The interpolation area is a weighted average of the scattering points, and the weight assigned to each scattering point decreases as the distance between the interpolation point and the scattering point increases. Values at unknown points are estimated as a weighted average of the values available at known points [52,53]. Summary statistics of IDW interpolation are presented in Table 2. For OK interpolation, the data were verified by a normality test (QQ Plot) to evaluate the distribution of the data, which confirms that the skewness value is greater than 1 (the results do not follow a normal distribution). These data were normalized using the logarithmic decimal method (Table 4). Subsequently, variograms of the Kriging method were calculated from the normalized data. The best theoretical optimal model used here for EC and SAR in both periods was the Exponential model to estimate the distribution of soil salinity (Table 5). The spatial dependence "nugget-to-sill" for soil salinity was in the class (<25%), which corresponds to high spatial dependence (Table 5). In addition, the range of the chosen model is 25.243 m (Figure 3).  Subsequently, variograms of the Kriging method were calculated from the normalized data. The best theoretical optimal model used here for EC and SAR in both periods was the Exponential model to estimate the distribution of soil salinity (Table 5). The spatial dependence "nugget-to-sill" for soil salinity was in the class (< 25%), which corresponds to high spatial dependence (Table 5). In addition, the range of the chosen model is 25.243 m (Figure 3).  This information, derived from the semivariograms, highlighted the reality of different spatial dependence for the soil properties collected in the field ( Table 5). The proportion of nuggets to the threshold mimics spatial autocorrelation.
The IDW interpolation method ( Figure 4) and Ok interpolation ( Figure 5) were used to perform a better spatial interpolation method for unsampled points, to obtain the soil salinity distribution pattern in the oasis. This information, derived from the semivariograms, highlighted the reality of different spatial dependence for the soil properties collected in the field ( Table 5). The proportion of nuggets to the threshold mimics spatial autocorrelation.
The IDW interpolation method ( Figure 4) and Ok interpolation ( Figure 5) were used to perform a better spatial interpolation method for unsampled points, to obtain the soil salinity distribution pattern in the oasis.   According to the maps of the spatial distribution of EC and SAR using IDW ( Figure 4) and OK interpolation ( Figure 5), the highest soil salinity (EC) was mainly observed in the southern parts of the study area in the first period, which is related to the direction of irrigation water flow observed during the survey. This could be a result of the stagnation of drainage water due to the presence of a depression nearby that can naturally drain irrigation water. Large quantities of water subjected to evaporation will eventually lead to an increase in salinity, whereas in the second period, the salts migrated to external parts of the study area, which is related to high evaporation due to excessively high temperature and reduced irrigation frequency.
A soil salinity map is needed to assess and monitor salt accumulation. This information can be used as an effective tool to help farmers manage soil salinity problems affecting the agroecosystem. Moreover, the simplicity of this approach, with its satisfactory accuracy, can greatly contribute to the prediction and mapping of soil salinity [67,75]. In this sense, farmers will have information to detect saline soils and improve their management in the next period of cultivation. According to the maps of the spatial distribution of EC and SAR using IDW ( Figure 4) and OK interpolation ( Figure 5), the highest soil salinity (EC) was mainly observed in the southern parts of the study area in the first period, which is related to the direction of irrigation water flow observed during the survey. This could be a result of the stagnation of drainage water due to the presence of a depression nearby that can naturally drain irrigation water. Large quantities of water subjected to evaporation will eventually lead to an increase in salinity, whereas in the second period, the salts migrated to external parts of the study area, which is related to high evaporation due to excessively high temperature and reduced irrigation frequency.
A soil salinity map is needed to assess and monitor salt accumulation. This information can be used as an effective tool to help farmers manage soil salinity problems affecting the agroecosystem. Moreover, the simplicity of this approach, with its satisfactory accuracy, can greatly contribute to the prediction and mapping of soil salinity [67,75]. In this sense, farmers will have information to detect saline soils and improve their management in the next period of cultivation.

Comparison of IDW and OK
The comparison of IDW and OK was performed using ME and RMSE values to define the best method of spatial interpolation suited to the region. Measures for evaluating model performance are summarized in Table 6. Based on a comparison between the results, the ME of both interpolation methods was satisfactory (values were close to zero) [68,76], but the RMSE value was lower when using IDW. According to these geostatistical results, there was a slight difference in the performance of IDW when predicting salinity, compared with OK. This may explain the good accuracy of the IDW method.
The OK map represented the spatial variability with less detail, demonstrating the importance of this interpolation to variability on long-distance scales, while the IDW map represented the spatial variability with more detail on short-distance scales, supporting fairly dense sampling.
Other studies that have been conducted using this approach, for example [77][78][79], confirm that IDW interpolation is more efficient than OK in predicting distribution, especially of the physico-chemical properties of soil, and this was probably due to the method of sampling and the similar distances of all samples. In contrast, other studies, for example [80][81][82], considered the performance of ordinary Kriging to be the best in comparison with other predicted data and with reference to actual data, concluding that it was also very reliable and more accurate.
We presume that the sampling interval that occupies the entire prediction scale can explain the performance of IDW in interpolating the soil property [65]. Other studies [83,84] provide maps of salinity distribution by using geostatistical methods and considered the co-Kriging method to be the most accurate method for predicting EC and SAR levels. In fact, the sampling method may determine the best geostatistics method.
Understanding and controlling the spatial distribution of soil salinity can provide a database for implementing an agricultural management policy capable of controlling and preventing soil salinization [85]. To explore the spatial variability of soil salinity, we applied a geostatistical method in an agricultural area with an arid climate in order to identify the influencing factors. The analysis of the maps obtained is consistent with the real situation of the terrain, which is characterized by high salinity in areas of stagnant drainage water.
This distribution of soil salinity is explained by the method of irrigation through submersion and the direction of flows used in this agricultural perimeter, coupled with severe climatic conditions, especially in the summer period when temperature and evaporation are excessively high. In addition to this is the poor management of natural resources and modification of soil structure caused by humans.
Our contribution focused, in particular, on oasis farmland in an arid region. The results and methodology could be extended to other regions that experience similar problems in arid and semi-arid areas. Agriculture is critical for many communities, and by understanding and mapping the problems it is possible to achieve sustainability and create opportunities for these communities associated with oasis systems. Therefore, our study will be able to open a path for reflection and research, allowing the implementation of projects dealing with soil salinity on a large scale, as well as a detailed investigation of factors influencing the sustainability of agricultural production in drylands.
It is clear that sustainability is not optional, and we need to know how the system works to improve it [86]. These results enable planners, decision-makers and engineers to undertake and guide agricultural development work, which consists of setting up an adequate drainage network capable of removing excess salt-laden water from the farming area and providing additional irrigation (leaching dose) to maintain drainage channels.
These arid regions are known for the scarcity of water and its high mineralization. In order to avoid aggravating the water crisis and prevent considerable food shortages, the productivity of water use must be increased. Moreover, the treatment and use of agricultural drainage water (ADW) is becoming an obligation for development needs in arid and semi-arid zones [86]. In the future, one option to be considered in these regions is the reuse of agricultural drainage water, which could be treated using recycling techniques [87], for instance, through the establishment of treatment plants and/or the desalination of water by reverse osmosis systems with solar energy supply. In this sense, keeping drainage systems in good conditions is necessary.

Classification of Soil Salinity
After detecting and mapping the salinity, it is important to determine the location of salt-affected soils and the grade of this affection. The soils (sampling points) were classified by comparing the EC, SAR and pH values with the threshold values set by Abegunrin et al. (2013) [88], considering salinity and sodium or alkalinity risks as the EC, SAR and pH values were below 4 dS/m, 13 and 8.5, respectively (Table 7). With the exception of a few samples, most of the soils could be considered non-saline and non-sodic soils ( Figure 6). It is clear that sustainability is not optional, and we need to know how the system works to improve it [86]. These results enable planners, decision-makers and engineers to undertake and guide agricultural development work, which consists of setting up an adequate drainage network capable of removing excess salt-laden water from the farming area and providing additional irrigation (leaching dose) to maintain drainage channels.
These arid regions are known for the scarcity of water and its high mineralization. In order to avoid aggravating the water crisis and prevent considerable food shortages, the productivity of water use must be increased. Moreover, the treatment and use of agricultural drainage water (ADW) is becoming an obligation for development needs in arid and semi-arid zones [86]. In the future, one option to be considered in these regions is the reuse of agricultural drainage water, which could be treated using recycling techniques [87], for instance, through the establishment of treatment plants and/or the desalination of water by reverse osmosis systems with solar energy supply. In this sense, keeping drainage systems in good conditions is necessary.

Classification of Soil Salinity
After detecting and mapping the salinity, it is important to determine the location of salt-affected soils and the grade of this affection. The soils (sampling points) were classified by comparing the EC, SAR and pH values with the threshold values set by Abegunrin et al. (2013) [88], considering salinity and sodium or alkalinity risks as the EC, SAR and pH values were below 4 dS/m, 13 and 8.5, respectively (Table 7). With the exception of a few samples, most of the soils could be considered nonsaline and non-sodic soils.   The Zelfana region has undergone a major change in agriculture compared to oasis agriculture, linked to the almost exclusive production of dates, followed by the production of vegetable crops throughout the year. Irrigation with groundwater can enrich soils with salts or cause their mobility in the soil profile, which is associated with problems of soil management and conservation. This will undoubtedly have negative repercussions on the soil (salinization and sodification leading to degradation of the soil structure) and on agricultural dynamics in Zelfana. The sensitivity of crops to salt stress results in reduced plant growth (reduced leaf length of fronds) and reduced productivity. Each crop has a specific tolerance to saline stress, which can reduce yield by about 50% for several cultivars such as pepper, lettuce and potato for an EC threshold of 5 dS.m −1 [89].
The date palm has a high tolerance for soil salinity [26,[90][91][92], but several publications consider that an increase in salinity corresponds to a reduction in yield < 25% [93,94]. According to Ramoliya et al. (2003) [95], this crop can withstand a relatively high salinity of 12.8 ds/m with no visible effect on the fruit, and other varieties of date palms could tolerate up to 9 ds/m of soil salinity [96]. In order to obtain a better quantitative and qualitative yield, it is necessary to consider the nature of the soil and its salinity as well as the crop's tolerance to salinity, to innovate an adequate strategy to ensure good, efficient and sustainable management capable of maintaining the maximum yield of crops for future generations. One of the best strategies for the rehabilitation of saline soils is the proper selection of crops, not only for economic purposes but also for the control of soil and saline water pollution [97]. Crops could extract salts and act as phytoremediators because of their ability to adapt to local saline environments. Due to the results obtained, it seems that acquisition of the data has given us a clear vision of the current state of soil salinization in this area, which is useful for managers and farmers.
The right choice of market gardening crops according to the nature of the soil requires an adequate selection of crops according to their resistance to saline stress, followed by good use of water resources through the installation of irrigation and drainage networks that will maintain the aquifer level and ensure the discharge of contaminated water below the critical level of salinization in order to guarantee the sustainability of agricultural systems.
In recent years, remote sensing, GIS and modelling have become preferred and inexpensive technological tools for mapping soil salinity due to the large areas covered, which is of utmost importance from both agricultural and environmental perspectives. This opens up great prospects for the use of this technique in the sustainable management of soils.

Conclusions
The main objective of this study was to estimate the soil salinity in an oasis by comparing two spatial interpolation techniques based on GIS, in order to highlight a better sustainable management strategy. It is crucial to select an appropriate method for estimating the spatial distribution of soil salinity with minimum error, for good management practices. In this study, IDW and OK were the  4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
The Zelfana region has undergone a major change in agriculture compared to oasis agriculture, linked to the almost exclusive production of dates, followed by the production of vegetable crops throughout the year. Irrigation with groundwater can enrich soils with salts or cause their mobility in the soil profile, which is associated with problems of soil management and conservation. This will undoubtedly have negative repercussions on the soil (salinization and sodification leading to degradation of the soil structure) and on agricultural dynamics in Zelfana. The sensitivity of crops to salt stress results in reduced plant growth (reduced leaf length of fronds) and reduced productivity. Each crop has a specific tolerance to saline stress, which can reduce yield by about 50% for several cultivars such as pepper, lettuce and potato for an EC threshold of 5 dS.m −1 [89].
The date palm has a high tolerance for soil salinity [26,[90][91][92], but several publications consider that an increase in salinity corresponds to a reduction in yield < 25% [93,94]. According to Ramoliya et al. (2003) [95], this crop can withstand a relatively high salinity of 12.8 ds/m with no visible effect on the fruit, and other varieties of date palms could tolerate up to 9 ds/m of soil salinity [96]. In order to obtain a better quantitative and qualitative yield, it is necessary to consider the nature of the soil and its salinity as well as the crop's tolerance to salinity, to innovate an adequate strategy to ensure good, efficient and sustainable management capable of maintaining the maximum yield of crops for future generations. One of the best strategies for the rehabilitation of saline soils is the proper selection of crops, not only for economic purposes but also for the control of soil and saline water pollution [97]. Crops could extract salts and act as phytoremediators because of their ability to adapt to local saline environments. Due to the results obtained, it seems that acquisition of the data has given us a clear vision of the current state of soil salinization in this area, which is useful for managers and farmers.
The right choice of market gardening crops according to the nature of the soil requires an adequate selection of crops according to their resistance to saline stress, followed by good use of water resources through the installation of irrigation and drainage networks that will maintain the aquifer level and ensure the discharge of contaminated water below the critical level of salinization in order to guarantee the sustainability of agricultural systems.
In recent years, remote sensing, GIS and modelling have become preferred and inexpensive technological tools for mapping soil salinity due to the large areas covered, which is of utmost importance from both agricultural and environmental perspectives. This opens up great prospects for the use of this technique in the sustainable management of soils.

Conclusions
The main objective of this study was to estimate the soil salinity in an oasis by comparing two spatial interpolation techniques based on GIS, in order to highlight a better sustainable management strategy. It is crucial to select an appropriate method for estimating the spatial distribution of soil salinity with minimum error, for good management practices. In this study, IDW and OK were the two methods assessed to study soil salinity in the Zelfana oasis, Ghardaïa. Hence, the IDW method was more accurate than the results obtained by the OK technique for the estimation of salt-affected areas. In our case, the relation between drainage system and soil salinity seems to be a key factor in managing the agrosystem, emphasizing the importance of sustainable water management in oasis systems.
The obtained results showed that the majority of soils belong to a salinity class that can be adequate for date palm tolerance but are not recommended for other crops, especially in the southern part of the study area. Based on the EC, SAR and soil salinity distribution maps, it seems that in this area, soils with high salinity do not support low tolerance market gardening. Thus, this model (IDW) for monitoring the spatial distribution of soil salinity is a good indicator and can be used by phoeniculturists in the Zelfana oasis and generalized by decision-makers in similar areas, so as to implement effective management programs and the sustainable use of agricultural systems.
To overcome the constraints of field data collection, an extensive data collection, which is regarded as time-consuming and costly, mapped and identified sites where salinity is detected, helping to reduce the number of sample points and determine the actions necessary to maintain the system and aid sustainability in the long-time. However, further research is necessary to investigate the ability to use remote sensing techniques if they can increase the accuracy of soil salinity modeling and mapping for areas dominated mainly by date palm and market gardening in arid environments.