Morphometric Analysis of Kakoi River Watershed for Study of Neotectonic Activity Using Geospatial Technology

Morphometric analysis is defined as the quantitative measurement of landscape shape. Morphometric studies of a watershed helps the researcher to compare different landform and calculate the geomorphic indices that may be useful for identifying a particular characteristic such as the level of neotectonic activity in a watershed. The present study has been carried out in Kakoi River Watershed which is a part of Lower Subansiri River Basin in Assam. The study area is geo-dynamically unstable region characterized by active faults, continuing crustal movements and complicated structural region which fall under the seismic zone V of India. In recent century, the Lower Suabsiri River Basin has gone tremendous morphological changes due to active tectonics activities. After the Assam Earthquake of 1950, some of the rivers like Subansiri and its tributaries Dirgha, Kadam and Kakoi suddenly changed their channels giving birth to new channels. The present study is to examine the neotectonic domains and the changes of river course. The geomorphic indices such as Hypsometric Integral (HI), Elongation ratio (Re), Asymmetry Factor (AF), Stream-length Gradient Index (SI), Mountain Front Sinuosity (Mfs), Basin shape index (Bs), Ratio of Valley floor width to Valley Height (Vf) and Channel sinuosity (S) etc. have been studied from remote sensing data. The results shows that the study area is under active tectonic area based on Vf (0.301), Rl (0.461), Bs (3.6), Sl (165.73) and other parameters such as S (1.38) and Mfs (2.02) indicated as moderate active tectonic region. Utilization of geospatial technology and remote sensing data in the present study becomes more reliable and helpful in analyzing, monitoring and understanding the landform changes in a watershed.


Introduction
Morphometry is the mathematical and geometrical representation of earth's surface and shape dimension of earth's landform [1] [2] [3].The study of morphometry analysis were initiated by Horton [4] [5] and Schumm [6] in the field of hydrology and are widely used by hydrologist and geomorphologist till date to study the spatial variability in the river basin.The morphometric analysis studies the basin processes and basin characteristic and describes the history of the river basin.It provides quantitative description of the river basin to understand initial slope, structural controls, geology and geomorphology [7].And thus it proved to be initial steps to understand the basic of basin dynamics.
Morphometric analysis on watershed of a river basin is most preferable and advantageous as it is the basic unit of hydrology rather than analyzing on the individual channel or inconsistent segment areas of the basin [8].Watershed is the appropriate unit to study several processes of the land surface because the major runoff is conveyed to single outlet [8].The morphometric analysis of a watershed reveals valuable information regarding to hydrological processes, geomorphic processes and development of the land surface [9] [10].Also it is helpful in predicting watershed maturity, erosion intensity, hydrological peak etc.
Morphometric analysis is very helpful to understand the underlying factors such as exogenesis and endogenesis factors which influence the hydrological behaviors [11].The endogenesis processes are the complex integration of the effects of the vertical and horizontal motion of earth crust and erosional processes [12].The morphometric analysis based on the geomorphic indices helps in establishing the endogenesis processes in a watershed [13] [14].Several geomorphic indices combine together to highlight the tectonic activity and provide relative assessment on the degree of tectonic activities in a watershed area [15].Tectonic plays an important role in evolution and modifying the drainage basin and is well reflected by structural, fluvial and morphometrics parameters [16].The geomorphic indices have been developed as a basic reconnaissance tools to identify areas experiencing rapid tectonic deformation [15] [17].
The Indo-Ganga-Brahmaputra fluvial plain is subjected to active tectonics processes from the time of its form [18].The active tectonic and climate change plays an important role in the development of geology and geomorphology of this region [19].The Himalayan belt from Kashmir to Arunachal Pradesh is being tectonically active with frequent earthquakes and neo-tectonic movement [20].As a result, the morphology of the river in this region is influences by tectonic activities.
The Brahmaputra River is one of the largest rivers in the world and is characterized by heavy floods and sediment load, flood damages and instability [20].
Brahmaputra basin is located in a geo-dynamically unstable region characterized by active faults, continuing crustal movements and complicated structural region which fall under the seismic zone V of India [21].According to plate tectonics, the Indian plate moving in the north-northeasterly direction is under The channels of the Brahmaputra underwent several major transformations as short as well as long terms of landforms.This includes either shifting of the river courses or rise in the river bed.Large and small tributaries from both north and south joined the Brahmaputra River while flowing through Assam, India.The northern tributaries such as Subansiri River, Ranganadi River, Jia Bhareli River, Sankosh River, Pagladia River, and the Manas River come from fragile Himalayas with highest rainfall and have steeper slopes [21].In general, the drainage anomalies characterized by migration of river course, paleo channels, etc. are common in the northern tributaries of Brahmaputra River [22] [23] [24] [25].
The Subansiri River is one of the important north-bank tributaries of Brahmaputra River.In the recent century, the lower Subansiri River basin has gone tremendous morphological changes due to active tectonics activities.After the Assam Earthquake of 1950, some of the rivers like Subansiri and its tributaries Dirgha, Kadam and Kakoi suddenly changed their channels giving birth to new channels destroying villages, agricultural fields or tea-plantation [26].
Geospatial technologies such as GIS and remote sensing data have been used successfully to evaluate the terrain and morphometric parameters of a watershed with easier and better accuracy [27].The analysis for the present study is carried out by GIS and remote sensing data to calculate the geomorphic indices of the watershed digitally and precisely.This helps to generate data to understand the watershed characteristics and changes of the landform.Thus it provides an insight into hydrologic condition of the watershed necessarily for developing watershed management strategies [27] [28].
The main objectives of the morphometric analysis of the Kakoi River watershed are: • To examine the active tectonic domains in the watershed; • To examine the river course changes due to the active tectonic activities using geospatial technology.

Study Area
The Kakoi River watershed extends from 27˚25'23.24"Nlatitude to 27˚14'12.83"Nlatitude and from 94˚5'45.79"E longitude to 94˚6'46.55"Elongitude.About 29.3% of the Kakoi River is located in the Papum Pare district of Arunachal Pradesh and rest is located in the Lakhimpur district of Assam (Figure 1).The total area of the watershed is about 140.11 sq.km.The source of Kakoi River is from the Siwalik range of lower Himalaya.The river travels about 38.17 km from the source till it merges with Subansiri River at Badhakarh Kuchiamari (Figure 2).
The Kakoi River flows about 11.17  above the mean sea level (Figure 2).
The study area can be divided into the three physiographical characteristics, i.e., the hilly tract, the foothill region and alluvial flood plain [29].The hilly tracts comprise Siwalik sediments of lesser Himalayas.The foothill region is characterized by older terrace deposit called Bhabar Zone which runs parallel to the foot hills of the lower Himalaya [30].This terrace deposits are characterized by undulating surface comprising boulders, pebbles of quartzitic and gneissic rocks with fine sand, silt and clay acts as metrix.The percentage of the boulders and pebbles diminish southwards from the foothills and the percentage of the sand and the silt increases towards the south [30].The third physiographical characteristics of the study area is alluvial flood plain consists of younger and older alluvial deposits (Figure 2).This region represents various sub-features such as Palaeochannel, swamps/marshy land, river terraces, flood plains, point bars, channels bars and river channel etc.The slope of the study area drops from the northern and eastern corner to the south.
The climatic condition of the study area is subtropical and humid with high rainfall.The study area receives SW monsoon rainfall from the month of April and continues up to September/October.The average annual rainfall is about 1300 mm [31].The highest rainfall area is located near the foothills of Arunachal summer temperature is about 29˚C to 36˚C [31].The relative humidity is about 74 to 89 percent with a mean of 81 percent [30].The winter is characterized by scanty rainfall and foggy morning and afternoon with minimum temperature of 10˚C [31].

Materials and Methodology
The present work is base on the morphometric analysis of Kakoi River Watershed to examine the neotectonic activity and river course changes.To access the morphometric parameters, the Cartosat DEM (2015) with resolution 1 arc second (30 m) has been used.The morphometric parameters have been generated in GIS environment using ArcGIS 9.3 software.The Cartosat DEM has been used for the present study because it gives higher value of stream parameters than other DEM [32] [33] [34].
The morphological characteristics of the watershed are calculated with ArcGIS 9.3 software and it provides main input for the present research (Table 1).The Basin area (A) is the entire area that is comprised by the drainage basin boundary.Basin perimeter (P) is the total length of the basin boundary and basin perimeter is the length of the main stream order, plus the length of the extension between the end of the main stream order until the longer point of the basin, parallel to the main stream order.
The stream number is counted according to Horton [5] and was developed by Strahler [35], which is most commonly used method today (Table 2).In this smallest headwater tributaries are called first order stream.When two first order streams meet, they together create second order Stream.And again when two second order streams meet, they create third order stream and so on.

Mountain Front Sinuosity (Smf)
Bull and Mc Fadden [17] proposed Mountain front sinuosity (Smf) to study the balance between the tendency of stream and slope processes produce an irregular mountain front and vertical active tectonics that tends to produce a prominent straight front.Mountain front sinuosity (Smf) is the ratio of the total length of the mountain front as measured along the foot of a mountain and the straight line length of the mountain front.
where, Smf = Mountain front sinuosity; Lmf = Total length of the mountain front; Ls = Straight line length of the mountain front.

Channel Sinuosity (S)
Channel sinuosity is the ratio of channel length and river valley length.The following formula is used to calculate channel sinuosity (S) where, S = Channel Sinuosity; SL = the stream length; VL = the valley length.

The Ratio of Valley Floor Width to Valley Height (Vf)
The valley floor width to valley height ratio is another index to access the area on the basis of the tectonic activity [15] [17].This index uses one vertical and one horizontal dimension at a given point along the stream in the erosion system.
The valley floor width and valley height ratio is defined as follows ( ) ( ) where, Vfw = the width of the valley floor; Esc = the elevation of the valley floor or stream channel; Eld = the elevation of the left valley; Erd = the elevation of the right valley.

Stream Length-Gradient Index (SL)
The stream length-gradient index (SL) is calculated along a river and used to evaluate the erosion resistance of the available rocks and relative intensity of active tectonics [37].The stream length-gradient index is proposed by Hack [38].
To calculate this index, the following formula is used For a Section of watershed For Total watershed where, ΔH = the change in elevation of the reach; ΔL = the length of the reach; L = the total length from midpoint of the reach of interest upstream to the highest point on the channel.
For the present study, total watershed formula has taken as consideration.

Elongation Ratio (Re)
Relative young basins in active tectonic areas are more tend to elongation in shape.Ratio of diameter of a circle (Dc) having same area as the basin to the maximum length (Lbm) of the basin [6].The calculation formula is: where, Dc = the perimeter of the basin; Lbm = basin length of the basin.
Values of the Rl index less than 0.5 indicate tectonically active regions, values between 0.5 to 0.75 moderate active regions and values larger than 0.75 inactive ones [39].

Drainage Basin Shape Index (Bs)
Drainage basin shape index is defined as the ratio of the length of the basin to the width of the basin at its widest point.
Bs Bl Bw = where, Bs = Basin shape; Bl = length of the basin; Bw = width of the basin.
High value index imply elongated basin indicating high tectonic processes and low value index imply circular basin indicate less active tectonic processes [17] [40].
Hypsometric analysis is the study of the distribution of horizontal cross-sectional area of the landmass with respect to elevation.HI is calculated by following formula mean elevation minimum elevation HI maximum elevation minimum elevation The Here, hypsometric curve is used to show the relationship between elevation of the basin and the area of the basin.Both hypsometric integral and hypsometric curve is important to show the conditions of the basin.

Asymmetry Factor (AF)
The asymmetry factor developed to detect tectonic tilting transverse to flow at a drainage basin.The equation of asymmetry factor is

Ar AF 100
At where, Ar = the area of the basin (facing downstream) of the trunk stream; At = the total area of the drainage basin.
AF index equal to 50% the drainage network formation is under stable condition.If AF > 50% indicate tilting left downstream and AF < 50% indicate titling right Downstream.Such tilting is result of active tectonics rather than lithological factor [41].The absolute difference (AF-50) is taken as consideration to evaluate the tectonic influence in the drainage basin.AF-50 value close to 0 indicates less active tectonic activity while value increase from 0 indicates higher tectonic activity [40].

Transverse Topographic Symmetry Factor (T-Vector)
Another quantitative index to evaluate basin symmetry is the transverse topographic symmetry following the basic technique presented by Cox [42].T-factor index is calculated with formula where, Da = the distance from the midline of the drainage basin to the midline of the active meander belt; Dd = the distance from the basin midline to the basin divide.
This above method is applied to the main stream of the basin.To calculate the T-vector altogether 102 points were taken considering at the confluence points of all the major tributaries of the main river.
The Analytical Hierarchy Process (AHP) was introduced by Saaty [43] and most widely spread and used theory for decision making.It is used to analyze complex problem with multi criteria.In the present study, Analytical Hierarchy Process (AHP) has been used for discriminating basin with active tectonic control.

Data Analysis and Results
The Kakoi river is 5 th order basin which traverse through high land, mid land and low land.It is observed that all the stream of the Kakoi watershed follow Horton's Law of stream number (Figure 3).According to Horton's law the number of streams in different orders in a given drainage basin tends to closely approximate an inverse geometric series in which the first term is unity.The drainage pattern is dendrite and the tributaries are oriented mostly in North-east direction.
It is observed that the stream length of the Kakoi River Watershed follows Horton's law of stream length.But an anomaly has observed in the stream length of the Kakoi river watershed.It is found that 1 st , 4 th and 5 th order stream shows some deviation from the straight line which according to Horton's law is due to structural control of the basin (Figure 4).
The morphometric indices that are taken as consideration to determine the        lithology plays important role.Hl and Vf has shown less influence of tectonic activity.Most of the basin area is low laying plain where lithology plays important role in the evolution of the basin.Erosion and deposition is the main activity which influence the landscape of the basin.

Discussion and Conclusions
The morphometric analysis of the Kakoi River basin has indicated that the basis has undergone dynamic changes from both lithology and active tectonic activity.Lastly we can conclude that the basin area is controlled by active tectonic as well as lithology.This present study has been carried out based on morphometric analysis (surface features anomalies).However, detailed subsurface lithological, structural, and geophysical studies are required to analyze the micro and mega scale changes of watershed, which are helpful for sustainable watershed development especially in tectonically active regions.

F
. Sangma, G. Balamurugan DOI: 10.4236/ijg.2017.8110811386 International Journal of Geosciences thrusting the Eurasian plate and is causing deformation and instability in the Brahmaputra basin.The flood plain of Brahmaputra, domesticated over the centuries suddenly face new challenge due to the earthquake of 1897 and 1950.
km in the hilly terrain of Arunachal Pradesh creating deep gorge and rest flows in the alluvial plains of Assam with shallow depth.The minimum height is about 21 m and the maximum is about 1728 m International Journal of Geosciences
results of HI indicate three stages of the basin: 1) Basin having HI> 0.5 representing with deep incision and slight erosion from active tectonics; 2) Basin having 0.4 < HI < 0.5 indicate approximate equilibrium and an intermediate stage of incision and erosion from recent active tectonic and; 3) Basins having HI < 0.4 characterized by low relief and severe erosion [37] [40].

Figure 3 .
Figure 3. Regression of logarithm of number of stream segments (Nu) versus stream order (u) for Kakoi River Watershed.

Figure 4 .
Figure 4. Regression of logarithm of mean stream length (Nu) versus stream order (u) for Kakoi River Watershed.

Figure 7 .
Figure 7. Migration of the river course.
The northern part of the basin has high relief with maximum 1728 m and part of lower Himalaya mountain range.These regions are continuously uplifted because of the Indian plate movement towards the Eurasia plate.As a result, thrust is created along the foot hills of the lower Himalaya range.The mountain front sinuosity (mfs) indicates the present of active tectonic in the north part of the basin.The valley floor (Vf) ratio also indicates that the presence of high tectonic activity at the northern part of the basin where V shape valley has formed.Highly dissected and structural hill is the physiographic characteristic of the northern part of the basin.Most of the area of the basin is low-lying alluvial plain with low relief in the southern part.The lithology plays an important role in modifying the landscape rather than active tectonics.Deposition and erosion are the main activities dur-F.Sangma, G. Balamurugan DOI: 10.4236/ijg.2017.8110811400 International Journal of Geosciences ing the monsoon season in this area.The deposition of the sediment along the river bed decreases the depth of the river.As a result, the overflow of water during the monsoon season is common in this area.The dipping of sedimentary layers during runoff may be the other reason for the tilting of the basin.The migration of the river course is more prominent in the alluvial plain of the basin than the hilly region in the north.The lithology plays an important role along with significant active tectonic activity in river course migration.Most part of the river course is migrated towards north-west and south-west direction.

Table 2 .
Stream order and stream length.

Table 4 .
T-index for Kakoi River Watershed.
dicates that the lithology and the tectonic plays important role while flowing through the alluvial plain.Elongation ratio (R e ) of the Kakoi watershed has indicated that the watershed has less influence from recent active tectonic rather F. Sangma, G. Balamurugan DOI: 10.4236/ijg.2017.8110811399 International Journal of Geosciences

Table 5 .
AHP comparison for different morphometric indices.

Table 6 .
AHP weight and rank for different geomorphic indices.