Abstract

Mountainous rivers are the most significant source of water supply in the Himalayan provinces of India. The drainage basin dynamics of these rivers are controlled by the tectonomorphic parameters, which include both surface and subsurface characteristics of a basin. To understand the drainage basin dynamics and their usefulness in watershed prioritisation and management in terms of soil erosion studies and groundwater potential assessment and flood hazard risk reduction in mountainous rivers, morphometric analysis of a Himalayan River (Supin River) basin has been taken as a case study. The entire Supin River basin has been subdivided into 27 subwatersheds and 36 morphometric parameters have been calculated under four broad categories: drainage network, basin geometry, drainage texture, and relief characteristics, each of which is further grouped into five different clusters having similar morphometric properties. The various morphometric parameters have been correlated with each other to understand their underlying relationship and control over the basin hydrogeomorphology. The result thus generated provides adequate knowledge base required for decision making during strategic planning and delineation of prioritised hazard management zones in mountainous terrains.

1. Introduction

Morphometry is the measurement and mathematical analysis of the configuration of the earth’s surface and of the shape and dimension of its landforms [1]. The form and structure of drainage basins and their associated drainage networks are described by their morphometric parameters. Morphometric properties of a drainage basin are quantitative attributes of the landscape that are derived from the terrain or elevation surface and drainage network within a drainage basin. Application of quantitative techniques in morphometric analysis of drainage basins was initially undertaken by Horton et al. [28] from topographic maps using manual methods.

Remote sensing and Geographical Information System (GIS) techniques are increasingly being used for morphometric analysis of drainage basins throughout the world [913]. Quantitative techniques have been applied to study the morphometric properties of different drainage basins in India [1423]. Several authors have studied morphometric properties of drainage basins as indicators of structural influence on drainage development and neotectonic activity [2427]. In many studies morphometric analysis has been used to assess the groundwater potentiality of the basins and to locate suitable sites for construction of check dams and artificial recharge structures [2832]. Watershed prioritisation based on morphometric characteristics has also been carried out and aids in the mapping of high flood potential and erosion prone zones [3337].

Present study bridges the connection between surface morphometry and subsurface geology of a drainage basin to produce effective information as a part of basin management. So the objective of the present research is to study the morphometric parameters of Supin River basin and to identify the influence of the underlying geology on the morphometric parameters of the basin and finally to generate a substantial knowledge base regarding the relationship between surface morphometry and subsurface lithology for integrated basin management.

2. Drainage Basin Setup

The Supin River basin is located in Uttarkashi district of Uttarakhand (Figure 1) and situated between 78°10′–78°38′E and 31°00′–31°04′N. The tributaries of Supin River, such as, Har Ki Dun Gad, Borasu Gad, and Ruinsara Gad, are fed by Jamdar glacier (Figure 2(a)). The basin covers an area of 565.41 km2 and has a perimeter of 190.47 km. The study area has three climatic zones: subtropical (1,500–1,700 m), temperate (1,700–3,500 m), and alpine (>3,500 m). The region receives heavy snowfall between November and March. The rainfall varies from 1,000 to 1,500 mm annually.

In Garhwal Himalaya, one of the major NW-SE trending thrusts is the Main Central Thrust (MCT). The MCT separates the Greater Himalayan Sequence (GHS) from Lesser Himalayan Sequence (LHS) [38]. The Main Central Thrust (MCT) separates the basin into two halves on its hanging wall (MCT sheet) and footwall sides which represent the Ramgarh Thrust (RT) Sheet [3944] (Figure 2(b)). The Munsiari Thrust (MT) separates the Wangtu gneiss and Rampur group of the Lesser Himalayan Sequence (LHS) within the RT Sheet [45] (Figure 2(b)). The basin is underlain by rocks belonging to three main geological formations: Martoli, Vaikrita, and Garhwal. The upstream portion of the Supin basin, which is mainly drained by its three major tributaries, namely, Har Ki Dun Gad, Borasu Gad, and Ruinsara Gad, is underlain by granite-gneisses and two mica schists belonging to the MCT sheet. The middle portion of the basin is underlain by rocks consisting of Greater Himalayan Gneisses (Augen gneisses and porphyritic granites) and phyllites, quartzites, and biotite grade schists separated by the Munsiari Thrust (MT), whereas near the mouth of the basin, it consists of leucogranite (Figure 2(c)) [4648]. The presence of MCT and MT fault zones across the drainage network of the Supin River basin influences the drainage and relief characteristics of the Supin River.

3. Methodology

The extraction of drainage network has been done from Shuttle Radar Topographic Mission (SRTM) Digital Elevation Model (DEM) (90 m) data using TNTmips software environment. The generation of depressionless DEM is always the preparatory step for morphometric analysis of drainage basin. Depressions are data errors or result from the averaging involved in assigning elevation values to cells (pixels) of finite area. These spurious depressions interfere with the correct routing of flow paths during the watershed analysis, especially in areas of low relief. The Watershed process solves this problem by first locating and filling the depressions. This depressionless DEM is used to compute the flow direction and flow accumulation raster. Further simulation of these two raster produces the standard flow paths and subwatersheds.

The Supin River basin has been classified into 27 subwatersheds. Only those watersheds have been considered for this study which includes streams of at least three different orders. Thereafter, 36 morphometric parameters (Table 1) have been computed for the entire Supin basin, as well as for each of the subwatersheds. The morphometric parameters have been evaluated from four different aspects—drainage network, basin geometry, drainage texture, and relief. The different parameters were then correlated to understand how they interact with and influence each other. The reliability of the correlation of determination has been tested with the help of Students’ t-test and calculated values. Hierarchical cluster analysis (Euclidean distance) has been used as the measure of association which enables the grouping of the subwatersheds into five major categories according to the four morphometric aspects. Hypsometric curve for the entire Supin River basin has been computed along with hypsometric integral (HI) values for all the subwatersheds.

4. Results and Discussion

The Supin River having a length of 40.08 km drains an area of 17.68 km2 (Table 2). The two major thrusts, that is, MCT and MT, cross the upstream and downstream sections of the basin, respectively. The MCT crosses the basin in northwest-southeast direction with subwatersheds like WS3, WS15, WS16, and WS17 lying on its hanging wall and WS18 and WS19 lying on its footwall. The MT crosses the basin in an east-west direction and WS6, WS8, WS20, WS21, WS26, and WS27 are along the hanging wall of MT whereas WS2, WS22, WS23, WS24, and WS27 are along the footwall of MT. In case of WS4, WS5, and WS7, the northern boundary of the subwatershed nearly coincides with MCT and the southern boundary nearly coincides with MT (Figure 2). The major 5th order tributary (Ruinsara Gad) of the Supin River follows the trend of MCT in the SE part of the drainage basin. Moreover, a few other 4th and 5th order tributaries also follow the trend of MCT and MT fault zones. These fault zones are usually zones of weak, fractured, and brecciated zones which are easily incised by streams. Along the stream profile of the Supin River, at a few places sudden change in slope (~15°) of the stream has been marked as knick points () which are controlled by the dip of MCT and MT fault zones (Figure 2(c)).

In the following section the various morphometric parameters have been discussed with regard to the derived cluster groups (Figure 3).

4.1. Drainage Network

Segmentation and hierarchical ordering of streams is necessary to address the hydrodynamic character of a drainage basin. Stream ordering has been done for Supin River basin following the hierarchical ranking proposed by Strahler [6]. Two 5th order streams of Har Ki Dun Gad and Ruinsara Gad combine to form Supin, the 6th order stream in the basin (Figure 2(a)). The total stream length of Supin basin is 1,439.204 km (Table 1) of which the 1st and 2nd order streams constitute 87.05%. The stream length ratio () varies from 1.05 to 2.06 and is high for 3rd and 5th order streams in the Supin basin (Table 2). With increasing stream order there is a decrease in stream number () and a simultaneous increase in mean stream length () (Table 2). The RHO coefficient () and bifurcation ratio () values for Supin basin range from 0.62–2.10 (Table 1) and 2 to 5.14 (Table 2), respectively.

The variation of between successive stream orders of Supin River basin is due to the greater number of streams belonging to lower orders indicating that the basin is still in its youthful stage of development (Table 2). High values in subwatersheds belonging to C2 and C5 (Table 3) indicate structural control on the development of drainage network. The value signifies the storage capacity of a basin and determines the relationship between drainage density and physiographic development of the basin. Subwatersheds belonging to C4 and C5 (Table 3; Figure 3(a)) having high values of are at a greater risk of being eroded by the excess discharge during flood.

4.2. Basin Geometry

Basin shape is controlled by structure, lithology, relief, and precipitation and varies from narrow elongated forms with irregular basin perimeter to circular or semicircular forms. Circularity ratio () of Supin basin ranges from 0.30 to 0.56 (Table 1) with high values in WS1, WS4, WS8, WS9, WS12, WS17, WS21, WS23, and WS24 and low values in WS7, WS10, WS11, WS16, and WS22 (Figure 4(a)). High values of form ratio () and elongation ratio () are found in WS3, WS5, WS17, WS19, WS21, WS23, and WS25, whereas low values are found in WS11, WS12, WS13, WS14, and WS16 (Figures 4(b) and 4(c)). The relative spacing of channels in a drainage basin is expressed by texture ratio () and drainage texture (). The and values of Supin basin range from 0.27 to 1.88 and 0.44 to 2.55, respectively (Table 1). High values of are found in subwatersheds located in the upper reaches of the basin, whereas low values are found in subwatersheds located near the mouth of the basin (Figure 4(d)).

bears a strong negative relationship with compactness coefficient (), whereas has a positive relationship with (Table 4) [28]. Low values of (C4 and C5) (Table 3; Figure 3(b)) are associated with high relief and steep slopes and imply the youthful nature of these subwatersheds. The subwatersheds belonging to C1, C2, and C3 (Figure 3(b)) having high values of , , , , and and low values of are more circular in shape than those belonging to C4 and C5 (Table 3). Although the circular subwatersheds are more efficient in the discharge of run-off [56], they are at greater risk from flood hazard because they have a very short lag time and high peak flows than the elongated basins [57]. The elongated subwatersheds, on the other hand, have low side flow for shorter duration and high main flow for longer duration and are less susceptible to flood hazard.

4.3. Drainage Texture Analysis

Texture indicates the amount of landscape dissection by a channel network and includes stream frequency (), drainage density (), constant of channel maintenance (), length of overland flow () and infiltration number (). High values of and are found in WS3, WS5, WS11 WS14, WS22, and WS25, whereas low values of and are found in WS1, WS6, WS18, WS23, and WS27 (Figures 5(a) and 5(b)). The value of Supin basin varies from 0.33 to 0.52 (Table 1). The value of Supin basin ranges from 0.16 to 0.26 and value ranges from 2.94 to 8.61 (Table 1). Low values and high values are found in WS3, WS5, WS11, WS12, WS14, and WS22 (Figures 5(c) and 5(d)) which have corresponding high and values (Figures 5(a) and 5(b)).

and provide a numerical measurement of landscape dissection and run-off potential [16] and bears a negative relationship with and (Table 4). Subwatersheds belonging to C3, C4, and C5 (Table 3; Figure 3(c)) with low values of and have developed numerous drainage lines on the surface. The high values of and in these subwatersheds maybe due to high relief, steep slopes, and also low permeability of the underlying rocks.

4.4. Relief Characteristics

The relief properties in morphometric analysis bring into consideration the influence of aspect and height over a large basin area. Relief ratio (), ruggedness number (), and dissection index () indicate the erosion potential of the processes operating within a drainage basin. The total basin relief (Z-z) of Supin basin is 4,296 m (Table 1). The value of Supin basin ranges from 0.07–0.44 with high values in WS2, WS5, WS17, WS19, WS21, WS23, WS24, WS25, WS26, and WS27 and low values in WS1, WS9, WS10, WS11, WS12, WS13, and WS14 (Figure 6(a)). High values of and are found in WS3, WS4, WS5, WS7, WS9, WS19, and WS20 which indicates the steepness of slope and high degree of dissection of these subwatersheds, whereas low values are found in WS10, WS11, WS12, WS13, and WS14 (Figures 6(b) and 6(c)). is related to the channel gradient and has a negative relationship with basin length (Table 4). Subwatersheds having high values of (C3 in Figure 3(b)) also have high values of (C1 and C5 in Figure 3(d)) since has a positive correlation with (Table 4). The low of Supin basin (0.27–0.61) indicates that the basin is moderately dissected (Table 1).

4.5. Hypsometry Analysis

Hypsometric curve of a drainage basin represents the relative proportion of watershed area below or above a given elevation. It is a measure of the erosional state or geomorphic age of a drainage basin as it represents the mass of drainage basin remaining above a basal plane of reference. Convex hypsometric curves indicate youthful stage, S-shaped curves indicate a mature stage, and concave curves indicate peneplain stage [5]. Here the hypsometric curve has been treated as a cumulative probability distribution and statistical moments have been used to describe the curve quantitatively [58]. The hypsometric curve has been derived by fitting a 5th order polynomial function which is as follows: The estimated coefficients are , , , , , and which can be obtained from regression of the hypsometric curve (Figure 7(a)). The statistical moments of hypsometric curve include skewness of the hypsometric curve (−0.474), kurtosis of the hypsometric curve (2.495), skewness of the hypsometric density function (−0.413), and kurtosis of the hypsometric density function (2.392). These are collectively known as the derived parameters of hypsometric curve. The area below the hypsometric curve is known as the hypsometric integral (HI) and is calculated using the following formula: where HI is hypsometric integral, is mean elevation, is maximum elevation, and is minimum elevation. HI has been calculated for all the subwatersheds of Supin River basin. The convex hypsometric curve (Figure 7(a)) of Supin River basin shows HI value of 0.58. The HI values for the subwatersheds range from 0.36 to 0.62 (Figure 7(b)).

The derived parameters of hypsometric curve are sensitive to subtle changes in overall basin slope and basin development as the mass is removed by erosion over a long geological time period [59]. Present study shows that head ward development of the main stream and its tributaries produce high values for hypsometric skewness [59]. On the other hand, the density function of the curve is closely related to rates of change in overall basin slope and tendency toward geomorphic equilibrium [59]. Density skewness interprets the behaviour of slope change in the basin. Positive value of density skewness is an indication of fluvial dominance over the landforms of this region. Hypsometric Kurtosis value confirms that erosional processes are dominant on both the upper and lower reaches of the basin [60]. Mid basin slope is moderate as the density Kurtosis value is platykurtic in nature. The HI value of Supin basin (0.58) indicates that the basin is at a youthful stage of development. [6]. A few subwatersheds showing low HI values (Figure 7(b)) have attained relatively mature stage and are in a state of reaching dynamic equilibrium. Most of the subwatersheds showing high HI are located on the hanging wall of MCT and MT (Figures 2(b) and 7(b)) which also indicate a subsurface control on the maturity of the Supin River basin.

5. Conclusion

There is a structural influence on drainage development with trellised pattern being the dominant drainage pattern. Channel extension in this mountainous region takes place mainly through headward erosion which is characteristic of basins in early mature stage of development. The area is well drained by the 1st and 2nd order streams having highest drainage area which produces a rugged terrain, moderately dissected by deep incised valleys. The high relief and moderate permeability of the surface in relatively circular subwatersheds with high circularity ratio generates high run-offs which make these subwatersheds more susceptible to floods, soil erosion, and debris flow. In a few elongated subwatersheds the management of flood flow is easier because of the low side flow for shorter duration and smaller peak flows for longer duration.

Fourth and 5th order streams flowing along regional fault zones may generate landslides due to increase in the level of erosion. In addition, an increase in the incision of the streams along the weak and fractured fault zones results in increase in the sediment load of the streams which in turn may trigger flash floods. Sudden topographic breaks along the 6th order stream profile of the basin at multiple places influence the tectonomorphic landforms developed along Supin River.

The significance of studying morphometric properties of Supin basin lies in the fact that it will help in future watershed management and hazard management studies. Due to the mountainous nature of the terrain the region suffers from frequent flash floods and landslides. Hence, such type of studies will provide knowledge and database for decision making for strategic planning and delineation of prioritised hazard management zones. Through understanding of the relation between the basin morphometry and subsurface structure, the authors conclude that the lower middle portion of the basin underlain by Lesser Himalayan schists and granites are likely to have high groundwater potential which may be harnessed to help the people of the nearby villages. Moreover, it may also help in assessing the groundwater potential of the region and delineating effective water harvesting sites. These morphometric techniques may be applied to other mountainous river basins around the globe. Morphometric analysis thus has a wider significance in watershed prioritisation and management, soil erosion studies, groundwater potential assessment, and flood hazard risk reduction in the Supin watershed that forms an important tributary of the Tons River basin.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Assistance from the postgraduate students of the Department of Geography, University of Calcutta, during field work was gratefully acknowledged. Infrastructural facilities were provided by the Department of Geology, University of Calcutta, and the Remote Sensing & GIS Department of the Vidyasagar University. Thanks are due to Professor Sunando Bandyopadhyay, Department of Geography, University of Calcutta. Comments from the anonymous reviewers are gratefully acknowledged.