Received 19 May 2016; accepted 25 June 2016; published 28 June 2016
One of the issues in tectonic geomorphology is Earthquake risk assessment in specific areas through the study of the problems associated with active faults and determination of seismic zone status in upper Pleistocene and Holocene  . The ongoing competition between tectonic processes that tend to create topography surface processes that tend to erode and disintegrate them represents a tectonic geomorphology  , Neotectonic studies phenomena including all the factors, processes and practices arising from new activities of the Earth (on a geological scale) and forms created by these activities, moreover, by this type of investigation, relatively stable and tectonically inactive regions can be identified. This relatively new branch of applied geomorphology can be used in solving the problems caused by tectonic geomorphology where landforms of roughness and evolution of landscapes caused by new tectonic activities are discussed  .
As Iran is located on one of the two earthquake-prone belts in the world and a lot of faults, earthquake in Iran tableland is normal. The main cause of earthquakes is Neotectonic movements which can be effective in other natural disasters such as range movements  . Due to the severe shortage of historical and seismic data and device in different parts of the world, particularly Iran, need for introducing more researchers and planners of Iran to tectonics and active tectonic areas is more visible. Qualitative measurement is not responsible for evaluating these activities. Quantitative and numeric assessment make it possible for geomorphologists to evaluate different landforms truly and tangibly; some tools are introduces for this aim, morphometric indexes are among them.
Assessment of research resources is one of the important tasks in any research work. Regarding tectonic activities using geomorphic indexes, a lot of studies have been conducted in Iran and other countries such as Ghanbari  who investigated earthquakes in relation to active faults, especially quaternary faults in Azerbaijan using mathematical relationship on fault of Tabriz and Salmas. Madadi, Rezayi Moghaddam and Rajayi  analyzed tectonic activities in northwestern slopes of Talesh using geomorphological methods and for this purpose, geomorphic indexes are used such as Mountain Front, ratio of valley width to its height, river gradient index and etc. Mokhtari  studied tectonic and expanding holes of Marand from sedimentary tectonics the perspective. In his study, role of tectonic factors in creating tectonic holes in Marand is explained by analysis of data related to sedimentary morphological factors. Roostae and Nayeri  analyzed effect of lithology and tectonic on longitudinal side view of Mahabad River and for this purpose, Geomorphic indexes were used such as SL index, standardized index SL/K, drainage density index, and concavity index.
Bull  studied Neotectonic activities and distinguished active regions and inactive regions by sinuosity of Mountain Front Methods (Smf), ratio of valley width to its height (Vf), and Drainage basin shape (Bs). Keller, Larry and Tierney  presented geomorphological indexes related to lateral growth of earth folds. There are six geomorphological indexes including: drainage density, cutting head of folds, reduction in topographic height of fold head in side view, reduction in side slope of the fold, deformation of young alluvial deposits, and growth of certain drainage patterns. Jamieson, Sinalair, Kirestein and Purves  in their article investigated the effect of tectonic factors on split Ratio, drainage density, hypsometric integral of rivers network in Himalayas in northern India. Singh and Jain  conducted a research entitled tectonic pressure on development front plates in edges of the stream in North West of Himalayas in India. Reyaz Ahmad, Shakil Ahmad, Rakesh and Ishtiaq  studied geomorphologic tectonics of Karwa basin in the valley of Kashmir using geomorphic indexes such as SL, Smf, RC, Rb, D and AF; they found that there is a high degree of active tectonics in the area.
Khafr basin is a narrow basin located between two mountain ranges namely Simakan and Meymand on one side, and White Mountains (Goshtasb) on the other side and Ghareaghaj river runs through it. Existence of this river caused a desirable climate a blooming environment in this area. Knowledge of local neotectonic conditions can be effective largely from earthquakes’ life and property damages and damages caused by natural disasters such as range motions.
Therefore, aims of this study are tectonic activity evaluation and investigation of tectonic activities in Khafr basin with regard to the validity geomorphic indexes through quantitative approach.
2. Studied Area
Khafr basin geographically is located is located on 52˚42'43" to 53˚30'3" longitudinal east and 28˚5'56" to 29˚8'20" northern latitude. Khafr basin in Jahrom city, Iran, is about 1182.5 square kilometers and one of the most populous and important part of Jahrom city. The maximum height of the basin is 3163 meters, its minimum height is 1165 meters, and its average height is 1735 meters. The main hydrographic network and the only major flows within the Ghareaghaj River is approximately 59 km (Figure 1).
Figure 1. Location of the studied area in a map of Iran and Fars.
3. Geological of Studied Area
Investigation of geological map of 1:250,000 Shiraz and Jahrom showed that the case area is consisted of several formations that are different in terms of resistance, these formations include: Bangestan (Bgp) is often of limestone, crack systems and karst phenomena usually are developed in them. Tarbur (Ktb): This formation is related to Upper Cretaceous, mainly is composed of massive reef limestone and its color is cream-colored, Gurpi (Gu): In the cutting of this formation is seen 320 m gray-blue marl and shale and secondarily layers of clay limestone, this formation has little resistance against erosion and has created the gentle topography. Pabdeh (Pd) This formation is related to Eocene and consists of clay sediment in terms of lithology. Sachun) Pesc( Mainly includes white to gray gypsum, yellow dolomites, silty lime stones and red marl. The upper part of this formation is composed of massive and marly gypsum with dolomite. Asmari Jahrom (OME as ja): Jahrom formation mainly composed of fossiliferous limestones with good layering and sometimes with layers of marly limestone, Asmari formation is mainly consists of cream to brown resistant limestone with shale layers. Age of this formation is considered Oligocene to lower Miocene. Razak (Mrz) is related to lower Miocene and is composed of red, green and gray Marly rocks, in some cases gypsum can be found. Aghajari (Mpl aj) Limestone it begins with brown to gray sandstone layers and a thickness of 2 to 5 meter , and continues with a thick layer of red marl and thin layers of sandstone and silt, sandstones of this formation are usually highlight mode but marl and silt have been deep erosion. Bakhtiari (Pl bk): age of Bakhtiari formation is the late Pliocene to early Pleistocene that is composed of conglomerate and sandy marly conglomerate and current alluvium. Sediments of Quaternary exist on the surface and different heights of each type of rock from any era and created several low steep lands which are separated by rough stone heights and formed big and small fields. Most residential and agricultural areas of Khafr are built on current alluvium which encompasses approximately 372.5 km2 that is 32% of Khafr section (Figure 2).
4. Methodology and Materials
In this study aerial photos with scale of 1/55,000, satellite images, geological maps 1/250,000 of Shiraz and Jahrom sheets, moreover Digital Elevation Model (DEM) is used to analysis and extraction of tectonic activity
Figure 2. Geological map of studied area.
indicators and drainage networks. Arc GIs 10.2, Excel and Global mapper soft wares were used to prepare different maps, finding required data and data analysis.
In this study stream length-gradient Index (SL), drainage basin asymmetry ( Asymmetric factor) (AF) , Index of drainage basin shape(BS), ratio of valley floor width to valley height index (VF), drain density index(D), river sinuosity index (S) and transverse topographic symmetry index (T) is used to investigate the effect of tectonic activities on valleys and rivers and as well as calculating the required indexes in Khafr basin, main stream, and subsidiary stream (By Strahler method).
4.1. Stream Length-Gradient Index (SL)
It is one of indexes for tectonic force assessment which is stated by the following correlation:
SL = Stream length gradient
∆H/∆L = Stream slope or gradient of a sheet
∆H = height difference of the case sheet
∆L = branch length of the case sheet
L = length of the river from central point of the measured sheet to the headstream  .
Theses parameters are measured by topographic maps and digital elevation model (DEM). The values of SL index were classified into three categories: 1 (SL ≥ 500), 2 (300 ≤ SL < 500) and 3 (SL < 300)  .
4.2. Ratio of the Basin Form Index (BS)
Tectonically active basins have elongated shapes. Form of the basin gradually becomes circular by a pause in amount of uplift over time  . This index is divided into three classes of activities: BS > 4 active, 3 < BS < 4 semi-active and BS < 3 inactive  which is computed by Equation (2):
Bi = the basin longitude: the distance between lowest height of the basin and its farthest point
Bw = width of the basin: the widest part of it
4.3. Drainage Basin Asymmetry Index (AF)
Asymmetry index for determining the tectonic tilting of the drainage basins and larger areas are as follows in Equation (3)  :
Ar: right area of the river, At: total area of the drainage basin.
Af-50 < 7 inactive, 15 > Af-50 > 7 semi-active and Af-50 > 15 is placed in active class.
4.4. Drainage Density Index (D)
Drainage density is the ratio length of the stream to area of the area as follows (Equation (4)):
D = Drainage density index
L = total length of streams in km
A = total area in square kilometers
This index indicates rupture of topography which can be related to tectonic, type of bedrock and amount of rainfall 
Low density is between 0.6 and 0.6, average between 0.6 and 0.9, high if it is above 0.9  .
4.5. The Ratio of Width to Depth of Valley (VF)
This index is calculated by Equation (5):
Vfw = width of the valley floor
Erd = average height of water dividing line on left and right side of valley
Esc = average height of the valley floor from the open waters 
This index reflects the difference between V-shaped and U-shaped valleys, this means that in V-shaped morphology of the valley Vfw is closer to zero which indicates destruction and that in U-shaped morphology of the valley, amount of Vfw increase which indicates widening and leveling of valley and low intensity of tectonic processes and sustainability of valley  . Based on Hamdooni’s classification  Vf less than 0.5 is active, between 0.5 to 1 is semi-active and greater than 1 is inactive.
In this study, to calculate Vf index, Global mapper software was used by eight cross-section on the main stream and sub-streams in the highlands of 1370.82, 1560.77, 1661.84, 1726.64, 1857.26, 2103.88, 2295.35 and 2338 that show the valley floor heights, they are drawn perpendicular to the stream, then the necessary calculations have been performed on them.
4.6. Sinuosity of River (S)
Sinuosity index is calculated by Equation (6):
S = Sinusitis or the twists and turns of the river
C = length of the river
V = length of the valley  .
Rivers that have a lot of twists and turns are very close to equilibrium, while the direct river route represents younger and existence of neotectonic activities in the area. The high values of this index shows that the river is closer to equilibrium while, low value of this index indicate tectonic activity in the case area  .
4.7. Transverse Topographic Symmetry Factor (T)
Index T is a vector with a certain orientation and different values between 0 - 1. In perfectly symmetrical Basins, the value of this index is zero. By reducing the symmetry of the basin, T index increases and gets closer to 1 which is computed by Equation (7):
Da = the distance between Midfielder of drainage basin and main river route
Dd = the distance between Midfielder of drainage basin and Water dividing line  .
4.8. Active Tectonics Index (Iat)
After calculating the geomorphic indices in the study area is estimated tectonic activities in the area using Iat index. Iat index is achieved from the average of different geomorphic indicator classes according to the following equation:
S: Total classes geomorphic indicators that is calculated
N: Number of indicators
The values of the index were divided into four classes to define the degree of active tectonics: 1―very high (1.0 ≤ Iatb1.5); 2―high (1.5 ≤ Iatb2.0); 3―moderate (2.0 ≤ Iatb2.5); and 4―low (2.5 ≤ Iat)  .
To calculate the amount of longitudinal slope of Khafr basin in three altitude on the main stream, 4 height areas were considered on sub-stream, and longitudinal side views were drawn (Figure 3) and SL values for each of them were calculated and the results are given in tables (Tables 1-4).
According to results of SL index value, the main stream can be classified according to Hamadooni  that is
Figure 3. Longitudinal profiles of SL.
Table 1. Values of the SL classes (stream-gradient index) in 3 elevation ranges in the main stream.
Table 2. Values of the SL classes (stream-gradient index) in 4 elevation ranges in Sub-stream.
Table 3. Values of the SL classes (stream-gradient index in 4 elevation ranges in Sub-stream 2.
tectonic tranquility in Class 3, sub-streams in Class 2 which means a moderate tectonic activity. In general, sub-streams at lower altitudes have more tectonic activity (Figure 4).
According to calculations (Table 6) AF index in Khafr basin is 57.75 that if it is deduced from 50, it will be is
Table 4. Values of the SL classes (stream-gradient index) in 4 elevation ranges in Sub-stream 3.
Table 5. Value and tectonic activity class of BS index in Khafr basin.
Table 6. Value and tectonic activity class of AF index in Khafr basin.
Figure 4. SL, AF, BS indexes along the drainage network.
equal to 7.75 that represents a tectonic function on the right bank of the main stream as well as greater length of sub-drainages on this side and it is classified in class 2 in terms of activities which shows a semi-active tectonic in the area (Figure 4).
To calculate Vf index, was used by eight cross-section on the main stream and sub-streams in that show the valley floor heights, they are drawn perpendicular to the stream, then the necessary calculations have been performed on them (Table 7). Based on the results, dominant tectonic activity in the area is active and semi-active. This indicates that tectonic activities did not let erosion to happen in the rivers of the region.
Regarding the fact that tectonic activity is directly related to increased Vf index from upstream to downstream. By taking this index into account, tectonic activity is active in heights; in other parts of the basin tectonic activity is semi-active which confirms Hamadooni’s classification, so according to this classification, the case basin is semi-active in terms of tectonic activity. Due to resistant lithological units in upstream, V-shaped valleys are often observed (Figure 5). According to Iat index, Indicators mentioned (SL, Af, Vf, Bs) are in class 3, that represents tectonic activity in the area is average (Table 8).
In this study, Khafr basin is divided into two ranges: northeast and southwest. Based on ranking the length of sub-streams were computed for each range separately and ultimately for the entire basin. High level of index reflects more tectonic activity. The following results were found: for the northeast (1.32), the Southwest (1.4) and the whole basin (1.3) (Table 9 and Table 10).
According to the index, the case basin in both the North East and South West ranges and in whole, have a high drainage density that is indicates an active tectonic in the area (Figure 4).
To estimate the required data for T index in the case area, values of Da and Dd were calculated in 10 separate sections and the results are presented in Table 11.
Given that the index value is zero which represents a symmetric basin and the results of the calculation is greater than zero, asymmetry in the basin proven, so it can be stated that the case area is tectonically semi-active.
Other indexes used to assess tectonic activities in Khafr basin, the sinuosity index of main stream. Thus, according to the obtained value from this index it can be claimed that the case area tectonically is active and did not reached equilibrium (Table 12).
Table 7. Values of the VF classes in Khafr basin.
Table 8. Classification of geomorphic indicators according to Iat index.
Table 9. Length of streams on the slopes of basin (Km) in Khafr basin.
Table 10. Drain density (D) in Khafr basin.
Table 11. Values of the T index in Khafr basin.
Table 12. Values of the S index in Khafr basin.
Figure 5. Latitudinal profiles to calculate VF index.
In this study, based on geomorphic indexes the roles of tectonic activities in Khafr basin area were investigated. With regard to the value of SL index, main stream is in Class 3 that is tectonic tranquility, sub-streams are in class 2 with the moderate tectonic activity. Generally, sub-streams at lower altitudes are tectonically more activity. that reflects the asymmetry of; according to calculations of AF index, tectonic in the case area is semi-active, and tectonic activity is much more on the right bank compared to the left side of the main stream which indicates longer sub-drainage on this side. The shape of the basin (BS) in the case area was 2.7 that proves inactivity of the basin. VF index is measured for the main stream and 3 sub-stream in various locations. Based on the results, dominant tectonic activity in the area is active and semi-active. Indicators mentioned are in class 3 based on Iat classification, that represents tectonic activity in the area is average
To calculations of drain density index, length of Khafr drainage basin in the northeast and southwest ranges was calculated separately and totally. Based on the results of investigations the northeast (1.23), the Southwest (1.4) and the whole basin (1.3) were found which indicates the amount of high-density and therefore high tectonic activity. To compute T index in the case area, 10 separate sections were measured and the results show that the case area is semi-active tectonically. Besides, based on the sinuosity index of main stream it can be claimed that the case area tectonically is active and did not reached equilibrium.
Generally, intensity and role of tectonic activity in various parts of the basin is different, the results of geomorphic indexes and tectonic evaluation prove these differences. Finally, the area can be considered as a moderately active region in terms of tectonic activity.