Received 5 December 2015; accepted 12 February 2016; published 16 February 2016
The Saravan fault is located in southeastern of Iran, 27˚ to 28˚ latitudinally and 61˚30' to 62˚45' longitudinally, that has been covered by Geological Map 1:250000 of Saravan. In Iran Geotectonic Division, Saravan fault and range of its performance are located in the suture zone and the east Flysch of Iran  . The most stones of this region are Flysch of Eocene. This fault where extended from southern range of eastern-northeastern Saravan Mountains and it is 270 km in length and it has been separated Ba Damo-Siyahan Stratigraphy from Saravan  .
The most part of the fault (all over Saravan plain) has been extended N50W generally and about 50 km of the east fault (Dehak village to the border between Pakistan) is also roughly along east-west  . It is found that the main process of these structures such as the main faults and the axis of folds and so on have been extended from northwestern to southeastern, i.e. N130 - 140.
Morphotectonic Analysis Method that used to interpret Geological Structures has been derived from Topographic and telemetric data. This analysis examines the spatial relationships between geomorphic (Landforms) and tectonic processes that caused them to be explored and understood by morphotectonic in line with the dynamic or kinetic analysis of the process.
2. The Study of Geomorphologic Indices
Some quantitative indices has been considered to evaluate the ratio of regional geotectonic activities, these indices indicate movement and the lifting-up one block than other block along the fault and the morphotectonic indications resulted from reversed and extensive movements.
In this study, it has been considered not only the new Neotectonic evidences but also parameters such as mountain front sinuosity, rate index (V), Hydraulic Gradient Index of Stream (SL) that has carried out by field sampling, using satellite images, topographic maps and SRTM data.
3. The Mountain Front Sinuosity
In fact, this index indicates the relationship between tendency of rivers to produce irregular front and vertical tectonic activity to produce a straight front on the mountain that surface processes lead to be irregularity of mountain front and morphotectonic processes lead to be these fronts smoother and linear. On the basis of resulted calculation, the range of changes is a value higher than 1 in Smf index  .
If the ratio of this index is equal to 1, it indicates the activity of geotectonic processes and the youth of mountains, and the increase of this value indicates lesser tectonic activities and the dominance of surface processes. In following relation, Lmf value is the length of the mountain meander and Ls is the length of the straight line of the mountain front.
According to the above mentioned features, the mountain front along Saravan fault is divided into 20 fronts to obtain Smf index (Figure 1 & Figure 2) and then it was calculated from the north to the south for different blocks.
By studying these mountain fronts parallel to the fault, it can be determined relative intensity of tectonic activities and also the regions where have the highest lifting-up activities. Then, by comparing it with the other morphometric indices, it can be obtained the final sum up of the condition of activities and active regions along the Saravan fault (Table 1).
As it is shown in Table 1, the calculated values of Smf index are small (lower than 2) for the most regions and it indicates that there is a relatively high activity in all over Saravan fault. In some regions, the considered mountain fronts have the least curve and sinuosity then their Smf index is the least.
Figure 1. Method of calculating Smf index.
Figure 2. Calculation of Smf index, measurement direction is from north to south.
Table 1. The derived numbers for Smf index in 24 blocks of Saravan fault.
Since it can’t be expressed any idea about regional tectonic based upon one morphotectonic index, thus other indices were also studied and finally the regional tectonic activities are determined according to all indices.
4. Truncated Percentage Index (Cleaved) of Mountain Front (Facet %)
Facet is defined as triangular or polygonal morphology which is created between the interfaces of two channels resulting from erosion (Figure 3). The percentage of Facet in mountain front or F% is another quantitative index to determine relative geotectonic activity of a mountain front which consists of the ratio of the total length of facets in a mountain front (Lmfd) to the length of straight line against it (Ls). A mountain front is considered truncated if it includes several definite and significant channels, and divided it into distinct facets and clearly causes to deviat the contour lines in the topographic maps. Some of the facets might have been eroded by channels which cause to be irregularity of its morphology that it is found in the maps and images. More active Mountain fronts have less truncation and appear as a continuous mountainous front with great and definite facets that have little internal truncation.
The percentage of truncation or F% has been approximated 100% in the active mountain fronts which it also reduces by decreasing the activity. This index is ranged 0% - 100% (Table 2). The mountain along Saravan fault was divided into 28 blocks to obtain faceting index (F%) and then this index was calculated for different blocks from north to south.
Figure 3. Method of calculating facet index  .
Table 2. The obtained numbers for index F% in 28 blocks of Saravan fault.
Most active regions based on these calculations are located on the middle part of Saravan Fault. It means that the least faceting of mountain front has been created by current streams in the region. But the southern and northern regions of Saravan Fault have lesser values of F% and according to this index; the level of relative lifting-up activity is less. These results are similar to the results of mountain front sinuosity significantly.
5. Length-Gradient Stream Index (SL)
Stream gradient index is one of the indices to evaluate active tectonic in relation to the shape of stream channels that defined with the following equation: 
where ΔH/ΔL is local stream gradient (ΔH is height differential of the range and ΔL is the range length) and (L) is channel length from the straight stream line to the center that its gradient is calculated or in other words, (L) is the total channel length from the point which index is calculated (Figure 4) to the highest channel point  . SL index in different parts of fault has been studied for 10 basins. The given index shows different SL values in different parts of the fault which are given in the rows 1 - 10 of Table 3. The obtained results represent that the fault is active.
Table 3 shows the obtained values of SL for the important streams that passes through Saravan fault, therefore are affected by its tectonic movements. As a whole, SL values have a different shifting range on the basis of structural geotectonic activity and its position in a structural state, topographic status, rate of flow and the lithologic type.
Figure 4. Method of calculating SL index.
Table 3. The obtained numbers for SL index.
As it can be shown in the above table, SL values have significant numerical differences as compared with each other. Because the value of SL difference can be unexpected in different streams regarding some parameters such as stream length. It can be also observed that streams 4 and 9 where located in the middle part of Saravan fault have high SL index.
Therefore, it can be said that geotectonic movements, especially lifting-up movements have increased the stream gradient. As well, an important factor such as lithologic matter can be also influenced on the SL value. Hence, more resistant lithology, more SL value and less resistant lithology, less SL value.
6. Ratio index V
V ratio index represents a relation between the width of valley bottom near to the mountain front and the mean height of its top chords that can be calculated the ratio of tectonic activities and regional erosion relatively and V ratio can be obtained by using this method (Figure 5) as following:  
This ratio gives us information about regional lifting-up so that if its value is about 1, the valley is U shape or if the value is very small, the valley is V shape that can indicate the vertical performance of active tectonics. The values more than 1 represent that the width of the valley is higher than its depth and regional tectonic activities has been approximately stopped.
The measured values of V index are shown for the case study of region in (Table 4).
According to the above table, the least value of V ratio is related to Saravan regions and zarawy and sheikh Rostam villages. In the regions with high V ratio, Meander streams are flowing and in the regions with low V ratio, the streams are flowing in strength with less width.
7. Study on Tectonic Earthquake
7.1. Historical Earthquakes
As history of past earthquakes with available experimental instructions can be used to access knowledge of his-
Figure 5. Indicates the method of calculating parameter of V ratio.
Table 4. The values of V index along Sravan fault.
torical earthquakes, it is clear that there are many problems in this field. The former historians have explained some earthquakes similarly then it causes to estimate the past earthquakes imprecisely. Studying on history of the past earthquakes (ancient, historical and 20th century) in an area is one of the important and basic data to estimate and to know historical data from descriptions of earthquake events. It is clear that the history of earthquakes in long intervals should be gathered and studied to access the specifications of tectonic earthquake. Historical earthquakes can be known on the basis of historical books but there isn’t any precise data about the exact focus of earthquakes because there weren’t any modern instruments and tools at the last time. Some historical earthquakes have been recorded in the studied area as following:
In this year, a strong earthquake occurred in Sistan that many houses were demolished. Hirmand River was dried twenty years after this earthquake and famine and pestilence were other dangerous problems in this time  . Also, Barbarians record in the first “Iran Earthquakes Catalogues” that Sistan was shaken with an earthquake in magnitude of 7 and 7.5 Richter for two times in 815.
In Qajar period (1838), a horrible earthquake also occurred in Sistan (along eastern boundaries of Loot Plain). It destroyed the deserts and covered about 150 km. on the westside, the lowlands were destroyed as caravan ways weren’t secure and safe till many years later. All villages at a distance of 75 km from this earthquake were destroyed but some people were killed and its aftershocks have been lasted for two years  .
7.2. Instrumental Earthquake
It is necessary to study on recorded instrumental earthquakes that occurred between 1900 and 2009 to determine earthquake formula of the studied area. We have taken the earthquakes in magnitude of more than 3 in local scale. It covers an area, 150 km, with a capital named Saravan with geographical coordinates (Longitude: 62˚33' and Latitude: 27˚36'). This information has been gathered with international catalogues (ISC, ISS, and USGS), earthquake databank and International Institute of Earthquake Engineering. Instrumental Earthquake Databank of this area has been presented in attachment 1 (Figure 6).
8. Maximum Calculation of Peak Ground Acceleration (PGA)
Peak Ground Acceleration (PGA) can be calculated by Formula 5 and its result has been presented in Table 5.
Figure 6. Focal location of instrumental earthquakes in Saravan between 1900 and 2014.
Table 5. Estimation of maximum calculation of peak ground acceleration on Saravan from surrounding fault.
According to these calculations, Maximum Calculation of Peak Ground Acceleration that imposed on Saravan is 0.92 of Gravitational Acceleration and it is a result from Saravan Fault Actions.
According to Satellite Images, Saravan Fault is a dextral reverse fault. There aren’t many instrumental earthquakes in this area and surficial focus is located on central district then it has more earthquakes than other areas. This fault passes through 12 km from Saravan that it is the most dangerous fault for this city. It can shake this area (270 km) with magnitude of 8.29 that its acceleration will be equal to 0.92 of gravitational acceleration. By combination of the whole morphotectonic and earthquakes data, it can be determined the tectonic activities of this region. The studying of regional morphotectonic units indicate that the whole region is active tectonically and regarding to the studies, Saravan fault has more activities in the middle part than its two northwestern and southeastern extremes. The mean calculated values of geotectonic indices of mountain front sinuosity, the ratio of the valley bottom on its height and the V ratio are 76.1, 0.98 and 77, respectively. These values indicate high tectonic activity in Saravan region.