JWSS - Journal of Water and Soil Science

---

The type and intensity of soil erosion in a region generally depend on climatic conditions, ups and downs, soil and land use. Of these, land use is most important. Using different systems of ploughing after unconscious and non-scientific change of land use affects soil physicochemical characteristics. This fact especially in marginal lands and mountainous regions is more visible. In order to investigate sensitivity to soil loss and erosion in various land uses of Aghajary deposits, part of Margha catchment with an area of 1609 hectares in Izeh city was selected. This was to determine the relationship between soil loss by rain simulator and some soil physicochemical characteristics like percentage of very fine sand, sand, clay, silt, pH, Ec, moisture, Calcium Carbonate and organic materials in different land uses. Then, sediment sampling in 7 points, three replicates and in various intensities of 0.75, 1 and 1.25 millimeters in minute in range, residential and agricultural land uses was done using rain simulator. In order to investigate effective factors in sediment production and erosion, samples of soil layers (in depth range of 0-20 cm meters) equal to the number of sediments were taken. For statistical analysis, EXCEL and SPSS 11.5 software were used. In total, the amount of runoff in residential land use was highest and in agriculture land use was lowest. The amount of sediment in agriculture land use was highest and in residential land use was lowest. Then, the most important factors in sediment yield were diagnosed by multi regression. The results showed that sediment yield and erodibility in land uses have meaningful differences in various intensities of precipitation. Regression models showed that in the production of sediment in various land uses, from among the measured factors, silt, sand very fine, lime, Ec, organic materials and pH had the greatest role. Sand percentage in the residential land use, and very fine sand and organic matter in agriculture land use had the most important role in sediment production. But in range land use, moisture percentage and pH had the biggest role in sediment production.

---

In order to evaluate the effect of concentration and time of foliar spraying of nitrogen fertilizer as supplemenfary to soil nutrition on grain yield and some important traits of rapeseed (Hyola401), an experiment was conducted in complete randomized block design with 16 treatments and three replications in in 2008-2009, 2009-2010 growing seasons, and in 2003-2005 growing seasons in paddy fields of Rice Research Institute of Iran (Rasht). In this experiment, nitrogen concentration (from urea source) at two levels (5 and 100.00) and application time at seven levels including 1-seedling stage: 6-8 leaves 2: beginning of stem elongation 3: prior to flowering 4: 6-8 leaves + stem elongation 5: 6-8 leaves + prior to flowering 6: stem elongation+ prior to flowering 7: 6-8 leaves + stem elongation+ prior to flowering with two control treatments including zero nitrogen fertilizer and basal fertilization as 1/3 at seed sowing, 1/3 at stem elongation and 1/3 prior to flowering stages of 180 kg pure N ha were considered. Results of combined analysis showed significant differences between the experimental treatments in most traits. Spray application of nitrogen (100.00) in stem elongation+prior to flowering stage produced maximum silique length (6.8 cm), number of secondary branch (7.1) and plant height (141.1 cm). Also, the spray application (100.00) in 6-8 leaves+stem elongation+ prior to flowering stage had maximum biological yield (10684.6 kg.ha-1), grain and oil yields (3686.2 and 1489.3 kg.ha-1, respectively) and the maturity period (202 days). Spray application (100.00) in stem elongation+ prior to flowering stage did not induce significant difference with spray application (100.00) in the three stages. Maximum oil content was observed in control treatment (zero nitrogen fertilizer) with an average of (44 %). According to the results of the present experiment, it seems that spray application of nitrogen (100.00) at final growth stage of rapeseed (stem elongation+ prior to flowering stage) increases grain and oil yields.

---

Soil erosion is one of the most important factors affecting soil quantity and quality and is environmental problems in developing countries like Iran. It can have deteriorating effects on ecosystems. This research was carried out in farm lands of the Chehel-Chai watershed, Golestan province to investigate the effect of soil properties on runoff and soil loss. Runoff and soil loss were measured in a completely randomized design in 36 plots with 10×10 m sizes in farm lands. For this reason, this study was conducted using rainfall simulator with 2 mm/min intensity and 15 min duration in 4 replicates. Soil samples were also taken in each plot. Sampling was conducted in October 2009. Results of the Pearson correlation showed that among soil properties, the contents of the lime, silt and fine sand had positive correlations with runoff at 1% confidence level. Also, soil surface resistance at 1% confidence level, the contents of the organic matter and nitrogen at 5% confidence level had negative correlations with soil loss. Finally, the results of multiple linear models showed that the content of lime is effective in estimating runoff and soil surface resistance, and organic matter is effective in estimating soil loss.

---

In this study, the influence of spatial heterogeneity of rainfall on flood hydrograph prediction in three mountainous catchments in south west of Iran was studied. Two interpolation techniques including Thiessen polygons method and Inverse Distance Weighting method were applied to compare the rainfall patterns of surrounding rain-gages in hydrograph simulation with rainfall patterns of nearest rain-gage from the catchment outlet. It was found that the best simulated hydrograph is obtained from rainfall pattern of the nearest rain gage. Moreover, the results did not show any relationship between spatial variation of rainfall and outlet hydrograph. Formation of different local rainfall patterns due to non-stationary rainfall field provoked by irregular topography and their effect on interpolation procedure caused important biases in interpolated rainfall hyetographs obtained by Thiessen and IDW methods. It seems that the observed biases in the response of the catchments are the result of inaccurate representation of spatially averaged rainfall rather than its spatial variability. Hence, in mountainous catchments with irregular topography, the lack of sufficient records caused by poor rain gage arrangement can be highlighted as the dominant source of uncertainty in modeling the spatial variations of rainfall.

---

The main scope of this research is evaluation of Soil Conservation Service Procedure in derivation of initial abstraction of precipitation in watershed scale. For this purpose Dalaki watershed which is located in south east of Iran was selected then by using hec-hms and GIS models and a number of observed rainfall runoff events some parameters like CN of watershed ,K and X of Muskingam method and initial abstraction of precipitation were calibrated through two different search algorithm of univariate and Nelder & Mead methods. The early results of this research indicated the superiority of Univariate search algorithm over the Nelder&Mead method both in calibration and also validation processes. Then using calibrated CN and Initial abstraction parameters which were derived through Univariate search algorithm, the factor between initial abstraction and potential retention of surface runoff (S) in each of sub basins were estimated. 0.13, 0.43 and 0.19 were derived as the above mentioned factor respectively for Minimum, Maximum and mean of the above mentioned factor in this step of the research which showed an acceptable compatibility to the offered factor of 0.2 by SCS. Then in rainfall runoff modeling process of this watershed SCS offers a reliable method of initial abstraction estimation.

---

Soil surface roughness (SSR) is one of the important factors in wind and water erosion studies and control. Several parameters such as surface rock fragments influence SSR. Main objectives of this paper are to study and compare (Allmarass) random roughness (RR) and tortuosity (Tb) indices in coarse, medium and fine grain plains of Yazd-Ardakan and investigate the relationship between indices and desert pavement coverage. Roughness data were obtained by pin roghness-meter and roller chain in the three mentioned plains. RR and Tb were measured in 90cm transect length with systematic sampling pattern after determining the boundary of three kinds of plains. In each plain, 30 transects were randomly sampled and the height of soil surface roughness fractions was recorded. Distance of pins in the applied roughness meter is about 2cm. Therefore, in each 90 cm transect the height of 46 points was measured. Desert pavement coverage was measured in 20×20cm2 plots. The result showed that desert pavement coverage in coarse, medium and fine grain plains were in the range of 55-100, 40-85 and less than 5%, respectively. The relationship between RR and desert pavement was significant. The result of ANOVA (Duncan) showed, RR and Tb were significantly different in coarse, medium and fine grain plains (p-value<0.01). In addition, by increasing desert pavement percentage RR and Tb increased exponentially in base of Neper number. The relation between RR and desert pavement coverage is stronger than Tb and desert pavement coverage. Correlation between the two measured indices was calculated and compared in different plains. The result revealed that about 54, 33 and 14% of the arability in Tb could be explained by RR in coarse and medium grain plains, respectively. The correlations of two indices were significant in coarse and medium plain and insignificant in fine grain plain. The trend of RR and Tb decreased from mountain to plain center. Although RR increased slightly in fine grain plain, the difference was not significant.

---

In this study, the effect of soil structure under saturated and unsaturated flow conditions on nonreactive bromide (BR) transport was investigated. The soil structure treatments consisted of undisturbed columns (prismatic and granular structures), and disturbed columns (single- grain structure). A constant concentration (C0= 0.005 M) of bromide was supplied on the surface of the columns in a steady-state flow condition. For the saturated flow condition, a flux equal to the highest saturated hydraulic conductivity (Ks) of the columns was applied on all of the columns. To create the unsaturated flow condition, a flux equal to the half of the lowest Ks of the columns was imposed on all of the columns. The leaching of the columns was followed for five pore volumes (5PV) and the bromide concentration of the effluent was measured at 0.2PV intervals using bromide selective electrode. The breakthrough curve (BTC) of single- grain structure was sigmoidal (S-shaped) and similar to piston-capillry flow form. In contrast, BTCs of the granular and prismatic structures had a steep initial part and later gradual tailed part. The preferential pathways caused the early appearance of bromide in the leachate of columns of these two structures. Tailing of the BTCs might be due to dispersion and diffusion between mobile and immobile water fractions. In saturated condition, the bromide plume appeared earlier than that in the unsaturated condition because of domination of mass flow and rapid macroporous stream. The results demonstrated the importance of soil structure, preferential pathways, and flow conditions in solute and pollutant transport.

---

The aim of this study was to evaluate soil properties changes in a soil toposequence series in Khuzestan province. Twelve soil profiles were dug in north to south direction in this province, according to the changes of the slope in Karoon’s alluvial deposits. The results showed that soil texture regularly changes to heavier form from north to south. It was slightly of redoxomorphic features in low land. Therefore, salic horizon was formed in this land. Despite high water table in this area, gray or mottle forms were few due to high soil salinity. The carbonates were eluviated from surface to subsurface horizon therefore, calsic horizon was formed due to illuviation. Secondary carbonate accumulation was maximum in the studied area. Both classic horizon and cambic horizon were formed in some pedons under long time cultivation in north regions. These developments were originated from organisms that grew under good drainage class. The smectite was seen in the south part of region. It may have been transported by river and deposited under high salinity and flocculated in contact to Persian Gulf’s saline water. The palygorskite was not seen in all patterns. This clay was removed in pedons under long time cultivation in the north area due to more leaching. High irrigation water can accelerate palygorskite transformation. Also, the XRD patterns showed that chlorite, illite, kaolinite and quartz were main clay minerals in the studied soils.

---

This experiment was carried out at the Agricultural Research Station of Bu-Ali Sina University to study the effect of different tillage methods on bulk density, yield, and yield components of barley (Hordeum vulgare) cultivars under rainfed conditions. Three levels of tillage (CT: conventional tillage, MT: minimum tillage and NT: no tillage) and five barely cultivars (V1: Local V2: AbidarV3: Valfagr V4: Bahman and V5: Makouei) were evaluated in a factorial arrangement in a completely randomized block design with three replications in the growing season of 2010-11. Traits of bulk density plant height, grain yield, yield components, biological yield, and harvest index (HI) were evaluated. The results showed that maximum value of bulk density (1.09, 1.26 and 1.29 g cm-3 for 0-10, 10-20 and 20-30 cm of sampling depth, respectively) was observed at MT treatment. But, there was no significant difference between MT and NT treatments for this trait. In comparison of cultivars, maximum bulk density was achieved with Valfagr cultivar. Using chisel plow (MT treatment) produced maximum plant height (70 cm) and yield components (513 numbers of spike m-2 and 19.2 grain spik-1). Also, among cultivars maximum plant height (72.7 cm) and grain number spike-1 (23.2 grain spike-1) belonged to Valfagr cultivar. Also, results showed that V3×MT treatment had the highest grain yield (3100 kg ha-1). Therefore, the findings of the study recommend using chisel plow (MT treatment) and Valfagr cultivar.

---

In this research, sediment production and delivery amount by Darabkola forest roads was estimated using the SEDMODL model. To evaluate the model results, the sedimentation rate in the above roads was directly measured using rainfall simulator. Also, the paired t-test, BIAS, RE and RMSE were used to assess the results. The analysis showed that the rate of sediment production from study roads' surface using the SEDMODL model and direct measurement under the rainfall simulation were 420.97 and 341.19 tons per year, respectively, and rate of sediment delivered to the stream with sediment delivery ratios of 42% and 51%, respectively, was about 177.58 and 174.02 tons per year. Also, results of the statistical methods of BIAS, RE and RMSE for the aforesaid model were 0.04, 17.59 and 0.71, respectively, and at 95% confidence level, no significant difference was obtained between the observed and estimated data. Therefore, the aforesaid model has the appropriate accuracy and efficiency to estimate the sedimentation rate of the Darabkola forest roads. It was also found that from among the input parameters of model, longitudinal slope of road, precipitation and sediment delivery factors were the most influential factors in the sediment production and transport, respectively.

---

Nowadays, due to the effective role of nitrogen fertilizer in growth, yield and crop quality, farmers apply large amount of chemical fertilizers. High application of nitrogen fertilizers has caused soil and water pollution and environmental dangers, higher nitrate accumulation in plant, and different disease risks in human and livestock. In order to investigate the effects of nitrogen and cycocel application effects on soil nitrate pollution and agronomic traits of rice, an experiment was arranged in split plot based on a completely randomized block design with three replicates at Sari region (north of Iran) in 2010. The main factor was nitrogen in four levels and cycocel was considered as a sub factor in three levels. Results showed that maximum and minimum plant height and fourth internodes bending moment were obtained in 0 and 150 kg N ha-1, respectively. Higher filled spikelet percentage per panicle, grain yield and harvest index were obtained in 100 kg N ha-1 application. Application of N up to 150 kg N ha-1 increased soil nitrate by 44.7 percent. As cycocel application decreased, the plant height and panicle length were reduced, but tiller number per plant, filled spikelet percentage per panicle, and grain yield were increased.

---

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.

---

Rainfall-runoff modeling and prediction of river discharge is one important parameter in flood control and management, hydraulic structure design, and drought management. The goal of this study is simulating the daily discharge in Kasilian watershed by using WetSpa model and adaptive neuro-fuzzy inference system (ANFIS). The WetSpa model is a distributed hydrological and physically based model, which is able to predict flood on the watershed scale with various time intervals. The ANFIS is a black box model which has attracted the attention of many researchers. The digital maps of topography, land use, and soil type are 3 base maps used in the model for the prediction of daily discharge while intelligent models use available hydrometric and meteorological stations' data. The results of WetSpa model showed that this model can simulate the river base flow with Nash- Sutcliff criteria of 64 percent in the validation period, but shows less accuracy with flooding discharges. The reason for this result can be the small and short Travel time noted. This model can simulate the water balance in Kasilian watershed as well. The sensitivity analysis showed that groundwater flow recession and rainfall degree-day parameters have the highest and lowest effect on the results, respectively. Also, ANFIS with the inputs of rainfall 1-day lag and evaporation 1-day lag, with Nash-Sutcliff criteria of 80, was superior to WetSpa model with Nash-Sutcliff criteria of 24 percent in the validation period.

---

Splash erosion is one of the most important water erosion types, causing initiation of other types of water erosion. The objective of this study is to model the splash erosion using fuzzy logic approach in part of northern Karoon basin. The major land usage in the area are irrigated farming, dry land farming, pasture and degraded pasture. For the purposes of this study, soil properties including organic matter; CaCO₃; surface shear strength (SSS); particle size distribution; mean weight diameter (MWD) and soil splash erosion were measured under four different slope conditions (S:%) and rainfall intensity (RI:mm.h-1): 5-50, 5-80, 15-50, 15-80, respectively, using multiple splash sets (MSS) at 80 different locations. Splash erosion was modeled based on combinational rule of inference under five conditions for selection of different operators. The efficiency of the models was evaluated using mean square error (MSE) between observed and estimated values. Results revealed that all models are capable of predicting splash erosion. Also slope, rainfall intensity, MWD, SSS, fine sand and coarse silt attributes were found to be appropriately and precisely using splash erosion.

---

Knowledge of variation in soil properties from each event to another is very important for the determination of critical periods during which soil is susceptible to erosion processes. This study was carried out to investigate soil loss in sequential rainfall events in Zanjan Province. Toward this, ten soil textures samples were taken and were transported to small plots (60 cm×80cm) with 20-cm depth) on a 8% slope land at three replications. The plots were exposed to ten simulated rainfalls with an intensity of 55 mm h^-1 for 30-min and 5-day intervals. A total of 300 simulated rainfall trials were carried out at the plots.  Results indicated that soil moisture, runoff production and soil loss were significantly affected by rainfall events (P< 0.001). Increasing soil moisture and consequently decreasing soil infiltration capacity were the most crucial element in increasing runoff production and soil loss in the sequential rainfall events, in a way that about 84% of soil loss variation in the rainfall events could be explained based on antecedent soil moisture. After the fifths rainfall event, no significant differences  was found in soil infiltration capacity as well as runoff production because of soil moisture reaching to the water-holding capacity. Nevertheless, an increasing trend was observed in soil loss after fifth event which could associate with presence of more erodible soil particles on the surface and consequently increasing the concentration of surface flows.

---

In this research, a comprehensive simulation model for water cycle and the nitrogen dynamics modeling including all the important processes involved in nitrogen transformations such as fertilizer dissolution, nitrification, denitrification, ammonium volatilization, mineralization, immobilization as well as all the important nitrogen transportation processes including nitrogen uptake by the plant, soil particles adsorption, upward flux, surface runoff losses and drain losses, was used for fertilizer management modeling in a sugarcane farmland in Imam Khomeini Agro-Industrial Company using a system dynamics approach. For evaluating the model the data collected from Imam Agro-Industrial Company equipped with a tile drainage system with shallow ground water and located in Khuzestan province, Iran, were used. The statistical analysis of the observed and simulated data showed that the RMSE for determining the accuracy of simulation of the nitrate and ammonium concentration in drainage water is 1.73 mg/L and 0.48 mg/L, respectively. The results indicated that there is good agreement between the observed and the simulated data. Nine scenarios of fertilization at different levels of urea fertilizer were modeled including one scenario of 400 kg/ha, two spilit scenarios of 350 kg/ha, two spilit scenarios of 325 kg/ha, two spilit scenarios of 300 kg/ha, one scenario of 280 kg/ha and one scenario of 210 kg/ha. Results of the modeling showed that the scenario of 210 kg/ha has the highest nitrogen use efficiency (52.3%) and the lowest nitrogen losses consisted of denitrification, ammonium volatilization and drainage losses (17.82, 7.16 and 92.59 kg/ha, respectively). The results revealed that increasing the consumption of urea fertilizer greater than 210 kg/ha increased the overall nitrogen losses and reduced the nitrogen use efficiency. Meanwhile, this model can be used for managing the fertilizer and controlling the nitrate and ammonium concentrations in the drainage water to prevent the environmental pollution. Also, the system dynamics approach was found as an effective technique for simulating the complex water-soil-plant-drainage system.

---

In agricultural development many factors such as weather conditions, soil, fertilizer, irrigation timing and amount are involved that are necessary to be considered by the plant growth simulation models. Therefore, in this study, the values of soil water content at different depths of soil profile, dry matter production and grain yield of winter wheat were simulated using AquaCrop and WSM models. The irrigation treatments were rain-fed, 0/5, 0/8, 1 and 1/2 times of full irrigation conducted in Agricultral College of Shiraz University during 2009-2010 and 2010-2011. The models were calibrated using measured data in the first year of experiment and validated by the second year data. The accuracy of soil water simulation was used to refer to the accuracy of simulated evapotranspiration. The accuracy of soil water content at different layers of root depth in the validation period was good for the WSM model (Normalized Root Mean Squer Error, NRMSE= 0/14). But the AquaCrop model showed less accuracy for soil water content (NRMSE=0/26). However, the values of predicted and measured crop evapotranspiration were close together at full irrigation treatment, the accuracy of AquaCop predictions was decreased with inceasing water stress. WSM model has had a good estimation of the dry matter and grain yield simulation with NRMSE of 0/15 and 0/18, respectively. However, they were simulated with less accuracy in the AquaCrop model with NRMSE of 0/19 and 0/39.

---

The status of soil K forms and its relationship with K content in wheat grain at harvest period may be important for nutrient fertility management. For this purpose, 40 surface (0-20 cm) and subsurface (20-40 cm) soil samples and also grain samples from wheat fields of Darab region were randomly sampled in 2014. Then physicochemical properties of the soils and the contents of soluble, exchangeable and non-exchangeable K and contents of K in wheat grains were determined. Results indicated that content of soluble, exchangeable and non-exchangeable K in the studied soils ranged 15-70, 91-443 and 396-1182 mg kg^-1, respectively. Significant relationships were obtained between soil K forms and clay, calcium carbonate and CEC. Although it is expected that content of easily available K (soluble and exchangeable K) was low at late stage of wheat growth, most soils had sufficient K content. The positive and significant relationship among different K forms was indicative of the K forms equilibration at the time of wheat harvesting. Although there is not a relationship between K content in wheat grain and soil K forms, it was shown that K content of wheat grain was correlated with soil exchangeable K and clay content in fine-textured soils (clay > 30%). It means that ammonium acetate may be a suitable extractant for estimation of soil K status and plant available K of heavy-textured soils in the studied region.

---

This research is aimed to determine the contribution of sub-basins flow to total watershed flood in Khersan river basin located in Kohkilooyeh and Boyer Ahmad province. To do this, the rainfall-runoff model HEC-HMS was used to simulate peak runoff values for 11 sub-basins. HEC-HMS input was constructed using GIS. The results suggest that the change in different return periods is accompanied by small change in prioritization of flood-potential of the sub-basins; so that for return periods of 2, 50 and 100 years, the most contributions came from sub-basins 1 through 11, respectively. With respect to area and flow rate, contribution of sub-basins to watershed total flow was different. The effect of area was between 0.31 to 1.03 percent; namely, sub-basin 6 showed the highest rank and basin 7 showed the lowest one. With respect to peak flow rate, the effect of individual exclusion of sub-basins, resulted in contribution between 51.2 to 1004.2 m³/s, that is, sub-basin 6 showed the lowest effect and the sub basin 11 showed the highest contribution.

---

This study aimed to investigate the dynamic changes of rill erosion and evaluate the ability of Hairsine-Rose model in estimation of sediment concentration. The experiments were carried out at the Soil and Water Conservation Research Station of Tehran University in Kuhin, Qazvin. Four flow rates were tested at three slopes in trapezoidal hand-made rills of 3 m long and 5 cm base width. Outflows were sampled periodically to determine changes in sediment concentration. The experiments were repeated in early autumn and mid spring. Sediment concentration showed a dynamic change with time which was affected by flow rate, slope and season. The sediment concentration was high in the first few minutes of the experiment but over time, dropped and finally reached a more or less stable state. The effects of flow rate and slope were more detectable on sediment concentration at the early unsteady conditions than at the final steady conditions. Though the temporal trends were similar, the sediment concentration was generally higher in autumn than in spring. Hairsine-rose model showed a better performance at lower flow stream powers but generally over predicted sediment concentration showing a systematic error probably due to model formulation.