فصلنامه پژوهش های ژئومورفولوژی کمی
سال هشتم شماره 1 (پیاپی 29، تابستان 1398)

Underground water is one of the most important global natural resources for communities as well as for agricultural and industrial activities. Considering the unique weather conditions in Iran, lack of rainfall, distribution and timing of inadequate space, shortage of water resources in terms of quality and quantity and the widespread expansion of Karst Formation throughout the country, the study of water resources of karst is very important. The main characteristics of feeding area in the karstic aquifers impact on the type of feeding, flow and the vulnerability of aquifers. Therefore, identification of the feeding areas in karstic aquifers plays a key role in assessment of the hydrodynamic and hydrochemical properties of aquifers, as well as their scientific and optimal management and utilization. Karstic water resources are the most important water supplies in the Kermanshah province. The Paraw-Bisetun aquifer is located within the carbonate rocks domain of the high Zagros mountains. Due to its high altitude and precipitation, it is susceptible to karst conditions, as witnessed by its particular geomorphology and the developed karst features such as sinkholes, uvalas, and caves. The karstic aquifer of Paraw-Bisetun has 15 important springs, supplying water for drinking and agriculture in the neighboring area. In this research, the KARSTLOP model was used to model the locations of feeding areas of the Paraw-Bisetun aquifier. The KARSTLOP model abbreviated to the first letters of the eight layers K: karstification process (Karst development), A: Weather conditions, R: Runoff, S: Slope, T: Tectonic, L: Geology, O: Covering layer and P: Vegetation. This model can be implemented in the Arc GIS software environment and the individual judgment of the researcher is involved in the selection of analytical parameters and definition of evaluation criteria and algorithm. To determine the source of water karsticity of the springs fed from the Paraw-Bisetun aquifer, to perform hydrological and hydrochemical calculations related to them, and to prepare their Piper diagrams, the RockWare Aq.Qa and the PHREEQC software were used to calculate the saturation ratio. The Fill Sink method was used to extract the sinkhole from a 10 meter DEM of the study area and field observations were performed to identify the surface karst geomorphology, to record the position of landforms and to verify the accuracy of the results of the Fill Sink method. Investigation of the hydrochemical characteristics of the outlet springs of the Paraw-Bisetun aquifer indicates the calcific facies and the carbonate type of these springs, relevant to their karstic origin. Also, the molar ratio determination indicates that low durability of water in the aquifer and reflects the presence of limestone in the feeding area of the springs. According to the final map of the KARSTLOP model, the percentage of feeding aquifer is between 44.2 and 88.7%, of which 70-80% of the area with the highest area is 65% of the area, corresponding to its central regions spanning from north-west to south-east. Subsequently, the 60-70 percent floor, with an area of 24 percent, is in line with the margins of the aquifer. The floor with a nutritional level of more than 80 percent is located in the upper high point and parts of the plain areas of the other mountains of the region, which is in accordance with the rock outcrop of the Bisetun carbonate formation with the highest density of faults. In general, it can be deduced that the aquifer of Paraw-Bisetun has a high nutritional potential due to carbonate rocks outcrop and without soil cover in most of the regions, so that about 96% of the aquifer area has a 60% higher nutritional potential. The high percentage of nutrition of the aquifer and the distribution of spatial areas of the feeding area, which indicates the dominance of centralized penetration by the sinkholes and vertical caves in the highlands, suggests that for the discussion of contamination and its prevention, it requires conservation and management plans. The high percentage of nutrition of the aquifer and the distribution of spatial areas of the feeding area, which indicates the dominance of centralized penetration by the sinkholes and vertical caves in the highlands, suggests that for the discussion of contamination and its prevention, it requires conservation and management plans. The high percentage of nutrition of the aquifer and the distribution of spatial areas of the feeding area, which indicates the dominance of centralized penetration by the sinkholes and vertical caves in the highlands, suggests that for the discussion of contamination and its prevention, it requires conservation and management plans. The high percentage of nutrition of the aquifer and the distribution of spatial areas of the feeding area, which indicates the dominance of centralized penetration by the sinkholes and vertical caves in the highlands, suggests that for the discussion of contamination and its prevention, it requires conservation and management plans.

Landslide susceptibility mapping (LSM) is a proper method to predicting landslide hazard risk in order to reducing its consequences. We prepared the LSM mapping of Kan basin by used of index of Entropy model and SVM-S. The validation of the produced maps is evaluated by used of the area under the curve ROC. Study area The study area is located on the north west of Tehran Province, between 35°46′ and 35°58′ N latitude and 51°10′ to 51°23′ E longitude. The area of the basin is 204.385 km2 (Figure. 2 (a)). Material and methods In the present study we used of 8 parameters consisting the distance from river, distance from fault, distance from road, land use, rainfall, aspect, slope, elevation, and lithology. Distance from the road was computed from the road at the interval 200 m using ArcGIS software (Fig.2f). The distance from road and distance from fault was calculated in same way (Figure. 2 d, h). The land use map has reclassified to 5 class (Figure.2g). The lithology parameter has been obtained by the reclassification of the geological map of Tehran at the scale of 1:100000 (Figure.2h). The digital elevation model (DEM) was extracted from the 1:50000 scale topographic maps. The parameters of slope degree (whit 5 classes), aspect layer were produced by used of the digital elevation model. (Figure.2c & d). The topography layer was reclassified into 5 class (Fig.2a). the introduced layers were used in this study according to the type of the models to produce of LSM maps. Entropy index Entropy is a measurement of the instability, imbalance, and uncertainty of a system (yang et al, 2010). The equations used to calculate the information coefficient dj representing the weight value for the parameter as a whole, are given as follows: pij=xij/(∑_(i=1)^m▒xij ) (1) (2) Ej= -k+∑_(i=1)^m▒pij ln⁡(pij) K=〖c(lnm)〗^(-1) (3) ((4 dj = EJ – 1 After calculating the total weight (wj) using Equation 5, the landslide risk of the case study is evaluated: (5) Hi = ∑_(i=1)^m▒xij In equation 5, H is the coefficient of landslide risk, wj is the final weigh of all the factors and Xij is the weight of each factor (zongji et al, 2010). The final landslide susceptibility map was prepared by the summation of weighted products of the secondarily parametric maps. H = (S×0/54) + (Df × 0/74) + (E × 0/082) + (Dr × 0/51) + (Dri × 0/51) + (A × 0/066) + (lu × 0/16) + (Lt × 0/0064) Support Vector Machine SVM algorithm as one of the most popular methods for solving regression problems has had significant results in landslide sensitivity zoning. consider a set of linear separable training vectors Xi (i = 1, 2, . . ., n). The training vectors consist of two classes, which are denoted as Yi = ±1. The goal of SVM is to search an n-dimensional hyperplane differentiating the two classes by their maximum gap. Mathematically, it can be expressed as: 1⁄2=∥w^2∥ (6) Y_i=((W.X_i )+)≥1 (7) A Lagrangian formulation is introduced to solve the problem (equation 8). Thus, the goal is now to minimize the Lagrangian L with to W and b and maximize with respect to λi. For this reason, we used of following equation: L=1/2∥w^2∥-∑_(i=1)^n▒Y^i ((W.X_i )+b)-1) (8) four types of SVM is existed: linear, polynomial, radial basis function (RBF) and sigmoid. The mathematical representation of each kernel (linear, polynomial, radial basis function, and sigmoid) is listed below, respectively: K (X_j 〖.X〗_i )= X_j^i.X_j K (X_j 〖.X〗_i )=(γ∙ X_j^i+r) ._ γ>0 K (X_j 〖.X〗_i )=e^(-γ〖(X_i-X_j)〗^2 ) ._ γ>0 tanh (γ .X_i^T. X_j+r) γ, d, and r are user-controlled parameters, as their correct definition significantly increases the accuracy of the SVM solution. In the present research we used of Sigmoid function. To measure the validation of the models, we used of a relative ROC by comparing the existing landslide location with the two landslide susceptibility maps. The success rate curves were obtained by used of the 70% training dataset (29 landslide locations). ROC curve (AUC) represents the quality of the probabilistic model (it is ability to predict the occurrence or nonoccurrence of an event). Result and discussion The area of the low, moderate, and high classes based on the SVM model were found to be 109.485 km2, 38.7 km2, and 56.2 km2, respectively, whereas based on landslide susceptibility map by used of index of entropy, the 118.175 of the study area has low susceptibility risk, and the moderate, and high susceptibility zones have the 41.2 km2, 45.02 km2 of the study area, respectively (Fig. 3). Based on the entropy model, the 8 numbers of the landslides points located on the high-risk zone and the 8 numbers of the landslide points located on the moderate risk zone and low risk zone have 10 of the landslide points. Based on the LSM map produced by the SVM-S model, the 13 numbers of the landslide points located on the high risk zone and the 5 number of the landslide point located on the moderate-risk zone. The ROC plot assessment reveals that the AUC in the susceptibility map based on the index of entropy model was 0.86 and the AUC in the susceptibility map based on the Logistic Regression model was 0.91 (Fig. 5). Conclusion The high-risk zone on the LSM map produced by the SVM model is located in the north east and the west and south of the basin and based on the LSM map produced by the Entropy model is located in the north east and the south of the basin. The LSM map has produced in a regional scale, so further study need be carried out at the site-specific level to determine the exact extent site of the slope instability. Keywords: LSM, Index of Entropy model, Kan basin, Support Vector Machine, Sigmoid function, SVM-SIGMOID.

Introduction One of the proper methods for utilizing water resources is the use of subsurface water in coarse gravel sediments by the construction of underground dams. Underground dam is a structure that blocks flow of subsurface water in alluvial deposits, which makes it easy to store water in porous layers and make the use of this water easier. Locating the site of underground dams requires extensive studies. The importance of locating due to impact of many criteria and factors such as hydrological and meteorological criteria, general geology, stratigraphic layering and permeability, topography, groundwater status and socio-economic factors. It takes a lot of time and money to study these factors using traditional methods. River grade, shape of the valley, volume of river runoff received, subsurface water quantity, depth of the bedrock to the surface of the earth, the high volume of alluvium in determining the location of the underground base are very important Locations with faults, spring and wells are inappropriate areas. The purpose of this research is identify suitable sites for construction of underground dams in the North-west of Kermanshah province as a region with appropriate rainfall in the wet season and abundant water requirement in dry season. In this research, with using decision support system (DSS) method, according to the factors and conditions in the study area, removal of step by step inappropriate places has been done. At the end, among appropriate sites, the suitable sites was determined. In this study, the thickness of sedimentary layers has been calculated using satellite images and trigonometric relationships with the least cost and credible accuracy. Using the drawing of altitudes and surface elevations buffer for streams, valleys with a standard width were identified. Methodology In the first stage of  this study, the digital elevation model (DEM) layer with 12M pixel size for study area was provide. Using it, areas with 0 to 5% slope classes were extracted and other areas were deleted. Using DEM and ARC INFO software, a layer of streams was created and overlaid with the slope layer. Streams with 0 to 5% bed slope were extracted. This new layer overlaid with the fault layer and the streams with fault in their beds were removed. It is economical to construct an underground dam in the local area of ​​the streams, with a maximum width and depth of the valley not more than 15 meters. For greater accuracy, a 60-meters surface buffer (30meters from the center of the stream to each side) was determined in the ARC GIS environment for the streams. Using ARC INFO software, a 15-meters altitude buffer was also drawn up for the streams. These two buffer were overlaid and the interleaved sites were extracted as appropriate valleys. For match the extraction points with the actual situation in the area, using the information layer of the villages, topographic maps and Google Earth software extraction points from the overlying surface of the elevation and surface were examined. In order to identify the sediment surface of the streams and their bedrock, extraction point sites overlaid with a 1: 1000000 geological layer of area and locations in inappropriate geological formations were removed for non-permeability. At the last stage, thickness of the sediments that were considered almost suitable at the field visit stage was calculated using the trigonometric relationships. Results and discussion Based on the initial result of overlaying of altitude and surface privacy (as an innovative operation), 257 locations were identified for the construction underground dams. Since this point of the sites is initial, based on local conditions, were examined. Using the topographic map of 1: 50000 area and villages of the region, points were evaluated. The point of a location located more than 500 meters from residential areas and agricultural lands was removed from the final map. Determination of the distance between the underground dams and residential areas depends on the topographic conditions. Forzieri et al. (2008) considered up to 15 km distance in Mali and Haji Azizi et al. (1390) 5 km distance from residential areas. At the end of this stage, number of suitable locations dropped from 257 to 33. With overlaying of the layers of the appropriate points with geological map, from a total of 33 points, 5 points were eliminated due to presence of dissolved formations (limestone) and the lack of coarse grained sediment for water storage. At the field visit, total of 33 locations were evaluated and verified. In terms of proper valley width to construct underground dams, the accuracy of the whole point of the area was confirmed. Regarding the suitability of seasonal flows for the construction underground dams, in the field survey and survey of hydrometric station data in the region, 10 other points due to being in permanent rivers were eliminated from the total points. Alluvium thickness has a direct and large impact on the volume of water stored in the reservoir (Mohebi e Tafreshi et al., 2014). According to studies, the thickness of 5 to 20 meters of sediment is suitable for construction underground in the study area. So, calculating the thickness of the sediment, 6 other site points were removed due to the unsuitable thickness of the sediments. Conclusion Based on the results of this study, 11 suitable spot locations were determined. By constructing underground dams in these areas, it is possible to store sub-surface water from the rain that goes unused, and in dehydrated seasons for drinking, agricultural irrigation, and the creation of new agricultural lands Used. In field studies, at the foot some of slopes, there were a lot of gravel. If, in one place, all conditions for the construction underground dam without adequate thicknesses are provided, it is possible to carry these deposits to the desired location to provide a suitable reservoir for storing water in that location.

Geomorphodiversity Investigation of Damavand volcano and its surroundings based on the GmI Index Extended Abstract Introduction Geomorphodiversity defined as “the critical and specific assessment of the geomorphological features of a territory, by comparing them in an extrinsic and in intrinsic way, taking into account the scale of investigation, the purpose of the research and the level of scientific quality”. Area with high geomorphodiversity can provide a wide range of services, including provisioning (food products), regulating (erosion regulation) and supporting (land and water as a platform for human activities). The geomorphodiversity of Damavand volcano and its surroundings require special attention to identification, management and conservation, prevent degradation and construction in vulnerable areas. Promoted recognition of Geomorphodiversity can provide useful information about the management of geomorphological heritage and how to better protect them against destructive human activities and environmental hazards. Also, the preparation of Geomorphodiversity maps, the identification of areas with high potential for regional development and geotourism, provide comprehensive information for designers, planners and responsible organizations in these areas and create solutions for the environmental complexities of the study area; therefore The main objective of the present study is identification and quantitative assessment of geomorphodiversity and the preparation of geomorphodiversity maps of the study area. Methodology Damavand volcano was inscribed as Iran's first natural property in the national heritage list at the Cultural Heritage, Handicrafts and Tourism Organization on July 3, 2008. The data of this research include geology, drainage network, slope, roughness and landforms, which are related to the evolution of the physical landscape of the study area. Landsat satellite imagery, thematic maps, DEM10 meter, etc. are used to achieve the goals. In this study, geomorphodiversity quantitative index has been used for analyzing the results and the results has been validated by comparing with geomorphology map and field control. Geomorphodiversity index (GMI) is presented by Melelli et al. (2017). The value of geomorphodiversity is calculated from the sum of five factors including geological diversity, drainage network density diversity, roughness diversity, slope position index diversity, landform category diversity. All of them are grids of different terrain parameters. Results and discussion In the geomorphodiversity of the study area, the V1 class is 1.05%, V2, 25.84%, V3, 37.17%, V4, 22.14% and V5, 4.84% of the total area of the region. Distribution of geomorphodiversity of the study area is mainly related to the topography condition; the lowest value of the geomorphodiversity (V1) is related to the Damavand volcano cone, especially the lava plain in the northwest. The highest value of geomorphodiversity (V5) extends along the valley of Haraz, where different types of geomorphological processes are meeting, such as fluvial, glacial, hillslope, igneous, tectonic and etc. In other words, increasing the value of geomorphodiversity along the Haraz valley can be interpreted as a result of the frequency of the geomorphological processes presence. Validation results also showed that there is a good correlation between the geomorphodiversity index value, the number of different types of landforms and the average number of landforms. The most extensive area with the lowest value of geomorphodiversity is located in the northwest of Damavand and in the Sardagh plain. Along the Haraz valley, which has the highest geomorphodiversity, there are different types of landforms, including river terraces, lake terraces, landslide, talus, basaltic prisms, deep valleys, travertine landforms, high channel density and so on. In general, in the study area, the high value of geomorphodiversity is mainly located at the meeting place of Damavand volcanic lava with rivers. Conclusion The geomorphodiversity of Damavand volcano and its surroundings are one of the most important fortunes in the region. Geomorphological landforms have been maximized in these areas and offer a unique geomorphological complex. The geomorphological landforms density in these areas are maximized and provide a unique geomorphological collection. Areas with high geomorphodiversity are rich and unique collection of different types of landforms and geomorphological processes that require special attention for geotourism, scientific-educational, national and world heritage, geopark, and so on. In addition, the determination of areas with high geomorphodiversity is another way to promoting geomorphological heritage; In other words, geomorphodiversity are the backbone for identifying and assessing geomorphological heritage for various purposes of geotourism, national and world heritage, geopark, conservation of geomorphological heritage. To ensure that the values of these areas can be preserved for current and future generations, managing and conserving them is very important. Nevertheless, studies of geomorphodiversity are at an early stage and need more help and reflection by geomorphologists in collaboration with other researchers in the Earth sciences. Keywords; Geomorphodiversity, Landform, Geomorphological heritage, Quantitative assessment, Damavand volcano

Extended Abstract Introduction Earthquakes and landslides are the most harmful natural disasters that have synchronization and time and space correlation with each other. This study was carried out in seismic zone affected by Arasbaran earthquake (Qaradagh), in which there were a lot of casualties and destruction. On 11 August 2012, the Varzeghan twin Earthquakes [Mw 6.4 and 6.2] struck Varzeghan, Ahar and Heris region in NW of Iran. It killed more than 306 people and a large number of people were injured and has caused increment in slope instability, and transformation of geomorphic processes into risk factors. Methodology In this research, landslide susceptibility was zonated by using two models of analytic network process (ANP) and logistic regression in Qoshadagh fault system and an appropriate model was introduced. Our analysis of geohazards distribution allowed evaluation of geomorphic, climate- hydrologic, human-land use and morphologic controls on earthquake-induced-land sliding, process mechanisms, and hazard process chains, particularly where they affected local populations. For this study, both the Sentinel-2A MSI and Landsat-8 OLI data (2017) were used. Fourteen effective factors of landslide (Lithology, Land use, Fault, DEM, Climate, Aspect, TWI, SPI, Slope, Road, River, Rain, LS & Soil) were generated in ArcGIS, and then weight of each factor was determined in Super Decisions software, and final zonation maps were regained from ArcGIS afterward. Results and discussion After preparing the standard maps and options, the network analysis method (ANP) and logistic regression were used to investigate the susceptibility of land slides. In order to implement the network analysis process, the related criteria layers were developed in ArcGIS software. Then, in the SuperDisign software environment, the main model of the network analysis process was designed based on completed questionnaires by the experts. Then, a logistic regression model was used to analyze the spatial relationship between the land slide event and the effective factors in this event. The maps related to the factors affecting the land slides of the study area, which are independent variables in the land slide event, were introduced into Idrisi-Selva software and processed for logistic regression modeling. Landslide distribution layers in the area were also converted to the binary map 0 and 1 by the Calculator Image function. This means that the slip pixels on the map are shown with the number 1 (slipping) of non-slip pixels with 0 (no slip). Finally, after entering the data into the logistic regression model, the coefficients of the model are extracted using the effective parameters in the Idrisi-Selva software. Distance from the fault and precipitation factors had the most, and Land Use Factor had the least effect on landslide occurrences. Occurrence of 62.2% and 71.1% of landslides in high and very high risk classes in ANP and logistic regression, respectively, indicated an acceptable accuracy for landslide prediction maps. Validation results of methods with ROC index showed that AUC of the maps in model was 85.52%, and in analytic network model it was 81.35%, with a standard error of 0.06; while both represented a very good predictive capability Conclusion The purpose of the present study was to study the ability of this model in comparison to other methods, in addition to modeling the sensitivity of Earthquake-Induced-Landslides using the network analysis method (ANP). Therefore, effective factors in the field of scaling were identified and in the Super Decision software environment, the weight of each factor was determined. Then, the weights received in ArcGIS software environment were converted to the final map of the zoning of the landslide. Accordingly, among the fourteen effective factors in the occurrence of landslides in the area, the distance from the fault has the most effective factor in the two models (ANP) and (LR) and has the highest coefficient of effect on land-slip occurrence while land use in Both models have the lowest coefficient. As expected, the main faults, in particular the fault and the fault lines, and interspersed with the twin Arasbaran earthquakes, as well as the northern and northeastern slopes, are more susceptible to instability. The map of the model implementation categorized the susceptibility of earth sculpting in 5 classes with very low, moderate, high and very sensitive sensitivity. Very high and high risk classes have been shown to be 6.6% and 13.4% of the area of the region, prone to hazardous landslides. These results show a high correlation and correlation with the model implemented in the logistic regression method, which is 5.8% and 17.7%, respectively. The occurrence of 62.2% and 71.1% of landslides in high and very high risk classes, respectively, in ANP and logistic regression, indicate the acceptable accuracy of predicted maps for landslide. The results indicate that ANP methods and logistic regression are accurate in the study of landslide in the area affected by the Arasbaran earthquake. Logistic Regression model had better results. Using these methods together and comparing them with regard to the dependencies of landslip issues can be very useful for identifying areas prone to landslide. As suggested, these methods have acceptable results in analyzing the sensitivity to landslides. The highest density of landslides, even old landslides, are not accidental in the two-axis seismic centers with magnitudes of 6.4 and 6.2 in 2012, indicating a history of seismicity and high tectonic activity in the area. Keywords: Analytical Hierarchy Process, Logistic Regression, Landslide hazard zonation, Arasbaran earthquake, Qoshadagh Fault Zone.

Introduction due to  the phenomenon of global warming that we are facing today, studies of glacier areas, and sub glacier and its effects are of particular importance, because by examining the  past climate,   future  theclimate change processes are better analyzed. One of the most important climatic heritage of the 4th era in Iran's highlands is glacier geomorphology and its glacier phenomenon. Climate change is one of the most characteristic features of the Quaternary period. Landforms have changed a lot during this period. Due to the warm and dry climatic conditions in central Iran and the presence of glacial landforms which shows the cold and wet climate in the past, we decided to study the current climatic conditions in the region with the aim of studying the past climate conditions. Location of study area   Shirkouh Highlands are located in 40 kilometers south-southwest of Yazd. Shirkouh is the most important infiltration mass of Central Iran. The mid-Jurassic Batulite of Shirkuh consists of three major units:  granodiorite, monzogranite and leukogranite. Shirkouh mountainous has created two different worlds of climate at a short distance in its environs. The location of Shirkouh elevations in the direction of northwest-southeast is like a long wall against the air masses and it is a temperature moderator in the region. Methodology   In the present study, first, by identifying the glacial landforms (circuses, moraines), the boundaries and extent of the last glacial period in the Shirkouh Mountains were determined and then by discovering their local-spatial relationship with climatic (temperature and precipitation) and morphological (elevation and direction and gradient) parameters, the temperature drop in the time of glaciers existence and the amount of ELA have been determined. In the first step, by using Digital Elevation Models (SRTM 30 meters from the American Geological Survey) and 10 meters DEM of the organization of mapping of Iran and the Google Earth image, the position of the circuses has been determined and then, their accuracy has been verified by field surveys. In the next step, the correlations between the quantitative characteristics of the glacial landforms with morphological (elevation, slope, slope direction) and climatic parameters (temperature and precipitation of14 weather stations of the region from 1985 to 2008) have been determined. Then, the amount of drop in temperature required to form a glacier in the last glacial period in the mountainous area has been calculated by comparing the current and the past ELA. Discussion and results Detection of glacial phenomenon of the Shirkouh area: Based on the circus reflected in topographic maps and satellite imagery (Google Earth) and field surveys in the Shirkouh mountains, 62 circuses were identified.  Due to the Wright method by passing 60% of the glacier circus location in the Shirkouh region, the Quaternary snow line is located at an altitude of 2491 m. This line represents the height of the zero-degree companion in the glacial period, by which we can reconstruct the past temperature conditions. Field observations in the Shirkouh region confirm the existence of glacial conditions such as circuses, glacier valleys and moraines. Based on the meteorological data and current climatic conditions, the current zero-degree temperature line is at 5438 m height, above the peak. According to the temperature equation, the past snow line is situated in the range of 22 centigrade degree of today’s temperature.  Regarding the paleo climate of the region: the distribution of circuses is different in Shirkouh altitudes. the highest number of circuses have northeast, north, west and east direction and the smallest of them have southeast and southern slope direction, respectively. The lowest circuses are situated in the range of 1800 to 1900 meters and the highest of them are located in the altitude range of 3800 to 3900 meters. Conclusion With the studies of the prior researches, including the Hagedorn, the glacier processes at this region were identified at the height of 1925 to 2240 meters in the Quaternary   and Moraines have also been seen at a height of 1600 meters. According to the fact that during the past decades we are facing with the global warming, due to the amount of annual snowfall and the lack of glaciers feeding, in contrary to the Quaternary period, today, we are no facing with the glacial processes in the Central parts of Iran. due to  the phenomenon of global warming that we are facing today, studies of glacier areas, and sub glacier and its effects are of particular importance, because by examining the  past climate,   future  theclimate change processes are better analyzed. One of the most important climatic heritage of the 4th era in Iran's highlands is glacier geomorphology and its glacier phenomenon.

Introduction Soil erosion and sedimentation are one of the most important environmental , agricultural issues and food production in the world. Soils that are eroded by runoff rain drops are transported to the rivers and released from the catchment areas. Therefore, regional analysis of the sediment load of rivers, especially in arid and semi-arid regions, and its relation to the characteristics of watersheds is very important in estimating the amount of erosion and sediment. Dez Dam Basin is one of the important catchments of the western part of the country in terms of water, energy and agricultural production. Dez river collecting runoff from vast areas of Lorestan, Isfahan, Khuzestan, Chaharmahal Bakhtiari and central provinces, and after watering it into the plains of its route to the Karoon River. The lack of sufficient number of sedimentation stations and the limitation of the number of statistical years as well as the low number of flood sampling during the outbreak of the Dez river indicate the importance of this study in modeling the suspended sediment estimation in the Dez dam basin. Therefore, the aim of this study is to model the relationship between the amount of suspended sediment load with hydrological characteristics and the environment of the dam dam. Matarials and methods The catchment area of the Dez River in the middle Zagros Mountains and between the geographic coordinates is 17 ° 31 and 15 ° 34, the northern latitudes are 17 ° 48 and 50 ° 50, respectively. The watershed of the Dez River is considered to be part of the Persian Gulf watershed for the general classification of the watersheds of Iran. The dominant climate in the study area is warm, dry and semi-arid soils. The average annual rainfall of the basin is 670.5 mm. The annual temperature varies from 6.7 to 24.4 degrees Celsius. The Dez River is more than 500 km long in the main drainage of the Dez Basin. For this purpose, 22 sedimentation stations were selected at the upstream of the Dez dam. 16 stations were used to develop suspended sediment estimation models. Then the developed models were validated with the information of 6 other sediment stations. Based on the suggestions of past studies, 15 different physiographic, climatic and hydrological features were considered as input parameters in modeling estimation of suspended sediments. Height, area and basin environment were obtained from maps of 1: 25000 scale with the use of ArcGis9.3 software. The average slope of the watercourses and basins was extracted from the DEM map. The average annual precipitation and average rainfall of rainy and floody days of the year included Azar, Di, Bahman, Esfand, Farvardin and May. Principal Component Analysis (PCA) was used to select the most important parameters. The stations were then divided into two homogeneous groups based on the selected parameters and cluster analysis. Finally, the models were developed step by step using the multivariate regression method. R2, error rate (RE) and mean error (ME) were used to evaluate the models. Result and discussion The results showed that the sub-basin of the Takht-Dareh valley with an average of 1256 tons per year is the lowest and the tele-alarm-dose basin has the highest suspended sediment with 8800919 tons per year. Based on the results of the analysis of the main components, four factors of annual discharge with a weight of 94, average annual rainfall with a weight of 88, Horton coefficient with a weight of 90/0 and average height of the basin with a weight of 0.93, the most important parameters Estimates of sediment were suspended. According to the effective factors, sub-basin groups were grouped using cluster analysis. According to the cluster composition at 25, sub-basins were classified into two homogeneous groups. Suspended sediment estimation models were developed based on four factors and using stepwise stepwise multivariable regression method. All developed models in the homogeneous region 1 have more R2 coefficients than those developed in homogeneous region 2. The results of the evaluation of the models showed that the developed model using the four factors with 8% = RE, t ME= 25375 t ya-1, and R2 = 0.93 had the best performance and the least percentage of estimated error. The results also showed that Dubai was the most influential parameter on suspended sediment estimation models annually. Conclusion The purpose of this study was to model the relationship between the environmental characteristics of the Dez dam basin and the amount of sediment suspended in this river basin using stepwise multivariate regression method. The results also showed that the regression model derived from all four factors of annual discharge, average height of the basin, horton coefficient and mean rainfall, with the highest determination coefficient (R2 = 0.93), the lowest error (RE = 8%) and the lowest mean error (ME=125375 t ya-1) has a very good performance in estimating suspended sediments. In general, it can be said that the combined use of principal component analysis, cluster analysis and multivariate regression There is an appropriate and acceptable efficiency step in estimating suspended sediments.

Extended Abstract: In general, satellite data classification with digital interpretation method is separation of similar spectral sets and classification of images into groups or classes placed in each class, spectrum or a value (not being separable statistically). In other words, when a classification is carried out on data sets or images, in fact spectral classification is conducted and in the new image, every class or category is indicator of a specific phenomenon with a unit spectral value. The base for data classification is comparison of spectral value of picture pixels with the samples introduced by its interpreter or with initial classes and   categories, constituted automatically in the time of digital interpretation. Hence, the pixels, values of which are not statistically significant, are placed in the same spectral group or class. So, preparation of land use maps by digital classification of remote data sensing depends on the methods used in classification, land coverage condition and environmental condition of the area. The aim of the present research study is to compare pixel-based classification algorithm (maximum likelihood) with object-based method (support vector machine) to prepare land use maps with the aid of SPOT5 satellite HDR image sensor and Landsat 8 satellite OLI image sensor in Ahar Chai basin from Ahar area to Varzegan. Methods and materials: In this research study a SPOT5 satellite HDR image sensor dated back to August 23, 2006 with four spectral bands and separation power of 10 m with path and passage 138-272, 138-273, 139-272 and 139-273 and Landsat 8 satellite OLI image sensor dated back to 2017/07/21 with route and passage 168-33, isolation power of 30m, increased to 15m through combination with band 8 were used to collected the necessary data. Topography map at scale 1/50000 and geological maps of 1/100000, land control data harvested by GPS, Garmin model and ArcGIS 10.4.1 and ENVI 5.3 software were used in this study. Information extraction from remote sensing data especially land cover can be obtained by digital classification. In practical some people are more comfortable using visual interpretation to retrieve land cover information. However, it is highly influenced by subjectivity and knowledge of interpreter, also takes time in the process. Digital classification can be done in several ways, depend on the defined mapping approach and assumptions on data distribution. The study compared several classifiers method for some data type at the same location. Maximum likelihood Classifiers (MLC) which use pixel-based and object-based classification using advanced supervised machine-learning algorithms such as Support Vector Machine (SVM).Therefore, the aim of the present research is to compare maximum likelihood in pixel-based processing and support vector machine algorithm in object-based processing in order to evaluate their performance in classification of LANDSAT 8 and SPOT 5   satellite images.  The necessary pre-processing including geometrical corrections and atmosphere corrections were conducted on the image. Discussion and Findings of the study indicated that both support vector machine algorithm in object-based classification and maximum likelihood in pixel-based classification enjoy higher and appropriate accuracy in classification and exploration of land use maps for Ahar-Varzegan area. However, considering the results obtained from the two algorithms used in this research, it can be concluded that   support vector machine algorithm in object-based classification of satellite images provide the conditions for higher accuracy compared to maximum likelihood in pixel-based method. One of the main reasons of achieving higher accuracy in support vector machine algorithm classification is that in this algorithm, in addition to spectral data the data related to contexture, shape, position and content are also used in classification process, hence classification accuracy is increased significantly. After doing classification procedure using both algorithms of pixel and object-based method we proceeded on evaluation of the results reliability in them The results showed that, Findings obtained from both classification algorithms indicated higher accuracy of object-based classification with total accuracy of 94.99 % and Kappa coefficient for exploration of land use of Ahar-Varzegan area, this is while pixel-based classification with total accuracy of 89.39% and Kappa coefficient of 0.89% indicated acceptable level of classification. So that Feizizadeh et.al (2009) in exploration of East Azerbaijan land use with the aid of the two pixel-based and object-based algorithms concluded that object-based classification with total accuracy of 95.10 % has higher efficiency in exploration of the province land use compared to pixel-based classification with total accuracy of 88.37%.

Estimation of erosion caused by ephemeral gullies with EGEM model in Lanbaran basin, Ahar Chay Abstract Gully erosion is one of a variety of geomorphic phenomena. Every year, this phenomenon causes lots of damage to agricultural land, rangeland and roads. Identifying the causes and modeling it play an important role in the management of catchment areas. One of the types of gully erosions is the ephemeral gully erosion. Currently, the only model that can assess this type of erosion is the EGEM model. Therefore, the study area has been observed since 2008 after rainfall. The data needed to implement the EGEM model was collected. SPSS and Excel softwares were used for quantitative analysis. According to the analysis of the measured parameters and the results of the model, the highest correlation coefficient exists among the estimated and measured depth, cross-sectional area, volume and width with the corresponding correlation coefficients of 0.849, 0.845, 0.818 and 814/0. Keywords: ephemeral gully – Lanbarn – Ahar Chai - EGEM model Introduction One of the most important types of erosions is gully erosion. This type of erosion causes significant changes in agricultural lands in arid and semi-arid regions of the world, and has harmful economic and social consequences (Esmaeiliy et al., 2015). (Kheir, et al., 2006; Marker, et al., 2008; Desmet, et al., 1999). Gully erosion is among the most important and complex types of erosion for the control of which the effective factors need to be identified (Ahmadi, 2011; 243). Gully growth is often considered as a major measure of environmental change, because its growth is not rapid and is not a common cause of erosion (Bayati Khatybi; 2004; 54). And because of the concentration of water flow in small grooves, a large volume of fine grains can easily be displaced in a short time due to ghost erosion compared to other erosion factors (Ryan, 2010; 288). According to Poesen, the gullies are channels with a cross-sectional area of more than 929 centimeters. Of course, this term was first used by Hughes in 1977 (poesen; 2003; 94). This type of erosion causes a sharp decrease in soil fertility and the production of abundant sediment loads. Sedimentation on the path to human structures and their destruction are among the effects of this type of erosion (Esmaeiliy et al., 2015). Recent field-based studies in Europe have evidenced annual soil losses due to EG erosion ranging between 2 to 90 m3 hectares. (Poesen et al., 1996; Casal´ı et al., 1999; Capra and Scicolone, 2002; aygarden, 2003; Valca´rcel et al., 2001). Material and Methods The sub-basin of Lanbarn is located on the eastern side of the Aharchay River. The area of this sub-basin is 20118 hectares and it is considered as the main stream of the Aharchay river. The basin was fielded between 1378 and 1395 after rainfall. There occur ephemeral gullies in this area in autumn, summer and spring, especially with heavy rainfall. Therefore, to study the ephemeral gullies, the required parameters were collected from 70 ephemeral gullies. To collect the data, field observations, GIS software and a GPS device were used. GIS, Excel and SPSS softwares were also used for processing. Results and discussion According to the results of the model, the highest correlation coefficient was observed among estimated and measured depth, estimated and measured cross-sectional area, estimated volume and measured volume, and estimated and measured widths. Therefore, the model predicts the cross-sectional area, the volume and the width of the ephemeral gullies efficiently. Also, regarding the effect of morphometric characteristics of the gully, the area of the basin, the width, watershed length and concentrated flow length of the ephemeral gully play the most part in soil erosion. The depth and cross section of temporary gullies have the least impact on erosion. Conclusion - Ephemeral gullies in this area are mostly formed in agricultural land. The intensity of rainfall, the expansion of communication routes and the burning of vegetation play an important role in the development of ephemeral gullies. - The standard version of the EGEM model can be used in the catchment area of Lanbarn. - To achieve optimal results, this model should be used in other climatic regions of Iran. - In this basin, the drainage area, width, watershed length and concentrated flow length of gully play the most part in soil erosion. - The depth and cross section of ephemeral gullies have the least impact on the eroded soils in this area. - - - - - - - - - - - - - - - - - - - -- -- --

Land use reflects the interactive characteristics of humans and the environment and describes how human exploitation works for one or more targets on the ground. Land use is usually defined on the basis of human use of the land, with an emphasis on the functional role of land in economic activities. In this research, a pixel-based method and object-oriented method for mapping the map of Maragheh city using ASTER sensor images in a 17-year time series were compared between 2000 and 2017. Also the effect of land use changes on erosion was studied. In order to compare the results, both methods used the same educational data for classification. Then, the most important methods for assessing accuracy including precision and kappa coefficient of classification were extracted, and it was determined that the result of the object-oriented classification was better than the 3% Offers. The increase in accuracy in a method based on object-oriented classification largely depends on choosing the appropriate parameters for classification, defining the rules, and applying the appropriate algorithm to obtain the degree of membership. After classifying and comparing the extracted maps, we attempted to reveal the changes that occurred during this period and it was determined, that the rangeland and Bayer classes have a decreasing trend and dense garden classes and water with increasing trend. According to the user-mapped maps of the two methods of maximum-likeness and object-oriented classification and comparison and adaptation of these maps with ground-based properties, the results of the object-oriented classification method were confirmed. According to the results of the study, with object-oriented method, in the Maragheh area, densely populated gardens, gardens, housing, housing, agriculture, industry, and communication in the object-oriented method have increased, and agronomic, pasture, dryland and Bayer uses Reduced area. This indicates the importance of agriculture and horticulture in this city. In 2000, 9.08 percent and 15.88 percent of study area are located in very high-risk and high-risk erodibility and these values for 2017 are 13.66 and 29.76 percent respectively. Also results showed that dense and low dense agricultural land, residential and industrial land use have been increased. The increase in accuracy in a method based on object-oriented classification largely depends on choosing the appropriate parameters for classification, defining the rules, and applying the appropriate algorithm to obtain the degree of membership. After classifying and comparing the extracted maps, we attempted to reveal the changes that occurred during this period and it was determined, that the rangeland and Bayer classes have a decreasing trend and dense garden classes and water with increasing trend. According to the user-mapped maps of the two methods of maximum-likeness and object-oriented classification and comparison and adaptation of these maps with ground-based properties, the results of the object-oriented classification method were confirmed. According to the results of the study, with object-oriented method, in the Maragheh area, densely populated gardens, gardens, housing, housing, agriculture, industry, and communication in the object-oriented method have increased, and agronomic, pasture, dryland and Bayer uses Reduced area. This indicates the importance of agriculture and horticulture in this city. In 2000, 9.08 percent and 15.88 percent of study area are located in very high-risk and high-risk erodibility and these values for 2017 are 13.66 and 29.76 percent respectively. Also results showed that dense and low dense agricultural land, residential and industrial land use have been increased. The increase in accuracy in a method based on object-oriented classification largely depends on choosing the appropriate parameters for classification, defining the rules, and applying the appropriate algorithm to obtain the degree of membership. After classifying and comparing the extracted maps, we attempted to reveal the changes that occurred during this period and it was determined, that the rangeland and Bayer classes have a decreasing trend and dense garden classes and water with increasing trend. According to the user-mapped maps of the two methods of maximum-likeness and object-oriented classification and comparison and adaptation of these maps with ground-based properties, the results of the object-oriented classification method were confirmed. According to the results of the study, with object-oriented method, in the Maragheh area, densely populated gardens, gardens, housing, housing, agriculture, industry, and communication in the object-oriented method have increased, and agronomic, pasture, dryland and Bayer uses Reduced area. This indicates the importance of agriculture and horticulture in this city. In 2000, 9.08 percent and 15.88 percent of study area are located in very high-risk and high-risk erodibility and these values for 2017 are 13.66 and 29.76 percent respectively. Also results showed that dense and low dense agricultural land, residential and industrial land use have been increased.

Today, along with morphometric and distance-based methods based on the use of satellite imagery and digital elevation model (DEM), several algorithms have been developed to evaluate tectonic activities in the form of land-based software applications. This includes TecDEM software. (Shahzad and Golghoun, 2011). The purpose of the present study was to investigate the morphotectonic analysis of the catchment areas of the Karkas Mountains by using TecDEM-based geomorphic components and analyzing the digital elevation model and field evidence from the analysis of morphotectonic data. In this research, we tried to study the geomorphometric methods and comparative analysis of geomorphic phenomena in catchment areas. The effect of each of these basins on the tectonic movements of Qom-Zefreh's fault system by quantifying landform phenomena using the features of the single program, Check and evaluate. Concerning the right-angled motion of the right-strike motion components along with the vertical components of the fault system (QZFS), is the study of the research problem, have the surface dynamics and catchment areas in the studied area been uniformly affected by the tectonic processes of this system? On the other hand, have the changes of the pattern tectonic processes have affected the hydrogeomorphological model of catchment areas? Therefore, in this research, using the basic principles of geomorphic tectonics, the adaptation of land surface changes to changes in their shaping processes has been investigated.

Tectonics play an important role in the evolution of large-scale gravitational phenomenon (Galadini, 2006), mainly through the formation of steep slopes. Competing tectonic and surface processes build and destroy topography in active orogens, hence, thrusting, crustal thickening and isostatic response result in rock uplift and relief production (Agliardi et al., 2013). In some cases, the faults play a primary role in increasing the local relief and their activity is an important geomorphic factor conditioning the gravitational movements (Galadini, 2006). We have studied this kind of gravitational movements and slope instabilities termed “Deep-Seated Gravitational Slope Deformation (DSGSD)”. This paper focuses on a study aimed at defining the role of structural setting, local uplift and morpho-structural evolution on the onset and development of a DSGSD that affects the western parts of the Siahcheshmeh pull-apart basin (SPAB) in a releasing bend of the Gailatu-Siah Cheshmeh-Khoy fault. DSGSD: DSGSDs are gravity-induced and large to extremely large mass movements generally affecting the entire length of high-relief slopes, extending up to 200–300 m in depth, which can frequently extend beyond the slope ridge and evolving over long periods of time. (Crosta et al., 2013). DSGSDs are not considered hazardous phenomena because they evolve very slowly. Despite their slow deformation rates, they can cause damage to surface and underground man-made structures (Soldati, 2013). The main feature that distinguishes DSGSDs from landslides is the absence of a continuous or well-defined sliding surface (Soldati, 2013) and discontinuous or poorly defined boundaries, both laterally and at their lower ends (Crosta et al., 2013). Gailatu-Siah Cheshmeh-Khoy fault: The 200 km long (Karakhanian et al., 2004; Berberian, 1977) Gailatu-Siah Cheshmeh-Khoy fault, with the same trend as the North Tabriz, Chaldiran, Nakhichevan and Pambak-Sevan-Syunik faults, is regarded as a part of the active strike-slip fault system in the Arabian and Eurasian collision zone, which extends from 42˚ E to 48˚ E with the Tutak and North-Tabriz faults in the west and east, respectively (Selçuk et al., 2016). This system includes a series of right-lateral strike-slip faults between the southern front of the Lesser Caucasus to the northeast and Bitlis-Zagros suture zone to the southwest. The available literature, fault plane solutions, offsets of various geomorphological and man-made features indicate the right-lateral strike-slip nature of the Gailatu-Siah Cheshmeh-Khoy fault. The trace of this fault is very obvious and displays a series of well-developed and preserved morphologic structures indicating recent activity of the fault, such as fault scarps and horizontal deflection in the Quaternary features, pull-apart basins, hot water springs and uplifted terrace deposits. Discussion and Our geological and structural survey along the Gailatu-Siah Cheshmeh-Khoy fault confirms the presence of a large slope instability in the western flank of SPAB. In order to understand the relationship between the nucleation and evolving DSGSDs with structural aspects of this fault, we focused on slip rate of this fault in two segments, the Gailatu-Siahcheshmeh (northwestern sector of Gailatu-Siah Cheshmeh-Khoy fault) and the Siahcheshmeh-Khoy fault segments (Southeastern sector of Gailatu-Siah Cheshmeh-Khoy fault), which overlap at a right step-over in the SPAB. Along the Gailatu-Siahcheshmeh fault segment, Quaternary lavas, known as Maku basalts, form a few ridges that are elongated parallel to the strike of the fault and displaced as a result of this fault activity by ~ 725±50. Using the about 400 kyr published age of these basalts (Pb206/U238 and Ar40/Ar39 dating methods, Allen et al., 2011; Lechmann et al., 2018), a mean slip rate has been calculated 1/65 ± 0.1 mm/yr. On the Siahcheshmeh-Khoy fault segment, we excavated a trench to determine the fault geometry and its rake, and to assess recent offsets. Radiocarbon dating of the youngest deposits in the stream wall that displaced right-lateral by 42±4 m, yield 6764±283 calBC, suggest a horizontal slip rate of 4.6±0.3 mm/yr. The northwestern and southeastern terminations of Siahcheshmeh-Khoy fault segments form the eastern and western flanks of SPAB, respectively. Hence, the higher slip rate of Gailatu-Siahcheshmeh fault compared to Siahcheshmeh-Khoy fault, causes uplift of the western SPAB sectors. This is accompanied by thrust faulting in a general northwest-southeast direction as a splay configuration at the termination of Siahcheshmeh-Khoy fault. Consequently, local uplift has been taken place in the western flank of SPAB that is readily obvious from higher altitude of this flank relative to the eastern side. Therefore, DSGSDs have been occurring almost along the entire slopes facing the pull-apart basin. On the other hand, decreasing altitude in the SPAB in the releasing bend and normal faults are additional controlling and intensifying factors for DSGSD. As a result, most of the expected structural features, especially splay strands of Siahcheshmeh-Khoy fault and normal faults at the margin of SPAB, have been covered by DSGSD phenomena. Therefore, except at a small part of the southwest of SPAB, we could not find exposure of normal faults along the western flank.