نشریه پژوهش های حفاظت آب و خاک
سال بیست و ششم شماره 3 (امرداد و شهرویور 1398)

Due to lack of organic matter in most soils of Iran, use of different organic fertilizers in soil is very important for improving soil fertility, increasing water and nutrient storage capacity in soil, increasing water infiltration, reducing soil erosion and improving plant growth and yield. However, these fertilizers alone cannot provide all the requirements of plants for different nutrients and are not sufficiently available to farmers everywhere. On the other hand, high consumptions of chemical fertilizers cause environmental pollution and increased agricultural production costs. So, in order to increase effectiveness of organic and chemical fertilizers and to achieve sustainable agriculture, integrated application of organic and chemical fertilizers is recommended. The deficiencies of micronutrients in the most agricultural soils often reduce the yield and quality of wheat. Therefore, the appropriate concentrations of these nutrients in wheat are important not only for the optimum yield of the wheat plant and its quality improvement but also for the health of humans and animals. In order to improve micronutrients nutrition of wheat and its quality the integrated management of soil fertility is essential.

This research was carried out to study the effects of urea combining with farmyard manure (FYM), municipal solid waste compost (MSWC) and municipal sewage sludge compost (MSSC) on concentrations of micronutrients and sodium in seed, leaf and stem of wheat. The experiment was done in a randomized complete blocks design with 15 treatments and three replications under field conditions at Khalatposhan Agricultural Research Station, University of Tabriz, Tabriz, Iran. The treatments included were: 1) control, 2) 150 kg urea/ha, 3) 300 kg urea/ha, 4) 30 ton MSWC/ha, 5) 30 ton MSWC/ha + 150 kg urea/ha, 6) 60 ton MSWC/ha, 7) 60 ton MSWC/ha + 150 kg urea/ha, 8) 30 ton MSSC/ha, 9) 30 ton MSSC/ha + 150 kg urea/ha, 10) 60 ton MSSC/ha, 11) 60 ton MSSC/ha + 150 kg urea/ha, 12) 30 ton FYM/ha, 13) 30 ton FYM/ha + 150 kg urea/ha, 14) 60 ton FYM/ha, 15) 60 ton FYM/ha + 150 kg urea/ha. The size of each plot was 2.0 m × 1.9 m. At the end of growth period, the plants were harvested and the concentrations of Fe, Cu, Zn, Mn and Na in their seeds, leaves and stems were measured.

The results showed that the effects of treatments on the concentrations of Fe, Zn, Cu and Mn in different sections of wheat were significant. The application of 300 kg urea/ha significantly increased Fe concentration in seed and stem, Cu concentration in leaf and stem, Mn concentration in leaf and Na concentration in stem as compared with the control. There was a significant positive correlation between Zn concentration in leaf and urea fertilizer level. The use of FYM, MSWC and MSSC significantly increased the concentrations of Fe and Zn in seed, concentrations of Fe, Zn and Mn in leaf and concentrations of Fe, Zn, Cu and Mn in stem relative to the control. The concentrations of Fe, Zn, Cu and Mn in different sections of wheat in integrated treatments especially in combined application of 150 kg urea and 60 ton FYM, MSWC and MSSC per hectare were greater than other treatments while the seed Na concentration in these treatments was lower than other treatments.

In general, in order to decrease chemical fertilizers use and achieving sustainable agriculture, improvement of Fe, Zn, Cu and Mn nutrition of wheat plant, improvement of wheat seed and straw quality, combined application of 150 kg urea and 30 or 60 ton FYM or MSWC or MSSC per hectare could be recommended at similar conditions.

Due to the location of Iran in a dry and semi-arid region, it is always affected by sloping instability and erosion as the most severe type of erosion, gully erosion. This erosion pattern occurred in different parts of Iran and continuously over many years, and during erosion and the transfer of high sediment volume, it destroyed roads, facilities, pastures, slopes, etc. This will require the identification of high risk areas and the development of sensitivity maps. In recent years, the processing of satellite imagery as an advanced method has been widely used by researchers to increase the accuracy and save time and money. The object-oriented analysis of images is one of the most important methods for extracting information from satellite imagery, which is based on spectral, form and spatial characteristics and using expert knowledge to identify complications.

In this research, the Lighwan watershed basin was studied as one of the most important sub-basins of Aji Chay in East Azarbaijan Province and the satellite images of Sentinel-2 (2016) with spatial resolution of 10, 20 and 60 meters for the processing and identification of gully Was used. The images were processed using the eCognition software and applied with different types of algorithms to design a semi-automatic model based on object-oriented analysis. Finally, in order to evaluate the accuracy of the model, the identified Gully were mapped out and calculated using ArcGIS software to match the ground reality map and the formation of the error matrix, manufacturer accuracy, user accuracy and kappa coefficient for each of the algorithms.

The results showed that the Density and Compactness algorithms had the highest and lowest accuracy of the manufacturer (manufacturer accuracy was 88 and 78 respectively). While based on Kappa coefficient, the asymmetry algorithm has the highest accuracy compared to other methods (kappa = 0.91). Then, the shape index and density algorithms with kappa coefficient equal to 0.89 and 0.85 provided acceptable accuracy for the classification and identification of the gully.

The use of object-oriented methods due to the increased accuracy of classifying and identifying surface effects and phenomena can be used as a suitable solution for future soil studies and natural phenomena. In the present study, semi-automatic semi-automatic model for ditch identification was presented using spectral and geometric properties of Sentinel-2 satellite images and object-oriented processing in eCognition software environment. Key words: Object Oriented Algorithms, Object-Oriented Processing, Segmentation, Classification

The present study was carried out in west of Iran including: Hamadan, Kurdistan and Kermanshah provinces in an area of 72336 square kilometers. The selected region is located 45 degrees, 20 minutes to 49 degrees and 36 minutes longitude, and 33 degrees, 37 minutes to 36 degrees and 30 minutes latitude. Regarding the studied period and the overlapping of recorded rainfall data in the meteorological stations, 13 synoptic stations including Hamedan, Toyserkan, Nahavand, Malayer, Kermanshah, Islamabad, Sarpol-e-Zahab, Kangavar, Sanandaj, Baneh, Bijar, Marivan and Saqez in three provinces was used. RegCM4 model has been run for 10 years by definite and constant boundary conditions for all four schemes including Grell, Kuo and Emanuel. Grell scheme itself has been divided into two different schemes of Arakawa-Schubert (AS) and Fritch-Chapel (FC). In order to investigate errors of model outcomes, simulated and observed rainfall For the 13 selected stations in Hamedan, Kermanshah and Kurdistan provinces were extracted using coding in the NCL environment.

In all three provinces, Emanuel Scheme has the lowest value for RMSE, and Grell A.S. Schema has the highest amount of RMSE in the provinces of Kermanshah (49.84 mm) and Kurdistan (55.35 mm), and for the Hamadan province, the Kue Scheme has the highest value (36.13 mm). In the case of Bias Error, in Hamedan and Kurdistan provinces the least amount of error is devoted to the Emanuel Scheme (about 18 mm in both provinces), but in Kermanshah province, KUE scheme has the lowest value of the Bias Error ( 24.96 mm). In winter, except the KUE scheme, the three other schemes simulated the rainfall values more than real values, and Emmanuel scheme had the highest over-estimated (about 63 mm). Also in this Season, Grell F.C. With an average of 13 mm, the difference from the actual values, showed the lowest error. In spring, Kue Scheme also simulated rainfall values for all stations less than actual values, but the other three schemes simulated the rainfall values more than real values, and the Emanuel Scheme with an average of 62 mm error was the best For the whole study area. In summer, all schemes simulated rainfall values more than real values, and KUE (1.2mm bias error) and Grell A.S. (122 mm bias error) showed the best and worst results for the whole region, respectively. In the fall, only the Kue scheme simulated the rainfall values less than actual values, and the best simulation was for the Malayer station with a 99.2 mm error. The other three schemes, like the other Season, simulated rainfall values more than real values, and Grell F.C. With an average 6 mm error in the region and Kue with an average of 143 mm bias error in the region, respectively, showed the best and worst results.

on an annual scale for all three provinces, Emanuel Scheme showed the least RMSE values (18.96, 28.95 and 20.69 mm for Hamedan, Kermanshah, and Kurdistan provinces respectively). In addition, Grell A.S. Scheme has highest RMSE and MBE errors for Kermanshah (49.84 mm) and Kurdistan (55.35 mm). On seasonal scale, Emanuel scheme shows the lowest RMSE for the autumn (68.76 mm) and winter (66.8 mm).

Many built structures in rivers have been destroyed or damaged due to lack of proper understanding of scouring. Therefore, predicting the maximum scour around structures in the river has a great importance. Groins widely are designed and constructed in order to river training in the country. This method, which is one of the most effective ways to stabilize the river banks are expanding day by day. In this research, bed deformation around permeable and close groins with different percentage of permeability was numerically simulated by using some of turbulence closure models. Then numerical simulation results were compared with experimental advantages.and at leastAnd, by using the best model of turbulence and calculating the transverse flow criteriaes created by the presence of a groin, the effect of the permeability of groine on the transverse flow power and the longitudinal distance in the downstream of which the continuity of the transverse flow continues is determined.

In other to compare the obtained results from numerical simulation with laboratory advantages, laboratory results around close and permeable groins with 30, 50 and 70 percent of permeability in a flume were used. The length of mentioned flume was 14m, the flume width and height’s was respectively 2 and 0.5m. The slope of its floor was near Zero. Discharge is 50 liters per second, and the average diameter of the bed materials is 0.26 mm and clear water conditions. For numerical modeling, Flow-3D software was used.

The comparison of results of numerical simulation and experimental advantages shows that numerical simulation can predict the bed topography around permeable groins with different percentage of permeability and close groins in channel by live bed as well. In the following, the effect of permeability on power of vortices created around the permeable groins was studied. For this purpose, some criteria that demonstrate the power of secondary flows were computed. Then the changes and continuity of these criteria before and after the permeable groins with different permeability was studied. In the following the same criteria were predicted by ANN method and the advantages of ANN method was compared with results of numerical simulation method.

numerical simulation and especially Flow-3D software can predict the location and maximum scouring around impermeable and permeable groins with different percentage as well. And also the changes of power of secondary flow at the upstream and downstream of groin have a relation with percent of permeability of groin. The ANN results for predicting the criteria of power of secondary flow in different cross section and for different percent of permeable groin shows that this method also has good performance. The results for the calculated cross-flow criteria also indicate that the deviation of the flow by a groin can be roughly ignored with a permeability percentage of more than fifty percent, and for impermeable groin, cross flow continue to extend over 1.4 times the length of the groin. Keywords: permeable and close groins, scour, secondary flow, vortex, numerical simulation, ANN

The prevailing rainfall in warm and dry areas is thunderstorm and cause the great floods. Destruction and sedimentation of the floods cause a lot of damage. considering that Iran is one of the hot and dry regions and the majority of rainfall is a thunderstorm, it is necessary to study the effects of floods caused by thunderstorm. In this study , the effect of flood waves on the bed load transfer in a laboratory slopble channel with unsteady flow was investigated.

In order to perform this study was used a channel that width, height and length were respectively equal to 0.3, 0.4 and 10 meters. The channel was sloppy that we used for of this study, two slopes of zero and 1.5% in the longitudinal. The walls of the canal are made of glass and it is possible to observe the passing flow. The sediments were of river sand in three types of uniform grains with a diameter of 0.85, 1 and 1.2 mm. We produced the triangular hydrograph with peak discharge of 20, 30 and 40 lit / s / m. It was also used to control the erosion of substrate materials and reduce the volume of transfer of the bed from the blade as a protective method against flow.

In this research, experiments were conducted to investigate the factors affecting sediment transport in seasonal rivers. Peak discharge, symmetric flow rates, bed slope and bed particle size play a decisive role in sediment transport. The results showed that with increasing gradient from 0 to 0.015, the cavity length caused by erosion was 1.92. By examining the effect of sediment particle diameter on sediment transport, the experiments showed that by decreasing the diameter of non-sticky particles, the eroded cavity length increased by 1.84 times. Also doubling the peak discharge rate from 20 to 40 liters per second, the length of the eroded cavity was 2.2 times, and the maximum scour depth increased by 56 percent. To reduce the transfer of bed loading to the bottom of the blades, they were zigzag, diagonal and perpendicular. The results showed that using the blade, the sediment transport to the downstream decreases, in which case the perpendicular vertex of the substrate had the best performance and reduced the length and depth of the eroded cavity by the inputs, which in Ultimately, the amount of sediment transfer volume decreased by 89% relative to the non-blown state.

The results showed that the increase of peak water discharge in load transfer was effective, and also the bed sloping and reduction of the particle diameter average size caused a greater increase in bed loading. Based on the experiments, the most critical mode was selected and the blade effect on the erosion of the bed was investigated. The results of the surveys showed that the greatest reduction in erosion was due to the vertical placement of the blades, the crowns of the blades being flattened.

The high use of pesticides in agriculture and the entry of some of these compounds in water and products has led to one of the main strategies in agriculture, increasing the health of the community using less pesticides and, consequently, producing a healthy product and in the future, isolation of pesticides in contaminated water before release in the environment. The poison in agriculture is a chemical substance designed to kill pests. The poison in the food chain is sometimes condensed up to many thousand times, and most of them, especially if they have an aromatic ring, cannot be easily decomposed. The purpose of the study was to remove methylene blue from aqueous solution by beta-cyclodextrin / Zinc oxide composite.

Beta-cyclodextrin / Zinc oxide nanocomposite was synthesized by sol-gel method in aqueous solution under nitrogen atmosphere and successfully identified using FTIR, SEM and XRD techniques. The maximum methylene blue absorbance wavelength was determined using a UV-VIS spectrophotometer and a range of 400 to 800 nm wavelength range at 665 nm. The methylene blue adsorption on β-cyclodextrin /Zinc oxide composite was evaluated discontinuously. The initial pH (1, 4, 7, 9 and 12), the initial amount of adsorbent (0.005, 0.01, 0.015 and 0.02 grams), the initial methylene blue concentration (5, 10, 20 and 30 mg/L), and contact time up to 40 minutes as well as the desorption process were studied. Langmuir, Freundlich, and Tempkin isotherm adsorption models were investigated. Experimental data were studied with different kinetic models. Thermodynamic parameters in surface adsorption including Gibbs free energy change (ΔG0), entropy change (ΔS0) and enthalpy change (ΔH0) were measured by examining the adsorption process at several different temperatures.

The lowest amount of adsorption was observed at pH 7 and highest in acidic pH=1. It seems that increasing methylene blue adsorption in acidic pH is due to the conversion of the pigment to anion, which results in more potent adsorption. Up to 15 minutes after the start of the process, adsorption is performed at more and then performed at a lower rate. The adsorption takes to equilibrium after 30 minutes. The maximum adsorption capacity occurs at a concentration of 5 milligrams per liter than the contaminant in the presence of 0.005 grams per liter of absorbent material at 20 degrees Celsius. The process of adsorption is accompanied by the reduction of entropy. The thermodynamic constants of ΔH0 and ΔS0 are 55.15 kJ / mol and 195.62 j/k.mol respectively. At 20 ° C, the amount of ΔG0 is 2131.65 j/k.mol and increases with increasing temperature up to 37 ° to 5457.19 j/k.mol. The Tempkin isotherm has a good correlation with experimental data with a correlation coefficient of 0.9818 and a constant KT of 0.493. The kinetic model of pseudo-second order with a correlation coefficient of 0.9578 and a constant speed of 0.07 min-1M-1 is a suitable kinetic model for describing absorption.

Forest degradation and land use change are among the factors affecting the changes of soil properties. Microbial characters and enzyme activity, as soil health indicators, are the most dynamic and sensitive soil properties, which play an important role in the nutrient cycle, long-term fertilization and energy flow in the soil. These characteristics provide useful and complete information on the biogeochemical cycles, since they react quickly to changes in the soil environment and provide comprehensive information on the soil physical, chemical and biological properties.

With the aim of studying and evaluating the effect of forest, rangeland and crop cover on soil microbial and enzymes activities, the mountainous habitat of Kodir was considered from the Kojur region in the south-east of the Noshahr city. In this study, four types of vegetation including natural forest (Carpinus orientalis - Quercus macrocarpa), rangeland dominated by Astragalus balearicus -Teucrium subspinosum, rangeland dominated by Stachys byzantina and agriculture field (culticated by Triticum aestivum) were selected. Following field trip, in each of the studied land uses, three transects (50 meters apart from each other) with 200 meters in length were considered. Soil samples (25 × 25 cm area) from a depth of 15 cm were taken at the first, middle and at the end of each transect. In total, 9 soil samples from land uses were transferred to the laboratory for analysis of soil physico-chemical, biological, microbial and enzyme activities.

The ANOVA indicate that the higher values of aggregate stability, clay and water contents, organic carbon, total nitrogen, particulate organic carbon, dissolved organic nitrogen, phosphorus, potassium, calcium, magnesium, fine root biomass, nitrate and nitrogen mineralization and the lower amounts of soil bulk density and carbon to nitrogen ratio were found in forest site. The greater amounts of sand content and the lower values of silt and dissolved organic carbon were observed in the agricultural field. Soil pH, electrical conductivity, particulate organic nitrogen and ammonium were not significantly different among the studied land uses. The highest values of basal respiration, substrate induced respiration, microbial biomass (carbon and nitrogen, phosphorus) and enzyme activities (i.e. urease, acid phosphatase, arylsulfatase and invertase) were found in forest ecosystem, while the studied microbial indices (i.e. qCO2, microbial ratio and carbon capability index) did not show significant statistical differences among the studied land uses. Principal component analysis (PCA) also showed higher values of soil microbial and enzyme activities, biological and fertility in the forest site with a completely different location on the axis.

In general, the results of this study showed that different soil properties in forest ecosystem have better condition than the other studied land uses, while deforestation and land use change is due to decreasing organic matter quality, soil microbial and biochemical activities.

Dams and diversion dams are structures that prevent of fish migration to upstream. Fishways provide access of fishes to their favorite habitat. Canoe_fishway channel is considered friendly environment because of having maximum compliance with it in many countries. In this fishway, synthetic vegetation grasses are used for energy dissipation, opposite of other usual fishways. In fact synthetic grasses are organized of vertical elements that play role of vegetation stems. In addition to creating a safe environment for fish passage, small boats can pass on it easily, because of having proper stems bend radius. Brushes are placed with different arrangements in the Canoe_fishway channel. Each of the arrangements has different effects on amount of dissipation. The main object of the research is experimental investigation of arrangements effect on water depth in the basins and velocity in the slot fish inlet. Also in terms of hydraulic condition is introduced the best arrangement.

Tests were done in a channel with length of 12 meters and width of 0.5 meters, with four arrangements of brushes that were designed with the number of fixed elements. For making of brushes, polyethylene pipes are used with the diameter of 6mm and 20cm length that placed on PVC plates with the thickness of 16 mm. by considering each arrangement, brushes were placed on the bed of the channel. Depth of Flow in the center of basins was measured by depth gauge and velocity in the slot fish inlet was measured by flowmeter in 1cm distances from the bed upto water surface. Also by dimensional analysis be determined effective dimensionless parameters on roughness coefficient (n).

Results show that submergence depth (y/h) has the most effect on roughness coefficient (n) in comparison with other parameters. By increasing the submergence depth (y/h), the Manning coefficient (n) increases and it decreases when submerged completely. The maximum coefficient (n) occurs at the base of submerging. In non-submerged condition in comparison with submerged condition, increase rate of roughness coefficient (n) is for, first arrangement 18%, second arrangement 5%, third arrangement 6% and forth arrangement 20%. This research shows process of velocity diagram (V) is increasing from the bed upto water surface. The velocity range is between 0.6_0.7 m/s. The velocity profile shows increasing process upto height of grasses is very slow. But in higher depth is so fast. Also the process of discharge_depth curve is always increasing that it increases with more slope in non-submerged conditions than submerged condition. According to dimensionless discharge (Q*) are found that roughness coefficient (n) range has difference in constant experimental condition for different arrangements. Roughness coefficient (n) increases 90% for forth arrangement.

According to the obtained results of this study, brushes are effective in roughness coefficient. In this research forth arrangement creates more roughness coefficient, because of having the most overlapping against water flow. So this arrangement in the canoe_fishway provide safe place for fishes, as well as canoes.

Prediction of runoff in order to effective operation of flood control reservoirs and earth flood walls is essential. Predictions also make possible emergency operation of reservoirs by estimating time of floods occurrence and expected damages. Predictions are based on recent meteorological and hydrological conditions of the basin and may include future conditions. Although, most of applications are for flood prediction, these may support water supply, hydroelectric requirements, environmental needs and other requirements for operation. Thus, various complex relations and models such as conceptual rainfall-runoff, linear time series, and hybrid models are developed, but in most cases the calculated values obtained by the models are significantly different. The objective of this research is assessment of precision of HMS-SMA conceptual model and bilinear model in prediction of daily runoff of Maroun basin located in Khouzestan province of Iran. The distinction between the current study and previous studies is that the comparison of the HEC-HMS and bilinear time series models has not been considered so far to predict daily runoff in Iran.

In the current study, daily discharge data of Maroun river within 17 years (1995-2011) at Idenak hydrometric station located in Maroun basin were analyzed. Maroun basin to upstream of Maroun dam according to topography and location of hydrometric stations was divided to four sub-basins and each of the sub-basins were introduced to HMS-SMA conceptual model independently. Amongst different rainfall-runoff models in HEC-HMS conceptual model, Clark’s unit hydrograph was used because of its applicability and acceptable performance in large basins. Furthermore, according to literature review and suggestions, SMA model of linear reservoir base flow was used for estimating of base flow. Flood routing in various reaches were performed by Muskingum method. Additionally, due to possibility of snowfall in the basin, for modeling of snow melt, index temperature method was used. Bilinear models were introduced by Granger and Anderson. In fact, bilinear models are extension of second order Taylor series. Considering to non-persistence of mean and variance of time series of discharge of Idenak station, the time series transformed to a persistent one by differencing, Box-Cox and square root transformation.

Considering to various hydrograph shapes in verification and calibration processes, it is evident that HMS-SMA model has a good precision in estimating low flows compare to high flows. The mean and variance of the time series transformed to a persistent one by second order differencing. Then various bilinear models with various orders were fitted to time series data. Verification of the fitted bilinear models to daily discharges of Idenak station was achieved by portemanteau statistic. Finally, bilinear model in the form of BL(2,2,1,1) with the least Akaike criterion was selected as the best model and was applied for comparing to the predicted daily discharges by the HMS-SMA model.

Summing up assessment criteria of the HMS-SMA and bilinear models show that bilinear model in the form of BL(2,2,1,1) with a coefficient of determination, sum of residual errors and mean square root of errors equal to 0.91, 8.9, and 17.8 respectively, is the best model with high precision in modeling and prediction of daily discharges of Maroun basin compare to HMS-SMA model. Furthermore, it is concluded that by increasing order of moving average in bilinear models, their ability for predicting daily discharges decreases.

The relation between water table depth and evaporation rate from bare soil is of great importance in arid and semi-arid areas. In this region, due to over irrigation, the water table is very close to the ground surface which leads to salinization of the soil. Evaporation from soil columns in the presence of a water table is a important that has received great attention for many decades. The soil-drying process has been observed to occur in three recognizable stages. The first stage is evaporation with constant intensity. Secondly, evaporation is in descending order. The third stage is the residual evaporation of low intensity, which begins after excessive drying of the surface layer of the soil and its effect on reducing the hydraulic conductivity of the soil. Precise determination of evaporation from bare ground surface in semi-arid or arid regions then becomes critically important. The purpose of this study was to determine the effect of water table depth on the evaporation from soil surface, as well as the determination of different evaporation stages.

The soil used in this experiment was loam with a bulk density of 1.32 gr/cm3. The experiment was conducted in a greenhouse and the duration was 74 days. Soils Were sived through a 2-mm mesh and then packed into the soil columns using soil funnel. The soil columns were cylindrical PVC tubes of 200 mm inside diameter. Water table was stabilized at depths of 400, 600 and 800 mm from the soil surface, and the experiment was repeated twice. For stabilizing the water table in different depths, Each soil columns contained a pipe from the botton, to supply water from bottles that maintained the water table constant. The water losses from the soil profile was measured at different depths and times using Delta-T Device Moisture Meter.

The results showed that water content between 0 and 160 mm in the soil column were decrease during the experiment and close the water table depth remained saturated. In the steady-state, the rate of water loss from bottle next to the soil column is equal to the rate of evaporation from the soil surface. In a non steady-state, the rate of evaporation equals the sum of the rate of water loss from bottle and the water depleted from the soil profile. The maximum evaporation from the the water table was related to a depth of 400 mm and equal 384.6 mm and the highest water loss from the soil profile was related to a depth of 800 mm and equal 51.3 mm. By increasing the water table depth from 400 to 800 mm (increased 100%), the evaporation rate from the water table 24 percent and the total evaporation from the soil surface decreased by 16.5 percent. The length of the first stage of the evaporation for the water table depth 400 mm was 2 days and for 800 mm less than 1 day.

The results of this study gave us information on the flow of water above the shallow water table. Water content changes in the soil surface soil are higher than close watertable. Movement of water close water table is liquid and close the soil surface is vapor. The duration of the first stage of evaporation has decreased with increase water table depth. In general, it can be concluded that evaporation from the water table can provide a large portion of the water requirement by the plants.

The hydrological simulation models represent a simplified representation of the real hydrologic system that helps to study the functioning of the basin in response to various inputs and better understanding of hydrological processes. These models are able to estimate the runoff values with the lowest possible time and costs using the simulation of rainfall-runoff process. Despite high efficiency, these models have uncertainty. One of the most important issues among researchers is the elimination of these uncertainties. Hence, the application of the combination technique as one of the most important approaches to improve results of simulations is main aim of this study. In this study, four models including the simple combination models (SMA), Weighted Average Method (WAM), Multi Model Super Ensemble (MMSE) and Modified Multi Model Super Ensemble (M3SE) were used analyze the hydrological simulations, in the Gharesou catchment, located in Kermanshah province.

The Gharesou catchment, with nearly 5354 km2 area, is located in the northwestern parts of Karkheh basin and the western parts of Iran. In this study, the baseline daily data including observed temperature, rainfall and runoff during the period of 1997- 2008 were gathered from selected stations in the study area. 70 percent of the data was used for the calibration period (1997- 2005) and the remaining 30 percent for validation (2006- 2008). To this end, the models in the RRL package such as Simhyd, AWBM, Sacramento and TANK and the SCS-Milc and Hymod models, which are coded in the Matlab programming language, were used. Then, four combination methods including SMA, WAM, MMSE and M3SE were used to improve the results. Finally, the performance of each method was evaluated using normalized root mean square error (NRMSE) and Nash-Sutcliff (NS).

In the present study, all simulated models provide acceptable results. The results of the combination methods showed that the application of these methods led to improve the simulation results. Also, the most improvement of results was achieved by M3SE and MMSE, respectively. For the M3SE method, the value of the NS and NRMSE evaluation criteria were 0.80 and 0.97 in the calibration period and 0.87 and 0.53 in the validation period, respectively.

As resultant it can be expressed that Multiple Combination techniques improved the results of simulated flow by each simulation model obviously. It also may be resulted that recent technique (M3SE) is more efficient than other due to incorporating the bias correction step. Finally it is observed that multimodel simulation generated by M3SE can be better at least comparable to the best-calibrated single-model simulations.

Hydraulic conductivity and retention are important characteristics of porous environment in relation to imposed contaminants. Retention curve and hydraulic conductivity are dependent on soil characteristics and soil solution. When Chlorinated pollutants enter to groundwater will contaminate it due to its compounds. Thus preventing them to enter soil and water is crucial to avoid any contamination. In order to investigate the hydraulic behavior of Trichloroethylene in soil, the retention of Trichloroethylene and water were determined. Determining the saturated hydraulic conductivity and estimating the unsaturated hydraulic conductivity by Mualem -van Genuchten, Mualem - Brooks- Corey and Mualem- Kosugi models in the two-phase NAPL-air systems were further objectives of this study.

In this study, the hydraulic behavior of both Trichloroethylene fluid and water were examined. In order to draw the water and Trichloroethylene retention curves in SiL soil the hanging water column method was used. The constant head method was employed to determine the saturated hydraulic conductivity. The Soil retention parameters for Trichloroethylene and water were obtained based on van Genuchten, Brooks-Corey and Kosugi retention models, using the RETC program. The unsaturated hydraulic conductivity, for both fluids as a function of Matric potential was obtained based on Mualem -Brooks-Corey, Mualem -van Genuchten and Mualem-Kosugi models. The performances of these models were assessed by some statistics including ME, RMSE, EF, CD and CRM.

The obtained results indicated that in a certain amount of liquid phase, Trichloroethylene has lower retention and larger hydraulic conductivity compared to water in soil. According to lower surface tension and viscosity of Trichloroethylene compared to water, the saturated hydraulic conductivity of Trichloroethylene and water were 136.75 and 94.5 cm/day, respectively. For water, the van Genuchten retention model demonstrated highest EF (0.93) and lowest RMSE (0.018) values compared to Trichloroethylene fluid. The other two models were also provided more efficiency for water than Trichloroethylene. In the case of Trichloroethylene, the van Genuchten and Brooks-Corey models showed highest efficiency. Generally, the accuracy of all three models for Trichloroethylene was less than water.

The validation results of hydraulic models shows that the van Genuchten model provides better prediction for retention in soil compared to Kosugi and Brooks-Corey models in the two-phase air-water and air-Trichloroethylene system. Trichloroethylenefluid according to its higher hydraulic conductivity and lower retention compared to water, if enters to soil would have faster movement to the groundwater as a consequence of high infiltration rate to aquifer.

Ber (Ziziphus spp.), is one of the most important of tropical fruits in hot and dry areas that cultures in southern parts of the Iran, including provinces of Khuzestan, Hormozgan, Bushehr, Sistan and Baluchistan and county of Jiroft and kahnuj. According to the climate conditions of these areas such as low rainfall and high evapotranspiration, water deficit is one of the most important environmental factors that influences on cultivation of plants in these areas and creates problems for farmers that find the solution is essential. On the other hand, one of strategies for challenge with the water deficit problem and reduce the adverse effects of drought stress on the growth and yield of crops and fruit trees in tropical regions is the application of the super absorbent materials.

In order to evaluate effects of superabsorbent polymer on survival and vegetative growth of Ber (Ziziphus spp.) grafted seedlings under drought stress, an experiment with 3 levels of irrigation as main plot including: 100, 80 and 60 % of cumulative evaporation from pan class A, 4 levels of superabsorbent polymer: 0, 40, 80 and 120g from A200 source as sub plot, carried out in split plat base on completely randomized block design with three replications on 36 Ber During two years in Date palm and tropical fruits research center in Ahwaz-Khuzestan Province. Before planting, super absorbent polymer was added to the soil of pots in deepness of root development. In this research, Irrigation was carried out by the under pressure irrigation system (Babler). Total water volume needed for each irrigation treatment was compete based on the suggested method of the Food Agriculture Organization (FAO). At first, based on the data of the meteorological data taken from the meteorological station adjacent to the experimental site, the amount of evapotranspiration was calculated, and then the amount of irrigation water base on 60, 80 and 100% of cumulative evaporation from pan class A was compete and accessed into seedlings through the pipeline. At the time of the experiment, the percentage of survival and vegetative growth characteristics such as seedling height, root stock diameter, scion diameter, mean number of leaves per stem length, mean number and length of internodes and chlorophyll content were measured. Statistical analysis of data using SPSS software and mean comparison with Duncan's multiple range test was performed.

Analysis of variance showed that the effects of treatments were significant on some of the studied traits. Based on the results of variance analysis of data, different levels of irrigation and superabsorbent had a significant effect on seedling survival at 5% probability level. As well as the interaction effects of two mentioned factors were significant. Comparing the mean of different levels of irrigation, treatments of 80 and 100% of water requirement were in the same group and did not show any significant difference, but showed significant difference with 60% water requirement. The results showed that under drought stress, application of super absorbent had a significant effect on survival and vegetative growth characteristics such as seedling height, root stock diameter, mean number of leaves per stem length and internodes number. However, on some vegetative growth characteristics such as scion diameter, length of internodes and chlorophyll content did not have a significant effect.

Base on the finding of this research, application 80g superabsorbent and irrigation with 80 % of water requirement for each Ber seedling, the while not significantly reducing the survival and vegetative growth, saved 20% in water use. So, consumption of 80g superabsorbent at the planting time of Ber (Ziziphus spp.) seedlings is recommended.

 Access to hourly temperature due to the more detailed analysis of plant growth processes is of fundamental importance in crop modeling and freezing studies. Also, hourly data is needed for more accurate analysis of climate change effect and atmospheric hazard phenomena on the growth and development of plants. The Bi-sinusoidal model is a precise method in daily temperature modeling which, while considering the sinusoidal nature of temperature variations, is very accurate in detecting the time of minimum and maximum temperature of the day. So far, the accuracy of this method has not been studied in different climates of Iran. The objective of this research is to evaluation of Bi-Sinusoidal model for estimation of hourly temperatures from maximum and minimum daily temperature in different climates of Iran.

For investigation of the efficiency of Bi-sinusoidal model, the data of meteorological stations in different climates from ultra-dry to very humid climates were used. For this purpose, daily and three-hour recorded temperatures at eight stations include: Ahwaz, Ardabil, Bushehr, Gorgan, Mashhad, Rasht, Tehran and Zahak in 2000 and 2005 were used. In this model, sunrise time is considered as the occurrence time of minimum temperature and maximum temperature occurrence is assumed after passing two-thirds of daytime length. These times can be accurately calculated with astronomical calculations. All the modelling calculations were performed in MATLAB software environment. To evaluate the mean error and mean bias of the model, RMSE and MBE indices were used, respectively.

The results showed that in arid and extra-arid stations, RMSE vary between 1.5 to 2 and in humid and sub-humid stations it’s close to 3 Celsius degrees. Also, it seems that the performance of this model is not related to the season, in other words, the hourly temperature error estimation in hot and cold months is not significantly different. The MBE showed that the model underestimate hourly temperature in warm months and overestimate in cold months. However, the bias error is negligible in most of the months and is less than 0.5 degrees. Variation of actual and modelled temperature showed that circadian fluctuation of temperature in dry region is more similar to sinusoidal changes than humid area.

Investigation of circadian temperature fluctuation showed that this model has some trouble in detection of occurrence time of minimum and maximum temperature in humid stations that this is one of the main sources of this model error. But this model simulates the sinus trend of temperature variations properly. According to development of this model based on ordinary circadian temperature fluctuation, in days with the meteorological phenomena such as warm and cold advection to the region, or in rainy days, the model accuracy in estimating hourly temperatures maybe reduced.

Monitoring and evaluation of hydrological drought has become a necessity due to its frequent occurrence in various basins of Iran, and in order to prevent and manage the resulting crisis. In most studies, drought indices have been used to investigate it. Drought indices like Streamflow Drought Index, using statistical calculations on river flow data, ultimately leads to a number that does not have the ability to completely investigate and analyze drought in the studied area. Low flow or minimum river flow is one of the most important drought indices that using different criteria can get a comprehensive understanding of the hydrological drought situation in the region. In most studies, deficit characteristics or low flow frequency analysis have been used. Using different low flow indices is necessary to recognize different aspects of drought and its management. Therefore, the aim of this research was to study hydrological drought of Armand river using all of the low flow indices including flow duration curve, deficit characteristics, low flow frequency analysis and base flow index, which can be used to determine characteristics such as dry and wet periods, amount of flow during drought, drought threshold, duration of drought, volume and intensity deficit and frequency of low flows. Also, in this study, the environmental flow of river was calculated and analyzed.

The Armand River basin with 9961 Km2 area located in Chaharmahal and Bakhtiari province was selected as the study area and the daily flow statistics of the Armand hydrometric station were used from 1957 to 2013. First, by plotting the flow duration curve, three low flow indices Q70, Q90 and Q95 also, the river status including wet, normal and drought periods was determined. Then, the minimum environmental flow was calculated based on the two indices Q75 and Q90, and the deficit characteristics using the threshold level Q70. In order to estimate the return period of low flow series (AM7), a suitable two-parameter gamma distribution function was considered appropriate. Finally, the base flow index was calculated using recursive digital filter method in different time scales.

Based on the results, the Armand River in 28, 22 and 14 years (from 57 years of statistical period) has a flow of less than Q70, Q90 and Q95, respectively (45, 32 and 27 m3/s). Due to the gentle slope of the flow duration curve of basin, also the high value of the BFI index (0.95), it was concluded that groundwater had a high participation in the river flow. The dry period in this river begins when the flow reaches less than 45 m3/s (the threshold level Q70). The results showed that the deficit volume and intensity increased over the past 50 years. Based on the frequency analysis of low flow, in the return periods of 2, 5, 10, 20, 50 and 100 years, the low flow was 34.26, 26.8, 23.1, 21.1, 19.1 and 17 m3/s. Comparison of these values with the minimum environmental flow showed that the river flow becomes less than its minimum environmental flow on average every two years.

Despite the sustained flow in the Armand River (due to the snow-rainy regime of the region, the ability of the basin to save water and drain it during the dry season), the dry periods with different intensities occurred over the past 57 years. Also, the severity of drought has increased due to the increase in the volume deficit of the river in recent years, and if this continues in the future, the severity and frequency of the drought occurrence will increase and the overall flow of the river will decrease.