Received 20 May 2016; accepted 5 August 2016; published 8 August 2016
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  . 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  . Exchangeable Sodium Percentage (ESP) index determines structural stability and class of soil sodification in a vast area  . Coke and the Muller  found out that ESP alone is not a satisfactory indicator for soil behavior.
According to Heidari  and Pazira  , many scholars have studied soil desalination. However, soil desodification was not studied by many scholars  . Mohsenifar et al.  studied soil desalination and desodification in Southeast region of Khuzestan Province.
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.
2. 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 TypicCalciorthids 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-so- dicity 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
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 (ESPi) and (ESPf) 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).
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:
ESPi, ESPf are exchangeable sodium percentage of soil saturation extract before and after leaching, Dlw the net depth of water for leaching, and Ds the depth of the soil layer (m). Actually, Dlw represents the depth of water that, after offsetting the moisture deficit in the related layer, leaves it by gravity. Deducting ESPeq 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.
3. Results and Discussion
According to analyses as well as fitting four mathematical models i.e. logarithmic, power, exponential and inverse 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:
Replacement of related variables results in:
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:
Replacement of related variables results in:
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 (ESPf) and required net depth of water for reclamation (D1w) 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 1. Soil desodification curve in Deylam soil series (treatment 1).
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.
 Shi, Z., Cheng, J.L., Huang, M.X. and Zhou, L.Q. (2006) Assessing Reclamation Levels of Coastal Saline Lands with Integrated Stepwise Discriminant Analysis and Laboratory Hyperspectral Data. Journal of Pedosphere, 16, 154-160.
 Cook, C.D. and Muller, W.J. (1997) Is Exchangeable Sodium Content a Better Index of Soil Sodicity than Exchangeable Sodium Percentage: A Reassessment of Published Data. Journal of Soil Science, 162, 342-349.
 Haidari, N. (1994) Studing of Leaching Models and Determine Leaching Efficiency at Saline and Sodic Soils. M.Sc. Thesis, Department of Irrigation and Reclamation Engineering, Faculty of Agriculture, University of Tehran, Karaj, 125 p. (In Persian)