# Evaluation of Different Soil Salinity Mapping Using Remote Sensing Techniques in Arid Ecosystems, Saudi Arabia

Authors: Mohamed Elhag
Journal: Journal of Sensors
 (2016)

---

## Abstract

Land covers in Saudi Arabia are generally described as salty soils with sand dunes and sand sheets. Waterlogging and higher soil salinity are major challenges to sustaining agricultural practices in Saudi Arabia principally within closed drainage basins. Agricultural practices in Saudi Arabia were flourishing in the last two decades. The newly reclaimed lands were added annually and distributed all over the country. Irrigation techniques are mostly modernized to fulfill water saving strategies. Nevertheless, water resources in Saudi Arabia are under stress and groundwater levels are depleted rapidly due to heavy abstraction that may exceed crop water requirements in most of the cases due to high evaporation rates. The excess use of irrigational water leads to severe soil salinity problems. Applications of remote sensing technique in agricultural practices became widely distinctive and cover multidisciplinary principal interests on both local and regional levels. The most important remote sensing applications in agricultural practices are vegetation indices which are related to vegetation and water especially in an arid environment. Soil salinity mapping in an arid ecosystem using remote sensing data is a demanding task. Several soil salinity indices were implemented and evaluated to detect soil salinity effectively and quantitatively. Thematic maps of soil salinity were satisfactorily produced and assessed.


---

## Body

## 1. Introduction

Saudi Arabia has practiced an eminent development in vital sectors over the last three decades. As a result, rapid escalation in terms of water demands in domestic, industrial, and agricultural activities have been countersigned. Agriculture is the keystone water consumption sector; it consumes about 85% of total national water use (Ministry of Water and Electricity) [1]. During the period of 1974–2006, Saudi Arabian government subsidized the agricultural sector to improve living standards in urban areas and to achieve sustainability. These subsidies lead to significant increases in agricultural areas, primarily for rigorous irrigated crops such as wheat, maize, and alfalfa; the irrigated areas increased from less than 400,000 ha in 1971 to about 1.62 million ha in 1992 [2, 3]. The extensive groundwater extrication subsequently leads to higher soil salinity rates in most of the agricultural zones [4–6].Excessive irrigation of less developed drained soil has led to high floating water tables, waterlogging, and soil salinization [7]. In the absence of well drainage network in shallow soil, water table depth intends to develop remarkably and, eventually, irrigable lands will be abandoned [8]. In the incidence of a higher soil salinity, crop production can be dramatically affected when salts accumulate in the soil surface [9–11]. Conservatively, soil salinity and waterlogging can be regulated through sustaining the salt concentration at the root zone below the crops constraint [12]. This can be realized by maintaining the equilibrium between soil recharge and discharge by draining the surplus irrigational water [13].Dwivedi [14] used the first four bands of MSS to designate salty soils. Chaturvedi et al. [15] and Singh and Srivastav [16] used brightness and thermal temperature layers to detect soil salinity. Menenti et al. [17] used 3 bands of TM (1, 5, and 7) to detect salt minerals. Moreover, Mulders and Epema [18] used three different bands of TM (3, 4, and 5) to detect calcareous surfaces as an indication of high salt content. Saha et al. [19] used the thermal infrared band of TM to classify inland waterlogged areas. Rao et al. [20] used the false color composition of TM bands (2, 3, and 4) to distinguish between moderate and high soil content soils. Naseri [21] used both unsupervised and supervised classification algorithms to delineate soil salinity on a regional scale using three bands of MSS and four bands of TM. Darvishsefat et al. [22] used principal component analysis to map soil salinity using four bands of ETM+ (3, 4, 5, and 7). Zhang et al. [23] used different bands of OLI-8 to develop a soil salinity detection algorithm based on the Soil Salinity Information Extraction (SSIE).Implementation of remote sensing data in soil salinity mapping is firm to employ since salinization is not a static process. Therefore, a number of derivatives and alternatives to normalized difference vegetation index have been anticipated to discourse soils salinization monitoring and mapping [24–26]. Different remote sensing indices such as the salinity index (SI), the normalized difference salinity index (NDSI), and the brightness index (BI) in addition to normalized difference vegetation index (NDVI) were implemented to investigate the means by which these indices work for soil salinity mapping in the arid environment [27, 28]. The salinity index (SI), which combines the blue and red bands, is sensitive to the surface reflectance of salt-affected land with sparse vegetation cover [27].The perspective of using remote sensing data and GIS practices has been demonstrated in several scholarly works to be the most resourceful [29–32]. The comprehensive use of satellite remote sensing and GIS has been recognized to be a cost-effective method for monitoring soils salinization in poorly drained basins [33, 34]. However, state-of-the-art outcome specifies that there is no ideal combination of data type and questioning method which can be functional with identical realization across approximately several environmental conditions [25–27].The aim of current research project is to differentiate the salinized soils from nonsalinized soils implementing different approaches of digital image classification and band math techniques. To reach that objective, several soil salinity indices were applied to the Landsat 8 data verified with ground truth data to validate the percentage of estimations accuracy based on different error’s assessment approach.


## 2. Materials and Methods

### 2.1. Study Area

The study area, Wadi Al Dawasir town, is located in the plateau of Najd at latitude of 44°43′ and longitude of 20°29′, about 300 km south of the capital city Riyadh (Figure1). This study area is comprised of gravelly tableland disconnected by insignificant sandy oases and isolated mountain bundles. Across the Arabian Peninsula as a whole, the tableland slopes toward the east from an elevation of 1,360 meters in the west to 750 meters at its easternmost limit. Wadi Al Dawasir and Najran regions are the major irrigation water abstraction from Al-Wajid Aquifer. Agriculture in Wadi Al Dawasir area consists of technically highly developed farm enterprises that operate modern pivot irrigation system. The size of center pivot ranges from 30 ha to 60 ha with farms managing hundreds of pivots with the corresponding number of wells. The main crop grown in winter is wheat and occasionally potatoes, tomatoes, or melons. All year fodder consists of alfalfa, which is cut up to 10 times a year for food. Typical summer crops for fodder are sorghum and Rhodes grass, which is perennial but dormant in winter. The shallow alluvial aquifers could not sustain the high groundwater abstraction rates for a long time and groundwater level declined dramatically in most areas. Meteorological features of the area are speckled. Five elements of meteorology are constantly recorded through fixed weather station located within the study area. Temperature varies from 6°C as minimum temperature to 43°C as maximum temperature. Relative humidity is mostly stable at 24%. Solar radiation of average sunrise duration is generally 11 hrs/day. Average wind speed is closer to 13 km/hr and may reach up to 46 km/hr in thunderstorm incidents. Finally, mean annual rainfall is about 37.6 mm [5].Figure 1
Location of the study area.


### 2.2. Methodological Framework

The current research work is based on assessing a regression correlation between ground truth data and its corresponding salinity index’s values conducted from satellite images. Therefore, accurate synchronization of ground truth data collection and satellite bypassing were exercised.


#### 2.2.1. Soil Sampling and Analysis

A total number of 150 random soil samples were collected from Wadi Al Dawasir area with a minimum distance of 100 m between the samples locations to avoid data clumping (Figure2). The standard procedure of determining soil extract salinity in terms of Electrical Conductivity (EC) was followed according to Shaw [35] under laboratory condition.Figure 2
Soil sample location on false-color Landsat 8 image.


#### 2.2.2. Remote Sensing Data Enhancement

Remote sensing data was obtained from Landsat Operational Land Imager (OLI-8) which was acquired on April 17th, 2014. Typical atmospheric and radiometric corrections and spatial resolution enhancement were implemented for each band individually. Furthermore, overcorrection was analyzed for the red and near infrared bands [36]. Atmospheric correction is basically considering blue band offset. Therefore, the largest atmospheric consequence is theoretical to be the most truthful value within the blue band. The calibrated radiance value of this offset allows deciding κ value. Therefore, λ-κ rule the offset values for the green, red, and near infrared band under investigation [37].


#### 2.2.3. Soil Salinity Indices

Soil salinity detection using remote sensing techniques developed through several scholarly works. Landsat Multispectral Scanning (MSS), TM (Thematic Mapper), and Enhanced Thematic Mapper (ETM+) in addition to Operational Land Imager (OLI-8) images were the most used remote sensing data in soil salinity mapping. Several algorithms of detecting and mapping soil salinity were developed because of the sensors and/or study area variability. There were many attempts to use the remote sensing data and its techniques to detect the salt-affected soil. Soil salinity indices are principally adjusted to detect salt mineral in soils based on the different responses of salty soils to various spectral bands. Table1 summarizes different indices used to map soil salinity in current study.Table 1
Soil salinity indices based on different band ratios of Landsat.
Salinity index
Band ratios
Reference
Normalized difference salinity index
N
D
S
I
=
R
-
N
I
R
R
+
N
I
R
[43]
Brightness index
B
I
=
R
2
+
N
I
R
2
[44]
Salinity index-1
S
I
=
B
×
R
[44]
Salinity index-2
S
I
=
G
×
R
[44]
Salinity index-3
S
I
=
G
2
+
R
2
+
N
I
R
2
[27]
Salinity index-4
S
I
=
G
2
+
R
2
[27]
Salinity index-5
S
I
=
B
R
[45]
Salinity index-6
S
I
=
B
-
R
B
+
R
[45]
Salinity index-7
S
I
=
G
×
R
B
[45]
Salinity index-8
S
I
=
B
×
R
G
[46]
Salinity index-9
S
I
=
N
I
R
×
R
G
[46]
B, blue band; G, green band; R, red band; NIR, near infrared band of Landsat 8 image.


### 2.3. Validation

Linear regression model was used to validate the correlation established between the ground truths data measured from the 150 soil samples collected and the corresponding salinity index’s values estimated using the aforementioned equations. Root square, root square adjusted, root mean square error, and mean of response are used to signify the inequality of variance and correlation of the linear regression model [38].


## 2.1. Study Area

The study area, Wadi Al Dawasir town, is located in the plateau of Najd at latitude of 44°43′ and longitude of 20°29′, about 300 km south of the capital city Riyadh (Figure1). This study area is comprised of gravelly tableland disconnected by insignificant sandy oases and isolated mountain bundles. Across the Arabian Peninsula as a whole, the tableland slopes toward the east from an elevation of 1,360 meters in the west to 750 meters at its easternmost limit. Wadi Al Dawasir and Najran regions are the major irrigation water abstraction from Al-Wajid Aquifer. Agriculture in Wadi Al Dawasir area consists of technically highly developed farm enterprises that operate modern pivot irrigation system. The size of center pivot ranges from 30 ha to 60 ha with farms managing hundreds of pivots with the corresponding number of wells. The main crop grown in winter is wheat and occasionally potatoes, tomatoes, or melons. All year fodder consists of alfalfa, which is cut up to 10 times a year for food. Typical summer crops for fodder are sorghum and Rhodes grass, which is perennial but dormant in winter. The shallow alluvial aquifers could not sustain the high groundwater abstraction rates for a long time and groundwater level declined dramatically in most areas. Meteorological features of the area are speckled. Five elements of meteorology are constantly recorded through fixed weather station located within the study area. Temperature varies from 6°C as minimum temperature to 43°C as maximum temperature. Relative humidity is mostly stable at 24%. Solar radiation of average sunrise duration is generally 11 hrs/day. Average wind speed is closer to 13 km/hr and may reach up to 46 km/hr in thunderstorm incidents. Finally, mean annual rainfall is about 37.6 mm [5].Figure 1
Location of the study area.


## 2.2. Methodological Framework

The current research work is based on assessing a regression correlation between ground truth data and its corresponding salinity index’s values conducted from satellite images. Therefore, accurate synchronization of ground truth data collection and satellite bypassing were exercised.


### 2.2.1. Soil Sampling and Analysis

A total number of 150 random soil samples were collected from Wadi Al Dawasir area with a minimum distance of 100 m between the samples locations to avoid data clumping (Figure2). The standard procedure of determining soil extract salinity in terms of Electrical Conductivity (EC) was followed according to Shaw [35] under laboratory condition.Figure 2
Soil sample location on false-color Landsat 8 image.


### 2.2.2. Remote Sensing Data Enhancement

Remote sensing data was obtained from Landsat Operational Land Imager (OLI-8) which was acquired on April 17th, 2014. Typical atmospheric and radiometric corrections and spatial resolution enhancement were implemented for each band individually. Furthermore, overcorrection was analyzed for the red and near infrared bands [36]. Atmospheric correction is basically considering blue band offset. Therefore, the largest atmospheric consequence is theoretical to be the most truthful value within the blue band. The calibrated radiance value of this offset allows deciding κ value. Therefore, λ-κ rule the offset values for the green, red, and near infrared band under investigation [37].


### 2.2.3. Soil Salinity Indices

Soil salinity detection using remote sensing techniques developed through several scholarly works. Landsat Multispectral Scanning (MSS), TM (Thematic Mapper), and Enhanced Thematic Mapper (ETM+) in addition to Operational Land Imager (OLI-8) images were the most used remote sensing data in soil salinity mapping. Several algorithms of detecting and mapping soil salinity were developed because of the sensors and/or study area variability. There were many attempts to use the remote sensing data and its techniques to detect the salt-affected soil. Soil salinity indices are principally adjusted to detect salt mineral in soils based on the different responses of salty soils to various spectral bands. Table1 summarizes different indices used to map soil salinity in current study.Table 1
Soil salinity indices based on different band ratios of Landsat.
Salinity index
Band ratios
Reference
Normalized difference salinity index
N
D
S
I
=
R
-
N
I
R
R
+
N
I
R
[43]
Brightness index
B
I
=
R
2
+
N
I
R
2
[44]
Salinity index-1
S
I
=
B
×
R
[44]
Salinity index-2
S
I
=
G
×
R
[44]
Salinity index-3
S
I
=
G
2
+
R
2
+
N
I
R
2
[27]
Salinity index-4
S
I
=
G
2
+
R
2
[27]
Salinity index-5
S
I
=
B
R
[45]
Salinity index-6
S
I
=
B
-
R
B
+
R
[45]
Salinity index-7
S
I
=
G
×
R
B
[45]
Salinity index-8
S
I
=
B
×
R
G
[46]
Salinity index-9
S
I
=
N
I
R
×
R
G
[46]
B, blue band; G, green band; R, red band; NIR, near infrared band of Landsat 8 image.


## 2.2.1. Soil Sampling and Analysis

A total number of 150 random soil samples were collected from Wadi Al Dawasir area with a minimum distance of 100 m between the samples locations to avoid data clumping (Figure2). The standard procedure of determining soil extract salinity in terms of Electrical Conductivity (EC) was followed according to Shaw [35] under laboratory condition.Figure 2
Soil sample location on false-color Landsat 8 image.


## 2.2.2. Remote Sensing Data Enhancement

Remote sensing data was obtained from Landsat Operational Land Imager (OLI-8) which was acquired on April 17th, 2014. Typical atmospheric and radiometric corrections and spatial resolution enhancement were implemented for each band individually. Furthermore, overcorrection was analyzed for the red and near infrared bands [36]. Atmospheric correction is basically considering blue band offset. Therefore, the largest atmospheric consequence is theoretical to be the most truthful value within the blue band. The calibrated radiance value of this offset allows deciding κ value. Therefore, λ-κ rule the offset values for the green, red, and near infrared band under investigation [37].


## 2.2.3. Soil Salinity Indices

Soil salinity detection using remote sensing techniques developed through several scholarly works. Landsat Multispectral Scanning (MSS), TM (Thematic Mapper), and Enhanced Thematic Mapper (ETM+) in addition to Operational Land Imager (OLI-8) images were the most used remote sensing data in soil salinity mapping. Several algorithms of detecting and mapping soil salinity were developed because of the sensors and/or study area variability. There were many attempts to use the remote sensing data and its techniques to detect the salt-affected soil. Soil salinity indices are principally adjusted to detect salt mineral in soils based on the different responses of salty soils to various spectral bands. Table1 summarizes different indices used to map soil salinity in current study.Table 1
Soil salinity indices based on different band ratios of Landsat.
Salinity index
Band ratios
Reference
Normalized difference salinity index
N
D
S
I
=
R
-
N
I
R
R
+
N
I
R
[43]
Brightness index
B
I
=
R
2
+
N
I
R
2
[44]
Salinity index-1
S
I
=
B
×
R
[44]
Salinity index-2
S
I
=
G
×
R
[44]
Salinity index-3
S
I
=
G
2
+
R
2
+
N
I
R
2
[27]
Salinity index-4
S
I
=
G
2
+
R
2
[27]
Salinity index-5
S
I
=
B
R
[45]
Salinity index-6
S
I
=
B
-
R
B
+
R
[45]
Salinity index-7
S
I
=
G
×
R
B
[45]
Salinity index-8
S
I
=
B
×
R
G
[46]
Salinity index-9
S
I
=
N
I
R
×
R
G
[46]
B, blue band; G, green band; R, red band; NIR, near infrared band of Landsat 8 image.


## 2.3. Validation

Linear regression model was used to validate the correlation established between the ground truths data measured from the 150 soil samples collected and the corresponding salinity index’s values estimated using the aforementioned equations. Root square, root square adjusted, root mean square error, and mean of response are used to signify the inequality of variance and correlation of the linear regression model [38].


## 3. Results and Discussion

The use of different soil salinity equations gives the results the differential preferences to achieve an adequate soil salinity estimation on a large scale using remote sensing data. The evaluation of different soil salinity indices was based on the regression model fit, effect significance, and model parameters. Table2 summarizes different regression significances.Table 2
Regression significances.
NDSI
BI
SI-1
SI-2
SI-3
SI-4
SI-5
SI-6
SI-7
SI-8
SI-9
R square
0.62
0.94
0.89
0.90
0.95
0.90
0.44
0.44
0.91
0.86
0.98
R square adjusted
0.62
0.94
0.89
0.90
0.95
0.90
0.44
0.44
0.91
0.86
0.98
Root mean square error
1.51
760.08
1000.38
1139.67
773.06
1638.87
0.29
0.29
1780.57
1266.35
603.14
Mean of response
0.01
30262
14749
15716
33634
22431
0.18
0.01
19408
14958
28103Based on different “errors” estimated by analysis of variance, SI-9 represents the most straightforward linear regression model. Figure3 shows the correlation between the estimated soil salinity values using soil salinity index-9 and the collected ground truth values.Figure 3
Linear regression between measured and estimated soil salinity values.Due to the fact that the soil salinity indices used in the current research are varied (11 in total), categorization of different indices using principal component analysis will help decide which is similar/different to which. Figure4 represents the grouping to different indices according to principal component analysis on covariances.Figure 4
Variance based principal component analysis.Generally, different soil salinity indices fell into two groups. SI-6, SI-9, SI-3, and BI are in one group and the rest are in the other group. SI-6 and NDSI failed to join their corresponding group and remain as group outlier [39]. Lack of correlation is the main reason of the insignificance of SI-6 and NDSI. High correlation significance between SI-3, SI-9, and BI is the best to represent soil salinity indices according to Table 2.The results proved the acceleration in the rate of vegetation death and soil salinization based on the high correlation between higher soil salinity values and negative normalized difference vegetation index withR2 values of 0.92 [40]. Moreover, higher soil salinity was related to the improper and intense drainage systems set up in addition to the lack of an operative water resource management plan [41].The dynamicity of the soil salinization process added further complications to designating salt-affected soils in a systematic uniform perspective. The use of different algorithms based on implementing different combinations and/or ratios of OLI-8 bands in the form of soil salinity indices evidenced to be more efficient to overcome soil dynamicity problems [23, 42].The selection of the sensitive bands adequate for proper soil salinity mapping is not systematically generalized [42, 47]. Spatial variability and land cover variations are the most limiting factors of band selection of a sensor [23, 48]. Therefore, some of the used equations may have better results in specific areas implementing specific band ratios than others taking into consideration the use of OLI-8 as a source of remote sensing data [23, 49].Estimated values of soil salinity index-9 were then used to map and classify soil salinity in the whole study area. Regular interpolation technique with six different smoothing Kernel functions is implemented and evaluated according to Krivoruchko et al. [50] as shown in Table 3 and Figure 5.Table 3
Evaluation of different interpolation kernels function.
Kernel function
Exponential
Gaussian
Quartic
Epanechnikov
Polynomial-5
Constant
Regression function
0.006
X
+
29339
-
0.001
X
+
29659
0.003
X
+
29477
-
0.005
X
+
29778
0.005
+
29368
-
0.001
X
+
29679
Mean
1383
1421
1378
1412
1355
1428
Root mean square
4251
4257
4250
4257
4253
4266
Mean standardized
0.0065
0.0067
0.0065
0.0066
0.0064
0.0067
Root mean square standardized
0.0201
0.0201
0.0202
0.0201
0.0202
0.0200
Average standard error
211114
211797
210334
211526
209641
212674Figure 5
Soil Salinity classes in Wadi Al Dawasir area.In Figure5, only 3% of Wadi Al Dawasir area is categorized to be with low soil salinity followed by 10% of moderate soil salinity class. Meanwhile, more than half of the study area is subjected to a higher concentration of soil salinity (52%). In addition, there are soils in Wadi Al Dawasir which cannot be used for the current agricultural practice due to extreme condition of soil salinity (35%).Kernel function of polynomial-5 was used to produce a thematic soil salinity map in the study area (Figure6). Based on Soil Survey approach in conjunction with natural breaks classification method was used to reclassify the prediction risk maps into four classes according to their toxicity levels into Low, Moderate, High, and Extreme [51, 52]. The role of the classifier is based on squared error minimization from a class’s means by creating internally homogenous groups but maintaining heterogeneity between classes.Figure 6
Soil salinity map in Wadi Al Dawasir.Figure7 shows the percentages of Wadi Al Dawasir agricultural area which were fallen into the four soil salinity classes based on NDVI estimation. The total agricultural area in Wadi Al Dawasir expressed low soil salinity contribution (6.5%). Those areas are the best to preserve the natural resources and the agricultural practice in Wadi Al Dawasir. About 36% of the total agricultural area constantly needs efforts to sustain the current agriculture practice. Both of the aforementioned classes are less than the half of the cropped area of Wadi Al Dawasir. Closer to 58% of the cultivated area is on improper soils due to its high salinity content.Figure 7
Contribution of different soil salinity classes to the total agricultural area.Soil salinity index-9, as a remotely sensed indicator, typically depends on the vegetation growing conditions which are helpful for mapping spatially the distribution of soil salinity [53]. The response of certain features of crop stress, dehydration, and senescence using vegetation indices have been examined to determine their response in different spectra [54, 55]. According to Carter [54], cell structure damage and chlorophyll reduction could be reliably measured using reduced near infrared reflectance (NIR) and increased visible reflectance (VIS) among several plant species in response to stress. Tilley et al. [53] also investigated salt stress responses in both of VIS and NIR. Normalized difference vegetation index (NDVI) as one of the most implemented vegetation indices has been inspected in abundant scholarly work by using the different vegetation reflectance to estimate soil salinity [40].


## 4. Conclusions

Crop production in Wadi Al Dawasir is completely dependent on underground water resources represented in artesian wells. The use of the underground water resources is annually overexploited due to the horizontal expansion of the agricultural land. The driving forces of the agricultural land expansion and the use of the underground water resources are controlled only by the Ministry of Water and Electricity but not by the local inhabitants. Such uncontrolled conditions put the water resources and the soils of Wadi Al Dawasir in jeopardy. The use of excess irrigational water to cope with the extremely arid conditions in the study area continuously creates more drainage problems as the agricultural areas in Wadi Al Dawasir have no drainage system by any means. Water logging and soil salinization are common problems in agricultural areas described by no drainage system. Consequently, the hydrological sceneries and the geomorphological patterns of arid lands are very important aspects to be considered in any rational development plan. Soil salinity indices are site/condition specific and soil salinity index-9 is the best to represent soil salinity values in the current study area. The current agricultural practices in Wadi Al Dawasir are obviously not the optimum water and/or soil resources management. More than half of the study area is characterized by either high or extreme soil salinity classes. Crop reallocation and rational water use may serve the goal of sustainability to keep the current agricultural practice in Wadi Al Dawasir. Monitoring and temporal soil salinity investigation are the keystone exercise required for adequate natural resources management plans in the near future.



---

*Source: 7596175-2016-01-06.xml*