A Practical Model for Desodification of Saline-Sodic Soils of Central Khuzestan Plains , Khuzestan Province

The most important task in leaching practices is assessment of water quantity required for leaching of saline and saline-sodic soils. Therefore, reliable estimation of the required leaching water quantity is vital for reducing soil salinity to a desirable level. The present study aimed to investigate desodification of saline and sodic soils in central area of Khuzestan Province. Consequently, a large area of 3216 ha with S4A3 salinity/sodicity class in Khuzestan, Iran, was selected to obtain the required data. This experiment was conducted with two treatments and tree replicates. In the first treatment, the experiment was conducted by applying just 100 cm water depth in four 25 cm intervals. In the second treatment, 5000 kg/ha Sulfuric Acid was applied prior to salt leaching together with leaching water. The intermittent ponding method was conducted with double rings in a rectangular array. The required physical and chemical analyses were performed on the collected data. The leaching water was supplied from Shotait River. Four mathematical models were applied to the collected experimental data to derive a suitable empirical model. The results for large scale applications indicated that the proposed logarithmic model can estimate the capital leaching requirement much than the previously proposed models.


Introduction
Unstable and degrading environmental issues are major challenges in application of natural resources (water and soil) in recent decades.Soil salinization and sodification are important processes of land degradation, especially in arid and semi-arid areas where the soil is modified by leaching of water soluble salts using high quality water along with addition of amendments as needed.Khuzestan Plain with an area of about 3 million hectares and access to a third of surface water resources in Iran is one of the most important poles of irrigated agriculture.Highly salinized and alkaline soils can be found in these areas due to improper use of soil and water resources and uncontrolled irrigation.These factors are major barriers to agriculture in this region.
Various studies have shown that increased number of measurement points reduced uncertainty of model inputs [1].The results of physicochemical analyses of soils with electrical conductivity values and low sand and low organic matter have shown that soil modification is time-consuming [2].Exchangeable Sodium Percentage (ESP) index determines structural stability and class of soil sodification in a vast area [3].Coke and the Muller [4] found out that ESP alone is not a satisfactory indicator for soil behavior.
The objectives of this study were to introduce an empirical model to account for reclamation water and to compare the obtained results with some available models.

Materials and Methods
The studied region is a part of Shoaybiyeh plains located in Khuzestan Province with an area of 15,800 ha.Since about 13,295 ha of all studied area have saline-sodic soil, from low to extreme, it can be argued that salinity and sodicity is the main qualitative limitation of more than 84% of the region's soils.
In order to perform required experiments, salinity and sodicity maps were studied and Deylam soil in Typic-Calciorthids group and Aridisals order with S4A2 salinity-sodicity class (extreme salinity and relative extreme sodicity) before leaching tests was selected as the target soil.In the study zone, the area of Deylam soil series is 6205 ha (more than 39% of the region's soil).Table 1 and Table 2 present the situation, series, and salinity-sodicity classes of the soils in the study area, and some physical characteristics of the various layers in the soil profiles before leaching.
To desotificate soils, leaching of soluble salts from soil profiles was carried out.The experiment had two treatments and was conducted with three replications.Leaching was performed using the intermittent ponding method and the required water was taken from the Shotait River.Table 3 lists results obtained from the chemical analysis of the water used in both treatments.
In the first treatment, no amendment material was used and the 1 meter depth of water was applied in four   0.25-meter intervals.However, in the second treatment, five tons of 95% sulfuric acid was applied and leaching was then carried out by using irrigation water.
The six double cylinders used in the experiment were arranged 5 meters apart on the perimeter of a circle with a radius of 5 meters.Soil samples were taken from depths of 0 -25, 25 -50, 50 -75, 75 -100, 100 -125, and 125 -150 centimeters of the soil profiles before leaching and after application of 25, 50, 75, and 100 centimeters of water.The samples were sent to the laboratory for analysis.In each experiment, the electrical conductivity, pH, cation exchange capacity, sodium adsorption ratio, exchangeable sodium percentage, saturation percentage, lime, gypsum, and cations and anions (Na, Ca, Mg, Cl, sulfate, carbonate, and bicarbonate) in the soil extract solution were measured.The salinity values (electrical conductivity of the soil saturation extract) before, during, and after each water application were determined for the desired horizons in the soil profiles (that is, at the depths of 0 -25, 0 -50, 0 -75, 0 -100, 0 -125, and 0 -150 centimeters) relative to the mean weight calculations for different depths of water applied.The results for both treatments are listed in Table 4.
Since all the water applied may not be used for leaching soluble salts out of soil profiles and some if it may be used to offset the soil moisture deficit, even applying large volumes of water will not result in a chemical equilibrium of the soil with the water used in leaching.To overcome this problem, desalinization values were determined based on the averages of salinity by weight as follows: ESP i , ESP f are exchangeable sodium percentage of soil saturation extract before and after leaching, D lw the net depth of water for leaching, and D s the depth of the soil layer (m).Actually, D lw represents the depth of water that, after offsetting the moisture deficit in the related layer, leaves it by gravity.Deducting ESP eq from the numerator and denominator of the mentioned fraction will cause the results obtained from the effects of external factors such as the amount of evaporation, the condition of the internal drainage in the soil, the quality of water used for leaching, and the conditions under which the experiment is conducted, to become independent.In fact, this will convert the function from the explicit to the implicit state.After obtaining all the values of the leaching experiments, the required analysis was performed using SPSS, Curve Expert, and Excel.The four mathematical models including the power, exponential, inverse, and logarithmic models were fitted to the desodification values, were analyzed using statistical criteria such as coefficient of determination and standard deviation at the one percent level of significance, and the most suitable Desodification model for the tested soils was determined.

Results and Discussion
According to analyses as well as fitting four mathematical models i.e. logarithmic, power, exponential and in-verse models to the values extracted from X and Y variables, derived from field Desodification experiments in Deylam soil series, exponential model with a determination coefficient of 0.73 and standard error of 0.115 in a significance level of 1% is the best model for treatment 1 which is shown as follows: 0.135Ln 0.222 Replacement of related variables results in: 0.135Ln 0.222 Having relation (), required depth for removing soluble salts from soil profile with a given thickness and a given final salinity of saturated soil extract is computed as follows: ( ) According to analyses as well as fitting different empirical models, logarithmic model with a definition coefficient of 0.65 and standard error of 0.113 in a significance level of 1% is the best model for treatment 2 which is shown as following: 0.110Ln 0.210 Replacement of related variables results in: 0.110Ln 0.222 According to Figure 1 and Figure 2, in the soil series removing salts is easier in treatment 1 compared with treatment 2. Comparison of the curves show that in the case of applying sulfuric acid as reclamation substance, more leaching water is required compared with the case of no use of sulfuric acid.
From the curves shown in Figure 1 and Figure 2, final exchangeable sodium percentage of soil (ESP f ) and required net depth of water for reclamation (D 1w ) can be estimated.It should be noted that if one wish to estimate total required depth of leaching water, soil water deficit in the considered soil layer, surface evaporation and precipitation amount should be taken into account in the calculations and planning of leaching process.

Figure 2 .
Figure 2. Soil desodification curve in Deylam soil series (treatment 2).However, the curves are valid only for the region's soil and within initial electrical conductivity limits of 35.55 dS/m to 94.46 dS/m and exchangeable sodium percentage (ESP) of 26.81 to 31.36.

Table 1 .
Series, salinity and sodicity class and some physical characteristics of the studied soil.

Table 2 .
Some physical characteristics of different layers of soil profile before leaching process in Deylam soil series.

Table 3 .
Chemical decomposition of the water used to leach out soluble salts from soil profile.

Table 4 .
Calculated values (weight average) for (ESP i ) and (ESP f ) of saturated soil extract before and after applying different amounts of leaching water (gross) in different soil layers of Deylam soil series (with and without reclamation substances).