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

Arsenic is one of the metalloid which has recently attract attention. Long-term exposure to As can lead to skin, bladder, lung, and prostate cancer. The presence of As in soil and water causes its transfer to the different parts of plants. So far the effects of hematite nanoparticles (α-Fe2O3) on the concentration of As in plants grown in contaminated soils have not been evaluated, therefore this study was conducted to evaluate the effects of hematite nanoparticles on the concentrations of As, P, Fe, Zn and Mn of corn plant growth in contaminated soils.
Materials and Methods For this purpose a factorial experiment was conducted using a completely randomized design and three replications. The experimental factors were the application rates of hematite nanoparticles (0, 0.05, 0.1 and 0.2%) and levels of soil As (0, 6, 12, 24, 48, and 96 mg/kg). Soil samples artificially contaminated by different levels of As using Na2HAsO4.7H2O salt and incubated for 5 month. Hematite nanoparticles were synthesized from iron nitrate, and their properties were investigated by XRD, SEM and TEM techniques. At the end of incubation time, the hematite nanoparticles were added to the As contaminated and control soils and after 1 month the concentrations of available As in the samples were measured, then planted to corn. After 75 days of sowing, plants were harvested the dry weights of aerial parts and roots, concentrations of As, P, Fe, Zn and Mn of these parts were measured.
Results and Discussion The results showed that the concentration of soil available arsenic and arsenic concentrations of root and aerial parts increased as the concentration of soil total arsenic increased. As accumulated in the roots more than the aerial parts. Contamination of soil by As increased the concentrations of P and Mn and decreased those of Fe and Zn in corn root. The concentrations of P, Fe and Zn in the aerial parts of corn decreased and that of Mn increased as a result of As application. As reduced P translocation from the roots to the aerial parts of corn plants. The application of hematite nanoparticles significantly decreased the concentrations of arsenic in soil and in root and aerial parts of corn. The application of hematite nanoparticles increased the dry weight of root by increasing the root concentrations of Fe and Zn and decreasing the accumulation of As, P and Mn in the plant root. It also increased the dry weight of aerial parts as a result of increased concentrations of P, Fe and Zn and reduced accumulation of As and Mn. Those, but they decreased the concentrations of P, Fe, Zn, Mn and dry weight of root and aerial parts of plants grown in uncontaminated soils. The application rate of 0.2% hematite nanoparticles gave the maximum effect. The reductions in the concentrations of root and aerial parts arsenic at the application rate of 0.2% of the hematite nanoparticles were 54.64 and 40.20 percent respectively for the soils with total As of 96 mg As/kg.
Conclusions As at the concentrations of less than 12 mg/kg increased root and aerial parts growth and dry weights but at the concentrations higher than 12 mg/kg, by creating toxic effects and disrupting the balance of nutrients, reduced the growth and dry weights of root and aerial parts of corn plant. These results suggest that hematite nanoparticles can be use for remediation of arsenic contaminated soils on the condition that further research is carried out on their environmental effects.

The increasing flooding trend in recent years suggests that most of the country's regions are vulnerable to invasions of periodic and destructive floods. In this aspect, many cities, villages, industrial and agricultural facilities and residential areas are prone to flood occurrence, as well. Therefore, the basic identification of flood process within the catchment area is one of the most important measures in flood control and the mitigation of damages. The main objective of this research is to investigate and identify the flood areas and the effect of watershed management on flood peak discharge in the outlet of the Bar watershed, Neyshabour, located in Razavi Khorasan province.
Material and For this purpose, the basin was divided into 20 subbasins and the physical properties of the whole basin and subbasins were determined using the geographical information system oodand in a digital format. Then, by using the HEC-HMS hydrologic model, the corresponding flow discharges were calculated for each subbasin. Then, by successively deleting subbasins at each model runtime, i.e. Single Successive Subwatershed Elimination method (SSSE), the whole basin water discharge was calculated after the flood routing in the main streams without the subbasin by using the kinematic wave routing approach, Thus the effect/share of each subbasin in the production of flood is identified. Also, the flood discharge of the basin was calculated in the basin area unit and the flood index (f) was the basis for the priority of the basin.
In calibration process, two parameters of curve number and manning coefficient were selected as the most effective parameters on flood discharge. The high Nash-Sutcliffe coefficient in flood events showed that calibration of the model in the watershed basin was appropriately done. The results showed that the subbasin B1 (in the northern part of the watershed) in the return periods of 50 and 100 years had the hieghst peak discharge of 38.9 and 44.1 cubic meters per second at the outlet of the subbasin, and the subbasins B11, B13 and B19 (in the western parts of the watershed) showed the minimum peak discharge. Also, according to the index (f), in flood plains with return periods of 50 and 100 years, the subbasins B4 and B3 (in the northern half of the watershed) ranked first and second, respectively and the subbasins B6, B11, B12, B13, B14 and B19 (in the southern part of the watershed and in the eastern and western parts of the watershed) showed the lowest priority in terms of their participation in basin flood. In subbasin B1, the highest level of peak discharge has been observed in the highest level of biological operations, ranging from 41.27 to 44.73 percent. On the other hand, the results showed that the higher the proportion of the biological activity to the subbasin area, the more obvious the role of these projects in reducing peak discharge. According to the study, the role of structural activities in reducing the flood peak is lower than biological activities, and increasing the number of structures along the river route will reduce the peak peak area of the subbasin.
By investigating the effect of biological activities and the construction of gabion check dams on the flood discharges, it can be said that the role of biological activities in reducing peak flow and flood volume is much more effective than structural activities (construction of gabion). Therefore, the CN factor is an effective and controllable factor for flood discharge of the basin, on reducing peak flow.

The occurrence of dust storms as a global threat poses annually a lot of financial and human losses to the country. Hence, according to the topographic-climatic conditions and the significant increment in the dust storm occurrence frequency, investigation and survey to mitigate the harmful effects is inevitable.
In this research, dust storms were evaluated based on the spatio-temporal frequency of 25 synoptic meteorological stations in western and southwestern Iran during the statistical period of 22 years (1995-2016). Daily Data were analyzed using eight observation points in the form of special codes using statistical methods. The average dust accumulation in the cluster and average monthly, annual, warm and cold periods was calculated and plotted using R software for all stations.
The results of the annual survey of dusty phenomenon indicated that Abadan station had the highest number of days in the whole region with the average frequency of 16 days in year. Bostan, Sarpol-E-Zahab and Ahwaz are the most frequent. 63.02% of the outsider dust storms events which distributed from 9:30 to 18:30 local time. 83.47% of the local dust storms events at 3:30 local time, the highest frequency of dust storms in June and July is the lowest in November and December. Also, the percentage of local and outsider dust storms frequency during the statistical period showed that 74.94% of the dust storm days were outsider sources and 25.06% of the local sources. The average frequency of dust storm in the under studied area is 36 days in warm and cold period of the year recpectively. In all studied stations, most of the dusty days were related to the warm period.
Whatever the disturbance of a dust wave originates from its original source, its intensity decreases. Therefore, it can say that the primary sources of dust phenomena majority of the under studied area, are such as Iraq, Saudi Arabia, Syria, and North Africa deserts. Investegition of spatial pattern showed Dust storms were increased from the north to the south. While along west towards the east, it is in the Quadratic polynomial form. As a result, dust storm occurrence frequency increased from west to the east, and then it decreas towards the east. Zagros Mountains positioning, which run northwest-southeast through western Iran prevents dust storms from entering those area like a dam. Also the offset of this phenomenon is more time-consuming compare to other regions.

Improvement and development of forecasts in ungaged catchments needs to understand the interaction between the catchments parameters and hydrological response of catchments in different climates. Flow duration curve is one of the best methods for showing the hydrological response of basins and have different application in the fields of hydrology and related sciences. The shape of flow duration curve also reflected the impact of climatic, geologic and physiographic parameters on river flow and is hydrological response of catchment. In the past decade, several studies have been done on the impact of geometric and hydrological parameters of basins on the shape of flow duration curve. Most of these studies are empirical and classified in graphical and statistical methods. Graphical method, focus on o the effect of physiographic and climatic characteristics on the shape of flow duration curve. Statistical methods on statistical distribution and correlation of the physical characteristics of the catchments are concentrated. Factor analysis of parameters affecting the flow duration curve lead to accurate modeling and interpretation of catchments hydrology. The aim of this study was to determine the most important hydro-climatic and geometric factors affecting on the shape of flow duration curve and investigation of relation between them in the different climates.
In this study, for factor analysis of effective parameter on the shape of flow duration curve, climate maps were prepared and intersect with the layer of fourth-order watersheds, and watersheds located in every climate were selected. Then at least 30 hydrometric stations with appropriate data and common period (1976-2001) in each climate zone were selected. Then 10 geometric and hydrological parameters affecting the flow duration curve including: average height, area of the watersheds, Gravelius coefficient, slope, main river length and hydro climatologically parameters including: annual rain fall, base flow index, curve number, permeability and the number of rainy days, were calculated for each basin. Flow duration curves plotted using daily flow data and the slope between Q33 to Q66 was computed as an indicator of the shape of flow duration curve, using coding in MATLAB programming environment. Factor analysis was performed and effective factors on the shape of curve were identified. The regression between the flow duration curve index and selected factors in different climate zones were extracted and analyzed.
The results showed that the parameters of first factor including: curve number, base flow index, recession index and the number of rainy days in all climatic zones were common except dry climate zone. Curve Number in all the climate zones has the highest weight of influence. The Weight influence of geometric parameters in all areas was higher than hydrological parameter, except in very humid zone. Parameters selected for factor analysis in the humid zone were found the highest variance explained with 88 percent and the lowest in the Mediterranean zone with 72 percent. Normally distributed errors, the coefficient of determination more than 0.90 and the coefficient of Durbin Watson between (1.5-2.5) reflects the confidence on the regression equations to estimate the slope of flow duration curve in untagged catchments in different climatic zones.
Overall Conclusions of factor analysis in different climatic zones became clear that some of the parameters classified in the first class including: curve number, Base flow index, and number of rainy days were common in all climatic zones. With the one exception that in dry zone the number of rainy days placed in a second class. The factor of curve number in all climate zones has more effective than other parameters.

Nowadays, water scarcity is a major challenge for our country, Iran. Therefore, storage and optimal operation of limited resources, including water stored in dam's reservoirs, is one of the issues of interest for researchers in the field of water resources. In this paper,optimization of single-reservoir operation problem is solved by using one of the newest heuristic algorithms, named charged system search algorithm. Generally, this algorithm is based onthe electrostatics laws to determine the quantity of resultant force. At the first time, Kaveh and Talataheri (2010) proposed this algorithm and examined its capabilities for solving engineering problems and sample functions. The results showed that the algorithm has a good performance. Therefore, its use for solving engineering optimization problems is recommended. However, a review of literature shows that the use of this algorithm is very limited in the field of water resource engineering.
In this paper, the simple and hydropower operation of Dez reservoir, over 5 and 20 yearly operation time period are solved using the proposed algorithm. In order to solve these problems, two different formulations are proposed considering water release or storage volume as decision variables of the problem in the first and second formulation, respectively, and the results are compared to other available methods.
Comparison of the results shows the capability of the proposed algorithm, in which the results of first formulation are better than the second one’s. In other words, the results of first formulation for solving simple operation problem over 5 and 20 years are reduced 11.29% and 16.69% in comparison with the results of second formulation and also using first formulation for solving hydropower problem the results are improved 20.06% and 37.66%. Furthermore, the results of proposed algorithm for solving simple operation problem over 5 and 20 years are reduced 33.64% and 74.97% in comparison with the results of particle swarm optimization and also using proposed algorithm for solving hydropower problem the results are improved 6.53% and 41.48%. In addition, the results of proposed algorithm for solving simple operation problem over 5 and 20 years are reduced 7.79% and 35.59% in comparison with the results of genetic algorithm and also using proposed algorithm for solving hydropower problem the results are improved 11.32% and 67.43%.
Investigating these results with the results obtained using other existing algorithms indicates a better performance of charged system search algorithm for solving the reservoir operation optimization problem. According to these results, the use of this algorithm is recommended for solving other problems in the field of water engineering.

for sustainable development, water recourses quality control has special importance. In this way, recognizing the stratification of reservoirs is very important. Also surveying water quality in the case for salinity of surface water resources of dams that figures as the main resource of drinking water of a country is one of the challenges of water resources management. For studying the stratification in reservoirs, turbulence dynamic and internal flows must be analyzed in reservoir. Because of many phenomenons that effect on stratification, studying this phenomenon is possible just with the using of dynamic models. In this way using numerical models beside yield measurements are unavoidable.
Because of the the Taham dam's water usage as for drinking water of Zanjan city, so studying the water quality of has more importance, so in this research for studying the status of the water quality of taham dam, modeling of temperature and TDS (Total dissolved solids) with CE-QUAL-W2 software in one year period from (December 21, 2011) until (December 20, 2012) was done. The value of observed temperature and TDS were compared with modeled data. With the results of modeled data, variation trend of said parameters and the stratification phenomena in taham dam was studied.
In this research, the value of the observed temperature and salinity were in good agreement with modeled data. The results showed that the stratification phenomena happens in reservoir and lasts 8 months in the year. This phenomena begins at the end of March and improves in August and September (with 20 degrees difference between temperature of upper and lower water layers). Turbulence is also seen from January until March. In fall, because of decreasing the difference between temperature of upper and lower water layers, stratification in the reservoir is decreased and in winter it disappears. Results also shows the salinity stratification happens simultaneous with temperature stratification in the reservoir.
From the view point of stratification, Taham dam is in the case of Thermal Monomictic reservoirs. In fall and summer seasons, the value of salinity for all of depths of dam is a constant value. This process is continued until May and after this time because of beginning the stratification phenomenon the difference between layers in depth increase and at summer the difference become maximum.

Rivers are source of any civilization formation and also any societies. As a result of many factors like human interferences, earthquake, flood and sand removing from river margins the dynamic balance of rivers changes. The consequences of these undesired phenomena's are erosion of river banks and lateral migration. On the other hand, these derbies could fill dam reservoirs. Decreasing useful life of dams and also producing operation difficulties are other two side effects of erosion. Hence bank erosion and lateral migration creates a lot of economic and social problems. Noticing above mentioned notes we conclude that investigation about hydraulic of rivers and also understanding dike hydraulic and characteristics is necessary and important. One of the common and simplest method of coastal protection is using spur dike that are useful structure to improve the flow pattern and used for control and reduce the coastal erosion of rivers.
Bandal like structure is new type of spur dike that is compound from a submersible and no submersible spur dike.
In this research we focus on Bandal-like dikes. Explicitly we investigate the effect of dike submersion on erosion and sediment formation around a 90 degree mild bend. Experimental tests are performed under various hydraulic characteristics. Three types of dike are provided. Heights of dike are 6.5,8.5 and 12 centimeters and the width of dike is 20 percent of flood width. Also three discharge 25,27,30 and 33 liters per second for two type (submersible and no submersible) of spur dikes were tested. In this study used a 90 degree bend flume with 0.7m width
The result of this study show that the maximum of scouring rate, occur for submersible theme test and the minimum is for Bandal – like with 12cm height. Also the result show that depth of scouring with Bandal like spur is reduce than no submersible theme structure and with increased the height of structure depth of scouring will be reducing. The result show that with increased the discharge, the Thalweg divert to middle of canal.
The results of this study can be summarized as: 1) The maximum of scouring is for theme submersible test.
2) The minimum of scouring is for Bandal like with 12cm height.
3) Between Bandal like spurs, the maximum of scouring is for structure with 6.5 cm height.
4) Increasing the height of structure reduce the depth of scouring.

Groundwater as one of valuable water resources has always been of interest to researchers. Preparing an interpolated water level zoning map is one of fields that is can be acquired via best interpolation methods among different available methods. The Kriging method, based on semivariogram analysis, is one of traditional geostatistical methods. Interpolation precision is depended on the suitable selection of variogram. Empirical Bayesian Kriging (EBK) method is developed to estimated semivariogram parameters during simulation process. The objective of this study is investigation of EBK to increase the precision of groundwater level interpolation zoning and its results compare to other deterministic interpolation method including inverse distance weighted, radial basis function and local and global polynomial functions.
Material and This study is based on annual mean groundwater level of 57 deep well in Ghareh-Sou aquifer located in Golestan province during 2005-2016 period. Different deterministic and geostatistical interpolation methods evaluated using cross validation technique to groundwater level zoning. The best semivariogram selected for Kriging and EBK methods and finally the model with minimum error is determined and its map is drawn.
Findings: The results of cross validation in study area showed better results for degree two local polynomial methods among deterministic method even it had less error in comparison with Kriging method. The EBK methods, with simulation of fitness of suitable variogram on groundwater level data, led to decrease in Kriging error (23 to 16 meter) and had close precision to local polynomial method.
Although the error of prepared maps based on EBK and local polynomial methods have not significant differences, but there are considerable discrepancies between these maps. The EBK basis map show smoother spectrum of groundwater level changes and the drawn pattern is proportionate with general slope of study area.

Background and objectives Plant cover as a major element of natural ecosystems has high influence from climatic indices. In this regard many models and indices were produced with the aim of clarifying the relationship between climate and plant cover. The air temperature variations have a major role in management of plant cover distribution especially in humid region. The most studies about Plant cover in Iran were based on statistical methods and remote sensing technology and only monitor vegetation changes. So presentation of climatic indices for prediction of natural plant cover is necessary. The aim of this study is evaluation and introduction of most suitable climatic thermal index and use of this for prediction the appropriate elevation for exotic species cultivation especially conifers in north of Iran forests.
Materials and methods In the first step, three indices of WI (Warmth Index), BWI (Biological Warmth Index) and EWI (Effective Warmth Index) were calculated based on monthly mean temperature data of 34 north strip synoptic stations of Iran in decade of (2005-2014). Variation of indices with elevation was achieved as altitudinal gradient relations and the variation map of them was drawn. In second step to investigate the efficiency of calculated indices in determination the plant cover distribution, the study region was smaller and limited to educational and research Kheirud forest. In this stage the forest type’s distribution maps of Kheirud was collected and compared with indices predicted plant cover distribution maps.
Results In this research, some new thermal indices with plant and tree cover distribution in Golestan, Mazandaran and Gilan provinces have been comparatively studied and predicted tree cover by these indices was compared with available map of Habitats forest types of Patom and Baharbon. With respect to recorded species of the region, all of the indices present good predictions for the plant cover of the region under study. The results of statistical analysis showed that effective warmth index (EWI) with highest correlation (R2=0.964) and a lowest error value (RMSE=3.87 ) is the best indicator for estimation the plant cover of the region. This index by considering both range and average of the annual temperature presents the most accurate prediction (better differentiation) for plant cover especially in region that their WI index is similar. As well, the changes of indices with elevation were studied and related regression equation was presented. Furthermore suitable height in terms of temperature conditions was proposed for exotic species that have thermal potential for importing to the north of Iran forests.
Conclusion Predicted plant cover by EWI index has more overlapping with north of Iran plant cover in compare to other two indices. Furthermore with lower error in estimation the elevation of boundary variation has more accurate prediction. With respect to the passage of 5 decades from entering coniferous species in northern forests and cultivation of them in various elevation and The failure of most forest cultivation that conducted on the average height of less than 900 meters above sea level, Based on this study proposed that in future forestry the best elevation will estimated with respect to more suitable predicted elevation by EWI index.

Iran and accordingly the study area due to placing on the dry belt region of the world are always faced with the drought phenomenon and shortage of rainfall. Accordingly, the efficient use of water and the management of these areas are integral component requirements of such areas. Precipitation pattern in arid areas occurs as storms with the high intensity usually in short time and causing huge runoff in the basin that will carry valuable surface soil and decreasing water infiltration opportunity in the soil. Different techniques of rainwater harvesting from natural lands, improved micro-catchments or agricultural micro-basins will lead to developing sustainable economical agriculture and store runoff into soils, dams, terraces, ponds, ditches, and reservoirs. The study area with 470 mm rainfall and runoff coefficient of 0.26 has the high potential of water harvesting at the basin scale. Previous studies show that lower gammas fuzzy logic values target macro-catchment systems of water harvesting well, but in micro-catchment systems with using fuzzy operators are less researched.
In this research for finding suitable areas for water harvesting, 3 systems including semi-circular bunds, contour furrow, and bench terrace with using GIS in Nazlu basin with 1507 km2 area, of west Azarbijan were investigated. Watershed factors such as climate, physiographic and hydrological factors, six important layers of precipitation, slope, land use, hydrologic soil groups, and drainage and soil depth were selected. For intersecting and making layers dimensionless, Fuzzy logic with trapezoidal membership function and Fuzzy logic with the analytical hierarchical process (AHP) were used for the quantitative and qualitative variables respectively. AND for combining of layers Boolean AND, Boolean OR, Fuzzy algebraic Product, Fuzzy algebraic sum, Fuzzy gamma of 0.1 to 0.9 methods were examined.
The “measure of similarity” is used for comparisons of filed observed and combination outputs. The results show for the crescent bankette and the Gamma-0.9 Fuzzy, this figure equals 0.824 and for counter furrow, equals 0.906. These figures indicate that fuzzy gamma 0.9 has better suitability than other method and also Boolean logic model do not show suitable results. Similarity index in fuzzy logic increases up to 0.9 gamma and returns below 0.9. High Gamma for site selection of Micro-catchment water harvesting systems shows better results. The 30% of the basin is suitable for semi-circular bunds, the 11.5% for contour furrow and 6.3% for bench terrace.
The results show that for determining the locations suitable to micro-catchment water harvesting methods, Boolean methods are less useful due to low precision and the usage of newer methods of gamma fuzzy logic should be considered, instead.

Garlic (Allium Sativum L.) is a herbaceous herbaceous plant of the Liliaceae family which has been ranked second in the world in terms of production of onions after onions. Due to lack of water resources, saving and optimizing water use is inevitable. Generally, the highest yield and water use efficiency are achieved in full irrigation and mild stress, respectively. And with the use of organic fertilizers such as vermicompost, the negative effects of water stresses are reduced. Considering the importance of applying water stress in irrigation programs and increasing the amount of effective material of plants by applying stress, this study investigated the effects of different levels of stress and vermicompost on yield, yield components and active compound of Allicin in garlic herb which is the second highest among the onion plants.
The experiment was carried out in a split plot design in a randomized complete block design with three replications at the Research Farm of Torbat Hadiriyah University during the cropping season of 2015-2016. Different irrigation levels including 100 (control), 120, 80 and 60% of water requirement were considered as the main factor and vermicompost fertilizer at four levels of 0, 5, 10 and 15 ton / ha as sub factor.
The results showed that simple and cross effects of drought stress and vermicompost on plant height, leaf length, fresh and dry weight of the plant, fresh and dry weight of the tuber, fresh and dry weight of root, leaf area, leaf dry weight and garlic yield were significant. The results show that with increasing drought stress, plant height, dry weight of bulb and leaf area decreased. Consumption of 5, 10 and 15 tons of vermicompost increased 25, 18 and 10 percent of garlic yield, respectively. The highest and lowest yield of garlic was obtained in 100% water requirement and 5 ton per hectare of vermicompost (with a yield of 10170 kg ha-1) and 60% water stress without vermicompost consumption (with a yield of 3048 kg ha-1), respectively. Drought stress reduced and vermicompost increased the amount of allicin.
Garlic is susceptible to drought stress and Water scarcity decreases and application of vermicompost fertilizer improves its performance. The application of vermicompost while improving the growth characteristics, yield and yield components of garlic and reducing the damage caused by drought stress during the cultivated season, the continuation of plant performance improvement in subsequent years, as well as the improvement of soil physical and chemical properties.

Every year, heavy rainfalls cause flood events and huge losses to life and property over the flood prone catchments of Iran. Heavy rainfall forecasting is an important step in development of a flood warning system. In recent decades, the Numerical Weather Prediction (NWP) models were widely used for weather prediction. Several operational centers, including the European Centre for Medium-Range Weather Forecasts (ECMWF), the National Centers for Environmental Prediction (NCEP), the Japan Meteorological Agency (JMA), and the United Kingdom Meteorological Office (UKMO) offer valuable operational numerical predictions at a global scale. Regional models were developed to admit the need for comprehensive and high-impact weather forecasting with higher spatial resolution. Weather Research and Forecasting (WRF) model is extensively applied for regional rainfall forecasting. The WRF modeling system is data assimilation system and a mesoscale forecast. It is aimed to progressive atmospheric research and operational forecasting. In this paper, the performance of the WRF model is evaluated for heavy rainfall simulation in Kan Watershed, Tehran.
Three domains were used in the implementation of WRF model. Horizontal resolution of domains are 27km, 9 km and 3 km respectively. The initial boundary which used to run the model has been downloaded from the National Centers for Environmental Prediction (NCEP) from Global Forecasting System (GFS). It is worth noting that the physics scheme of model was selected using the results of the previous research on the selection of the best physics for WRF model. The evaluation was conducted on the short-term forecasting. For this purpose, three heavy rainfall events occurred over the study area have been simulated using the WRF model. Precipitation forecasts were also downloaded from NCEP's Internet web site. Then, the heavy rainfall simulated by WRF model and presented by NCEP were compared to the observed rainfall.
The results showed that rainfall amount has been underestimated by NCEP forecasts and the time of precipitation events has not been correctly predicted. It is also observed that the WRF model is able to capture the heavy rainfall events, So that the error indexes (RMSE and MAE) significantly reduced compared to global model.
The WRF model increased the accuracy of precipitation forecasting compared to the global model. Thus, it is recommended to use the WRF model coupled with a hydrological model to development a flood warning systems in in the flash flood-prone watersheds.

Population growth as well as global fresh water shortage problem have necessitated the use of unconventional water resources in agricultural sector, as the largest consumer of fresh water, especially in arid and semi-arid regions. Salinity is a major environmental stress, limiting plant growth and yield. Determining the extent of salinity tolerance in agricultural and horticultural crops helps overcome the problem. The goal of this study Investigation the use of Irrigation by Caspian Seawater on Sweet Pepper (Capsicum annum) cv. California wonder under Greenhouse conditions.
In this study, five mixing use of seawater and fresh water is evaluated on yield and yield components of Sweet Pepper (CV. California wonder). the research was done based in completely randomized design including 3 replications as pot planting under greenhouse conditions in Gorgan University of Agricultural Sciences and Natural Resources during 2016. Research Station is located in north of Iran at 36° 51' N latitude and 54° 16' E longitude at the south-east corner of Caspian Sea and its height from sea level is 13.3 meters. Soil texture is silty clay. In this study salinity factors consist of 5 levels (0.5, 3.0, 5.5, 8.0 and 10.5 dSm-1). Data obtained (Fruit Fresh and dry weights, Fruit number, Yield, Water use efficiency) were analyzed using statistical software SAS. Ver. 9.0 and the means were compared using LSD range test at 5 % percent.
The results showed that effect of irrigation water salinity on fresh and dry weight of fruit, number of fruit per plant, yield and water use efficiency were highly significant at 1 percent levels (p>0.01). The result showed that the 10 percent mixture of seawater and tap water irrigation regimes compared to other regimes after control regimes had the highest Fruit Fresh and dry weights, number of fruits, yield and water use efficiency, respectively.
The result showed that increasing salinity decreased all characteristic measured, increasing of irrigation water salinity from 0.5 to 3.3 ds/m (10 percent mixing of sea water tap water) has resulted to redaction of fruit fresh weight and yield to 50.5 and 48.4 percent. The results indicate that sweet pepper plants are sensitive to salinity, mixing of sea water and tap water to irrigate the plants in this greenhouse is not applicable. The results indicate that sweet pepper plants are sensitive to salinity, mixing of sea water and tap water to irrigate the plants in this greenhouse is not applicable.

Nowadays, extensive attention has been paid to removal and remediation of heavy metals from polluted sites. Biosorption is a cost-effective method for removing heavy metals from aqueous solutions. In this study, biosorption of Cd has been examined by shrimp shell.
The adsorption of Cd on shrimp shell was studied as a function of pH (Cd concentration: 8.0 mg L-1) in the range of 3.0-8.0, and as a function of metal concentration (Cd concentration in the range of 15-100.0 mg L-1) using a 24h batch equilibration experiments. The amounts of Cd sorbed on the adsorbent were calculated from the difference between the metal concentration in the blank and equilibrium concentration of Cd in the solutions. Chemical species of Cd(II) in the solutions was also predicted using Visual MINTEQ, a computer code developed to simulate equilibrium processes in aqueous systems. Infrared spectrums of untreated and Cd-loaded adsorbent were obtained using a Fourier transform infrared spectrometer.
The results of pH dependent experiments showed that the percent removal of Cd on the sorbent increased with the increase of suspension pH. The experimental data showed that shrimp shell used in this study removed more than 95% Cd from solution at pH=7.5. The Freundlich and Langmuir isotherms were described satisfactorily equilibrium data. Maximum sorption capacity (qmax) of the sorbent was 5.1 mg g-1. Pseudo second order kinetic model provided a perfect fit for the kinetic data of Cd sorption onto shrimp shell. Biosorption kinetics of Cd consisted of two steps: An initial rapid step where adsorption was very fast. The second phase is the gradual biosorption stage before the Cd uptake reached equilibrium. The equilibrium was achieved in 30 min for initial metal concentration of 5 mg L-1 and in 240 min for initial metal concentration of 50 mg L-1 where 91 % and 92 % of Cd were adsorbed, respectively. Based on the estimates obtained by Visual MINTEQ, the Cd2 and Cd(OH) species were the dominant species present in solutions in pH ≤ 8.5 and pH≥ 8.5, respectively. Infrared spectrum analysis of adsorbent before and after adsorption showed that the N atoms of the amine functional groups played a major role in shrimp shell for adsorption of Cd.
Based on the equilibrium sorption and kinetic information obtained for Cd removal by this sorbent, it may be concluded that shrimp shell is an effective sorbent for Cd and may be considered for environmental remediation in future studies.

Land suitability evaluation is inevitable for soil protection. Selection of an appropriate land suitabiliyty method seems to be necessary due to the limitation of the FAO method as a conventional method. In this regard, multi-criteria desision methods (MCDM) such as electrothermics, can be helpful due the considering the weights of effective factors and experts oponions in land suitability evaluation.
In this research, Square Root Method and ELECTRE Tri multiple criteria method have bean also used for suitability evaluation in an area about of 5000 ha of Chalderan region located between 47° 19′- 47° 23′ East longitude and 39° 01′-39° 08′ North latitude for potato cultivation. In general 60 soil profiles were dug, described and sampled and after soil samples analysis, soil texture, gravels, calsium carbonate, gypsum, EC, ESP, slope, flooding, drainage and region climate in different profile were produced by weighted coefficients. Criteria weights via pair-wise comparison matrix method and ELECTRE Tri profiles and thresholds limits determined using requirements table, soil and land properties data. In two selected land suitability method, land production potential was calculated by multiplication of the soil indices and radiation thermal production potential.
soils of study region has been classified on Entisols, Inceptisols and Aridisols based on 12th Keys to soil taxonomy. Statistics comparison of allocated grades mean to suitability classes’ reveald that results of two methods were significantly different and ELECTRE Tri method had the best grade. Also mean comparing of soil and land properties amount via Duncan test showed significant difference between suitability classes. The ELECTRE Tri methods was the best model from the viewpoint of statistical indices including R2 (0.87) and RMSE (241.2). Also, Root mean square method has a R2 (0.79) and RMSD (197.6). Predicted potential production by Root mean squre method and ELECTRE Tri method was 9330 to 30200 and 10140 to 32350, respectively. The GMER index also indicated lowestimation of two selected land evaluation method. ELECTRE Tri suitability evaluation that implies high accuracy of this method due to determination of transition suitable limits, appropriate values of thresholds and criteria weights as well as usage of fuzzy logic principals on calculation processes.
ELECTRE Tri method, integrated into expert opinion, is better suited to the problem of land suitability for agriculture. The time saving during the development of the land suitability evaluation for the agriculture of the potato was considerable. The coherence of the obtained rsults with observed yield confirms the system effectiveness.

Pressurized flushing is one of the common methods in sediment flushing of dam reservoirs. In this method, sediments settled in the reservoir are released out of the dam together with water flow using the bottom outlet. For increasing the performance of pressurized flushing, it is attempted to fall down water level to minimum height of operation. Also, in the previous researches, the performance of different hydraulic and structural methods was evaluated in order to increase the pressurized flushing. Researchers have always tried to recognize the process of flushing and effective parameters and therefore, propose efficient procedures to increase the performance of flushing. Aim of this study is to investigate experimentally the effect of the number of rows, height and arrangement of submerged vanes in the enhancement of flushing reservoir.
To carry out the experiments of this research, a physical model including sediment reservoir, settling basins and water supplied tank was used. Submerged vanes were set up in front of bottom outlet with two heights (HV/D>1 HV/D Results showed that with decreasing the raw number of vanes in parallel arrangement, the performance of sediment flushing was increased. On the other hand, one-raw submerged vanes (double vanes) with relative height HV/D=0.5 and 2 have the best flushing performance with 173 and 231 percentages respectively compared with control model. Also, submerged vanes with relative height HV/D=2 have better performance than ones with relative height HV/D=0.5. The radial arrangement of vanes increased effectively the flushing performance. So that the flushing performance of submerged vanes relative height HV/D=0.5 and 2 increased about 11.33 and 4.1 times compared with the control model. In the radial arrangement, the flushing performance of vanes with relative height HV/D= 2 were less than ones with relative height HV/D=0.5 because of focusing the approaching inflow on the back of the side vanes.
Radial arrangement could increase the flushing intensity (especially with increasing the water depth of reservoirs). This condition is a relative advantage for flushing. In the other words, the flushing operation has more efficiency in the full condition of dam reservoirs and it is not necessary to fall down the water level of reservoir.

The hydraulic jump is one of the applicable phenomenons in energy dissipation of rapid flows. The hydraulic jump is one of the rapid varied flows that with converting from the super-critical state of flow to the sub-critical state in a short distance results in sensible energy dissipation and increase of the flow depth. In this experimental study the hydraulic jump is formed on the smooth bed and the bed with attached floating objects. A set of sphere objects with a density lower than the density of the water were attached to the bed, where the bearing length had flexibility. The goal of this study is to investigate the sphere objects bearing length effect on the hydraulic jump length, the relative secondary depth, and relative energy dissipation of hydraulic jump and finally comparing the obtained results with other related work.
Material and The experiments were performed in a flume with transparent walls, 8 m length, 35 cm width and 40 cm heights. In order to form the hydraulic jump, the height of the walls were extended up to 80 cm in the beginning part of the flume and a shut with 30 degree angle and the height of 40 cm was set up. Then, for modeling the floating objects the size of the applicable diameter of the objects was set to 4 cm. To analyze the effective variables, the dimensional analysis using the Buckingham л-theorem was applied. In total 30 experiments were performed with bearing length and Froude number as variables, where 0 cm, 1.5 cm, 2.5 cm, 3.5 cm, and 4.5 cm were considered as the values of bearing length and 5.1-8.3 was considered as the value interval for the Froude number.
Based on the performed analysis, one of the main achieved results is that the flexibility and the oscillation in the applicable objects as the fixator of the jump results in the increase of energy dissipation while increasing the bearing length, results in the more energy dissipation. However, this effect is more significant for the Froude number between 5.1 - 7.5.
The obtained results show that the energy dissipation using the floating objects in average is around 69% that is approximately 10.2% greater than smooth bed. Additionally, other characteristics of the hydraulic jump were a part of the obtained results based on which, the length of the hydraulic jump and the relative secondary depth, respectively, in average are 35.5% and 19.3% lower than the classic stilling basin. It was observed that the bearing length has a significant effect on characteristics of the hydraulic jump, where as a result of changing the bearing length between the Froude number interval, 5.1-8.3, the jump length and relative secondary depth, respectively, decrease in maximum 19.1% and 15.2% compared to the bearing length of zero that does not have oscillation and in maximum the amount of energy dissipation difference is 10.6% greater than the bearing length of zero.

Evapotranspiration (ET) is an important parameter for climatological and hydrological studies, as well as for irrigation planning and management. Evaluation of simple reference evapotranspiration (ETo) methods has received considerable attention in developing countries where the weather data needed to estimate ETo by the Penman–Monteith FAO 56 (PMF-56) model are often incomplete and/or not available. So the goal of this study, was to Evaluation and Comparison 30 different methods according to the standard method of reference evapotranspiration PMF-56 equation and determine the best equation, as a Replace to PMF-56, in a warm arid climate is the city and station of Zahedan.
In this study, 30 commonly used ETo equations that belonged to four groups: (1) pan evaporation-based methods, (2) temperature-based methods,(3) radiation-based methods, and (4) mass transfer-based methods were evaluated against the PMF-56 standard model; and the best and worst equations of each category were determined from the Zahedan station located in a warm arid climate.The assessed methods were: Cuenca, Allen and Pruitt, Snyder, Modified Snyder, Pereira, Orang, Raghuwanshi and Wallender, FAO-56 (pan evaporation-based), Schendel, four new types of Hargreaves equation reported by Droogers and Allen and Trajkovic, Blaney–Criddle (temperature-based), Makking, Turk, Jensen–Haise, Priestley- Taylor, FAO-24 (radiation-based), Dalton, Trabert, Meyer, Rohwer, Penman, Albrecht, Romanenko, Brockamp and Wenner, WMO and Mahringer (mass transfer-based).
Initially, in each group the monthly mean values ETo (using 30 model) were determined and compared their relative performance with respect to PMF-56 ETo estimates in the study.
The ETo calculated by the Snyder equation (R2=0.85¡ NRMSE=0.08¡ PE=8.30), Hargreaves- 4 (R2=0.93¡ NRMSE=0.04¡ PE=8.33), Turk 24 (R2=0.96¡ NRMSE=0.03¡ PE=6.52) and Dalton (R2=0.81¡ NRMSE=0.18 PE=12.89), pan evaporation-based, temperature-based, radiation-based and mass transfer-based respectively best matched the ETo estimates by the PMF-56 equation with the lowest errors rates.
In general, the comparative results showed that the mass transfer-based equations had the worst performances, while the radiation-based and temperature-based models (as Turk, Jensen–Haise, Hargreaves- 4 and Blaney–Criddle) were the best-suited equations for estimating ETo in this Warm arid climate (Zahedan). Considering the unavailability of full weather data for applying the PMF-56 model for estimation of ETo in many regions of the world and in different parts of the Iran, especially in developing countries, the results will be useful for choosing the simpler ETo methods in Warm arid climates. It is better to achieve real results for other places with similar climatic conditions to be assessed.

Soil moisture characteristic curve is one of the most important soil hydraulic properties, and shows the relation between soil suction and soil moisture content. Direct measurement of this curve is costly and time-consuming. Therefore, an alternative to measurement is to estimate this curve using the more easily available soil properties. On the other hand, many models have been presented for fitting to measued data of this curve. However, their fitting ability on different kinds of soil have been rarely investigated.
In this study, 30 measured surface soil samples with high variaty texture from the depth of 0 to 30 cm in different regions of Fars province have been used. Selected soils were classified into three textures of fine, medium and coarse. Then, 10 useful soil moisture characteristic curve models (Gardner, Brooks and Corey, Campbell, van Genuchten, Tani, Boltzman, Feremi, Fredlund and Xing, Groenevelt and Grant and Dexter et al.) were fitted on the data, and the parameters of each medel in each soil sample were determined, and then the most appropriate fitted model in each soil group was determined. To determine the best model, the combination of statistical equations and linear correlation was used, and then the value of mean square deviation (MSD) was calculated to evaluate the different models for fitting soil moisture characteristic curve.
The results indicated that the fitted models were better in medium texture group than the fine and coarse texture groups. The results showed that in fine texture group the models of Brooks and Corey, Campbell, van Genuchten, Fredlund and Xing, Groenevelt and Grant and Dexter et al., in medium texture group the models of Gardner, Groenevelt and Grant and Dexter et al., and in coarse texture group the models of Gardner, van Genuchten, Groenevelt and Grant and Dexter et al. were appropriate. Consequently, the results according to the MSD values indicated that the model of Dexter et al. was the most appropriate model for fitting to the measured data of soil moisture characteristic curve. Furthermore, the results of this study showed that the famous model of van Genuchten was good for futting to the measured soil moisture characteristic curve data. On the other hand, the models of Tani, Boltzman and Feremi were not appropriate.
The results of this study investigated that the model of Dexter et al. was better than the other models for three textural group (fine, medium and coarse) for 30 soil samples of Fars province. Also. similar results have been reported by other investigators for the soil samples in Gilan province.