Morphometric Analysis of Tectonically Active Bhimtal-Naukuchiatal Lake Basin, Kumaun Himalaya

.


Introduction
Throughout the Quaternary Period, active tectonics along thrusts and faults has been reshaping the Himalayan topography, causing geomorphic changes.The evolution and change of varied landscapes in tectonically effective locations are also linked to tectonic elevate, deterioration, rapid erosion, and denudation mechanisms that combine tectonics and climate. 1,2,3,4,5,6Thrusts such as the HFT, MBT, MCT, and analogous secondary thrusts/faults such as RT, SAT, NAT, BT, and THF characterize the Kumaun Himalaya 7,8,9,10,11,12 and the tectonic activity along these has resulted in uneven slopes, asymmetrical topography and irregular landscapes. 13,14,15,16,17,1819,20,21In certain regions, rivers are more responsive to provide the geomorphic elements that make up the fluvial landscape architecture. 22,23,24,25,26Since neotectonic activity can have a substantial impact on fluvial systems, qualitative and quantitative investigations of drainage pattern are valuable for evaluating how it affects geomorphic processes and landscape evolution 22,27 Some of the notable characteristics of the tectonically formed landforms in vicinity of thrusts/faults zones of the Himalaya in general and Uttarakhand Himalaya in particular are the active and settled landslides, formation of palaeolakes, raised paired or separated terraces, incision, entrenched meandering, change of valley channels, occurrence of deep valley, waterfalls and displacement of previous faults by the later ones . 28,29,30,31,32,33,34,35,36,37In tectonically active regions, morphometric tools are widely used to understand how tectonics affects the landscape. 38,39,40,41,42,43,44,45 number of studies have shown that the mountainfront in the area is affected by neotectonic. 46 number of features such as landslides, asymmetrical river terraces, river incision, the formation of deep gorges, triangular facets, and tilting of recent deposits, suggest that the Gaula basin is structurally regulated and neotectonically active. 47,48,49Analyzing linear, areal, relief and morphometric parameter of basin play important role to understand the structural/lithologcal control of the basin. 49,50e youthfulness and spectacular exposure of Himalayan terrain makes it ideal orogeny to study variety of geological process. 51,52,53,54,55,56,57,58he Himalayan orogeny and the lithospheric deformation have aggravate dintense research in recent years. 59,60,61,62The landscape modification in the Himalayan thrust zones is result of tectonics and surface geological processes viz, erosion, weathering, upliftment and frangible lithology.The traces of these modification can be evaluated with the help of longitudinal profile modeling of river. 63Degree of tectonic activities attribute to a range of geomorphic and structural features which leads to landscape evolution.Over the years evaluation of tectonic activity with the help of geomorphic and structural expressions has attain its own scientific platform to infer tectonic activities. 64,65Apart from this toposheets and satellite data along with multidisciplinary approach act as add on for morphometric indices analysis for a comprehensive evaluation of tectonics. 66,67,68 combination of satellite based tectonics, advance geospatial technology and field observations have been used largely for better understanding of landscape evolution by mapping, evaluating and analyzing the morphometric parameters. 48These studies aid in identifying active tectonic deformations, regions that are vulnerable to natural disasters, and land use planning. 69In response to tectonics and differential uplift over a range of timescales from thousands to millions of years morphometric analysis detects silent changes in a river.
The geomorphic indices identify endogenic geomorphic processes when they identify basin anomalies brought by tectonic activity. 50e Vf assumes that every area susceptible to up-lift has carved ravine with tapered floors and deeper levels of V-shaped ravine. 70Greater Vf values are typically linked with sluggish tectonic activity and wide, flat ravine floors, on the other hand less Vf values are reflective of comparatively quick uplift in V-shaped ravine and tributaries that have recently been incised and are connected with grater uplift rates. 70It also shows the developmental stage of a river valley and is an excellent indicator of the extent of incision. 71Furthermore, the appearance of the valley wall is determined by the lithology of the foundation rock and the river's erosive capacity. 72he hypsometric profile depicts variation in altitudes throughout a geographical area, from one watershed to the next.This curve is ratio of total basin elevation and total basin area.The basin's overall area (A) is the addition of the areas between each adjoining contour.In a basin surface area (a) is area above a certain line of altitude (h). 70,73The relative area (a/A) fluctuates from 1.0 to 0.0 and value 1.0 indicate minimum level of basin and 0.0 indicate highest level of the basin.The SL gradient is important aspect for identifying irregularities in natural concaveness within longitudinally.The steepness index allows for the normalization of slope measurements as well as for the recognition of aberrant locations in each part of the river from the upstream to the terminal.The SL index is extremely reactive to variations in river gradient, allowing for the estimation of correlations between probable tectonic forces, rock resistance, and topography.Along a graded stream, the SL index stays roughly stable, and fluctuations in this index tend to be due to tectonic and lithological restrictions. 74anching design of a channel network is strongly connected to the bifurcation ratio.The high R b shows an early hydrograph maximum along the possibility of sudden flooding during stormy incidents in places where these stream orders occur. 50Drainage density is a helpful mathematical expression of topography segmentation and washout capacity since it influences the time travel by water.A more drain density suggests a well-drained watershed with a quick hydrological reaction to precipitation, whereas less Dd suggests a weakly drained watershed with a moderate hydrological reaction. 50mith (1950) categorized the drainage density into texture divisions viz, Grater than 2 very course, 2-4 course, 4-6 moderate 6-8 fine and greater than 8 very fine. 75Less Dd indicate areas with resistant rocks with extremely permeable subsurface and when elevation ratio is minimal whereas more Dd shows strong or impermeable underlying material, scarce vegetation, and steep topography. 50Stream order, or classifying streams according to the quantity and kind of tributary junctions, is a simple, practical, and all-encompassing measure of stream size, discharge, and drainage area.Horton (1945) developed a system for ordering streams that was later refined by Strahler (1964). 76,77One of the most popular proxies, the SPIM based model evaluates the variability in rock upliftment in different catchments and is related to the steepness index (Ksn). 78,79,80In steady state conditions, the variability in uplift of longitudinal river channel can be conventionally evaluated using this model.The direction of movement of the drainage division with regard to the steady state line, where the drainage is either losing (victim) or gaining (aggressor), in inequilibrium state, is usefully understood using the Chi (χ) parameter. 81,82One can forecast the migration of the drainage divide when all the tributaries display varied values of Chi and depart from the ideal steady state since the migration always moves from low to high chi values.Asymmetric Factor is often used to estimate basin's maturity and the stage of equilibrium of the basin. 3When a river flows in a homogeneous terrain, which have not been affected by tectonic activity it anticipate the symmetric drainage network and dendritic drainage network. 70On the other hand asymmetric drainage network indicate upliftment of rock which tends to tilt the basin lithology and can be seen in the field as the uplifted part shows more length of tributaries as compared to other half. 83hen the AF value is less than or more than, 50 it suggest that the basin is tilted and developed under an active tectonic phase.When AF is more or less 50 indicate tilting towards left and right respectively. 70The longitudinal profile (Figure 5) is a curve derived from association between a channel's height (H) and distance downstream (L).It depends on numerous determining factors, viz., initial relief, rock type, palaeo-geographical evolution stage and positive or negative tectonic movements. 70The form of the profile, especially its degree of concavity is the result of the several factors, upwardly concave profiles could indicate protracted basin and channel deterioration due to longer time since base of ravine sinking.Upwardly convex outline indicates smaller down cutting, and a shorter time since ravine base decline. 70,75Since the concept of drainage basin evolution is not utterly understood yet in Kumaun Himalayan area, here present exercise is an attempt to evaluate variety of morphometric indices to assess relative tectonic activities for better understanding of tectono-geomorphic evolution.The objective of the present exercise is to calculate morphometric indices in the Bhimtal-Naukuchiatal lake basin with respect to the neotectonics activity along the thrust/ fault zones.

Study Area
Our study area extends (79°30'E to 79°36'E:29°17'N to 29°24'N) with Survey of India toposheet no.53O/11 (Figure 1).We studied a number of subsidiary thrusts and tectonic zones along which the mass movements, fault scarps and smaller palaeolakes have been developed, which shows that the area is affected by tectonics.The Bhimtal-Naukuchiatal area which comes under Kumaun Lesser Himalayan region shows rough topography which includes large valleys, low to high hills, escarpments, gorges, sag ponds, streams and rivers.The large valley formed by the NW-SE extended median part has a number of lakes and intervening terrace-like flats. 33The hills with quartzite are flanked by volcanic rocks on either side of the series of lakes.The denuding agents may have exposed the volcanic rocks that make up the core of the anticline by forming a broad erosional valley that slopes towards Naukuchiatal.The sharp, elongated ridges that run NW-SE throughout the area slope gently towards the NE but climb more steeply in the SW slopes.The high hills with summits 1,863m, 1,927m, and 1,786m are located to the NE of the chain of lakes.These ridges dip NE to create a river basin that runs from NW to SE, beyond which the terrain again rises.30TheBhowali town is located in the northwest corner at elevation of 1,715 m.The terrain to the southwest of the Bhowali-Bhimtal border is extremely rough, with winding ridges that fall sharply to the SW.The peaks 1,779m, 1,773m and 1,605m mark the highest points.These lofty peaks descend down into the Gaula river valley in the south.

Geological and Geomorphological set-up
Located at 1,356 m.s.l., the Bhimtal lake is 1.8 km in length, 0.454 km in width with a catchment area of 4.5sqkm.The shape is roughly triangular, the western limb being the largest and curved inward. 31ccording to Thomas (1952), the lake was formed EBRAHIMI & ASL, Curr.World Environ., Vol.18(2) 855-875 (2023) by choking of river through landslips. 46However, a detailed study around Bhimtal lake suggested that it was formed due to neotectonic activity and the palaeolake extended for about 11km. 30he stratigraphy around the region is displayed in Table 1.Lithology (Figure 2) around the lake mainly belongs to Bhimtal volcanics, consisting mainly of metabasites in association with meta-sedimentary formations, the later predominantly consist of quartzites with subordinate grits and phyllites and are associated with prehnite-pumpellyite metagreywacke facies and zeolite facies, grading into greenschist facies. 33The main geomorphological features are, debris flow type fans, old and new landslide,depositsalluvial fans, valley scarps, triangular fault facets, talus cones, fluvial, colluvial and alluvial fans, palaeolake deposits along the Nagarigaon-Naukuchiatal fault (Figure 3) and waterfalls, entrenched meandering, shear zones as well as a deep gorge through which the palaeolake was breached. 30In the current exercise, various morphometric characteristics are examined, including linear parameters (stream order, bifurcation ratio), aerial or basin parameters (drainage density), relief parameters (relief), and morphometric parameters (vally floor width-height ratio, HI, stream length gradient index, steepness index, Chi, and asymmetry factor). 50The hierarchical link between the many stream segments that make up a drainage network is expressed by stream order (Sµ) which is initial stage for the whole excercise. 50The number of streams in particular order (Nu) divided by the number in the next lower order is expressed by the dimensionless parameter known as the bifurcation ratio (Rb). 48,50,70Areal parameters determine how much erosion is occurring as a result of exogenic geomorphic processes. 70The total stream length in a basin divided by basin's area is known as the drainage density.It specifies the relative position of drainage channels and measures the average length of streams per unit drainage area. 47,50The relief parameters include basin relief (H), where H is elevation difference and is defined as the roughness of the topography. 48Based morphological condition of the basins, the morphometric parameters describe the processes that define the landscape and active tectonics. 70The ratio of the width and height of the valley is a value of the valley area that characterizes the uplift of an area. 50,76HI describes the relative area distribution at various elevations.Higher HI values of the basin suggest that majority of the land surface is elevated compared to the mean, viz., plane highland terrain dissected by extremely sculpted channel. 70edium to less integral values are connected with uniformly distributed catchment area. 50By examining potential tectonic tilting of the basins, it is feasible to determine the tectonic influence in a watershed. 48f depicts the flow-transverse direction tilting. 9he Stream Length-Gradient (SL) index efficiently delineates the topographic break zones by defining the change in stream slope along a longitudinal river profile. 84,85The use of chimaps to show the migration of divides has grown in popularity.The variable chi at the divides indicates whether a drainage basin is a victim or aggressor, or whether it will lose or gain its drainage region in the future.chi (X) is a horizontal transformed coordinate taken from a drainage basin's outlet. 84,85Ksn can be used to evaluate the longitudinal river channel's variability in uplift under steady state conditions.among the most widely used proxies, the steepness index (Ksn)-related SPIM driven model assesses the variation in rock upliftment in various catchments78.We selected a reference value of 0.45 to prevent anomaly changes. 82Distance upstream serves as the X axis, while Ksn (on the dexter side) and elevation (on the sinister side) serve as Y axis indicators.As Ks is normalized to a reference value, different values of for the longitudinal profile result in varied Ks values.The concavity is caused by the stream power incision law (SPL), which states that along-river slopes S are related to upward area a multiplied by the negative mn-ratio, and continuous uplift hypotheses.

S ~ A -m/n
According to the SPL, the amount of energy needed to incise into the stream bed depends on the stream gradient (S), the area of the upslope (A), as well as several factors (m and n).

Ratio of Valley Floor Width-Height (V f )
In present study, we took V f values (Figure 4) from five different locations with values 0.49, 0.28, 0.19, 0.37 and 0.40.Such low values specify that the area has been controlled by recent ground movements and upliftment.Because the valley is completely encircled by peaks, the valley walls are practically vertical, indicating that the area is significantly more prone to erosion.

Longitudinal Profiles
Study area displays an upwardly concave contour, indicating extended basin and stream deterioration linked to longer periods of time.It also suggests the less channel down cutting.Figure 5 is longitudinal curve shown by river from Nagrigaon (upstream) to Amritpur (downstream).

Hypsometric curve and Hypsometric Integral (HI)
The basin contain the Nagrigaon-Naukuchiatal fault, the Bilaspur fault, and the Sattal fault have medium HI values (0.59), and sigmoidal in shape hypsometric curves, indicating that the development process is at a mature level (Figure 6).These regions are least influenced by tectonic activity and are primarily eroded.

Stream-Length (SL) Gradient Index
A total 16 km long Bhimtal stream was taken for analysis of the SL curve (Figure 7) which demonstrates that the greatest values correspond to the knickpoints regression.The abnormalities are related to the structure that divided the river portions from higher to lower stream.In our area, SL values fluctuate from 0 to 1384, which indicate possible tectonic activity.Along the longitudinal profile, major knickpoints are found at 6.5km and 12.7km from the source of the stream, which indicate occurrence of thrust and faults in vicinity.3).originates along the slopes of Tani peak as Tandi Gad, meets Kalsa river (that flows from E to W) and finally flows due SE meeting Malwatal area.The Kalsa gad and Champi gad show phenomenon of the incised meandering.All these streams except Gaula, contain little water during the summer months.

Relief
The difference in elevation of the earth's surface or comparative vertical disparity of the land surface is termed as relief.Relative relief approaches can efficiently depict three-dimensional maps without taking into account sea level, as well as variances in altitude. 50It is generally, the "ruggedness" of the topography, particularly, the highest elevation minus the lowest elevation in a given area.The study area occupies highly variable relief and the relief values vary from minimum 509 to maximum 1,839 m and are classified into five groups (Figure 9).

Steepness index (Ksn)
The variable Ksn values shown by major streams are shown in Figure 10a.The southern part of the basin showing high values of K sn as compared to northern part which indicate the southern block of Nagrigaon Nakuchiatal fault is uplifted and northern block is down faulted block.The values of ksn is around 550 shown in Figure 10b with thick red horizontal line represents the average K sn value, and the green and blue curving areas represent positive and negative Ksn anomalies, respectively (Figure 10b).The values of ksn are around 550 which suggest that the basin is experiencing rapid tectonic upliftment.

Chi Integral (χ)
The Chi (χ) plot for different values of mn is shown in Figure 11. Figure 12 shows the lateral distribution of Chi along the longitudinal course of the main trunk and tributary stream.From the Chi plot of Bhimtal-Naukuchiatal, it is clear that towards the upper reaches, the erosion activity (elevation range 1300 to 1500) is very high.Towards the middle riches of the river, the Chi values indicate lesser erosion activity and, in this portion, the rivers are actively incising due to tectonic uplift.Further, towards the lower reaches, the gradient change is caused by tectonically enhanced upliftment processes.

Stream Ordering (Sµ) and Drainage Pattern
In the study area, we observed five orders of stream (Figure 3) in the drainage network.Sinuosity is a characteristic of a river to travel from one place to another across its flood plain.The stream may leave scars of where the river channel once was as it meanders through the flood plain.A river that does not meander at all has sinuosity of 1. 70 The more meander in a river, closer the sinuosity value will get to 0. Due to the ruggedness of the terrain, no clear picture of the drainage pattern is observed but, in general, it is noticed that various small streams flowing down the slopes of the high ridges, join one or the other streams, which ultimately give rise to a few south-easterly bigger streams which are more braided in nature.There are three big streams of which the most prominent one is the Gaula river in the extreme south, just outside the investigation area.

Discussion
The drainage system in the area follows a dendritic drainage pattern in whole Bhimtal Naukuchiatal basin.Lithology, structure, and topography all influence the drainage pattern.By observing the allotment of relative relief it is concluded that the basin consists mainly of mountainous terrain.
The results show that most of the river streams are controlled by the intense pattern of NE-SW trending faults that can also be seen in the longitudinal profile.The faults in and around the study area are responsible for recent mass movement and also make the more vulnerable for landslide disaster.Construction should be banned in vicinity of fragile rocks, thrust/fault zone as there is a limit to how much weight these brittle, worn, steep angle slopes can support.Plantation should be done in and around the active landslide zones and awareness programs related to landslide hazard should be organized among the local population for better understanding of their habitat environment.

Fig. 3 :
Fig. 3: Geomorphological map around the study area Methodology In this exercise, SOI toposheet No. 53O/11, ASTER-GDEM were processed in ArcGIS 10.4 version to calculate morphometric parameters, viz., V f , HI, SL, Rb, Dd, Sµ and average Sinuosity (S), Relief (H).Drainage basins are defined using the Hydrology tools in GIS after additional processing of the DEM data.The extraction procedure incorporates void filling, flow direction, flow accumulation quantification, and drainage basin extract.Using programme MATLAB, we computed the 30 m SRTM DEM and Topo Toolbox (WGS84) to analyze the dynamics of drainage pattern with Bhimtal-Naukuchiatal basin. 84,85Topo Toolbox allows for the identification of ksn, and Chi (χ) to understand the dynamics of drainage basins.

Fig. 4 :
Fig. 4: Overall drainage in the study area (different stream order are shown by different color code)

Fig. 8 :
Fig. 8: SL index in the study area

Table 4 : D d and average sinuosity Drainage density Average sinuosity (km/km 2 )
It flows due west and receives water from the various south-west flowing smaller streams that drain the slopes of Barchuna Dhar.Another striking stream, Nauti gad which originates near Bhowali and meets the Bhimtal lake.It flows out of the lake due south and meets the river Gaula.
This stream is known as Bhimtal Gadhera.Another important stream is Champi gad in the NE that EBRAHIMI & ASL, Curr.World Environ., Vol.18(2) 855-875 (2023) The structure of the valley floor is examined in relation to the valley's tectonic background.By calculating V f value for different location of the study area, low values of V f indicating high upliftment rate and more incision.Active tectonic modification can alter river system geometry which affects aggradational and degradational processes.This distortion is more obvious as streams and terraces are displaced across faults.Upwardly concave profile is obtained in the study area, which indicates the prolonged basin and the channel degradation associated with longer periods.As a result, morphometric variables and field geomorphic traits demonstrate that the region is structurally active.Differential upliftment, tilting in the Bhimtal Naukuchiatal basin all along and across the Nagrigaon Naukuchiatal fault, sattal fault and salari thrust is observed which is shown in the Figure2and the degree of tilting which is represented by AF values.Similarly we estimated channel sinuosity for all the tributary streams of main channel and observed that the sinuosity of trunk and tributary streams are varies from 1.06 to 1.29which indicate most of the streams are braided in nature.The observed values of sinuosity indicate youth stage of channel where channel migration and subsequent incision is governed by active deformation.The steepness values and chi graph also indicate the Bhimtal Naukuchiatal basin is experiencing differential upliftment.Vf, HI, SL, Bifurcation ratio (Rb), Drainage density (Dd), Stream Order (Sµ) and average Sinuosity (S), Relief (H), steepness (K sn ), Chi integral (χ) And Asymmetric Factor (AF).The low value of V f indicates high upliftment rate, more incision and less channel down cutting.The research area has intermediate HI values as well as a sigmoidal-shaped hypsometric curve, both of which signify a mature stage of river development.The SL values range from 0 to 1384, which suggests potential tectonic activity.The average bifurcation ratio as 4.46 indicates potential structural control.The Ksn values are approximately 550, suggesting that the basin is undergoing rapid tectonic upliftment.From the Chi plot, it is evident that erosion activity is particularly active in the higher reaches of channel (elevation range: 1300 to 1500), whereas the intermediate and down portions of the river show less erosion activity and active tectonic uplift.The calculated value of AF for Bhimtal Naukuchiatal is 40 which indicate the tilting of the basin is towards right.