جمشید یاراحمدی

Dust is one of the environmental issues that affect air quality and pollution, human health, soil fertility, and many social and environmental aspects of society. Wind erosion and sand dune movement have adverse environmental consequences. Identifying factors affecting sand mobility and predicting their condition in the future is necessary to control the wind erosion process. The purpose of this study is to identify and monitor the sources of dust and quicksand in the Jolfa, ShushQom and Gildir areas.

Analyzing wind climate data, identifying the sources of fine dust production, measuring and monitoring the production of dust and wind deposits were done using the RDD index. Also, the threshold speed and intensity of wind erosion were measured using a wind tunnel device. In order to investigate the effect of climate change, especially wind and rain, on the activity of sands and to predict the possibility of movement of sand dunes and sand dunes and the dust caused by them, Lancaster's global method was used.

The results of the Lancaster index showed that 10% of the statistical years were in an active state, which corresponded to the years 1996, 2000 and 2001. Whereas, the immobility of sand dunes (inactive) according to the Lancaster index has accounted for about 26% of the statistical years; 65% of the time in 2020, sand dunes were in an active or sand-moving state at the top of the dunes. The lands of the Jolfa micropollution center are mainly formed by the type of gravelly and sloping fan-shaped debris, and the soils of ShushQom and Gildir micropollution center are mainly a combination of mountains, hilly plains and flood plains. The results of wind tunnel tests on soil samples collected from the land surface of the Jolfa region showed that the erosion threshold of the collected sample is 8 m s-1.

According to the results of this research, abandoned rain fields, poor pastures, dry river beds and seasonal stream, etc., in the presence of environmental factors and accompanied by adverse climatic factors, are potential sources of dust and fine dust production, but the said lands in East Azarbaijan province are scattered and in certain conditions and times, they can produce a significant amount of fine dust.

The extended periods of drought, decrease in vegetation, and enlargement of the salt playa at Lake Urmia have resulted in the formation of local dust hotspots on the periphery of Lake Urmia in East Azarbaijan province. The objective of this study was to examine the spatial and temporal distribution of wind erosion in East Azarbaijan province using the DSI index. To achieve this, a set of weather data recorded over the past 30 years at specific codes in meteorological stations was compiled. Subsequently, the Mann-Kendall test was employed to assess the presence of a trend in the data series. Following this, the frequency of dust storms was calculated using the DSI index, and a spatial distribution map was generated in the GIS. The results of the Mann-Kendall test indicated that the slope of the trend for both local and extra-local dust storms in the province is on the rise. However, only the trend in extra-local dust storms is statistically significant at the 95% confidence level. The findings of the research revealed that the majority of the observed dust storms originate from outside the local area in East Azarbaijan province. Nevertheless, there has been a relative increase in the frequency of local dust storms in recent years, which can be attributed to the prolonged droughts resulting in reduced vegetation and the expansion of the playa of Lake Urmia. In this context, the highest frequency of dust storms has been observed at the Tabriz, Sahand, and Maragheh meteorological stations near Lake Urmia.

Dunes are one of the most important and sensitive landforms sensitive to wind erosion. Their activity level varies depending on the intensity of various climatic and terrestrial factors. Dust caused by wind erosion in the area of dunes is one of the most important dangers and serious challenges in arid and semi-arid regions of the world, which has been aggravated by climate changes in recent years. Due to the location of IRAN in the dry and semi-arid belt, it has been exposed to numerous local and extra-local dust systems. Dust storms have been a serious crisis in the country's provinces for the past few years and have brought dangerous consequences for the residents of these provinces in terms of environment, health, and economy. The prolonged recent droughts in East Azerbaijan and the drying up of Lake Uremia have made the threat of sand and salt storms in the coastal areas of Lake Uremia possible. The purpose of this research was to introduce the dunes in the Uremia Lake basin in East Azerbaijan province, where the mobility of dunes was determined using the Lancaster index, and the level of their future activity was predicted by carrying out a sensitivity analysis.

The dunes in the Lake Uremia basin in East Azarbaijan province were identified through the literature review and field visits, and then, using Google Earth images, their area was verified and then their spatial distribution map was prepared in the ArcGIS software environment. In the present project, the sand dunes of Shush Qom located in the village of Maqsoodlou and Gildir in the political area of Harris County are reported for the first time. Here, the mobility of sand dunes in the study area is calculated based on the Lancaster index. To calculate this index, meteorological data from the last 30 years of synoptic weather stations near the studied sand dunes have been used. Blaney Kridel method was used to calculate the annual potential evapotranspiration, and its calculations were done in the R program environment separately for each meteorological station. Also, the frequency of winds with a speed higher than the erosion threshold speed (6 m/s) at a height of 10 meters from the ground is considered. After calculating the amount of sand dunes mobility, the sensitivity analysis method has been used to reveal the possible effects of climatic change elements on the mobility of sand dunes.

The results showed that the sand dunes of Qom Tepeh, located on the edge of Lake Urmia and west of Tabriz metropolis, were active in 81% of the last 30 years, and in other words, they were not active in only 19% of the statistical period. Also, the results of the Lancaster index indicate the intensification of the mobility of the sands in the Qom Tepeh dunes in recent years. This fact can be justified in the form of the prolonged of droughts and the expansion of the dry areas of the Uremia lake bed; Also, field observations indicate the instability of the surface of these dunes due to the indiscriminate extraction of sand from these hills for use in the sand brick production plant and other agricultural and industrial uses. The dominant wind direction of Shush Qom region is south-southwest and winds with intensity greater than (>6m/sec) constitute about 12% of all observed winds at the Harris synoptic station. Based on the results of the Lancaster index, the dunes of the Shush Qom region have been active in half of the studied period; Of course, the short statistical period of the Harris station (last 12 years) has prevented a detailed analysis of the past behavior of this sandy area. In both study sites, the activity of the sands was mainly limited to the top of the dunes. Despite observing a direct relationship between the Lancaster index with the frequency of winds above the threshold of erosion and evapotranspiration and an inverse relationship with precipitation, the sensitivity of the Lancaster index to the decreasing or increasing changes of climate variables is not the same in the mentioned sites.

Sand dunes are one of the famous geomorphological forms in arid and semi-arid regions of the world, and their activity has intensified with climate changes in recent years. Therefore, investigating and monitoring the behavior of this phenomenon in terms of its impact on the emergence and escalation of the dust phenomenon is an undeniable necessity to reduce the effects of environmental crises. Even though today the Lancaster index is widely used in discussions related to the monitoring of sand dune activity in different weather conditions around the world, it is necessary to mention the lack of proper spatial distribution of weather stations and the short period of statistical observation. Also, the use of data from meteorological stations that are located further away from the sandy areas, makes the accuracy and correctness of the calculations of the Lancaster index doubtful, and the results should be evaluated with field observations and other valid scientific documents. On the other hand, in addition to analyzing the sensitivity of the Lancaster index, it is necessary to use the results of different climate change scenarios and appropriate downscaling methods in revealing the climate data changes affecting the Lancaster index in predicting the future behavior of sand dunes.

Groundwater depletion is one of the scenarios of climate change, which has become a global crisis today. Land use changes and land cover and rapid urbanization are the main causes. The decline of groundwater in recent years has caused land subsidence in many alluvial plains of the country. Uncontrolled and uncontrolled consumption of surface and groundwater resources, reduction of precipitation, concentration of consumption in some places (Imbalance between demand and water supply potential), inadequate cultivation pattern and lack of proper irrigation and digging of multiple wells and operation Irregularity of them in recent decades has caused a critical situation of groundwater resources in most plains of the country. To Groundwater level in most of the country's aquifers is constantly declining and the average annual decline over It has been 12 meters for the last 15 years. Decreasing the groundwater level of the plains increases the cost of water extraction and Increased energy consumption, reduced water quality and the emergence of land subsidence. As side effects, with decreasing volume The phenomenon of subsidence occurs in different spatial and temporal situations with different characteristics such as groundwater content, morphology, pedology and geology, which have different destructive effects. Shabestar-Sufian plain is one of the plains in the northwest of the country, which is facing many restrictions in terms of access to surface water and groundwater and is in a critical situation. In recent decades, the increasing trend of irrigated agricultural lands and overexploitation of groundwater resources has put the plain at risk of subsidence. Groundwater level drop is one of the most common causes of subsidence. Recognizing and examining ways to control it and prevent accelerating the process and expansion of subsidence phenomenon prevents the occurrence of hazards such as groundwater pollution, land slope change, well wall collapse, gradual sinking of masts and structures, changing the slope of rivers and roads Be. Risks that cause significant financial and human losses and on the other hand cause a waste of time and manpower, and pollution and destruction of the environment. In Shabestar-Sufian plain, there may be signs of subsidence in the form of sloping pipes, piping, shallow cavities, which is a warning for this fertile plain. If proper planning is not done regarding the exploitation of groundwater, we will soon see the destruction of lands, houses, water pollution, etc. In this regard,. The aim of this study was to identify areas prone to subsidence, for which the Fuzzy-AHP model was used. Information layers including land slope, altitude, geomorphological units, distance from the river, land use, lithology, groundwater level decline have been used as effective criteria. Expert Choice software was used to perform the calculations and after obtaining the final weights of each criterion, the obtained weight was applied on the layers according to the expert opinions in GIS software. Also, for the fuzzy model, fuzzy maps were combined using the fuzzy gamma operator. To modulate the very high sensitivity of the multiplication fuzzy operator and also the very low sensitivity of the sum fuzzy, the gamma fuzzy operator has been used and finally a map of subsidence prone areas has been prepared. Landslide risk areas are classified into five categories: very low (10.8%), low (14.2%), medium (7.22%), high (7.28%) and very high (2.23%). have became. About 52% of the southern, middle and western regions of the region are at very high risk of landslides due to excessive use of groundwater. In terms of spatial distribution, the highest risk of land subsidence is in the villages of Aq Kahriz, Dizaj Khalil, Shanglabad, Bagherabad, Kafi Malek, Kushk, Hallaji, Yousefabad, Ali Biglou, Haftcheshmeh and Qom Tappeh. Complete and common data were used between 2001 and 2017. Water loss map was prepared by idw interpolation method in gis environment. The results of groundwater status assessments in the region indicate that the study area in recent years has been faced with an increasing trend of groundwater level decline, so that the rate of water level decline has reached 4.7 meters in 18 years, a sharp drop in land subsidence in this region Has caused. . The calculated results are confirmed by examining the time series of water level drop in the wells of the region. Analysis of geomorphological parameters of the region has also shown that the same occurrence of subsidence has been in areas with a slope between zero to 5 valleys, have experienced less subsidence in the slope of more than 5 degrees or no subsidence has occurred in these slopes at all. It has a high subsidence in places with a height of less than 1500 meters, and with salt cover and alluvial materials, as well as areas with agricultural use. Areas close to the river were less prone to subsidence. Flat and smooth slope directions have been the main area of occurrence of this phenomenon. It can be concluded from the results that more than half of the 5% (51%) of the study are prone to high to very high subsidence, this means the critical situation in this plain.

Dust is a weather problem that can cause serious environmental hazards and cause adverse effects on agriculture, industry, traffic and people's daily life. In this study, sources of dust and fine dust production in East Azarbaijan province (eastern edge of Lake Urmia) were identified and investigated.

Progressive research from the point of view of the method in the category of quantitative research and data analysis using, identifying the sources of fine dust production using 6-year statistics (2017-2022), measuring and monitoring the production of dust and wind deposits using the RDD index, as well as measuring the threshold speed and intensity of wind erosion using the wind tunnel device was done.

The share of land use in micro-district centers of the province is 48.8% of agricultural lands, 19% of dried wetlands or lakes, 32.2% of degraded pastures, etc. From the point of view of morphology, micro-rounded centers of the province are about 2.53% of hilly plains and 97.47% of flood plains, etc. The micro dust production center on the eastern edge of Urmia Lake is the most active center in the province, which causes many problems for local residents and neighboring areas. This center of fine dust mainly includes rainfed agricultural lands, abandoned agricultural lands, dried up wetlands or lakes, and degraded pastures.

In East Azarbaijan province, some lands such as plowed agricultural lands, abandoned wet fields, poor pastures, the dry bed of Urmia Lake shore caused by water receding, the dry bed of rivers and seasonal waterways, etc. In case of adverse climatic factors and availability the presence of other environmental conditions can potentially act as sources of fine dust production, but the above-mentioned lands are scattered in East Azerbaijan province and in certain cases and times they are capable of producing fine dust in a considerable amount.

Nowadays, natural resources are exploited for the purpose of economical development in developing countries. Expansion of agricultural lands, supply of charcoal and fuel wood and wood production play an important role in forest degradation which affects biodiversity, soil conservation, the quantity and quality of water and the global climate conducted to the importance of forest conservation and reforestation. Therefore quantitative assessment of forests is required for conservation programs and forest monitoring is defined as a tool for sustainable forest management.Today, remote sensing techniques and satellite images can widely provide functional information in environmental studies. In this work, sentinel-2 satellite images with high spatial, temporal and spectral resolution were applied to determine the area, distribution and density of the Arasbaran forestsas well as other land use classes in the area.

Arasbaran area is located in Qaradagh mountainous region inthe north of East Azerbaijan province between 38°25′ 59 " N- 39°20′ 7.7 " N latitude and 46°09′ 18 " E- 47°16′ 5.3 " E longitude which covers an area of 551211 hectares and the deciduous forests of this area are known as the 11th Biosphere Reserved in Iran. The altitude varies from ca. 256 m tomore than 2000 m. the importance of the area is in having a rich flora (about 1334plant species) and unique vegetation among the vegetation of the country.
For the first time, the Sentinel-2 images with a combination of high spatial and temporal resolution were used to classify the land use of the area.The best band combination was found for bands 2, 3, 6, 12 and NDVI index. Land use classification included dense, semi -dense, sparse and very sparse forests as well as rangeland, agriculture, residential area-bare soil, garden and water was implemented using 9 different algorithms in a pilot area to find the best algorithm. 280 training sample points were collected from all different land use classes in the area.Consequently, supervised classification technique and Maximum Likelihood algorithm with the Kappa coefficient of 0.886 and anoverall accuracy of 89.6% was identified as the best classification method for the Arasbaran area.Accuracy assessment of the final map was done using ground control points and Google earth images with a total accuracy of 95%.Finally creating an error matrix with 880 ground reference test pixels revealed the accuracy indices.

The final land use map of the Arasbaran area based onthe Supervisedclassification technique and Maximum Likelihood algorithm was created.Based on the results, the accuracy assessment of the final map showed that the Kappa coefficient and the overall accuracy of the classified map were 0.88 and 89.8% respectively.The forest distribution and canopy cover density map were extracted from the land use area map. The total area of forests with a canopy cover of more than 5%, obtained 131019 ha consisting of 39% dense forest, 36% semi -dense forest, 17% sparse forest and 8% very sparse forest. In addition, the largest type of land use accounted for rangeland with 270000, forest with 131019, agriculture with 101974, residential area-bare soil with 30028, garden with 15434 and water with 2756 hectares respectively. Based on the error matrix table and correct classified points as well as total ground control points, the highest user’s and producer’s accuracy belonged to the densed forest class as well as the lowest user’s accuracy and lowest producer’s accuracy belonged to garden and agriculture classes respectively.

The results conducted supervised pixel-based image classification based on the Maximum Likelihood algorithms an acceptable method. It can be because of well -distributed training sample points, the high spatial resolution of sentinel-2 images or Environmental heterogeneity of the area. According to the results, dense forests declined(from 56910 to 50628 ha)however semi -dense and sparse forests have increased (from 35280 to 47930 ha)with respect tothe last forest survey project in the Arasbaran area in 2003.In addition, the results revealed an overlap between agriculture and garden as well as rangeland and residential area-bare soil classes because of multi culture of crops and fruit trees together as well as dried or low vegetation cover of rangelands in the area. These results can provide useful information for decision- making and sustainable forest management for reducing forest degradation and it seems to be an important next step to manage these forests based on conservation policies.

The aim of this study was to evaluate the accuracy of the climate suitability map of pistachio planting in the northwestern provinces of the country. Accordingly, the climate suitability map of pistachio planting was prepared based on long-term meteorological data using weighted overlap modeling in GIS environment. The accuracy of the mentioned map was evaluated by comparing it with the pistachio orchards map that located in this area using the Receiver Operating Curve (ROC). The results showed that 33, 33.4, 60 and 12.4% of the study area in East and West Azerbaijan, Kurdistan, and Ardabil provinces are suitable for pistachio cultivation respectively. The area under the curve (AUC) for the mentioned provinces were calculated as 0.9, 0.8, 0.9, and 0.62, respectively, which in general indicates the high accuracy of the results. In Ardabil province, insufficient attention to the limitation of high relative humidity in different stages of pollination, and fruit growth in the construction of pistachio orchards has caused a low area AUC. Citing the results of the climate suitability of pistachio cultivation map can be a good guide for site selection of new orchards and therefore, the reduction of crop yield due to climatic hazards will be minimized.

Changes in erosion and sedimentation of the basin are one of the most important factors that affect different parts of human life and natural life. it is very necessary to receive these changes quantitatively, which mainly take place under temperature fluctuations and climate changes in different regions, in order to be more prepared to deal with its negative consequences. In this research, erosion and sedimentation changes in Hajiler watershed were investigated and predicted using GeoWEPP and SWAT models. Based on this, first, by using the data of the current situation of the Ahar synoptic station and using the SDSM model, the changes of the statistical period2020-2040 in three scenarios RCP2.6-RCP4.5-RCP8.5 were investigated, then simulation and prediction of erosion changes was carried out. and sedimentation was done under the influence of climate change by using popular models. The output of the SDSM model indicates an increase in temperature and a decrease in rainfall for the basin until 2040.And the analysis of the simulation results of the sedimentation rate of the models showed that in the studied basin, the GeoWEPP with the selection of the domain method has a suitable level in estimating the sedimentation rate compared to observational statistics. The final model was chosen to predict the amount of sediment in the mentioned period of the basin. Using the downscaled results of the atmospheric general circulation model, the sediment changes in the statistical period of 2020-2040 under the above mentioned three scenarios were estimated as -1.97, 4.45, and 2.98, respectively.

This research aims to introduce the water supply method required for pistachio orchards in semi-arid climatic regions using the rainwater harvesting technique. For this purpose, a rain catchment system with an area of 1000 m2 was constructed at the Tasuj Station site. Based on the mean annual rainfall of the region, runoff threshold, runoff coefficient for insulation surfaces and the area of the catchment system, water storage pool was designed and constructed with a nominal capacity of 280 m3. In addition, a geomembrane cover is used to insulate the surfaces of the rain catchment system and the water storage pool. To reduce agricultural water consumption and increase productivity, the drip irrigation method with gravel filters was used to irrigate pistachio trees from harvested rainwater. The results showed that from the total rainfall of 267 mm, about 227 m3 of water was harvested from rainfall in the year 2021. According to the findings, the per capita water consumption of each fruitful pistachio tree has decreased from 8.65 m3 in the usual drip method to one m3 in the drip irrigation method with gravel filters; without having a negative effect on production yield. The simultaneous use of rainwater harvesting techniques and drip irrigation methods along with gravel filters will reduce the pressure of groundwater resources exploitation and increase the irrigation water efficiency. The results of this research were introduced to the farmers of the region in the form of educational-promotional workshops, which have been well accepted.

The issue of global warming, climate change and drought are among the major challenges facing the world today, which can cause widespread fluctuations in the Earth's climate. Areas with Mediterranean and semi-arid climates are highly dependent on temperature and precipitation and as a result are affected by climate change. Iran is one of the countries in the semi-arid and arid regions of the world that is more sensitive to climate change. Global warming due to increasing greenhouse gas concentrations and land use change has caused obvious changes in Iran's climatic parameters, including increasing temperature, decreasing rainfall and increasing the incidence of destructive atmospheric-climatic phenomena in the country. This research has been carried out with the aim of detecting future changes in temperature and precipitation in Liqvan watershed. For this purpose, statistical Downscaling of SDSM model and CanESM2 model under RCP2.6, RCP4.5 and RCP8.5 climate change scenarios on long-term data of Tabriz Synoptic Station (1951-2020) has been used. The results show that in four periods of 20 years (2020-2100) and based on the three existing scenarios, the temperature will increase and the precipitation will decrease. This temperature increase for the minimum temperature will sometimes be up to 14.35 ° C (January with RCP8.5 scenario and time period (2100-2081). Examination of rainfall in four periods and three scenarios shows that the amount of rainfall in October, November, December and April will increase and in the remaining months will decrease.In the present study, climatic parameters of temperature and precipitation were simulated using multiple linear model of SDSM and general circulation models of barley using data of Tabriz city for Liqvan watershed. In this research, the output of canESM2 model under scenarios RCP8.5, RCP4.5, RCP2.6 has been used for future periods.The results showed that temperature data correlated better with observational data (compared to rainfall data), this is because temperature variability is less than rainfall and temperature is a parameter with a normal possible distribution. One of the reasons for the decrease in rainfall correlation is that different factors affect rainfall and on the other hand, rainfall is a discrete variable. These results are consistent with the results of Sajjad Khan et al. (2006), Sarvar et al. (2010) and Nouri and Alam (2014). Therefore, solving the correlation problem in the development of future climate change models should be considered.Also, the results of this research are in line with the results of most researchers such as trainee et al. (2009), Sajjad Khan et al. (2006), Sarvar et al. (2010), Mino et al. (2012), Chima et al. (2013), Dehghanipour et al. (2011). And Shakeri et al. (2021) agree that the SDSM model has a good ability to small-scale temperature and precipitation data.Climate change can cause temporal and spatial changes in climate variables. The characteristics of these variables can have detrimental effects on ecosystem components. According to the results, it was found that during the 21st century, temperature is increasing and precipitation is decreasing.In Tabriz station, in general, precipitation will decrease in the three scenarios studied and in only one scenario and for the period 2100-2080, there will be an increase in precipitation. Also, rainfall will generally increase in winter and the rest of the seasons will decrease. These results are consistent with the results of Golmohammadi and Masah Bavani (2011) which introduced the period 2069-2040 as a period with increased rainfall in Qarasu basin but with the results of Rajabi and Shabanloo (2012) using the SDSM model in Kermanshah region in Western Iran has considered the period 2070-2070 to be drier, there is a difference, so that according to the results of the present study at Tabriz station, from May to September, we will have a decrease in rainfall for all periods and scenarios. Figures 8 and 9 show the changes in monthly and annual rainfall in the RCP scenarios and the four periods compared to the base period, which show a decrease in precipitation in most months except October, November and December compared to the base period.Changes in the average minimum temperature of Tabriz station in all months except October, November and December will increase in future periods. Figures 11 and 12 show the monthly and annual minimum temperature changes in the RCP scenarios and the four periods compared to the base period, which show an increase in the minimum temperature compared to the observation period, which is more than the other scenarios in RCP 8.5. .

Floods and their consequences, with the intensification of human exploitation of nature in the early twentieth century, have had negative effects on vital ecosystems. Also, adverse effects of erosion, while destroying the harvest site, lead to reduced production capacity and degradation of physical and chemical properties of soil in lands. To manage this phenomenon, the factors of production and flood must be identified and then areas with high potential in flood production must be identified to enable the possibility of executive and corrective operations at smaller and risky levels. The level of flood areas in the country is estimated at 91 million hectares, of which about 42 million hectares have moderate to very high flood intensity. Therefore, knowing the flood situation of the regions is a necessity to prepare strategic plans for sustainable management of basins. The purpose of this study is to estimate flood potential and determine the priority of flood areas by physical factors to combat erosion in Hajilar watershed by using a combination of data and information based on the CN-SCS method and also to determine critical flood areas.

In this study, in order to investigate the potential and zoning of flood risk in the basin, at first, by studying and examining the foundations and theoretical background of the subject, the physical factors affecting the occurrence of floods were identified.Then, the required information was collected and layers of each of the proposed factors were prepared in Arc GIS 10.7 software and Arc-Hydro and Arc CN-Runoff extensions. In this regard, the information layers of the waterway network, level lines, and elevation classes were prepared using a digital elevation model with a scale of 1: 120,000. Lithological information layers were obtained using the geological maps of Siah Rud, Tabriz, belonging to the Geological Survey of Iran. The precipitation map of the basin was prepared using data from meteorological stations within the study area and also adjacent stations using IDW interpolation method. Land use layers of the area were obtained using the area land use map and monitoring using Sentinel 2 satellite imagery. The soil map of the region has been prepared using studies of natural resources of Arasbaran basin prepared by the Forests, Rangelands and Watershed Management Organization of the country and controlled with lithological conditions and other environmental factors. Vegetation of the area was prepared through NDVI index. Then, by extracting the number of curves (CN) and the amount of penetration (S), the layers were combined and the runoff height of the basin was estimated by the runoff curve number (CN-SCS) method, the units were classified by SPSS software and the priority of areas in terms of floods in the basin was determined. In order to determine the peak discharge of flood through the obtained data, first based on the proposed Schwab relation, the water accumulation time of all sub-basins was prepared separately and, finally, the maximum peak discharge obtained from this runoff was obtained.

The results of the SCS-CN model which was intended to determine areas with different flood potential in the basin and compare it with environmental factors such as slope, lithology, land use, rainfall, etc in the basin indicate that this model is highly capable of estimating the flood potential in different areas of the basin. Results of changes in the number of curves based on effective environmental factors were also presented in Table 1.The potential for runoff production was high in high altitudes with poor pasture landuse or dryland agriculture with poor permeability soil, as well as in residential areas where the city surface consists of impermeable or low permeability surfaces.Also, there was a high flood potential in the  dense and medium pastures with hydrological group B. Irrigated and rainfed agriculture with hydrological group A has the lowest runoff potential and thus has the lowest amount of flooding in the basin.By obtaining the runoff height in different parts of the basin and determining the levels of changes in the values through quartering, it was found that mainly the central and lower parts of the basin are in the first priorities in terms of flood potential. Poor coverage, low permeability, and low rainfall in this area are some of the factors that increase flood potential. The southern and central regions are areas marked by low sensitivity to flooding. By prioritizing areas in terms of floods and mapping them, the results can be used in watershed management operations at the level of high-sensitivity units to reduce erosion and damage. Taking into account the values of runoff height relative to the sub-basins and water accumulation time obtained based on the Schwab estimation method, the results showed that the flood volume and maximum peak discharge in the basins have a good relationship with each other and the highest maximum peak discharge is related to H33 and H18 sub-basins with volumes of 51.44 and 48.67 and the lowest is related to H12 and H29 sub-basins with volumes of 3.77 and 3.86 m3/s.

Floods are among the most important environmental hazards that cause human and financial losses every year. In the meantime, using the SCS-CN experimental model, as a method for estimating floods in basins with different environmental conditions, and the inter-environmental approach in it can be a useful solution in watershed management studies. Considering that human intervention has caused an increase in floods in all areas, providing methods for accurate flood estimation is one of the basic needs of the relevant responsible organizations. Accordingly, the present study was conducted to determine the priority of areas in terms of flood potential and the results showed that 9% and 51% of the basin areas are in the very dangerous and high-risk flood categories, respectively. According to the final map obtained, Areas with very high risk and high risk are mainly located in residential areas and in the lower areas of the basin, which are the first priority in programs related to water resources, especially flood control and watershed management in the upstream units. Also, according to the results, areas with medium risk potential occupy 23% of the area and 17% of the area has low risk potential. The results of the study confirmed the high potential of the studied area in terms of flood risk, so lands with very high and high risk are lands that should be protected and appropriate watershed management measures must be conducted to control the speed of floods and reduce soil erosion.

Land subsidence is one of the natural hazards that occurs as a result of natural factors or a combination of these two factors, in the form of flooding. Monitoring land surface changes requires the use of precise techniques. In this regard, the radar interferometry technique with permanent dispersers with wide spatial coverage as well as high temporal and spatial resolution is one of the most accurate and cost-effective remote sensing techniques for monitoring and measuring the amount of displacement on the ground. In recent years, following climate change and successive droughts on the one hand, and the improper management of water resources, improper abstraction of groundwater and increasing growth of population, on the other hand, has led to land displacement - and subsidence, in particular - on the Shabestar-Sufi plain. Therefore, this study has evaluated and monitored the rate of land movement in this region by using the mentioned method and by processing Sentinel-1 images between 2016 and 2020. To study the changes in groundwater level in the region and to prepare a water drop map as well as the hydrograph of the observation wells in the region, water level data between 2001 and 2017 have been used. The results of the permanent dispersion method indicate the annual land displacement rate for the region from 2016 to 2020 in five periods - 8.54, -9.47, 8-9, -9.7 and -9.02 cm, respectively. Due to the lack of geodynamic and GPS stations in the region, validation of the results of the radar interferometry method was assessed by using piezometric well data and the predominance of groundwater level drop (4.7 m) in the region and hydrograph analysis of the aquifer unit, and by comparing the results of the previous research, and the accuracy of the results of the present study was confirmed.

To prevent the occurrence of floods, it is necessary to identify areas having high potential for this phenomenon. The aim of this research was to estimate runoff coefficient, understanding the effective factors in flooding, and the flood hazard of Arasbaran forests (Sotanchai basin). For this purpose, first, the required data, including the maximum daily precipitation, soil, slope, land use, vegetation density and drainage density layers, were prepared in geographic information system (GIS) environment. In the following using the analytic hierarchy process (AHP) method, an effective weight was assigned to each of the six factors influencing the flood hazard of the basin. Through applying the weight of each layer and combining it, the flood hazard zonation map of the basin was prepared with different return periods (10, 25, 50, and 100 years). The results showed that 15 - 19% of the basin area had a medium flood potential during the period of 10 to 100 years and 1-13% had a high flood hazard. During the 10- and 25-year return periods, no area of the basin was in the flood hazard category, but in the 50- and 100-year return periods, 0.32 and 0.96% of the basin area were in the very high flood hazard category, respectively. Generally, it can be concluded that due to good forest density and vegetation, flood hazard was low in the region, instead areas with high flood hazard were located in areas, where forests have been destroyed.

Today, the use of underground dam technology is considered as one of the efficient methods of rainwater harvesting systems and flood management in order to eliminate the shortage of water resources in arid and semi-arid regions. Underground dams block the subsurface flow of water across the river and stores it in the riverbed. The construction of an underground dam, while helping the environment, will compensate for surface water and expand cultivated lands. The content of water extracted from these dams depends on the hydraulic and hydrodynamic of alluvial sediments. The aim of the present study was to determine the hydrodynamics of the riverbed deposits of the Meshnaqeh underground dam in the north of Lake Urmia. Feasibility study and location of the dam axis was done using ArcGIS system. Three options were determined and finalized after a geophysical evaluation of the location of the main alluvium with a thickness of 14 m. Engineering studies showed that the sedimentation of bed was in gravel and sandy soils in UNIFIDE division with more than 78% of sand, and according ASTM method, were in the GC and GM. Lefranc experiment showed that alluvial sediments had a permeability of 0.0012 cm/s and a porosity above 35%. According to the special discharge, 163,800 m3 of water can be harvested in each flooding period. Drinking water and sanitation of the region is provided hundreds of kilometers away through Zarrinehroud water transmission line. The construction of this dam with a realistic view of the real value of water and its role in human societies, will have constructive environmental effects requiring independent study.

Introduction Climate change is one of the most important phenomena of the present century, which has created many problems and challenges both globally and regionally and nationally. In the second half of the twentieth century, global warming relative to The first half of this century has increased and it is predicted that this increase in temperature will continue in future periods, resulting in changes in the level of climatic conditions in different parts of the world. Due to the lack of atmospheric precipitation, due to the increase in temperature, the rate of evaporation has increased significantly and can greatly affect the aggravation of water shortage conditions in surface currents, especially evaporation from the surface of lakes behind dams. Percentage by evaporation leads the country to higher values ​​(Farajzadeh and Ghasemifar, 1398). Regarding the changes in Iran's water resources in the horizon of 2100, few studies have been done and most have been case studies (Fahmi, 1393). Although the results of these studies, based on the climatic models and different scenarios used, sometimes show contradictions, so it is necessary to do more studies in this field. 
  Methodology The present research has been done in three specific sections and the output of each section has been used as the input of the next section. In the first part, climate change in the form of precipitation variables in the study area is detected and subsequently, rainwater runoff in the Daryan catchment is simulated. Then, while identifying the characteristics of hydrological drought periods in the basin, the probability of occurrence, intensity and duration of hydrological drought periods are calculated based on the fit of different statistical distributions for different return periods in the third section. 
 
Climate change is one of the most important environmental problems of this century. Thus, evaluating the phenomenon of climate change and reducing its effects on both global and regional scales has attracted the attention of many researchers, planners and legislators (Yohe et al., 2007). Proper assessment of these effects requires the existence of climatic information with appropriate spatial distribution and long-term time series, as well as a thorough understanding of its future trends at the regional and local scale. Despite the fact that today the output of public circulation models (GCM) are the main sources of future climate data production. One of the most important consequences of climate change includes changes in the hydrological cycle and river flow regime of watersheds. Therefore, the present study aimed to investigate the possible effects of climate change on rainfall and runoff in the Daryan catchment area north of Lake Urmia. In this study, statistical method (SDSM) and data of CanESM2 Canadian climate model in the form of three scenarios RCP2.6, RCP4.5 and RCP8.5 in order to micro-scale the precipitation data of five synoptic stations adjacent to the sea basin and changes Its future is used. Here, the basic period (1961-2005) and future periods (2049-2020), (2079-2050) and (2080-2100) were selected. In this research, the threshold level method has been used to identify hydrological drought periods and extract its characteristics. The results of the analysis of the last 35 years of hydrological droughts in the Daryan Basin showed that 44 drought events occurred in this basin, which in total, led to a reduction in surface flow volume of about 140 million cubic meters in this basin.
Conclusion The simulation results of SWAT model showed that the annual average runoff of the sea basin in the first period (2020-49) in all three scenarios increases by 3.7 and 6%, respectively, compared to the base period. While in the rest of the periods of all scenarios, runoff reduction is predicted compared to the base period. Accordingly, a decrease in surface runoff compared to the base period is predicted for five months of the year (April to August) and an increase in the remaining months. Future changes in precipitation at Tabriz station, which is the basis for modeling runoff in the Daryan basin, are not very noticeable compared to the base period, and only in the period (2049-2020) all three scenarios are predicted to increase by 5, 2 and 8%, respectively, compared to the base period. In the other periods, in all three scenarios, a decrease in rainfall is predicted compared to the base period. Results of evaluating the effects of climate change on rainfall and surface runoff in the Daryan Basin with the results of other researchers in the catchment area of ​​Lake Urmia, including: Goodarzi and Fatehifar (2010) in the Azarshahrchai Basin, Qaderpour et al. (2016), Dariane et al. (2019) ), Sobhani et al. (2015), Goodarzi et al. (2015) and Salehpour and Malekian (2019) are consistent.

Since the most consumed water is in the agricultural sector, available water resources, especially in arid and semi-arid regions, are the main limiting factor for agricultural production. Therefore, by rain water harvesting systems and its proper management, it is possible to compensate part of water shortage in the mentioned region. The purpose of this study was to select the most suitable rain water harvesting techniques based on multi-criteria decision-making system using Analytic Hierarchy Process (AHP) and investigation of its effectiveness in supplying the water demand of pistachio gardens. The results of the AHP analysis indicated that the environmental impact criterion with a weight of 0.342 and permanency criterion with a weight of 0.284 were the most effective criteria for selecting of the best rain water harvesting system. Finally, the rain water harvesting system with geomembrane insulator with a weight of 0.371 was selected as the most suitable insulator cover, and two other rain water harvesting systems with plastic and isogam insulator were weighted with the weight of 0.350 and 0.279, respectively in the following priorities. Also, the results showed that from the insulated surface with an area of 800 m2 and an average annual rainfall of 300 mm, it is possible to harvest at least 200 m3 rainwater per year. With this volume of harvested water and the choice of drip irrigation with gravel filter, the annual irrigation water demand of 200 fruitful pistachio trees will be supplied during a crop season. Therefore, due to rainwater harvesting techniques, some parts of irrigation water demands of gardens will be provided in arid and semi-arid regions, as well as maintaining sustainability of production and increasing gardeners' revenues, the double pressure on surface water and underground water will also be reduced.

Landslides are one of the most common and dangerous threats in the world that generate considerable damage and economic losses. Landslides detection and monitoring are two important research aspects of landslides analysis. There are different geodetic and non-geodetic methods to measure slope instability. Geodetic methods includes ground observations via GPS, Total Station and laser scanners. Observation based on mentioned methods provide accurate and continuous measurements at limited points in unstable regions, but none of them has the ability to determine the extent and pattern of spatial unstable regions. On the other hand, the repetition of each of these methods, especially when it is intended for a large area, is very costly and time-consuming. The existence of such constraints has always been one of the fundamental challenges faced by researchers in relation to precise measurements and spatial monitoring of land surface changes. Due to the fact that landslides directly affect the surface of the earth, the use of remote sensing techniques in instability studies of the slopes seems to be very suitable. Among the remote sensing techniques, radar interferometry (InSAR), capable of working in all weather conditions and the duration of night and day, is one of the most effective and efficient techniques for detecting and monitoring the steady change of the earth's surface. Of course, the lack of spatial and temporal correlations in available radar images limits the use of conventional radar interferometry to monitor the displacement of ground level. Persistent Scatter Interferometry (PSI) is the advanced InSAR technique which has significantly improved upon traditional InSAR methods by increasing the accuracy of results (millimeter scale precision). The main objective of this project was to detection and monitoring of Ahar Chai sub basin landslides based on PSI technique.

This landslide with an area of about 42 hectares is located near the village of Dupiagh, located in Ahar city and 22 km on the Ahar-Varzaghan axis, in East Azarbaijan province, Iran. The geographical coordinates of the study area is N38 29´ and E46 49´. For detection and monitoring of the landslide, the PSIn-SAR method was implemented on 22 ASAR images (with descending orbital modes) recorded between October, 20030502 and 20100709 by ENVISAT satellite in VV Polarization. For detection and monitoring of Dupiagh landslide, the PSInSAR method was implemented on mentioned ASAR images. SARscape5.2 software on the platform of ENVI 5.3 was employed to process the radar images and to extract the persistent scatterers. The PSInSAR processing technique was performed within five stages including: (1) selection of master image or connection graph selection. (2) Interferometry or interferometric workflow stage consisting of co-registration and differential interferogram generation parts. (3) First step inversion for selection of the candidate persistent scatterer points. (4) Inversion: second step including phase unwrapping and filter implementation stages. (5) Geocoding or conversion of the phase into displacement. Finally, The PSInSAR method results compared with dual-frequency -GPS measurement.

Based on the ASAR images processing by PSInSAR method, results showed that some parts of Dupiagh landslide are still active and displacement rate of the this landslide was obtained 12.4 mm/year duration of 2003 to 2010. While, the ground surface displacement velocity was estimated between 58 to -22.5 mm/year in the Ahar Chaiy sub basin duration the mentioned period. The results of GPS measurement indicated that the Dupiagh landslide was inactive duration of observation baseline. Differential RADAR Interferometry (DInSAR) method in detection and monitoring of landslides has been reported in numerous studies. Also, the study suggests that PSInSAR is a powerful technique to determine displacement and spatial pattern of landslides.

Droughts are intrinsically regional. The critical conditions in the dry periods occur when there is a severe shortage of water resources for a long period on a wide area. Therefore, an investigation of the regional and not of a local drought offers better understanding of this phenomenon. This study aims at a regional analysis of the hydrological drought in East Azarbaijan Province for a period of 45 years (1972– 2017) using the data collected at 25 hydrometric stations with proper spatial distribution across the province. Homogeneous hydrological areas was identified based on hierarchical analysis method in the SPSS software using the following independent variables: Area, Perimeter, Length of the main basin, Main channel gradient, Drainage density, Basin shape factor, Total hydrographic network, Weighted average slope of the basin, Maximum, Minimum and Median height of the basin, and its length. Frequency analysis and probability distribution of the dependent variables (duration and intensity of hydrological drought spells) was performed on the basis of the two year return period for each hydrometric station. The factor analysis was performed using the PCA method in the SPSS software. The suitability of the data for the factor analysis was evaluated using the KMO coefficient. The value of this coefficient was estimated as 0.546, which indicates that the data were suitable for factor analysis. The result of factor analysis showed that the contribution of perimeter, height and shape of the basin were 40.13, 26.42 and 13.12% (Sum 79.68 % ) respectively in explaining the variance of the variables (i.e. the most important factors in zoning). Results of regional analysis using multivariable regression (stepwise method) indicated that perimeter, avrage height and the basin shape coefficient have a decisive role in estimating the drought severity and its duration with two-year return period in East Azarbijan. In the scarcity of hydrometric stations in most of the watersheds of the East Azarbaijan province, the regional equation presented in this study can be used to estimate hydrological drought characteristics of watersheds without observed surface runoff data.

Hydrological droughts occur after meteorological and agricultural droughts. In other words, this kind of drought is the final stage of the drought cycle and is affected by global climate change. Nowadays, many studies on hydrological droughts are based on the Threshold Level method, which is based on the Run Theory. According to this view, a drought is described as the length of a period in which hydrological variables, or discharges, are below the specified threshold level. To the researchers' knowledge, there have never been a major research on hydrological droughts in Iran and studies in this field are very limited in comparison with other types of droughts. Most of the researches in hydrological droughts have been done in the monthly and annual time intervals. However, the present study was conducted to investigate the hydrological droughts on a daily basis with the objectives of identifying and monitoring them and determining their occurrence and severity in the studied area.

In this research, the Threshold Level Method was used to identify the hydrological droughts in the DaryanChai sub-basin in the northern part of Lake Urmia. Using the computer program NIZOWKA2003 based on the partial time series (PDS), the droughts in the hydrometric station of the Daryan during a thirty one year period (1982-2014) were calculated. Hydrological drought characteristics including time of occurrence, duration, severity, and minimum observed flow during dry period were calculated based on the NIZOWKA2003 software. In addition, frequency analysis, probability of the number of drought occurrences, duration, and volume deficiency in different return periods were calculated using the mentioned computer program. Here, the threshold level was chosen based on the flow duration curve (FDC), which was based on the daily flow rates of the selected hydrometric station and prepared using the NIZOWKA2003 program. The IC method was used to integrate the minor and mutual dependency periods of hydrological drought. Based on the characteristics of dry periods, the probabilistic distribution and frequency analysis of dry periods, the probability of the occurrence of drought and characteristics of dry periods (duration and volume deficit), fitting different types of statistical distributions and selecting the most suitable ones based on the fitting values (chi Score) were calculated using the NIZOWKA2003 computer program.

The results of the calculations showed that a total of 38 periods of a hydrological drought occurred in this river. Approximately about 20 events from all detected drought occurrences lasted more than 200 days. The longest period of drought with 577 days lasted between the end of June 1988 and March 1990. Less than two months later, the second largest hydrological period with 365 days lasted between the beginning of May 1990 and the end of June 1991. Regardless of these two months, the largest hydrological drought period with 950 days was between 1988 and 1991. The severest hydrological periods were respectively 8987 and 6133 cubic meters, coincided with the largest hydrological drought periods.  Moreover, the calculations showed that this river's water volume has fallen by 117 million cubic meters over the period of 31 years.   Based on the results of Chi-square test and Akaya's criterion, Poisson distribution was the most appropriate distribution of the probability of the occurrence of drought occurrence. Among the various probabilistic distributions, Weibull's distribution was the most appropriate distribution of the duration of the hydrological drought events. Accordingly, the probability of the occurrence of a dry period of at least 13 day duration is likely to occur in this river. Additionally, the maximum observation dryness event of 578 days was estimated with a probability of non-exceeding 0.99 % at the Daryan hydro-station. The probable distribution of Weibull was selected as the most suitable distribution of the hydrological drought severity. Based on the Weibull method, the maximum water deficit volume was estimated at 13730 cubic meters with a low probability (0.01%). In contrast, the highest probability of the severity of drought (about 70%) in the Daryan Chai was estimated at 253,000 cubic meters. The results of the calculation of the return periods showed that in half of the cases, the probability of the occurrence of a dry period was with a duration of 170 days and an intensity of 2276 cubic meters with a 2-year return period.

Until the last decade, little research has been conducted on the hydrological droughts in Iran.In very few studies, the analysis of the low flows from the point of view of frequency analysis has been addressed.The lack or shortage of hydrometric data in many parts of the world has been the main limitation in hydrological drought studies. In this research, the Threshold Level Method was used to extract thehydrological dry periods and their characteristics. Despite the fact that this method has been extensively used in drought studies, its selection is still one of the controversial issues in such studies. Using the daily time series can provide complete and accurate information from the start and end dates of drought events. However, the use of this time basis in drought studies, especially hydrological droughts, is not common in the world and in Iran.The only studies in Iran, based on the daily time series, were conducted by Bayzidi and Saghafian (2011), Mesbahzadeh et al. (2017), and Mostafa Zadeh et al. (2018).  Almost all other studies have been based on monthly or annual basis. However, the use of a daily time series is often associated with the problem of having minor andmutual dependency periods. There are several ways to overcome this problem. The Interevent Critria. (IC) is one of these methods which was used in this research.The efficiency of this method has been previously reported by Zelenhasićand Salvai(1987), Hisdal et al. (2003), Tallaksen et al. (2004), and Baiziidi and Bayzidi and Saghafian (2011).

Watershed resident's benefit increasing and soil erosion controlling are the main goals of agricultural sustainable development. Landuse optimization operation is an effective solution in order to maximizing the benefits, as well as minimizing soil erosion damages. Current research was executed in southern parts of Ahar Chaey Watershed in East Azerbijan Province. In this research, Multi-Objective linear programing based on Simplex method applied for optimizing landuse in form of three different scenarios: current landuses, implementation of management into current landuses and performance of the standard landuses by WinQSB software. Results showed that due to landuse optimization, irrigated orchards could be increased from 1.736 to 6.629 km2. In addition, based on suitable soils and adequate rainfalls, it is expected to expand the rainfed orchards up to 9.72 km2. Also, the weighted average of soil erosion reduces from 16.3 to 14.8 and from 10.11 to 9.5 tha-1year-1 within the first and second scenarios, respectively. This rate will reduce up to 6.8 tha-1year-1 in the standard landuse situation. On the other hand, maximizing results showed that the annual net benefits increases from 6.6 to 12.9 and 8.46 to 14 million Rails per unit area in the first and second scenarios, respectively. This rate increases up to 15.67 million Rails per unit area for standard landuse situation. Sensitivity analysis results showed that irrigated and rainfed orchards are of high sensitivity in maximizing profit functions due to high profitability of these landuse in unit area. In contrast, the rangelands, forests and drylands are of high sensitivity in soil erosion minimization, so that with increasing their area, soil erosion will increase significantly in the region.

Recent permanent trend of drought in East Azerbaijan province has caused of the reduction in available water resources, a sharp drop in the groundwater levels as well as drying of the Uremia Lake and finally increasing saline lands surrounding the Uremia Lake. Also, the orchards and farm lands which strangely depended on groundwater have been indiscriminately and disproportionately developed related to the current climate condition. Therefore, modification of agricultural patterns is one of the best solutions for efficient water usage in agricultural section and is considered as the only practical solution in relation to the sustainable development of agriculture. This study was conducted in order to determine the suitability for growing pistachio orchards and replacing it with second and third degree gardens. Suitable areas were determined by modeling in ArcGIS software environment based on various affective climate parameters in various stages of plant growth. Results show that about 25 percent of the area is completely suitable, 27 percent is relatively suitable and 48.4 percent is completely unsuitable for pistachio plantation. Also, marginal areas of the Uremia Lake are quite suitable which indicates that the replacement of the orchards of pistachio with high water consuming orchards in the area is fully practical and feasible. The results could provide more reliable basis for agricultural decision-making and planning to provide a new model to replace the pistachio orchards with high water consuming orchards.

The lack of stored soil moisture due to the reduction of rainfall and runoff during the growing season as well as inappropriate distribution of the rainfall is an important factor that affects the establishment of vegetation cover and the production of the plant species. Hence, the application of integrated methods for appropriate and optimized use of the rainfall provides suitable conditions to irrigate dry land gardens, especially in cases that the intensity and duration of the rainfall do not result in suitable runoff. The study of vegetation cover and dry land gardens on the southern hillside of Oun Ibn Ali showed that precipitation and runoff reduction are one of the most important factors preventing the development of vegetation cover and dry land gardens in this area. This present research intends to provide the soil moisture content required for plant species in the region via the optimized use of the rainfall and help to solve this problem. The statistical model of this research was performed in the form of a randomized complete block design with three treatments. The first treatment was control treatment (A), which consists of Buckthorn tree (as a resistant plant to water scarcity and widely used plant in the mountain environment), soil of pit and rotting animal manure as 25 % of the whole of the pit. The second treatment (B) was water harvesting insulator surface system accompanied with Pomis at the end of the system. Water harvesting insulator surface system without Pomis at the end of the system was the third treatment (C). The research was performed during 2011 to 2013 in the Oun Ibn Ali region. The dimensions and shape of plots were 7*5 and rectangle. In addition, they were built on a slope about 30%. The height, collar diameter and leaf area were taken from Oleaster (Elaeagnus angustifolia) trees and soil moisture was also measured using a TDR device during the execution of the research. In this plan, the selection of Elaeagnus angustifolia is due to its resistance to drought and the use of this plant in the vegetation cover of Oun Ibn Ali. The results showed that the use of these treatments increased the soil moisture statistically significant at1 and 5 percent compared to control treatment. Comparing the volumetric moisture average of two years revealed that water harvesting insulator surface system accompanied with Pomis at the end of the system treatment had the maximum stored moisture and without Pomis and control treatments were observed in the next order, respectively. Moreover, the effect of the experimental treatments on the height, collar diameter and leaf area of Oleaster tree was significant at 1 %. The maximum and minimum of the height, collar diameter and leaf area belong to water harvesting insulator surface system accompanied with Pomis at the end of the system and control treatments, respectively. Therefore, it is recommended to use this method to create green space in this area.

Prolonged recently drought in East Azerbaijan has significantly been caused decreasing in the availability of water resources. In this time, gardens and cultivated area has been increased as irregular and disproportional in terms of climate changes. As the highest rates of water consumption is in the agricultural sector, it is very necessary to investigate the alternative water resources and optimization methods in water supplying. Now day, Rainwater harvesting is one of the practical methods for this mean. This research project was executed(in1ha area) with aims to apply gravel filter to infiltration and optimization as well as increasing water storage of harvested rain water into plant root zone. Research area is located in southwest hillslope of Own-Eb neh-ali mountain of Tabriz city. This project was conducted as Randomized Complete Block with five treatments and four replications. Treatments included: Treatment A, control, Treatment B, system by collecting pebbles and vegetation cover from its area with gravel filter, Treatment C, system by collecting pebbles and vegetation cover from its area without gravel filter, Treatment D, system by semi isolated area and gravel filter and Treatment E, system by semi isolated area without gravel filter. The effect of Gravel filter on the penetration rate and increasing of the water storage, were measured in 20 and 50 cm depths using TDR after each rainfall event. Results showed that mean soil moisture rate in Treatment D (semi-isolated system with gravel filter) in both 20 and 50 cm depth are 35.94 and 30.56 respectively and significantly higher than other treatments (p