Monitoring and Assessment of Salinity and Chemicals in Agricultural Lands by a Remote Sensing Technique and Soil Moisture with Chemical Index Models

: Agricultural land in the south of Iraq provides habitat for several types of living creatures. This land has a signiﬁcant impact on the ecosystem. The agricultural land of Al-Hawizeh marsh covers Cr, and Zn) in one image. This approach was repeated for each season in this study. In conclusion, the developed systematic and generic approach may constitute a basis for determining soil quality parameters in agricultural land worldwide.

GIS and multi-temporal remote sensing data, such as Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper (ETM) satellite images, have been used to analyze vegetation stress to develop soil salinity maps and determine the content of total dissolved solids (TDS), total suspended solids (TSS), Ca, Mg, Na, K, CO 3 , H CO 3 , Cl, NO 3 , and SO 4 in irrigated areas of Syrdarya province, Uzbekistan, with a higher degree of spatial accuracy and in a cost effective way [36]. Another study estimated soil salinity and conducted salinity mapping with considerable accuracy using multi-spectral remote sensing data with an advanced land imager (ALI) sensor and partial least-square regression model (PLSR), together with soil salinity data collected from the Yellow River Delta in China [37]. Another study used different satellite sensors (e.g., TM, MSS, ETM+, SPOT, Terra-ASTER, LISS-III, and IKONOS) and methods used to detect, monitor and map of soil salinity in arid and semi-arid regions [38]. Other studies assessed the metal contamination (As, Cr, Cd, Co, Fe, Cu, Zn, Hg, Pb, and Mn) of sediments in the Shadegan and Hawr Al-Azim wetlands in southwestern Iran based on GIS and remote sensing with mathematical models [39][40][41].
This article critically evaluates the existing literature that indicate more studies are needed to bridge the gaps identified in previous studies. A comprehensive survey of the literature revealed that only a few techniques have been developed for the detection and monitoring of physical and chemical soil quality parameters in the marshland area of Iraq. More systematic studies are required to support findings from previous studies and improve the accuracy of findings on the soil of agricultural land. The aim of this research was to examine whether the integration between mathematical models, remote sensing Landsat-8 OLI data and GIS could provide a powerful tool for predicting, assessing, monitoring, managing, and mapping salinity and chemical parameters such as iron (Fe), lead (Pb), chromium (Cr), and copper (Cu) in the soils of agricultural land in the Al-Hawizeh marsh in southern Iraq during the four seasons of 2017.

Study Area
The study area was the agricultural land that surrounds Al-Hawizeh marsh in southern Iraq, covering an area of 3500-4000 km 2 , as shown in Figure 1. It is considered an enriched resource for the production of several crops that have a significant impact on the ecosystem, as well as a habitat for biodiversity [13,14,18]. The neglect and wars of recent years are factors that have led to an increase of salinity and heavy chemical pollution in this area. Al-Hawizeh marsh is the major water source for these lands [15,16,18,20].
Geosciences 2020, 10, x FOR PEER REVIEW 3 of 20 suspended solids (TSS), Ca, Mg, Na, K, CO 3 , H CO 3 , Cl, NO 3 , and SO 4 in irrigated areas of Syrdarya province, Uzbekistan, with a higher degree of spatial accuracy and in a cost effective way [36].
Another study estimated soil salinity and conducted salinity mapping with considerable accuracy using multi-spectral remote sensing data with an advanced land imager (ALI) sensor and partial least-square regression model (PLSR), together with soil salinity data collected from the Yellow River Delta in China [37]. Another study used different satellite sensors (e.g., TM, MSS, ETM+, SPOT, Terra-ASTER, LISS-III, and IKONOS) and methods used to detect, monitor and map of soil salinity in arid and semi-arid regions [38]. Other studies assessed the metal contamination (As, Cr, Cd, Co, Fe, Cu, Zn, Hg, Pb, and Mn) of sediments in the Shadegan and Hawr Al-Azim wetlands in southwestern Iran based on GIS and remote sensing with mathematical models [39][40][41]. This article critically evaluates the existing literature that indicate more studies are needed to bridge the gaps identified in previous studies. A comprehensive survey of the literature revealed that only a few techniques have been developed for the detection and monitoring of physical and chemical soil quality parameters in the marshland area of Iraq. More systematic studies are required to support findings from previous studies and improve the accuracy of findings on the soil of agricultural land. The aim of this research was to examine whether the integration between mathematical models, remote sensing Landsat-8 OLI data and GIS could provide a powerful tool for predicting, assessing, monitoring, managing, and mapping salinity and chemical parameters such as iron (Fe), lead (Pb), chromium (Cr), and copper (Cu) in the soils of agricultural land in the Al-Hawizeh marsh in southern Iraq during the four seasons of 2017.

Study Area
The study area was the agricultural land that surrounds Al-Hawizeh marsh in southern Iraq, covering an area of 3500-4000 km 2 , as shown in Figure 1. It is considered an enriched resource for the production of several crops that have a significant impact on the ecosystem, as well as a habitat for biodiversity [13,14,18]. The neglect and wars of recent years are factors that have led to an increase of salinity and heavy chemical pollution in this area. Al-Hawizeh marsh is the major water source for these lands [15,16,18,20].

Methodology
The methodology adopted for this article is presented as a flowchart in Figure 2. It includes details regarding satellite images, image pre-processing, soil samples, soil analysis, mathematical models, data validation, image classification, and data visualization. Landsat-8 satellite images are represented as basic remote sensing data collected during 2017. Processed images from Landsat-8 satellite data and ancillary data were synergistically integrated with the mathematical models to predict the salinity and concentrations of chemicals (Fe, Pb, Cu, Cr, and Zn) in the soil for the agricultural land of the study area. The creation of mathematical models was informed by the nature of the satellite image processing and auxiliary data. The mathematical models had two equations: salinity equation (SE) and chemical equation (CE). The SE was based on the soil moisture index (SMI), which depends on the B6 and B11 bands. The CE was based on the clay chemical indices (CCIs), which depends on the B6 and B7 bands. The outputs of both salinity and chemical equations were the salinity and chemical concentration data obtained from the spatial distribution of satellite images. These data were categorized using a decision tree (DT) of image classification to classify the salinity and chemical parameters on satellite images over four seasons. Data visualization was used to check the salinity, Fe, Pb, Cu, Cr, and Zn classification results on the maps with the appropriate specifications as output results of data validations.
depends on the B6 and B7 bands. The outputs of both salinity and chemical equations were the salinity and chemical concentration data obtained from the spatial distribution of satellite images. These data were categorized using a decision tree (DT) of image classification to classify the salinity and chemical parameters on satellite images over four seasons. Data visualization was used to check the salinity, Fe, Pb, Cu, Cr, and Zn classification results on the maps with the appropriate specifications as output results of data validations.

Satellite Images and Pre-Processsing
The remote sensing data used satellite images from Landsat-8 OLI at longitude 47 • 32 -47 • 45 E and latitude 31 • 30 -31 • 42 N with a resolution of 30 m. The satellite images of the study area were synchronized with soil samples collected from the agricultural land during 2017. The mathematical models of this study were based on the B6, B7 and B11 bands of the satellite image, as shown in Table 1. Accurate pre-processing of Landsat-8 satellite images, including atmospheric correction, is the first important step of the research methodology. This process involves radiometric calibration, dark subtraction, and internal average reflectance (IAR) calibration, layer stacking, georeferencing, image enhancement, and region of interest (ROI) determination. The output of this process provides suitable images. These data were combined with mathematical models for assessing and monitoring the spatial distribution of salinity and chemicals in the study area.

Soil Sample Collection
Soil samples were collected at 15 cm depth using a clean shovel. These samples were collected from 60 stations positioned across different regions, covering the whole of the agricultural land by using the Global Positioning System (GPS). These soil sampling stations were numbered S1-S60, as illustrated in Figure 3. All soil samples underwent laboratory testing using the atomic absorption spectrophotometer (AAS) technique to measure the concentration values for salinity and chemicals (Fe, Pb, Cu, Cr, and Zn). These data were used to determine the status of soil quality.
All data were tested using the atomic absorption spectrophotometer (AAS) technique. AAS is a technique used to measure trace elements present in soil samples, in which free gaseous atoms absorb electromagnetic radiation at a specific wavelength to produce a measurable signal. The absorption
All data were tested using the atomic absorption spectrophotometer (AAS) technique. AAS is a technique used to measure trace elements present in soil samples, in which free gaseous atoms absorb electromagnetic radiation at a specific wavelength to produce a measurable signal. The absorption signal is proportional to the concentration of those free absorbing atoms in the optical path. Therefore, for AAS measurements, the sample must be converted into gaseous atoms, usually by the application of heat to a cell using an atomizer. The main type of atomizer used in AAS-based analytical techniques is the flame AAS (FAAS), which provides analytical signals in a continuous fashion [48,49].
The samples were grinded in a motor and sieved. Using 1 mm, the sample (1 g) was digested with concentrated nitric acid in a sand bath until dry. Then, deionized water was added and the sample was filtered then complete the filtrate to 25 mL with deionizer water. AAS was used to determine the chemical concentrations at specific wavelengths along with a hallow cathode lamp (HCL). Chemical speciation was determined by AAS after extraction of the metal with different chemical reagents.
The soil analysis results for the salinity and chemicals represent the real data, demonstrating the status of the soil quality in the agricultural land. The concentration values for the salinity and the chemicals during four seasons were lowest in the summer and highest in the winter, as shown in Figures 4 and 5. The minimum, maximum, and average of salinity and chemicals during the four seasons are shown in Table 2. The mathematical models used these data to combine with satellite images and assess the salinity and chemical values for the four seasons of 2017.
Geosciences 2020, 10, x FOR PEER REVIEW 6 of 20 seasons are shown in Table 2. The mathematical models used these data to combine with satellite images and assess the salinity and chemical values for the four seasons of 2017.

Mathematical Models
The developed mathematical model can retrieve the soil quality parameters (salinity, Fe, Pb, Cu, Cr, and Zn) from Landsat-8 data. Soil quality parameter estimations are based on the integration between ground data, satellite images after processing and mathematical models. The developed model is based on the soil moisture index (SMI) and clay chemical indices (CCIs). The input of this model is based on the ground data and post-processed satellite image. These models were introduced by [38,50].
The model consists of four parts: The first one depends on the B6 and B11 bands of Landsat-8 images, which calculates the SMI; the second part is the salinity equation (SE), which depends on the output of SMI to retrieve the salinity values from Landsat-8 images. The third part depends on the B6 and B7 bands of Landsat-8 images, which calculates the CCI. The fourth part is the chemical equation (CE), which depends on the output of CCI to retrieve the chemical values (Fe, Pb, Cu, Cr, and Zn) from Landsat-8 images. These four components are described below.

Soil Moisture Index (SMI)
The SMI equation was used to calculate soil moisture, depending on the relation between the B6 and B11 bands of the Landsat-8 images, as shown in Equation (1). The B6 (short wave infrared (SWIR)) band discriminates vegetation and moisture content of the soil by penetrating thin clouds with a wavelength range of 1.57-1.65 μm. Meanwhile, the B11 (TIRS-2) band estimates the soil moisture and the thermal mapping in a wavelength range of 11.5-12.51 μm [42,43,51].

Salinity Equation (SE)
The SE calculates the salinity value from the Landsat-8 images, depending on the relation between the output values of the SMI from the Landsat-8 images and the salinity values of the ground data, as illustrated in Equations (2) and (3), and shown in Figure 6, where the y-axis represents the salinity values from the ground measurements while the x-axis represents the values of the SMI from the satellite images.

Mathematical Models
The developed mathematical model can retrieve the soil quality parameters (salinity, Fe, Pb, Cu, Cr, and Zn) from Landsat-8 data. Soil quality parameter estimations are based on the integration between ground data, satellite images after processing and mathematical models. The developed model is based on the soil moisture index (SMI) and clay chemical indices (CCIs). The input of this model is based on the ground data and post-processed satellite image. These models were introduced by [38,50].
The model consists of four parts: The first one depends on the B6 and B11 bands of Landsat-8 images, which calculates the SMI; the second part is the salinity equation (SE), which depends on the output of SMI to retrieve the salinity values from Landsat-8 images. The third part depends on the B6 and B7 bands of Landsat-8 images, which calculates the CCI. The fourth part is the chemical equation (CE), which depends on the output of CCI to retrieve the chemical values (Fe, Pb, Cu, Cr, and Zn) from Landsat-8 images. These four components are described below.

Soil Moisture Index (SMI)
The SMI equation was used to calculate soil moisture, depending on the relation between the B6 and B11 bands of the Landsat-8 images, as shown in Equation (1). The B6 (short wave infrared (SWIR)) band discriminates vegetation and moisture content of the soil by penetrating thin clouds with a wavelength range of 1.57-1.65 µm. Meanwhile, the B11 (TIRS-2) band estimates the soil moisture and the thermal mapping in a wavelength range of 11.5-12.51 µm [42,43,51].

Salinity Equation (SE)
The SE calculates the salinity value from the Landsat-8 images, depending on the relation between the output values of the SMI from the Landsat-8 images and the salinity values of the ground data, as illustrated in Equations (2) and (3), and shown in Figure 6, where the y-axis represents the salinity values from the ground measurements while the x-axis represents the values of the SMI from the satellite images.
Geosciences 2020, 10, x FOR PEER REVIEW 8 of 20 Here, SE represents the required values of salinity from the satellite images; SMI represents the values of the SMI from the satellite images, based on the output results of Equation (1); and f, g, and i are constants that depend on both the values of the SMI from the satellite images and the salinity values of the ground data. The base e of the exponential function is equal to 2.718.  Here, SE represents the required values of salinity from the satellite images; SMI represents the values of the SMI from the satellite images, based on the output results of Equation (1); and f, g, and i are constants that depend on both the values of the SMI from the satellite images and the salinity values of the ground data. The base e of the exponential function is equal to 2.718.

Clay Chemical Indices (CCIs)
Clay chemical indices compose all of the component chemicals-such as Pb, Fe, Cu, Cr, and Zn-that are found in the clay soil of agricultural land. CCIs calculate the content of clay chemicals in the soil based on the relationship between bands B6 and B7 of the Landsat-8 images, as shown in Equation (4). The B6 (SWIR-1) band discriminates the moisture content of soil and vegetation by penetrating thin clouds with a wavelength ranging from 1.57 to 1.65 µm. The B7 (SWIR-2) band improves the moisture content detection of soil and vegetation through thin cloud penetration with a wavelength ranging from 2.11 to 2.29 µm [42,43,51].

Chemical Equation (CE)
The CE retrieves values of chemicals (Pb, Fe, Cu, Cr, and Zn) based on the relationship between the output results of the CCIs from the Landsat-8 images and the ground data, as illustrated in Equations (5) and (6) and shown in Figure 7, where the y-axis represents the chemical values of the ground measurements and the x-axis represents the CCIs values from the satellite images.
where CE represents the concentration values of the chemicals (Pb, Fe, Cu, Cr, and Zn) acquired from the satellite images; CCIs represents the clay chemical indices from the satellite images based on the output results of Equation (4); p and h are constants that depend on the chemical values of the ground data and the indices of the clay chemicals from the satellite images, respectively; and the power n is an integer. All of the previous steps of the model must be repeated for each season of study.
Geosciences 2020, 10, x FOR PEER REVIEW 9 of 20 where CE represents the concentration values of the chemicals (Pb, Fe, Cu, Cr, and Zn) acquired from the satellite images; CCIs represents the clay chemical indices from the satellite images based on the output results of Equation (4); p and h are constants that depend on the chemical values of the ground data and the indices of the clay chemicals from the satellite images, respectively; and the power n is an integer. All of the previous steps of the model must be repeated for each season of study.

Image Classification
Several image classification techniques can be used to classify the satellite images, such as maximum likelihood (ML), minimum distance (MD), artificial neural network (ANN), support vector machine (SVM), and decision tree (DT) [4,24]. The present study classified soil quality parameters using images of agricultural land in Al-Hawizeh marsh during four seasons of 2017 depending on

Image Classification
Several image classification techniques can be used to classify the satellite images, such as maximum likelihood (ML), minimum distance (MD), artificial neural network (ANN), support vector machine (SVM), and decision tree (DT) [4,24]. The present study classified soil quality parameters using images of agricultural land in Al-Hawizeh marsh during four seasons of 2017 depending on the output results of SE and CE equations. DT classification was used in this study. It is an advanced approach for image classification, which depends on the layered or stratified approach to solve the problems of distinguishing between spectral classes [4,24]. DT uses data from many sources to make a single decision tree classifier. This performs multistage classification using a series of binary decisions to place pixels into classes. Input data can be obtained from various sources and data types. The results of the decisions are classes that can be saved as trees to apply them for other datasets [52]. This classification refers to all parameters simultaneously in one stage. The output of DT classification results can display all soil quality parameters in one image. Thus, DT in this case is considered as the optimal approach. This approach was repeated for each season in this study.

Data Visualization
The output results are presented as DT classifications for both SE and CE equations after performing validation. The output of the DT results are presented as one image for salinity, Fe, Pb, Cu, Cr, and Zn. This procedure was performed for each season.

Results and Discussion
The mathematical models were combined with the Landsat-8 satellite images and synchronized with the ground data to retrieve soil quality parameters for agricultural land at Al-Hawizeh marsh. The values of soil quality parameters (salinity, Fe, Pb, Cu, Cr, and Zn) are evaluated from Landsat-8 data based on the mathematical models. These values indicate the status of soil quality for agricultural lands of Al-Hawizeh marsh from the satellite images during four seasons of 2017. Table 3 summarizes the evaluated soil quality parameters.  Figure 8 shows the seasonal variation of soil salinity obtained using the SMI model, which is observed to be lowest during summer and highest during winter. The average salinity concentrations during autumn, summer, spring, and winter are 1175, 1010, 1105, and 1789 mg/dm 3 , respectively, determined by the soil moisture index. This indicated a reduction of SMI in the summer and rise in the winter. The SMI values are decided by the spectral signatures of Landsat-8 Short Wave Infrared (SWIR-1) band B6 and Thermal Infrared Sensor (TIRS-2) band B11. These images (bands) are considered as suitable and sensitive as far as the soil moisture indexes are concerned. Salinity concentrations displayed a positive relationship with the SMI values [18,42,43,51].  Figure 9 shows the seasonal variation of iron content obtained using the CCIs model, which is observed to be lowest during summer and highest during winter. The average Fe concentrations during autumn, summer, spring, and winter are 813, 784, 842, and 1106 mg/dm 3 , respectively, which was determined by the clay chemical indices. This indicated a reduction of the clay chemicals in the summer and a rise during winter. The CCIs values are decided by the spectral signatures of Landsat-8 Short Wave Infrared (SWIR-1) band B6 having and Short Wave Infrared (SWIR-2) band B7 with wavelength range 2.11-2.29 μm. These images (bands) are considered as suitable and sensitive as far as the CCIs are concerned. Fe concentrations displayed a positive relationship among the clay chemicals indices values [18,42,43,51] Figure 9 shows the seasonal variation of iron content obtained using the CCIs model, which is observed to be lowest during summer and highest during winter. The average Fe concentrations during autumn, summer, spring, and winter are 813, 784, 842, and 1106 mg/dm 3 , respectively, which was determined by the clay chemical indices. This indicated a reduction of the clay chemicals in the summer and a rise during winter. The CCIs values are decided by the spectral signatures of Landsat-8 Short Wave Infrared (SWIR-1) band B6 having and Short Wave Infrared (SWIR-2) band B7 with wavelength range 2.11-2.29 µm. These images (bands) are considered as suitable and sensitive as far as the CCIs are concerned. Fe concentrations displayed a positive relationship among the clay chemicals indices values [18,42,43,51].

Iron Values
8 Short Wave Infrared (SWIR-1) band B6 having and Short Wave Infrared (SWIR-2) band B7 with wavelength range 2.11-2.29 μm. These images (bands) are considered as suitable and sensitive as far as the CCIs are concerned. Fe concentrations displayed a positive relationship among the clay chemicals indices values [18,42,43,51].

Lead Values
The Pb concentrations that were evaluated using CCIs model during the four seasons were lowest in summer and highest in winter ( Figure 10). The average Pb concentrations during autumn, summer, spring, and winter were 4.85, 3.79, 4.74, and 7.2 mg/dm 3 , respectively, determined by CCIs. This indicated a reduction the indices of the clay chemicals in the summer and a rise during winter. Pb concentrations displayed a positive relationship with the CCIs values [18,42,43,51].

Lead Values
The Pb concentrations that were evaluated using CCIs model during the four seasons were lowest in summer and highest in winter ( Figure 10). The average Pb concentrations during autumn, summer, spring, and winter were 4.85, 3.79, 4.74, and 7.2 mg/dm 3 , respectively, determined by CCIs. This indicated a reduction the indices of the clay chemicals in the summer and a rise during winter. Pb concentrations displayed a positive relationship with the CCIs values [18,42,43,51].   Figure 11 shows the seasonal variation of Cu content obtained using the CCIs model, which was observed to be lowest during summer and highest during winter. The average Cu concentrations during autumn, summer, spring and winter were 3.9, 3.1, 4.45, and 7.5 mg/dm 3 , respectively, which was determined by the CCIs. This indicated a reduction in the CCIs in the summer and a rise in the winter. Cu concentrations displayed a positive relationship with CCIs values [18,42,43,51].   Figure 11 shows the seasonal variation of Cu content obtained using the CCIs model, which was observed to be lowest during summer and highest during winter. The average Cu concentrations during autumn, summer, spring and winter were 3.9, 3.1, 4.45, and 7.5 mg/dm 3 , respectively, which was determined by the CCIs. This indicated a reduction in the CCIs in the summer and a rise in the winter. Cu concentrations displayed a positive relationship with CCIs values [18,42,43,51]. Figure 11 shows the seasonal variation of Cu content obtained using the CCIs model, which was observed to be lowest during summer and highest during winter. The average Cu concentrations during autumn, summer, spring and winter were 3.9, 3.1, 4.45, and 7.5 mg/dm 3 , respectively, which was determined by the CCIs. This indicated a reduction in the CCIs in the summer and a rise in the winter. Cu concentrations displayed a positive relationship with CCIs values [18,42,43,51].

Chromium Values
The Cr concentrations were evaluated using CCIs and were lowest during summer and highest during winter (Figure 12). The average Cr concentrations during autumn, summer, spring, and winter were 1.28, 0.73, 1.03, and 2.91 mg/dm 3 , respectively, determined by CCIs. This indicated a reduction of CCIs in the summer and a rise in the winter. Cr concentrations displayed a positive relationship with the clay chemical index values [18,42,43,51].

Chromium Values
The Cr concentrations were evaluated using CCIs and were lowest during summer and highest during winter ( Figure 12). The average Cr concentrations during autumn, summer, spring, and winter were 1.28, 0.73, 1.03, and 2.91 mg/dm 3 , respectively, determined by CCIs. This indicated a reduction of CCIs in the summer and a rise in the winter. Cr concentrations displayed a positive relationship with the clay chemical index values [18,42,43,51].   Figure 13 demonstrates the seasonal variation of Zn contents obtained using CCIs model, which was observed to be lowest during summer and highest during winter. The average Zn concentrations during autumn, summer, spring, and winter were 8.25, 6, 7.05, and 12 mg/dm 3 , respectively, which is decided by the clay chemical indices. This indicated a reduction of CCIs in the summer and a rise in the winter. Zn concentrations displayed a positive relationship with CCIs values [18,42,43,51].   Figure 13 demonstrates the seasonal variation of Zn contents obtained using CCIs model, which was observed to be lowest during summer and highest during winter. The average Zn concentrations during autumn, summer, spring, and winter were 8.25, 6, 7.05, and 12 mg/dm 3 , respectively, which is decided by the clay chemical indices. This indicated a reduction of CCIs in the summer and a rise in the winter. Zn concentrations displayed a positive relationship with CCIs values [18,42,43,51]. was observed to be lowest during summer and highest during winter. The average Zn concentrations during autumn, summer, spring, and winter were 8.25, 6, 7.05, and 12 mg/dm 3 , respectively, which is decided by the clay chemical indices. This indicated a reduction of CCIs in the summer and a rise in the winter. Zn concentrations displayed a positive relationship with CCIs values [18,42,43,51].   Figure 14 presents the minimum and the maximum spatial distribution of salinity and chemical (Fe, Pb, Cu, Cr, and Zn) contents for soils in the agricultural lands of Al-Hawizah marsh obtained using the mathematical models during all four seasons. Figure 14A represents the spatial distribution with the minimum and maximum concentration values for salinity, Fe, Pb, Cu, Cr, and Zn in the autumn, which were determined to be (196-9235), (274-35,531), (0. , (0.42-26), (0.03-10.7), and (1.6-2.95) mg/dm 3 , respectively. Figure 14B represents the spatial distribution with the minimum and maximum concentration values for salinity, Fe, Pb, Cu, Cr, and Zn in the summer, which were (992-18,214), (56-4604), (0.47-9), (0.5-18.4), (0.08-4.3), and (4.9-32) mg/dm 3 , respectively. Figure 14C represents the spatial distribution of the minimum and maximum concentration of salinity, Fe, Pb, Cu, Cr, and Zn in the spring, which were (813-7944), (290-11234), (0.37-13), (0.05-26), (0.015-5.9), and (2.7-4.6) gm/dm 3 , respectively. Figure 14D represents the spatial distribution of the minimum and maximum salinity, Fe, Pb, Cu, Cr, and Zn in the winter, which were (83-4628), (1013-3818), (0.29-25), (0.26-40), (0.017-3.73), and (2.74-4028) mg/dm 3 , respectively. Figure 15 illustrates two-dimensional image layers for salinity and chemicals distribution during all four seasons based on the proposed models.

Image Classifications
The salinity values were observed to be lowest during the summer and highest during the winter, which were decided by the SMI, indicating a reduction of the soil moisture indexes during summer and a rise in the winter. The soil moisture index values were decided by the spectral signatures of the Landsat-8 Short Wave Infrared (SWIR-1) band B6 and the Thermal Infrared Sensor (TIRS-2) band B11, which are considered suitable and sensitive as far as SMI is concerned. These salinity concentrations displayed a positive relationship among the soil moisture index values. [42,43,51].
All chemical values were determined using the CCIs and indicated a reduction of the clay CCIs in the summer and a rise during winter. The CCIs were determined by the spectral signatures of the Landsat-8 SWIR-1 band B6 and SWIR-2 band B7, which are considered suitable and sensitive for this purpose. The chemical concentrations displayed a positive relationship with CCI values [19].

Conclusions
The agricultural land surrounding Al-Hawizeh marsh in southern Iraq is considered a habitat provider for biodiversity and an enriched resource for the production of several crops. The neglect and wars of recent years are factors that may have led to an increase of salinity, iron, lead, copper, chromium, and zinc in the soil. Techniques have been developed for the detection and monitoring of physical and chemical parameters in soils.
AAS is a technique that can be used to measure the salinity and chemicals in soils. The mathematical models of this study were based on the B6, B7, and B11 bands of Lasndsat-8 satellite imagery. These bands were suitable to combine with the mathematical models for assessing and monitoring the spatial distribution of salinity and chemicals from satellite imagery. It is concluded that the integration of mathematical models with remote sensing data and GIS can provide a powerful tool for predicting, assessing, monitoring, managing, and mapping salinity and chemicals components.
The seasonal variation of average salinity and chemicals was obtained using the SE and CE equations. The salinity values were evaluated by using the SE during all four seasons. The average salinity concentrations during autumn, summer, spring, and winter were 1175, 1010, 1105, and 1789 mg/dm 3 , respectively. These values were determined by the SMI. This indicated a reduction of the SMI values in the summer and a rise in the winter. These salinity concentrations displayed a positive relationship with SMI values.

Conclusions
The agricultural land surrounding Al-Hawizeh marsh in southern Iraq is considered a habitat provider for biodiversity and an enriched resource for the production of several crops. The neglect and wars of recent years are factors that may have led to an increase of salinity, iron, lead, copper, chromium, and zinc in the soil. Techniques have been developed for the detection and monitoring of physical and chemical parameters in soils.
AAS is a technique that can be used to measure the salinity and chemicals in soils. The mathematical models of this study were based on the B6, B7, and B11 bands of Lasndsat-8 satellite imagery. These bands were suitable to combine with the mathematical models for assessing and monitoring the spatial distribution of salinity and chemicals from satellite imagery. It is concluded that the integration of mathematical models with remote sensing data and GIS can provide a powerful tool for predicting, assessing, monitoring, managing, and mapping salinity and chemicals components.
The seasonal variation of average salinity and chemicals was obtained using the SE and CE equations. The salinity values were evaluated by using the SE during all four seasons. The average salinity concentrations during autumn, summer, spring, and winter were 1175, 1010, 1105, and 1789 mg/dm 3 , respectively. These values were determined by the SMI. This indicated a reduction of the SMI values in the summer and a rise in the winter. These salinity concentrations displayed a positive relationship with SMI values.
The chemical concentrations were obtained using the CE equation. These values were observed to be lowest during the summer and highest during the winter. The average values for (Fe) during the autumn, summer, spring, and winter were 813, 784, 842, and 1106 mg/dm 3 , respectively. The average values for (Pb) during the autumn, summer, spring, and winter were 4.85, 3.79, 4.74, and 7.2 mg/dm 3 , respectively. The average values for (Cu) during the autumn, summer, spring and winter were found to be 3.9, 3.1, 4.45, and 7.5 mg/dm 3 , respectively. The average values for (Cr) during the autumn, summer, spring, and winter were found to be 1.28, 0.73, 1.03, and 2.91 mg/dm 3 , respectively. The average values for (Zn) during the autumn, summer, spring, and winter were found to be 8.25, 6, 7.05, and 12 mg/dm 3 , respectively. These values were determined using the CCIs in the soil. This indicated a reduction of the CCIs in the summer and a rise in the winter season. These chemicals concentrations displayed a positive relationship with values of the CCIs.
To conclude on the salinity and chemical values obtained based on the SMI and CCIs results, the SMI and CCIs were dependent on the amount of water coming from Al-Hawizeh marsh during the four seasons. This indicated a reduction in water during summer and an increase in winter. This was exacerbated by the high temperatures and dry soil in the summer. The salinity and chemical concentrations were above the maximum allowable concentration in agricultural land based on Iraq's standards for agricultural land pollution.
The DT classification scheme provided good results for this study. The DT classification depended on the output results of salinity and chemicals for both the SE and CE equations. This classification refers to all the parameters simultaneously in one stage. The output of the DT classification can display all the soil quality parameters (salinity, Fe, Pb, Cu, Cr, and Zn) in one image. This approach was repeated for each season. In conclusion, the developed systematic and generic approach may constitute a basis for determining the soil quality parameters in agricultural land worldwide.