JWSS - Journal of Water and Soil Science

---

The aim of the study was to investigate the effects of land use on soil quality parameters using multivariate statistical analysis. Soil samples (0-25 and 25-50 cm depths) were taken from three land uses in forest area of Marivan including forest, rangeland, and cultivated land. Soil characteristics of pH, EC, sand, silt, clay and CaCO3 content, water-stable aggregates and their organic carbon content were measured. Principal component, cluster and discriminant analyses were used to evaluate the soil quality. Principal component analysis classified soil properties into five factors. The most important factors were soil aggregates organic carbon content and aggregate stability indices. Schematic distribution of factors and also cluster analysis showed the same pattern. Soil aggregates organic carbon content, water-stable aggregates and aggregate stability indices were the most sensitive factors to land use changes. These soil properties and factors had the same pattern in forest and rangeland, but significantly reduced in the cultivated land use. Land use change from forest to cultivated land resulted in significant decrease of aggregates organic carbon content, water-stable aggregates and also an increase of pH. The results showed the usefulness of multivariate statistical methods for integration of the soil properties and determination of different soil quality indices.

---

Soil aggregate stability is considered as a key indicator of soil quality and health assessments in rangelands. Many factors and properties such as soil texture, organic carbon, calcium carbonate, sodium adsorption ratio, and electrical conductivity might affect soil aggregate stability. The effects of these factors on aggregate stability of 71 soil samples collected from 4 rangeland sites (2 in semi-arid and 2 in arid lands) in Isfahan province were investigated. Aggregate stability was measured using the wet-sieving method. To optimize the trial conditions for the investigated soils, three shaking times (5, 10 and 15 minutes) were used to impose different hydromechanical stresses on the aggregates of ten soils selected out of the studied soils. The structural stability was assessed using mean weight diameter (MWD) and geometric mean diameter (GMD) of the water-stable aggregates. Significant differences of MWD were observed between the shaking times. The 10-min shaking was selected as best for structural stability assessment in the studied regions because it resulted in better differentiation of soils on the basis of structural stability. Among the intrinsic properties, soil organic carbon content had the most important role in aggregate stability in all zones. However, electrical conductivity (in addition to organic carbon content) had an important role in aggregate stability in the arid rangelands. Log-normal distribution and GMD could represent better the aggregate size distribution when compared with normal distribution and MWD in the studied regions. Overall, wet-sieving method with shaking time of 10 min is suggested to assess the soil structural stability in rangelands of Isfahan province. Therefore, soil aggregate stability and the factors affecting this vital indicator can be used efficiently for assessing and monitoring management effectiveness and rangeland functionality trend.

---

Topography and land use are among the most important factors affecting the soil formation. Chemical forms of Fe and magnetic susceptibility (χ) are widely used for the evaluation of soil development. This study was conducted in order to determine the effect of these factors on χ. A toposequence was selected in Madvan Plain, Northern Yasouj. Nine soil profiles (paddy and dryland soils) were dug and sampled from diagnostic horizons. Magnetic susceptibility was measured by Bartington Dual Frequency, MS2 Meter at frequencies of 0.46 and 4.6 KHz. Results indicated that less amounts of χ, frequency dependence of χ (χfd%) and CBD extractable Fe (Fed) (3.1, 2.6 and 2.7 times, respectively), and more quantities of oxalate extractable Fe (Feo) and Feo/Fed ratio (5 and 7.2 times, respectively) were measured in paddy soils. The highest value of χ was observed in pedons located on plateau and piedmont plains, and the lowest belonged to those located on river terraces with aquic conditions. Compared to paddy soils, χ enhancement at soil surface was greater (17%) in dryland soils. A positive correlation existed between χ and some soil characteristics such as Fed, clay content and χfd%.

---

CEC of the soil is the exchange sites of organic and inorganic soil colloids. Modeling and Estimation of CEC is a useful indicator for fertility. The new alternative approaches for estimating CEC are indirect methods based on intelligent models. In this research in order to estimate CEC, 485 soil samples were prepared from two regions, chaparsar (Mazandaran in northern Iran) and Bostanabad (North of West Azarbaijan, Iran).In this paper introduces the application of genetic programming. Input parameters that are percent Clay, Organic Carbon and Silt, evaluate using genetic programming, neural network and Neural Inference Systems-Fuzzy models. The results indicate a good ability to intelligent models for CEC Estimation According to indices used in this study. Genetic programming model with a root mean square error of 1.78 and coefficient of determination 0.95 compared to other models have been more efficient and is able to provide satisfactory results, Also are the explicit solutions that reflect the relationship between input an output variable, was presented based on genetic programming. This preferred the genetic programming model adds the other models. Stepwise regression analysis to determine the contribution of each of the parameters indicated in the CEC that organic materials having Most coefficient of variation of 84% is justified CEC and clay and silt, respectively, with a correlation coefficient of 10% and 6% respectively.

---

In the present research, the performances of six empirical models, i.e., simple threshold exceedance, fixed proportion exceedance, quadratic function of storage, power function of storage, cubic function of storage, and exponential function of storage were investigated for estimation of groundwater potential recharge in a semi-arid region. First, the FAO Dual Crop procedure was used to calibrate evaporation from bare soil during the occurrence of potential recharge period. Then, the empirical models were calibrated utilizing soil moisture and potential recharge data. For validation of empirical models, soil moisture and potential recharge were simultaneously estimated for an independent event. Results indicated that 5 of the six models (except for the simple threshold exceedance model) were able to estimate potential recharge with a reasonable accuracy, showing the maximum computed value of NRMSE (Normalized Root Mean Square Errors) of 24.4 percent. According to validation results, exponential, cubic, and power function models provided better estimation of potential recharge in comparison with the linear models. Also, all of the applied empirical models were able to simulate soil moisture during the recharge period with an acceptable accuracy. Finally, the exponential model with minimum NRMSE value for soil water simulation and also acceptable performance of potential recharge estimation was recommended for estimation of potential recharge in the study area.

---

Texture fractal dimension is a physical index to describe soil particle size distribution having a variety of applications. Fractal dimension may be calculated from three relations of mass-time, mass-diameter and modified mass-diameter (Kravchenko-Zhang) with two linear and nonlinear options for fittings. The aim of the present study was to compare methods and select an appropriate one and fitting option for determining the fractal dimension using hydrometer data. Sixty soil samples were collected from four fields of Taqanak, near Shahrekord. After removal of organic matter and other initial treatments, hydrometer readings were obtained at 0.67, 1, 2, 5, 15, 30, 60, 120, 180, 1440 and 2880 minutes and were converted to mass-time or mass-diameter data. Nonlinear fitting of the Kravchenko-Zhang mass-diameter relation was selected as the most appropriate method of calculating the fractal dimension of solid particles, due to its highest coefficient of determination and smallest mean square error and lowest Akaike Information Criteria. Error analysis also confirmed this conclusion. There was a significant, though not very strong, relationship between the fractal dimension obtained by linear and nonlinear fitting of mass- diameter and Kravchenko-Zhang mass-diameter methods. These relationships can be used to correct the fractal dimension determined by other methods and fitting options.

---

Evaluation of chemical and biological indicators of soil in different land uses could be helpful in sustainable range management, preventing degradation of soil quality trend. This study was conducted in Friedan in Isfahan province in 2010 to compare chemical and biological indicators in three land uses (rangeland, degraded dry land and dry land), during two growing seasons (May and September) in three slopes (0-10, 10-20, 20-30 %). Nitrogen, phosphorus, potassium, organic matter, cation exchange capacity and microbial soil respiration were measured. Results showed that all measured characteristics except potassium decreased over an increase in the slope. Maximum values of phosphorus, organic matters, cation exchange capacity and soil respiration were obtained in pasture (28.4 mg/kg, 0.62%, 20.38 cmol/kg, 33.2 mgC/day, respectively)but potassium maximum rate was seen in dry land form (406.8 mg/kg).The effect of season on all measured parameters was significant except for N, while the highest amounts of phosphorus, potassium, cation exchange capacity and soil respiration (28.7 mg/kg, 377.3 mg/kg, 19.6 cmol/kg and 25.9 mgC/day, respectively) were seen in May and the highest organic matter rate (0.68%)in September. The results of this study showed that an increase in the slope, poor range management, and the end of the growing season could be major factors degrading the soil quality indices and soil productivity.

---

This study aimed to evaluate the performance of TM satellite data acquired in June 2009 to map soil salinity in southeast of Isfahan province. Ground salinity data (EC) was collected within 9 pixels, covering an area of approximately 8100 m2 using stratified random sampling technique at 53 sample sites. Spectral indices including TM bands, BI, SI1, SI2 and SI3, PC1, PC2, PC3 and also multiple linear regression modeling and maximum likelihood classification techniques were applied to the geometrically corrected image. Results of regression analysis showed that the TM band 4 had the strongest relationship with EC data (R2=0.48) and also the relationship of the modeling image using TM 3, TM 4, TM5 and PC3 was significant at the 99% confidence level. The accuracy assessment of the stratified TM4 and modeling image into five classes including 0-4, 4-20, 20-60, 60-100 and EC>100 ds/m indicated that there was more than 86% agreement with the field measurements of EC data. Therefore, it can be concluded that the discretely classified salinity maps have higher accuracy than regression methods for identifying broad areas of saline soils, and can be used as appropriate tools to manage and combat soil salinization.

---

Management of organic and inorganic treatments may have positive or negative effects on soil quality, plant growth and human nutrition. The objectives of this study were to determine the effects of organic and inorganic zinc fertilizer application on soil quality indicators and wheat yield. This research was conducted at Agricultural Research Station Roudasht, Isfahan, Iran. Sewage sludge and cow manure (5 and 10 t/ha), ash rubber (1 t/ha), powder rubber (200 kg/ha), ZnSO4 (40 kg/ha) were applied and wheat was cultivated. Soil samples were collected at tilling and harvest stages. After taking samples and measurements of the soil parameters, we determined the critical limits for each category and class rating for the each soil parameters, and the soil quality index was calculated. The results showed sewage sludge and rubber ash were significantly effective in increasing soil bioavailable Zn compared to other treatments. Application of sewage sludge and cow manure at 10 ton/ha improved soil quality. The expanded soil quality index can help better understand the effect of fertilizers on soil. A positive and significant relationship between soil quality indicators and Zn uptake and wheat yields was also observed. Our results indicate that addition of 10 t/ha sewage sludge as fertilizer can significantly improve soil quality, supplying the necessary amount of Zn for wheat growth.

---

Rill erosion is the detachment and transport of soil particles by concentrated flow of runoff. It is the most common form of water erosion in the hill slopes. Rill erodibility is the rate at which soil particle is detached and transported by shear force of the concentrated flow. The study was conducted to determine the rill erodibility in different soil textures in Zanjan province using a rainfall simulator. To this end, samples of eight soil textures consisting of clay, clay loam, silty loam, sandy clay loam, sandy loam, loamy sand, and sandy were collected from land surface and transported to small plots (120 cm  100 cm) on a sloped uniform land (10%). The plots were exposed to five simulated rainfalls with a constant intensity of 60 mm h-1 for one hour. Based on the results, there was a significant difference among the soil textures in the rill erodibility (p< 0.01). Rill erodibility of the soils significantly correlated with mineral fraction (sand, clay, gravel) and exchangeable sodium percentage (ESP). With an increase in sand and gravel percentage, soil infiltration rate strongly increased and consequently production of the concentrated flow steadily decreased. Multiple regression analysis indicated that the rill erodibility in the soils was remarkably related to ESP (R2= 0.85, P< 0.01). Clay soil showed to have the highest rill erodibility among the soil textures due to higher exchangeable sodium percentage (ESP= 13).

---

Soil is one of the most important components in forests and distinguishing soil types and soil capability are first steps in forest management. The main aim of this study was to determine relationship between slope aspect and position, and chemical properties of the soil. Soil sampling was done in Tang-e-Dalab in Ilam province which is a part of southern Zagros. Samples were collected in both northern and southern slopes of oak stand (Quercus brantii). In each slope, three transects 50m apart were sampled. Overall number of samples was 60. After data normalization, the means were compared by Duncan test. Slope aspects influenced organic carbon and total nitrogen of soil. These parameters were higher in northern slope than southern one. Slope position showed a significant effect on C, N and P. Also, concentration of C, N and P were increased by moving down the position. Most amounts of C, N and P were 5.84%, 0.58% and 108.19 mg/kg in bottom, middle and bottom of northern aspect, respectively. The least amounts of C, N and P were 3.31%, 0.24% and 37.83 mg/kg in bottom, middle of southern aspect and top of northern aspect, respectively. The results of this study confirmed that nutrient concentration in northern slope was more than southern slope and nutrient concentration in soil was increased by moving downward.

---

In this research, the hydraulic behavior of two kinds of envelopes including synthetic envelope, PP450 and gravel envelope with USBR standard in two soil tank models with silty loam texture was investigated. Three water heads including 55, 75 and 105 cm (water logging) from drain level were used. The discharge of pipe drain in the steady state condition for gravel envelope and at 55, 75 and 105 cm water heads was 188.9, 172.0 and 897.0% more than those in PP450, respectively. Envelope hydraulic conductivity rates at gravel envelope for 55, 75 and 105 cm water heads were 24.6, 14.0 and 21.2 times higher than those in PP450, respectively, and gradient ratios in these water heads for gravel envelope were 14.5%, 2.8% and 14.2% lower than those for synthetic envelope. There were also different behaviors in the two kinds of envelopes for hydraulic conductivity and entrance resistance of pipe and envelope in 55 and 75 cm water heads relative to 105 cm. In general, according to the measured parameters in this research, gravel envelope showed a better performance.

---

With the advent of advanced geographical informational systems (GIS) and remote sensing technologies in recent years, topographic (elevation, slope, and aspect) and vegetation attributes are routinely available from digital elevation models (DEMs) and normalized difference vegetation index (NDVI) at different spatial (watershed, regional) scales. This study explores the use of topographic and vegetation attributes in addition to soil attributes to develop pedotransfer functions (PTFs) for estimating soil saturated hydraulic conductivity in the rangeland of central Zagros. We investigated the use of artificial neural networks (ANNs) in estimating soil saturated hydraulic conductivity from measured particle size distribution, bulk density, topographic attributes, normalized difference vegetation index (NDVI), soil organic carbon (SOC), and CaCo3 in topsoil and subsoil horizon. Three neural networks structures were used and compared with conventional multiple linear regression analysis. The performances of the models were evaluated using spearman’s correlation coefficient (r) based on the observed and the estimated values and normalized mean square error (NMSE). Topographic and vegetation attributes were found to be the most sensitive variables to estimate soil saturated hydraulic conductivity in the rangeland of central Zagros. Improvements were achieved with neural network (r=0.87) models compared with the conventional multiple linear regression (MLR) model (r=0.69).

---

Soil pollution and accumulation of heavy metals in soils and crops are the most important bioenvironmental problems that threaten the life of plants, animals and humans. This study was conducted to explore contamination of heavy metals in soils of Hamadan province. A total of 286 composite surface soil samples (0-20 cm) were collected thoroughout the province. After preparation of the samples, the total contents of Zn, Pb, Cu, and Ni in soil samples were extracted using HNO3. Total contents of heavy metals were measured by ICP. Contamination factor results showed that most samples were moderately polluted and contamination factor for lead was highly polluted. Interpolated distribution map of contamination factors (CF) and pollution load index (PLI) of the heavy metals were prepared using GIS. The overlap of CF and PLI maps with geology and land use maps indicated that the concentrations of Ni, Pb, Zn, and Cu have been controlled by natural factors such as parent material, but agricultural activities according to excessive consumption of animal manure and chemical fertilizers can increase most of these elements in soil.

---

Considering soil compaction problem in sugarcane fields due to using heavy harvester and haulout equipment under unsuitable moisture conditions, this research aims to assess soil compaction in sugarcane fields located in Da'balKhazaei Plantation unitofSugarcane Development and By-product Company, Ahvaz. Undisturbed soil samples from the furrow (wheel tracks) were collected for measuring soil water content and bulk density. Considering the changes in soil texture of sugarcane fields, for expressing the degree of soilcompactness, in addition to soil bulk density (BD), relative bulk density (BD divided by reference BD) was also determined. The change in soil mechanical resistance with depth was determined by a cone penetrometer. Results showed that most of soil BD values measured in the sugarcane fields were in the range of small root development scale (high limitation). Comparingthe calculated RBD values with optimum value (0.85), it was observed that most of the values were higher than the optimum values recommended for root growth. This shows excessivesoil compaction in the sugarcane fields. The values of cone indices measured in soil profiles indicated that most of the values were higher than either limiting (2 MPa) or critical (3 MPa) values for root growth. Therefore, for improving soil physical fertility and achieving sustainability in crop production, management of farm machinery traffic in sugarcane fields, especially at the harvest time, needs to be reconsidered.

---

Transporting borrow materials for proper infrastructure of water channels to bear the load of such structures is important in the development plans. Therefore, in this research clay nanocomposite material with a weight ratio of %1 was added to the soil. Soil sample was taken from the bed of the C25 canal (distributary of GanjAfrooz diversion dam within Alborz project area) at various intervals and the depth of 1 meter. Unconfined compression strength and consolidation tests were performed on the selected soil. The results showed that the addition of nanoclay to the soil increased the rate of shear resistance, cohesion property and compressibility of soil, respectively, equal to 14.13, 14.13and 82.76 percent. Also, angle of failure and ultimate void ratio decreased. As a result, the addition of nanoclay to the soil makes soil strength and stability greater and there are no problems caused by bed erosion and transporting of borrow material for infrastructure of channel.

---

Consequences of heavy metal accumulation in soils are of great concern. One way of decontaminating heavy metals from soils is using chelating agents, particularly EDTA. In this research, three contaminated soils (with total concentration of these metals of 10.5, 55.8 and 80.6 mmol kg-1) were collected from the surface layer of the lands surrounding a zinc-lead smelting plant in Zanjan province. The extent of Zn, Pb and Cd release by Na2H2EDTA (100 mmol kg-1 of dry soil) from these soils in column leaching experiments (both continuous and pulse addition methods) assembled into half of saturated hydraulic conductivity was assessed. In preliminary experiments, the leaching was stopped due to a drop in hydraulic conductivity. Therefore, the continuous addition method was performed with calcium nitrate as the background solution and the pulse addition method was conducted using this background solution coupled with pH adjustment to 8. Based on the results, the percentage removal of Cd as well as Pb was relatively the same for the two addition methods while the removal of Zn was 13% on average higher in the continuous addition method than in the pulse addition method. For both methods, the removal efficiencies followed the order of complex stability constants (as Pb>Zn>Cd) in a limited concentration range of EDTA to complex heavy metals. Furthermore, in contrast to Cd and Pb, a direct linear relationship was found between the percentage removal of Pb and its total amount in the soils. Surprisingly, the Pb concentration was on average only about one-twentieth of the Zn concentration. The breakthrough curves of both methods showed the mobility order of Cd>Zn>Pb. In general, it seems that the removal pattern of soil heavy metals is dependent not only on the soil type but also on the removal method.

---

Transport of fluoride and consumption of groundwater with excess fluoride concentrations poses a health threat to millions of people around the world. The objective of this study was to simulate transport of fluoride (F) using HYDRUS-1D model. The study was conducted in lysimeters at Lavark research station site in Isfahan. The treatments consisted of two concentrations of F (157 and 315 mg kg-1). The duration of the study was 125 days. Some of soil physical and chemical properties, soluble F and total F concentration were determined during the study. The results showed the transport of F in calcareous soil profiles. This may be due to the high pH and desorption of F ion as a result of repulsion by the more negatively charged soil surfaces. The highest concentration of total F and water soluble F were observed in the 10 cm surface soil layer. The concentration of F decreased with increased soil depth. The correlation coefficient was significant between the water soluble fluoride and the total fluoride (1% level). Also, the difference between the observed t- value and a critical value on the t distribution is statistically insignificant. It showed that the model simulated successfully water soluble F concentration in the soil profile.

---

Lack of knowledge on soil geotechnical properties can cause many problems in the construction and maintenance of irrigation and drainage networks. In general, all of unconventional soils such as gypsiferous soils can cause some problems to irrigation canals. Some studies have been conducted on a variety of problematic soils, but still there is a need for more research activities and field studies. This research was conducted to study the impact of adding perlite and pumice (5%, 10% and 15%) and micro silica (1%, 5% and 10%) on some mechanical properties of soil including shear stress, bearing capacity and Atterberg limits. Statistical analysis was done to compare their averages (P<0.01). Results showed that micro silica had the most effect on shear, bearing and condensation parameters and Atterberg limits of gypsiferous soil, and it improved these parameters of soil. Pumice improved shear, bearing and condensation properties of gypsiferous soil. Perlite reduced the shear, bearing and condensation properties of gypsiferous soil.

---

Determining the spatial distribution of different contaminants in soil is essential for the pollution assessment and risk control. Interpolation methods are widely used to estimate the concentrations of the heavy metals in the unstudied sites. In this study, the performances of interpolation methods (inverse distance weighting, local polynomials and ordinary Kriging and radial basis functions) were evaluated to estimate the topsoil contamination with copper and nickel in Golestan Province. 216 surface soil samples were collected from Golestan province, and their Cu and Ni concentrations were measured. Soil contamination was determined using different interpolation methods. Cross validation was applied to compare the methods and estimate their accuracy. The results showed that all the tested interpolation methods have an acceptable prediction accuracy of the mean content for soil heavy metals. RBF-IMQ and IDW1 methods had the lowest RMSE, whereas RBF-TPS method with the largest RMSE estimated a larger size for the polluted area. The greater the weighting power, the larger the polluted area estimated by IDW. Compared with the ‘‘sample ratio over the pollution limits” method, the polluted areas of Cu and Ni were reduced by 8.38% and 6.14%, respectively.