SEDIMENT DYNAMICS AND STABILITY STATUS OF THE KARRA KHOLA, HETAUDA DUN VALLEY, CENTRAL NEPAL SUB-HIMALAYA

The Karra Khola in Hetauda, Siwalik region, originates within the Dun Valley, and contributes the Rapati Nadi in Hetauda Metropolitan city. The stability status of the river is of main concern because of rapid growing of the river corridor and peripheral land areas into suburban city. The river was surveyed for hydraulic parameters, sediment characteristics, and sediment loads. Rate of sediment transport and sediment yields were computed, and competency of the river was evaluated using Shield’s parameters and Reynolds numbers. The results show that the river sediments are sandy gravel to gravelly sands, and are moderately to very poorly sorted. The total sediment yield of the whole basin near the outlet is around 2% of the maximum total sediment yield. The


INTRODUCTION
The Siwalik Group is comprised of Late Miocene to Early Pleistocene sedimentary sequence deposited in the foreland basin of the uplifting Himalaya (Tokuoka et al. 1986, Kizaki 1994, Critelli and Ingersoll 1994, DeCelles et al. 1998, Gautam and Rosler 1999). Galay (1987) explained that the steep slopes of the Siwalik Hills and weakly consolidated rocks with thin soil cover generally tends to promote severe surface erosion in spite of the thick vegetation. Likewise, Shrestha et al. (2005) also suggested that the apparent reason for the huge sedimentation within the Siwalik basins could be because of fluvial characteristics. Tamrakar and Karki (2019) concluded that there were variation in specific sediment yields and sediment delivery ratio among the sub-basins of the Karra Khola, and reported that specific sediment yields were lower in most of the southern sub-basins compared to the northern subbasins of the Karra Khola, due to topography and streambed characteristics.
Erosion, transportation and deposition of sediments depend upon hydraulic action, which is again governed by topography and climate (Grantham and Velbel 1988), tectonics and weathering (Johnson 1990), composition of source area (Dickinson et al. 1983, Dickinson et al. 1986) and riverbed materials (Cavazza et al. 1993, Cox and Lowe 1995). Deposition occurs when the critical shear stress required to entrain particles of a specific size exceeds boundary shear stress, if the case is opposite, then the sediment gets entrained and transported. Several researchers have formulated the sediment transport using the concept of bed shear stress or unit flow discharge based on field and laboratory experiments (Du boys 1879, Schoklitsch 1934, Shields 1936, Meyer-Peter and Muller 1948, Einstein 1950, Bagnold 1966, Yang 1972, Van Rijn 1984. Hassanzadeh (2007) presented a dimensionless semi-empirical formula on the bed load, and showed that the values obtained from the formula were closer to the measured data than the common ones on sediment hydraulics.
The behavior of river can be understood by its sediment transport capacity and flow competence. The sediment transport capacity is the river's capacity to transport the total amount of its load across the specific site of the river bed (Bagnold 1960(Bagnold , 1966. It can be expressed in terms of sediment total load carried or in terms of river's energy required to transport sediment per unit bed surface area, which is expressed as the unit stream power or specific stream power (Bagnold 1966). Several researchers expressed usefulness of unit stream power in interpreting risk of erosion of the stream, and channel sensitivity to erosion and deposition (Gilbert 1914, Rubey 1933, and Velikanov 1955, Brookes 1987, Miller 1990, Magilligan 1992, Rathburn 1993 Stream flow competence is stream's ability to entrain and move the available bed materials of desired size grades that are often median grain size (Wilcock 1971). Generally, incipience happens when the boundary shear stress of the stream overcomes the critical shear stress of the bed material (Shields 1936, Knighton 1998 . Commonly, rivers with high capacity but low competence tend to deposit sediments, those with high capacity, high competence tend to transport enormous amount of sediment, and those with low capacity but with high competence tend to degrade river to achieve a dynamic equilibrium. The present study aims to establish sediment flow dynamics, dealing first with the sediment transport capacity and its bearing on potential risk of erosion, and secondly with the competency of the Karra Khola that is an axial river originated and flowing in the Dun Valley of the Siwaliks, central Nepal (Fig. 1). The status of capacity and competency of the Karra Khola leads to evaluation of stability of the river.

GEOLOGICAL SETTING OF THE STUDY AREA
The Karra Khola Basin (KKB) spans for 6-9 km in width and 17-18 km in length, and comprises of the Siwalik Hills, which are composed of middle Miocene to early Pleistocene molasse sediments of the Siwalik Group (West andMunthe 1981, Tokuoka et al. 1986). The two belts of the Siwalik Group are separated by the Main Dun Thrust (MDT) along which the Lower Siwalik Subgroup of the northern belt overthrusts the Upper Siwalik Subgroup of the southern belt (Fig. 2). The northern belt Siwaliks folded as the Ratomate Syncline and the Karra Khola Anticline (Schelling et. al. 1991, Kimura 1994, Kimura 1995, Tamrakar and Karki 2019, whereas the southern belt Siwaliks are gently north dipping homocline. The KKB covers portions of both the northern and the southern belts of the Siwalik Group. It was created due to the thrusting along the MDT, and then progressive upliftment of the Chure Range as a consequence of the thrusting along the Main Frontal Thrust, located further south of the basin. The southern part of the KKB consists of thick conglomerates with few intercalations of mudrock and sandstone of the Upper Siwalik Subgroup that forms the major part of the Churiya Hills and divides the southern boundary of the KKB. The northern portion of the KKB consists of intercalation of mudrock and sandstone of the Lower Siwalik Subgroup followed by thick salt and pepper sandstones with alternation of silt-shale, mudstone, and clay-shale of the Middle Siwalik Subgroup. Both the Lower and the Middle Siwalik Subgroups cropping out in the northern portion of the KKB show NW-SE strike and steep NE dipping (Fig. 2).
The unconsolidated Quaternary Alluvial sediments, which are the basin-fill sediments, have been deposited since the Mid-Pleistocene to Recent (Kimura 1994), and unconformably overlie the Siwalik rocks of both belts, mostly covering the wide extension of the central portion of the KKB. These Quaternary Alluvial sediments consist of coarse-to fine-grained sediments originated within the KKB, and can be distinguished into three morpholithological zones; Upper, Middle and Lower (Tamrakar and Karki 2019). The Upper zone comprises of bouldercobble to pebble-cobble gravel, the Middle zone of cobble-pebble gravel and sand, and the Lower zone of pebbly sand, sand and mud.
The emerging point of the main stem Karra Khola is located in the Siwalik Hills, at the upstream head of the Tamaune Khola, somewhere at Makawanpur Gadhi (Fig. 2). Upto the confluence of the Tamaune Khola and the Barauli Khola, the gradient of the main stem Karra Khola is quite steep, and further downstream the gradient lowers abruptly, and then diminishes gradually (Fig. 3). The main stem Karra Khola is the 5th order river and its major tributaries are the Gahate Khola, Trishuli Khola, Chisapani Khahare, Gauritar Khahare, Agarathe Khahare, and Khahar Khahare, which originate from the southern part of the basin (Fig. 4), and all are intermittent streams that contribute towards the left banks of the Karra River. The Bhedaha and Gardoi Kholas are remarkable tributaries from the northern portion and contribute to the right banks of the Karra Khola. Therefore, the main stem Karra Khola appears as an axial drainage in the basin.

Methodology
The main stem Karra Khola was surveyed for recording the nature of sediment distribution along the river channel and banks. A sediment distribution map was prepared in 1:10,000 scale. Altogether 19 main stem transects, which were continuously from the transect number T1 to T19, were selected for the cross sectional and longitudinal surveys. Cross-sections of all transects were prepared to obtain hydraulic parameters of the river. Along each of transects, channel sediments were assessed applying the Wolman pebble counting method (Wolman 1954, Rosgen 1996. The grain size distribution of the channel sediments were processed through graphic statistical measures after Folk and Ward (1957).
The Karra Khola Basin lacks gauging station, therefore, velocity (V) and discharge (Q) of the Karra Khola were estimated using the Manning's Formula and the continuity equation (Chow 1959): Where, V= velocity (m/s), R=hydraulic radius (m), a ratio of cross-sectional area (A) to wetted perimeter (P), S=average channel slope (mm-1), and n=Manning's roughness coefficient.
The Manning's roughness coefficient, n at each transect was estimated by firstly calculating the basin roughness value using the equation of Limerinos (1970) and then inserting this value along with estimated adjustment values in the equation of Cowan (1956) as below: Where, R = hydraulic radius (m), d84 = intermediate particle diameter (in meter), that equals or exceeds that of 84 percent of the particles, n0 = base value, n1 = additive value to account for the effect of cross-section irregularity, n2 = additive value to account for the variations in size and shape of the channel, n3 = additive value to account for the effect of obstructions, n4 = additive value to account for the type and density of vegetation, and m = adjustment factor for the degree of channel meandering that is the ratio of channel meander length to the straight channel length.
The water samples collected from each of the transects more or less at the bankfull stage were filtered in the laboratory using filter of uniform pores (0.45 µm) to calculate the suspended sediment load carried by the stream. The suspended sediment concentration (SSC) was calculated as the ratio of the dry mass of suspended sediment retained to the volume of water sample taken to filtration.
The stream power (Ω) is the rate of potential energy expenditure against beds and banks of channel per unit length, and can be expressed using the following expression (Bagnold 1960(Bagnold , 1966: Where, Ω = stream power (Ns -1 ), ρ = density of water (kgm -3 ), g = acceleration due to gravity, (ms -2 ), Q = hydraulic discharge of the stream (m 3 s -1 ) and S = stream slope (m/m).
The unit or specific stream power (ω) is defined as the total stream power (Ω) per unit width of the channel (W) (Bagnold 1966), and is expressed in W/m 2 , where W stands for Watt.
The competence of the Karra Khola was evaluated by using shear stress criterion of Shields, and in this the dimensionless boundary shear stress (*τ b ) was calculated applying Shields criterion as below: Where, τ b =boundary shear stress (Nm -2 ), γ = specific weight of water (Nm -2 ), R h = Hydraulic radius (m), S= Average stream slope (m m -1 ), γ s = specific weight of sediment = 26500 Nm -3 , and d 50 = median grain size. Then the dimensionless critical shear stress (τ c *) was obtained after computation of Reynolds number and extrapolating it using the Shields curve. The boundary Reynolds number was obtained using the following expression: Where, u* = shear velocity (ms -1 ) = (√(τ b /ρ)), τ b = boundary shear stress (Nm -2 ), ρ = density of water (1000 kg m -3 ), d 50 = median grain size (m), and ν = hydraulic viscosity (1.004 *10 -6 m 2 s -1 ). The stream transect site was taken as competent when the dimensionless boundary shear stress exceeded the dimensionless critical shear stress, otherwise it was regarded as incompetent. Then, both transport capacity and flow competence were considered with their thresholds to evaluate stability status of the river along its course.

CHANNEL BANK DEPOSIT
As the stream descends from the higher relief areas to the flat plain downstream, the channel bank deposits are extensively developed (Fig. 5). The sediments in such depositional elements are laminated fine sand, sandy mud, muddy sand and gravelly mud. The sediments are relatively coarser in the banks near to the source regions compared to the banks of the downstream regions.

CHANNEL SEDIMENTS
There are many convex bends in the Karra Khola corridor where the point bar depositional features are frequently developed along the Karra Khola. The point bars of both sandy gravel and gravelly sands are observed (Fig. 5). The point bar deposits show thick layers of fine-to medium-grained pebbly sand fining up section into finegrained sand.
The recent sediments along the Karra Khola is composed mainly of gravelly sand (gS), gravelly mud (gM), muddy sand (mS), sandy gravel (sG), and sandy mud (sM). The boulder sized sediments are present at the upstream channel, where coarse sediments are shed from mass wasting from the valley slopes. The sediment size decreases moving downstream and becomes sandy gravel to gravelly sand on both banks of the stream channel, but clast size fluctuates due to mixing of gravels from the intermitted contributing streams from lateral direction towards the confluence with the Karra Khola.
Grain size characteristics of the channel sediments of the Karra Khola were assessed ( Table 1). The particles range from clay to boulder in most of transects of the Karra Khola. The mean size (Mz) varies between 57.947 and 0.293 mm. The median size ranges from 2 to 73 mm. Sorting () of the Karra Khola sediments ranges from moderately well sorted to very poorly sorted. Sediments of the upstream stretch are more or less coarsely skewed, whereas those of mid to downstream stretches show nearly symmetrical to finely skewed, however, majority of transects (T3, T4, T9, T10, T11, T12, and T18) indicate finely skewed sediments. The peakedness or graphic kurtosis (KG) varies from very platykurtic (at T12 and T13) to leptokurtic (at T2 and T6), and the majority of transects are being characterized by the mesokurtic grain size distribution showing that the contrast to the median size with the background size remains moderate. Along the river course, the median size tends to increase except for some sites T3, T4 and T12 where the size diminishes or the size gets diluted (Table 1). (mk) *Inclusive Graphic Standard Deviation: vps-very poorly sorted, ps-poorly sorted, ms-moderately sorted, mwsmoderately well sorted **Inclusive Graphic Skewness: ns-near symmetrical, cs-coarse skewed, fs-fine skewed, sfs-strongly fine skewed, scs-strongly coarse skewed ***Graphic Kurtosis: lk-leptokurtic, pk-platykurtic, mk-mesokurtic, vlk-very leptokurtic, vpk-very platykurtic The channel sediments possess moderate to poor sorting in the upstream stretch, whereas from the midstream towards the downstream stretch, the sorting becomes notably poor to very poor (Table 1). Both skewness and peakedness fluctuate along the downstream course of the Karra Khola suggesting that the sediments vary within ranges of fine to coarse skewness, and moderately poorly to well sorting. Sediment Dynamics and Stability Status of The Karra Khola, Hetauda Dun Valley, Central Nepal Sub-Himalaya

HYDRAULIC PARAMETERS
The morpho-hydraulic parameters ( Table 2) were calculated during field survey. The channel cross-sectional area (A) at bankfull ranges from 1.12 m 2 (T1) to 18.02 m 2 (T19). The hydraulic radius (Rh) varies between between 0.27 m and 0.85 m. The bankfull cross-sectional area in the midstream stretch grows roughly by ten times and that in the downstream stretch by 20 times compared to the bankfull cross-sectional area at the upstream stretch. Manning's roughness coefficients vary from .03 (T3) to .099 (T10) for different transects. The upstream stretch (between T1 and T5) possess slightly entrenched channels with good riparian vegetation, and heterogeneous, coarse skewed sediments (e.g., Figs. 6a and b). The midstream stretch between T5 and T16 remains slightly entrentched to non-entrenched with wide channel and sparse riparian vegetation (e.g., Figs. 6C and d). The downstream stretch (downstream from T5) lacks good riparian vegetation, and channels are wide non-entrentched with abundant sediments deposited in the channel (e.g., Figs. 6e and f). The velocity at bankfull ranges from 0.50 ms -1 (at T3) to 2.18 ms -1 (at T9), and exhibits fluctuating trend along downstream course. Discharge at bankfull varies between 0.81 m 3 s -1 (at T3) and 22.65 m 3 s -1 (at T11) ( Table 2), and tends to increase downstream (Fig. 7a), with notable higher discharge in the midstream course of the Karra Khola. Discharge of the Karra Khola tend to increase, as the basin area increases (Fig. 7b).

SEDIMENT TRANSPORT CAPACITY
The maximum and minimum bed load transport (Table 3), are respectively, 1,063,206 t yr -1 (at T13) and 3,515 t yr -1 (at T4). The suspended sediment load varies from 2,095 t yr -1 (at T3) to 109,550 t yr -1 (at T9) ( Table 3). Except for transects T2 to T5, where the suspended sediment load is less than 10,000 t yr -1 , and transect T9 where it is the highest, the suspended load tends to vary roughly between 10,000 t yr -1 and 60,000 t yr -1 . The total sediment load fluctuates along the main stem course and varies from 8073 t yr -1 to 1077220 t yr -1 . Transects 4 and 5 reveal lower total sediment load compared to the other transects. The sediment yield of the river ranges from 561 to 148,015 t km -2 yr -1 . The sediment yield from the outlet of the Karra Khola is 3199 t km -2 yr -1 . The sediment yield abruptly decreases from T1 towards downstream, showing that the sediment yield tends to be high at source region, where the gradient of the river is high, and it tends to diminish along the river course ( Table 3). The rate of sediment yield is not uniform, and varies along the Karra Khola, showing that deposition of sediments has taken place intermittently along its course.
The stream power (Ω) is the rate of potential energy expenditure against beds and banks of channel per unit length. The quantity of sediment flowing together with the water is highly dependent on energy of flowing water that is a stream power or also known as transport capacity (Bagnold 1966). The total stream power in maximum is 22.17 KNs -1 at T9 and in minimum is 0.11 KNs -1 at T3 (Table 4). The graphical representation (Fig. 8a) of the above data helps to visualize the downstream variation of the stream power of the river. The trend of variation is quite fluctuating between 5 and 25 KNs -1 , and the few lowest stream power in the upstream segments. The midstream segments have fluctuating trends due to variation in discharge. The specific stream powers provides sense of risk of erosion (Bizzi and Lerner (2015), and thus can be used to interpret the potential of stream erosion. Highly fluctuating trends of specific stream power has been obtained (Fig. 8b). Except for few locations, (T3 and T4; Fig. 8b), transects T2 and T5 represent stream power varying between 60 and 300 Wm -2 , therefore showing moderate risk of erosion, whereas the majority of the locations bear specific stream power exceeding 300 Wm -2 , and indicate high risk of erosion and transportation of sediments.

COMPETENCE OF THE KARRA KHOLA
Dimensionless boundary shear stress and dimensionless critical shear stress from different transects of the Karra Khola were calculated to analyse flow competence of the river. Both dimensionless shear stresses were compared to define flow competence. The Karra Khola is competent enough to transport the median sized sediments in most of the river transects, but it is still not competent enough to transport the median sized sediments at its lowermost sections T16, T17, T18, and T19 (Table 5). When competence of river was correlated with many variables, only the ratio of boundary shear stress to median grain size is found to give a good positive correlation. When this ratio, τ b /d 50 increases the competence of river tends to increase (Fig. 9). Therefore, at the river segments downstream from T16 transects ( Table 5), most of median sized and coarser sediments are deposited on the river channels and bars due to great reduction in boundary shear stress and increased median grain size, and are supposed to be carried downstream only during flash floods.

DISCUSSIONS
The boulder-sized sediments which are resulted by mass wasting are found in the source region of the Karra Khola. From sampled sites towards downstream course of the river, grain size varies from clay to pebble, and there is a weak tendency of increase in median grain size with the distance downstream. Sediments are moderately well sorted to very poorly sorted, and often show textural inversion. The textural inversion is attributed to the mixing from the lateral tributaries of the Karra Khola. Slope and discharge of the river also play important roles for the transportation and deposition of sediments, as these parameters are among the governing factors in transport capacity and competence.
The discharge of the Karra Khola increases as the basin area increases, as in common hydraulic geometry of most of the rivers. (Singh 2003, Stewardson 2005, Navratil and Albert 2010). The number of contributing tributaries to the main stem increases as the basin area increases which is the main reason for increase in the discharge of the river downstream.
Sediment transport is critical to understand as it is influenced by many factors among which three fundamental factors are sediment supply, capacity, and competence. The capacity of rivers to transport sediments in suspension is proportional to their discharge and thus to rainfall and to the hydrologic properties of their drainage basin (Tucker and Slingerland 1997). The larger the river the higher the rate of sediment transportation. According to Mathur and Da Cunha (2001), the Hwang He River in China, which is also called Yellow River, carries 796 million tons per year and the Mississippi River of United States carries 406 million tons per year.
The sediment transport of the Karra Khola is dominated by bed load sediment, and the total sediment transport in the midstream is about a million tons per sq. km, but it reduces at the extreme downstream transect and becomes three hundred thousand tons per year. Downstream reduction in slope and boundary shear stress and increase in d 50 size are the main reasons for decreasing rate of bed load transport at the downstream stretch.
The sediment yield of the Kathmandu basin computed by Lama and Tamrakar (2016) for the Manahara River, Bagmati River and the Bishnumati River are respectively, 1,936,008 t km -2 yr -1 , 70,040 t km -2 yr -1 and 44,530 t km -2 yr -1 . The sediment yield of the Karra Khola, which is 3,199 t km -2 yr -1 , is very low compared to the sediment yields of the rivers of the Kathmandu basin. The sediment yields of the Karra Khola tend to diminish with downstream distance of the river.
Stream power is the function of slope and hydraulic radius. With fluctuation of hydraulic radius, stream power or transport capacity also fluctuates, therefore, transport capacity of the Karra Khola also varies with seasons. The maximum transport capacity calculated for the Karra Khola is 22.17 KNs -1 . The other major reason of fluctuations in the stream power is owing to the change in slopes along the river course. Bagnold (1966) found better correlation of stream power with sediment transport compared to that with shear stress.
The specific stream power of the Karra Khola ranges from 20.98 Wm -2 to 2866.34 Wm -2 . There are various cases in which threshold specific stream powers vary in creating channel pattern changes and increasing in risk of erosion. Sediment Dynamics and Stability Status of The Karra Khola, Hetauda Dun Valley, Central Nepal Sub-Himalaya Brookes et al. (1983) studied in Britain and found that small channels having specific stream powers < 35 Wm -2 were stable while channels with greater values were unstable. Miller (1990) studied stream power in regard to historic floods with extreme discharges in the central Appalachian Mountain and reported widespread erosion including scouring and channel widening had occurred at specific stream power between 200 and 500 Wm -2 , and sever erosion was recorded when the specific stream power was above 300 Wm -2 . Similarly, Magilliagan (1992) suggested that for the catastrophic modification of stream channel, the specific stream power must exceed 300 Wm -2 . Yochum et al.
(2017) studied the rivers of Colorado Front Range after the flooding in 2013 and found for channel slopes less than 3%, that there were substantial channel widening at SSP > 230 Wm -2 , and high potential for avulsion and pattern change at SSP >480 Wm -2 , and also reported that there were possibility of infrequent to numerous eroded banks with SSP exceeding 700 Wm -2 . Considering, various researches the specific stream power of 300 Wm -2 , can be considered as the threshold for defining the boundary between low to high potential of river erosion and modification.
Stream power is capacity to transport sediments from the river channel, whereas dimensionless shear stress ratio is the measure of flow competence. The sediment starts to move only if the stress provided on the bed by the flow exceeds a critical shear stress value. That means when the dimensionless boundary shear stress exceeds the dimensionless critical shear stress of substrate sediment, the flow is competent enough to entrain and transport the river bed materials. As such situation applies to all the studied transects, except for the four downstream transects, the Karra Khola is found to be competent river up to the transect T15. Costa and O'Conner (1995) showed that shear stress and stream power are indicators of the geomorphic work done by floods. In fact, the stability condition of the river is governed by a balance between river's transport capacity and competency. Three conditions arise after considering thresholds of specific stream power of 300 Wm -2 and of flow competence of τ b* /τ c* = 1: (a) low capacity but with high competency, (b) high capacity with high competency, and (c) low capacity and low competency or high capacity but with low competency. In the first case, transport capacity is below the threshold in which severe erosion and channel modification is unlikely, but the river possesses high competency. In this situation, the river gets starved of sediments and may tend to erode its channel bed, thus undergoes degradation (e.g. T2 to T5, and T8). In the second case, where both transport capacity and flow competency are high, the sediments are transported without significant deposition because the river is competent enough to move the river bed materials. Such situation arises in T1, T6, T7, and T9 through T15. In the third case, whatever amount that the river transports, sediments tend to deposit due to incompetency of river to entrain and erode, therefore, allowing much of the sediments for deposition (e.g., T16, .T17, T18 and T19).

CONCLUSIONS
The channel sediments of the Karra Khola are sandy gravel, gravelly to muddy sand, and gravelly to sand mud. The sediment distribution of the Karra Khola ranges from fine clay to boulder size, with median grain size ranging between 2 to 73 mm, and are moderately well sorted to very poorly sorted. They are coarsely skewed in the upstream stretch and becomes neary symmetrical to finely skewed in the mid to downstream stretches. Overall median grain size tends to diminish downstream, while other measures varies without a distinct trends. River discharge tends to increase with increase in downstream distance and with increase in basin area.
The total sediment load is 8073-1077220 t yr -1 , in which bed load dominates over suspended load except in few cases. The total sediment yield varies from 561 to 148015 t km -2 yr -1 , and that of the whole basin is only 3199 t km -2 yr -1 , which is only around 2% of the maximum total sediment yield. The sediment yield does not follow its abundance with the increased basin area. The stream powers too vary along the distance downstream without a distinct trend of variation.
The Karra Khola possesses enough flow competence along its course except for the lower four transects located at the downstream segment. The nature of sediment transport capacity and flow competence of the Karra Khola is not uniform throughout its course, and it changes as the factors responsible for sediment transport changes. Both transport capacity and flow competency are taken to evaluate stability status of the Karra Khola. As a result, the upstream stretch of the river represents mainly degradation, the midstream stretch represents high rates of transport with lateral instability, and the downstream stretch represents aggradation. The midstream stretch of the river, which has high capacity, and competence of flow are laterally more erosive compared to the other stretches.

SOURCES OF FUNDING
This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

CONFLICT OF INTEREST
The author have declared that no competing interests exist.

ACKNOWLEDGMENT
Authors are thankful to the Central Department of Geology for providing necessary facilities for this research. Authors also thank Miss Manju Subedi, Mr. Ananta Joshi and Mr. Rythum Rai for their assistance during field work.