Journal of Spatial Analysis

---

Has been stated in various sources, soil as one of the most important natural resources has a major, role on the lives of humans. Today soil erosion and sediment production, a problem that is increasing day to day process and loss of surface soils and sediment accumulation in dam reservoirs, canals and also sedimentation damage to the country's economy. One of the most important types of water erosion, gully erosion or (Galli formation). This type of erosion and loss of soil due to sediment production and enormous damages to land, roads and infrastructures, is of great importance. Soil erosion is one of the most important factors that threatens large areas of Iran annually and decreases or eliminates the quality of agricultural lands and rangelands. Due to highlands of Iran in comparison with the grounds and surrounding plains (mean elevation of 1250 m), it has been affected by water erosion. So it is very important to study erosion and present management strategies to reduce the impacts of erosion in basins of Iran. This study to Target morphometry gullies and the influencing factors on gully erosion in the South West sub-basin of ILAM (Cham Fusel).

The average height in the area between 50 and 1,200 meters above sea level and has a dry climate with an average annual rainfall of less than 200 mm and the average annual temperature is above 25 degrees C°. In terms of the main geological formation of this area is affected by gully erosion, Aghajari formation (red mudstone, siltstone and sandstone) are related to Miocene period. Basin area of over 150 hectares affected by gully erosion and slope threshold for gully erosion in the area between 8-2 percent.

The research method in this study is field research, library study and morphometric measurements of gully erosion forms. Besides, on the basis of the geological maps of 1/100000 and topographic 1/50000 and DEM 30 meter area, the software of GIS Arc has been used to make the required maps. Furthermore, the GPS device and a camera has been used in the field to harvest. o

The average height in the area between 50 and 1,200 meters above sea level and has a dry climate with an average annual rainfall of less than 200 mm and the average annual temperature is above 25 degrees C°. In terms of the main geological formation of this area is affected by gully erosion, Aghajari formation (red mudstone, siltstone and sandstone) are related to Miocene period. Basin area of over 150 hectares affected by gully erosion and slope threshold for gully erosion in the area between 8-2 percent. The study, which the sub basin in Cham Fusel ILAM province is located in the South West, with the aim of gully morphometric including deep, height, length, height from sea level, and geographical location as well as classification criteria in three sample gully gully of area, In order to identify factors affecting the development of erosion in the basin plain Cham Fusel was evaluated. The results of morphometric and field visits to the region, factors such as climate factors, slope and aspect, geology (Land genus), tectonic activity and faulting, land use, overgrazing of amount plants, including the most important factors recognized were identified in the development of gully. Which layers to extract some of these factors in Arc GIS software designed and based on creation of maps and data to analyze each of these parameters and their impact on Gully Erosion percent has been paid.

According to the morphometric data, it can be concluded that the gullies basin largely toothed and clawed, permanent, of medium to large gullies, which have expanded continuously. Below you can see images of morphometric gullies sample was collected during the field visit of the area. The results show that due to the widespread earthquake in the region in recent years and local faults can be concluded that one of the most important factors in the formation Galli uplift of the Earth. Also overgrazing by cattle ranchers in the area Tuesday Abdanan city, murmury, Dahlgren and ranchers outside the province of Hamedan and Kermanshah provinces such as ranchers who use the area as Qishlaq, Another important factor in the development of gully erosion in the area. Due to the climatic amount of rain, sleet, snow, ice, temperature and wind could be the climate factors in the rise of water erosion. As the geographical location and local climatological data suggests, the study area is influenced by the Mediterranean winter rains. Which are imported to Country from the West in winter and more in the form of showers and hail rainstorm conditions caused severe erosion gully in the area easy. Lack of growth of vegetation and bare Land is another contributing factor is considered in the development of gully erosion in the area. That is because the rainy season in the winter when the soil in the area is almost devoid of vegetation. The lack of vegetation in the region as one of the important factors, has caused soil erosion, its most destructive erosion of the gullies show. Also according to the map slope and aspect of the region was marked within the range of between 2-8 percent and Create Gully and aspect with the West and the South West and North West are more affected by gully erosion. Finally, all factors except factors of land use, all other factors were named among the important factors affecting the development of gully erosion are the region. It was also found other factors in the evolution of land use contrary to gully erosion region. That is because of the lack of culture in this area and dropping land by farmers as wasteland, and also because of the recent drought in the region has caused more and more extend of gully erosion and land with the ability crop damage and destroy all. The results obtained showed that the study area in terms of of gully erosion in critical condition.

---

 Quantitative assessment of landslide sedimentation in the ILAM dam Basin

Information on the accurate volume of landslides and sedimentation in landslides is a research necessity, with the assumption that the bulk of sediment accumulated in the ILAM Dam (located between , E and , N) is related to the surface landslides of the basin. Although the role of landslides in erosion, sediment transport and sedimentation of slippery basins is confirmed and different experts understand and determine the relationship between the fluctuation of slopes and the fluctuation system in many respects more important than other areas. Because according to the results they can assess the widespread environmental changes, but comprehensive research on the scale of catchment basins has done very little (Harvey 2002). So far, the study of wet landscapes in Iran has been more sensitive to the factors, their sensitivity and their hazards, and there has been no study on the sedimentation of landslides.

Data and Method
First, using a geomorphologic system methodology with topographic maps of 1: 50000, geological map of 1: 100000, aerial photography1: 20000, Landsat TM1988 ETM2002,2013 satellite imagery, and Google Earth in the GIS environment in the following sub-basins and landslide events at the following levels The basin was drawn. The discharge data of the water and sediment flow of three hydrometric stations GOLGOL,CHAVIZ and MALEKSHAHI Station were provided from the waters of the ILAM province. Two models of estimated MPSIAC and EPM models have been used to estimate soil erosion and subsoil sedimentation. The Moran spatial correlation model was used to introduce the spatial pattern of landslides, and the fuzzy logic model was used to determine the relationship between the dependent landslide to the independent variables and the potential risk of landslide hazard in the basin. In order to elucidate the quantitative results of landslide sedimentation, empirical models of estimation of sediment erosion, hydrological model of discharge curve and sediment, observational statistics of sediment during statistical period, landfall time occurrence in compliance with the hydrometric station sediment peak during the statistical period of computation Estimated a small amount of sedimentation of the landslides of the ILAM dam basin.
Result and Discussion
The spatial correlation model of Moran showed that the data have spatial correlation and cluster pattern. The average total sediment production in the MPSIAC model in the GOLGOL basin was estimated to be 13.3 tons per hectare per year under the CHAVIZ basin of 10.3 tons per hectare for one year and 4.00 tons per hectare in the sub-basin MALEKSHAHI. Using hydrological model of discharge-sediment curve, the mean sediment was calculated during the statistical period at the hydrometric station of the sub-basin of GOLGOL 18.8 ton per hectare, the station CHAVIZ 10.4 tons and the station MALEKSHAHI 0.9 tons of sediment per hectare per year was calculated. According to the results of the research methodology, the observation of the sediment in the two stations of GOLGOL and CHAVIZ compared to estimated sediment is related to the events occurring in these two sub-basins.
The data of 16 active landslides were recorded. Under the GOLGOL basin, 9 landslide events occurred, and in the CHAVIZ basin, 7 landslide events, the time of landfall occurrence recorded with sedimentary peaks, the length of the statistical period, the precipitate in the sub-basins was almost synchronized. Average relationship between suspended period of the statistical period - average of the peak delayed flight time of the statistical period coinciding with the occurrence of landslide = the amount of suspended load of landfall occurrence in the basin.
The amount of suspended land slip under the GOLGOL 75088.19 - 315.85=74772.34
Landing slope under the Chavez Basin 19907.30 - 20.24=19887
The area of the sub-basin is about 29,000 hectares and the active landslide area is about 100 hectares. According to the calculations, 77772.34 tons of suspended sediment is a sedimentary passage passing at the GOLGOL hydrometric station, which shows with a coefficient of 1.4 times the suspended sediment load of approximately 104681 tons of landslide sedimentation in this sub-basin, which shows the amount of sediment yield 100 hectares of landslide, so each landslide hectare had an average of 1046. 81 tons of sediment deposited at the GOLGOL hydrometric station. The area under the Chavez Basin is about 14000 hectares and the active landslide area of this sub-basin is about 65 hectares. According to the data of the discharge data, the sedimentation of the Chavez hydrometric station is 19887 tons of suspended sediment load, which shows a 1.4 equivalent of 27842 tons of landslide sedimentation in this sub-basin, equivalent to a slope of 428.33 tons per hectare.
Conclusion
According to the calculations, it is concluded that in the sub-basin of flowering GOLGOL, 37.35% is equivalent to 4.9 tons per hectare per year, the increase of sediment is related to landslide events. As a result, 28.2 tons of sediment per hectare were introduced in one year Dam reaches ILAM. The results showed that in the CHAVIZ sub basin, 38.2 percent is equivalent to 4.6 tons per hectare per year for the increase of sediment related to landslide events. As a result, an amount of 14.5 tons of sediment per hectare has entered ILAM dam in one year. In the sub-basin MALEKSHAHI, there was no increase in sediment during the period without active landslide occurrence. A total of 1237314 tons of landslide deposition have entered the ILAM Dam. To control this threat, the appropriate action by the executive office for sustainable development should be applied.

---

Deforestation or vegetation degradation is one of the main drivers of global earth changes, which has significant consequences in terms of ecosystem performance and biodiversity conservation. One of the ways for studying vegetation changes as the most important indicator of land degradation is remote sensing. In this study, in order to monitor the vegetation degradation trend in Ilam Province.After obtaining and preparing the required data (410 downloaded images) in the ArcGIS and Surfer software, the multiplication, mosaic and georeferencing operations are made. Converting format of images into ASCII is the next stage of the study. By converting this format, the total number of 953552 pixels is studied within the range; after removing the lost and negative values, 328042 pixels are analyzed. Besides, using parametric statistical method of the classical linear regression and programming in R software, the trend of slope variations and significance of slope variations of vegetations are obtained for the 17-year period (2000-2016). Results of this study show that the focus of the highest trend of declining slope variations (trend of negative slop variations) is in the NDVI index across the western half of the studied area and the focus of the highest trend of increasing slope variations (trend of positive slop variations) is in the NDVI index in the center and east. Significance of the trend of slope variations also approves this claim. Thus, the focus of the highest trend of slope variations (negative) in the west and southwest of the studied area along with the highest trend of slope variations (positive) in the center and east is significant at the probable level of 0.05
 

---

Developing urbanization and changing hydrological conditions of natural streams increases the flooding risk. This study tries to do flood hazard zoning in the Ilam city and determine the critical area of the urban regions against flooding by using AHP method and GIS environment. For this purpose, the parameters of the curve number, height, distance from the river, geology, land use, population, slope, soil, building density, worn texture buildings and accumulated flow as effective parameters in flooding hazard in Ilam city selected and of these parameters weighted by using Expert Choice software. The result of the Expert Choice software is transferred to the environment of GIS software and flood hazard map of study area prepared. Results of the study and flood hazard map show that areas with very low-risk, low risk, intermediate-risk, high-risk and very high-risk form the 0.8%, 8.5%, 49.6%, 32.54% and 8.56% of the of Ilam city area, respectively. Also, the central area of the city has the highest risk and the probability of occurrence of the flood due to the high density of population and residential areas in this area and its proximity to the seasonal rivers and old part of the city. Therefore, by examining the results of Expert Choice software, it is possible to identify the most effective factors in the occurrence of flood risk and prioritize them to address management solutions to eliminate or mitigate the effects of these factors.

---

.
 Introduction
Old trees are important and key elements of forest sites and are of great value in terms of forest management, reforestation, silviculture and ecology. Although old trees constitute a small percentage of forest trees, they account for a large share of forest carbon reserve and play a vital role in carbon storage. Understanding the how geographical and site distribution of these trees across the forest is essential to obtain information for forest restoration management. Therefore, this study was carried out to investigate the geographical and site characteristics of old trees of Wing nut, Ash, Hackberry, Sycamore, Elm, Olive, Cypress and Fig in Ilam province.
 
Materials and methods
After querying the villagers and local people and conducting numerous forest surveys, the old trees were identified and selected on the basis of the diameter of the breast. Then their geographical characteristics including city, district, village, geographical coordinates and site conditions including slope, aspect, altitude, soil depth, climate and proximity to water source were measured or recorded.

Results and discussion
The results showed that in terms of geographically distribution, the identified old trees have located in Ilam, Mehran, Malekshahi, Badreh and Dehloran cities. Topographically, the old trees of Wing nut, Elm, Ash and Fig were located in the 0-10% slope class, Hackberry and sycamore in the 0-10% and 10-30% slope classes, olive in the 10-30% slope class and Cypress in the 40-70% slope class. The old trees of Wing nut, elm, Ash and Hackberry were located in the north aspect, fig, sycamore and Cypress in the south aspect and olive in the west and south aspects. The old trees of Wing nut, elm, Ash, Hackberry, Sycamore and Cypress were dispersed at altitude class of 1100–1250 m and olive and fig old trees were at altitude class of 1250–1400 m above sea level. Climatically, the old trees of Wing nut, elm, Ash and Hackberry were located in the very cold Mediterranean climate, Cypress trees and some sycamore trees in the cold Mediterranean climate, and fig, olive and some plantain trees were in the semiarid cold climate. In terms of access to water resources, old trees of Wing nut, elm, Ash, Sycamore, Hackberry and Fig were located on the bed or margin of river, old Cypress trees had no access to water resources and some olive trees were close to water resources. In terms of soil subsidence, old trees of Wing nut, elm, Hackberry, olive, and fig were mostly in soils with medium depths. Old ash and sycamore trees were present in shallow to medium depths and old cypress trees were present in shallower soils. Although the identified old trees were present in limited sites, their long-term and sustained presence in these sites indicates that sites conditions are favorable for their survival.

Conclusion
Therefore, it can be concluded that the presence of low slopes, suitable soil bed and access to water resources were desirable characteristics for stability and survival of the studied old trees in these sites. Due to the above-mentioned characteristics, ecologically similar sites can be found in the forests of the province and can be restored by seed of old and resistant trees.
Keywords: Site, Topography, Climate, Geography, Old trees, Ilam

 

---

Introducing the appropriate model of oak forest and drought relations

in Ilam province

Introduction
The forest ecosystems of the Zagros vegetation region have a very long history of exploitation in various
forms. The material of the Zagros vegetation region is Iranian oak. In recent years, a significant
proportion of oak forests have dried up or have experienced drought. Although the main cause of drought
in these forests has not been determined yet, in the preliminary studies, factors such as climate change,
increasing dust, increasing drought periods, pest infestation and disease, high user changes have been
cited as reasons for drought in the Zagros forests. (Hosseinzadeh and Pourhashemi, 1396). Iran's location
in the arid and semi-arid zone of the world (sub-tropical region) has often been associated with
fluctuations in climatic and atmospheric elements and under the influence of atmospheric currents,
synoptic patterns, irregularities in precipitation and temperature patterns (Rahmati, 2016; 1383).
Comparison of the effect of climatic variables on healthy masses and affected by the decline of oak in
Khorramabad city based on rainfall and temperature data using Pearson correlation coefficient, on annual
growth rings of oak trees Effectiveness of drought of oak trees from both series In general, healthy trees
have been more affected by monthly and seasonal temperatures and have shown the highest coefficient of
correlation with the temperature of the region (Naseri Karimvand et al., 2016). , And the Standardized
Rainfall Index (SPI) and the correlation between these two indicators in assessing and monitoring drought
in different areas of Isfahan province, the results showed that NDVI plant index can be a good alternative
to climatic indicators in drought assessment and monitoring) with the conclusion and colleagues, 2011:
79).
Data and Method
So from SPI and NDVI indicators and Moran index and statistical regression statistics and satellite
images of Modis and Landsat have used to analyze the relationship between dieback of Ilam forests and
happened drought in the region. The precipitation data of 93 rain gauge stations were analyzed during the
statistical period and according to the dry coefficients of SPI index, drought zoning layers of Ilam
province were prepared for two time series of 2000 to 2009 and 2010 to 2019. Greenery's raster layers
were prepared from Modis satellite imagery for the mentioned time series. The results of analysis of
Moran's statistical showed a significant correlation between the SPI index and the NDVI index in spatial
dimensions. By a simple random method, 143 points of oak dieback with dimension of 30 m 2 , which each
point was equivalent to a pixel-size, were recorded with a GPS device, and by simulating in satellite
imagery, the droplet layer of oak dieback was extracted.
Result and Discussion
What is debatable about the results of the implementation of methods for obtaining drought ranges and its
relationship with oak drought points or masses is that the results of the models show a statistically close
and direct relationship between drought and oak drought. . The general trend of oak drought and drought
in these two decades has been from the southeast to the northwest of Ilam province, with increasing
temperature and decrease in rainfall in the southern and eastern regions of the province and increasing
rainfall and decrease in temperature in the central and northwestern regions of Ilam province. The data of
the synoptic stations are consistent. Analysis of satellite imagery and declining greenery in the models
although the study was aimed at meteorological drought and precipitation fluctuations, but spatial
changes of arid points and masses in the province were adapted to field visits and human intervention,
especially in the southeast with agriculture. Under the rubble and the remnants of the dried trees, the ax
has been placed on the roots of these trees, and this trend is spreading in almost other parts of the arid
areas of the province. Therefore, due to the irregularity in the pattern of precipitation and temperature of
the research country (Rahmati, 2016; Zandi Army, 2004) and the effect of monthly and seasonal

نشریه تحلیل فضایی مخاطرات محیطی، سال نهم، شماره 1، بهار 1401 2
temperature on the growth and decline of oak trees in the study (Naseri Karimvand et al., 2016) and other
related research and The flooding situation in the basins of Ilam province, the rainfall, the impermeability
of the soil and their erosion, and finally the lack of moisture in the months before the oak trees grow in
the soil and the increase in temperature in the dry season, which leads to reduced humidity and eventually
greenery. Variables affecting oak drying in linear regression are not responsive, but more accurate results
will be obtained in multivariate regression, although regression analyzes are spatially empty, and X and Y
represent a one-way, quantitative analysis based on the number of dried trees with pixel counts. Drought
range is measured which this defect in SPI method despite its spatial and statistical analysis using Moran
statistical index due to non-compliance in the coefficients of this index with the range of changes in
Moran statistical analysis in statistical analysis is a more appropriate explanatory coefficient than The
regression models showed but at a lower level than the NDVI method it placed. The advantage of NDVI
method with Moran statistical analysis is the relationship between pixel and pixel, ie in spatial analysis,
all pixels that have green changes have been analyzed in the same domain of spatial changes with oak
trees. High results and higher statistical explanation coefficient were obtained than other models.
Conclusion
Although linear regression between extracted oak dieback points with SPI and Moran statistical indicators
was significant, but the relationship between NDVI index and Moran statistic has the effect of
independent variable of drought trend in spatial and temporal dimensions on the dependent variable
process of oak drought with spatial analysis. And nonlinear regression has a more appropriate and
accurate statistical significance and explanation. So this method as desirable method has been introduced
for analyzing of drought and oak dieback.
Keywords: Ilam province, oak forest drought, drought, Moran index

---

Gully erosion is one of the advanced forms of soil erosion, which needs to be analyzed and identified in order to protect the soil. In this research, according to the complex system of factors influencing the creation of ditch erosion, 23 factors were analyzed in the two famous Dempster-Schiffer models and the entropy model, and using Google Earth images and field visits, 331 ditch points were identified, recorded, and a ditch distribution map was prepared. Spatial data of gully erosion distribution were divided into two random training (70%) and experimental (30%) groups. In this research, two indicators of tolerance coefficient and variance inflation factor were used to check the collinearity test, and as a result, two indicators of waterway density and relative humidity index were removed and 21 factors were used in the modeling process. The output results of the layers, weighting and classification and integration in two Dempster-Schiffer and entropy models are the extraction of the zoning map of the gully's erodibility sensitivity. and 30% of the calibration and validation of the models, the area under the ROC system performance characteristic curve and the area under the AUC diagram of the Dempster-Schiffer model with an explanation factor of 0.934 and the maximum entropy model with an explanation factor of 0.936, both models have an acceptable percentage of the area under the curve were that this issue shows the high performance of both models in the region. Among other results of statistical analysis, the prioritization of the impact of 21 factors in causing ditch erosion in the region was determined. The scientific results of the research can be promoted and taught, and from the practical point of view, the relevant executive body to control ditch erosion can take the necessary measures using the results of this research.