Abstract

West Lidder River, in the Northwest Greater-Himalayan mountain range, is the major source of irrigation and drinking water supplies for the Kashmir Valley with a population of seven million people. The major source of water for the whole Lidder River is snow and icemelt emanating from its two subcatchments East Lidder and West Lidder. Snowmelt significantly contributes to the evolution of drainage patterns in this area. Quantitative morphometry plays a vital role in routing the snowmelt and other hydrological processes. Morphometric analysis of the West Lidder River catchment was carried out using geospatial technique. The outcome revealed that the entire study area has uniform lithology and is structurally permeable. The high drainage density of all subwatersheds indicate more surface runoff. The morphometric analysis also indicates that the area is more prone to weathering due to very-coarse to coarse drainage texture. All the subwatersheds showed dendritic to subdendritic drainage pattern. An immense control of structure on the drainage in some subwatersheds is indicated by their high bifurcation ratios. Circulatory and elongation ratios show that the subwatersheds have elongated to circular shapes. From the integrated analysis of the morphometric parameters, important hydrologic behaviour of 17 subwatersheds could be inferred.

1. Introduction

The measurement and mathematical analysis of the configuration of the earth’s surface, shape, and dimension of its landforms is called morphometry [13]. To understand the evolution and behaviour of drainage patterns, several quantitative methods have been developed [4, 5]. In hydrology, basin drainage characteristics are fundamental in understanding various hydrological processes. Since watershed is the basic unit in hydrology; therefore, morphometric analysis at watershed scale is advantageous and preferable rather carry it out on individual channel or inconsistent segment areas. Watershed is an area of surface whose major runoff is conveyed to the single outlet and is the appropriate unit to study several processes of the land surface. For example, watershed is considered a fundamental erosional landscape element, wherein conspicuous interaction of land and water resources takes place. Being fundamental units of fluvial terrain, considerable research focal point has been done on watershed geometric characterization such as stream network topology and quantitative narration of shape, pattern, and drainage texture [5]. Hydrologic and geomorphic processes occur within the watershed, and morphometric characterization at the watershed scale reveals information regarding formation and development of land surface processes [6, 7] and thus provides a holistic insight into the hydrologic behaviour of a watershed. Basin travel time, time to hydrograph peak (basin lag time), and intensity of erosional processes operating at watershed scale can be predicted with better insight and accuracy from morphometric evaluation of a watershed. For unguaged watersheds wherein, in addition to hydrology, information regarding soil, geology, geomorphology, and so forth is also scarce, morphometric analysis provides a very good alternative to understand the underlying factors controlling the hydrological behaviour [8]. Furthermore, there are a myriad of practical applications of quantitative morphometric analysis such as river basin evaluation, watershed prioritization for soil and water conservation, and management of natural resources. Watershed morphometric analysis dispenses a quantitative description of the drainage system and thus enabling a better characterization of watersheds [9]. The role of landform processes, soil physical properties, and erosional characteristics in shaping the idiosyncrasy of different watersheds can be best evaluated through juxtaposing their morphometric parameters [7]. Using conventional techniques, morphometric characterization of many river basins and subbasins in different parts of the globe has been carried out [911]. With the advancement in geospatial and computer technology, assessment of the drainage basin morphometry has been more accurate and precise. Nowadays, using Geographical Information System (GIS) technique, various terrain and morphometric parameters of the drainage basins are evaluated with more ease and better accuracy. Satellite data and GIS tools have been successfully employed to generate data on the spatial deviations in drainage characteristics thus providing an insight into hydrologic conditions necessary for developing watershed management strategies [1214]. GIS, being a powerful tool for the manipulation and analysis of spatial information, provides a flexible environment for morphometric analysis. In the present study, morphometric characterization of 17 subwatersheds of the West Lidder River catchment was performed in GIS environment. The importance of morphometric analysis for this region lies in the fact that Lidder River forms the back bone of irrigation, drinking water supplies, and economy of the Kashmir valley. Snowmelt runoff contributes significantly to the drainage evolution of the area. Comparison of the quantitative morphometric parameters of subwatersheds of a river helps in understanding the geomorphological effects on the spatiotemporal variation of the hydrological functions. Moreover, some of the morphometric parameters such as circularity ratio and bifurcation ratio are input parameters in the hydrograph analysis [1518] and evaluation of surface water potential of an area [19]. In this context, this work represents a better understanding of hydrologic behaviour of the study area and its impact on various socioeconomic aspects of the valley of Kashmir. Moreover, a more rational approach of drainage morphometric analysis has been employed in this study compared to other automatically DEM-based drainage generation methods [20, 21]. It uses natural stream network and lake information to condition the DEM prior to the calculation of flow direction and flow accumulation grids. This greatly reduces error introduction in the form of spurious drainage.

2. Materials and Methods

2.1. Materials
2.1.1. Study Area

The West Lidder catchment occupies the Southeastern part of the Kashmir valley and is situated between 34°09′01′′N–34°14′00′′N and 75°06′00′′ E–75°23′29′′ E (Figure 1). Forming a part of the middle Himalayas and lying between the Pir Panjal range in the South and Southeast, the north Kashmir range in the Northeast, and Zanskar range in the Southwest, the study area has a unique geoenvironmental setting making it a suitable niche to many of the important glaciers such as Kolahoi. The West Lidder catchment has been engraved out by West Lidder River, a left bank tributary of River Lidder which finally drains as a right bank tributary of River Jhelum. It has a catchment area of 350.6 km², which constitute about 27.80 per cent of the total catchment area of river Lidder [29, 30]. The valley begins from the base of the snow field called the Kolahoi giving rise to West Lidder catchment. The latter joins the East Lidder which originate from another snow field called Sheshnag near the famous tourist town of Pahalgam. From this point onwards, it is referred as River Lidder. The latter joins the Jhelum (upper stream of Indus River) at Gur Village after travelling a course of 70 km [30]. Combined action of glaciers and rivers has bestowed West Lidder catchment a marked topography. The valley possesses distinctive climatic characteristics because of its high altitude location and its geophysical setting, being enclosed on all sides by high mountain ranges. The study area receives precipitation in the form of both rain and snow. The mean annual precipitation at Pahalgam rain gauge station is 1267.2 mm. The highest rainfall of 1629 mm was recorded during the year 1994 and the lowest of 899.9 mm during the year 2000. The area receives the highest rainfall in the month of March (208.8 mm) and the lowest in the month of October (45.9 mm). Such a climatic regime is peculiar characteristic of sub-Mediterranean type of climate. Moreover, 70 per cent of its annual precipitation concentrated in winter and spring months [31].

2.1.2. Data Sets Used

The data sets used in this study for accomplishing the respective objectives are as follows.(i) Quantitative morphometry: ASTER DEM 30 m.(ii) Survey of India (SOI) topographical maps at 1 : 50,000; Dated: 1969: Natural drainage delineation.

2.2. Methods
2.2.1. Drainage Generation

Drainage was generated from Advanced Space borne Thermal Emission and Reflection Radiometer (ASTER) 30 m digital elevation model (DEM) in GIS environment. The Arc Hydro approach was used for drainage generation which is more logical and consistent when compared to a manual approach [32]. There is no distinction between water filled stream channels and land surface in a DEM. In order to incorporate already existing streams and lakes into the DEM for drainage analysis, a process termed as DEM manipulation is followed. Natural drainage present in SOI toposheets was digitized and used to manipulate the DEM. For proper determination of flow direction and flow accumulation, DEM sinks were identified and filled. Boundaries of the 17 subwatersheds were derived by defining pour point for each subwatershed (Figure 1). The pour point is the location where water drained from whole of the watershed flows into the main river. Based on the cumulative number of the upstream cells draining to each cell, stream network in both the watersheds was defined. We used a critical threshold 0.06 which represents the 6% area in all 17 subwatersheds for defining streams. The critical threshold is the minimum upstream drainage required to initiate a stream. Area of the watershed was evaluated by calculating the geometry of the derived watershed polygons, and length of the watershed was calculated by summing the length of the main stream channel and the distance from the top of the main channel to the watershed boundary. By summing the lengths of all stream segments in each subwatershed the total stream length was calculated.

2.2.2. Quantitative Morphometry

In the present study, morphometric analysis of the parameters, namely, stream order, stream length, bifurcation ratio, relief ratio, drainage density, drainage intensity, drainage texture, stream frequency, drainage texture, form factor, length of overland flow, constant channel maintenance, circulatory and elongation ratio, area, perimeter, and length of all the 17 subwatersheds have been carried out using the standard mathematical formulae given in the Table 1 [9, 10, 19, 2228]. The values of various basin characteristics required for calculating morphometric parameters are shown in Table 2.

3. Results and Discussion

3.1. Stream Order ( )

The highest stream order among the 17 subwatersheds is five and is shown by four subwatersheds: WS7, WS8, WS9 and WS10. The lowest stream order is two and is shown by WS17. A perusal of Table 3 and Figure 1 indicates that the 17 subwatersheds draining into West Lidder River are contributing surface runoff and sediment loads differentially due to variations in their physical characteristics. There are only 7 subwatersheds in the East side of the river among which five watersheds have and only two watersheds have . In the West side of the river, there are 10 watersheds out of which four watersheds have , two watersheds have , three watersheds have , and only one watershed has . Drainage maps of all the 17 subwatersheds are shown below (Figures 4(a), 4(b), and 4(c)). Higher stream order is associated with greater discharge, and higher velocity [33]. West side of the river clearly contributes more to discharge and since higher velocity enhances the erosion rates, therefore, this side also contributes higher sediment loads into the West Lidder River. Further, the total number of stream segments decrease with stream order. This is referred to as Horton’s law of stream numbers. Any deviation indicates that the terrain is typified with high relief and/or moderately steep slopes, underlain by varying lithology and probable uplift across the basin [34]. In practice, when logarithms of the number of streams of a given order, are plotted against the order, the points lie on a straight line [22]. Similar geometric relationship was also found to operate between stream order and stream numbers in all subwatersheds of this study area. It indicates that the whole area has uniform underlying lithology, and geologically, there has been no probable uplift in the basin (Figure 2).

3.2. Stream Length (Lu)/Mean Stream Length (Lsm)

Analysis of the results shown in Tables 3 and 4 show that the total length of stream segments is the maximum in case of first order streams. It decreases as order increases in all the subwatersheds (Tables 3 and 4; Figure 4). The results reaffirm the fact that the area is underlain with uniform lithology with no probable basin upliftment. The observation demonstrates that the area depends only on the drainage characteristics for movement of water. Further, since there are more number of watersheds in the West side of the river and the total stream length of all the orders from this side is greater than the East side, therefore, the longer travel times [35] makes this side hydrologically very active. From this observation, it is concluded that during a same intensity storm event over whole of the East or West side of the river [36], the East side will show short basin lag times [37] compared to West side under similar soil moisture [38, 39] and vegetation cover [39]. Moreover, Table 4 indicates that Lsm in these subwatersheds range from a minimum of 0.21 km for stream order 2 of WS5 to a maximum of 6.69 km for the order 3 of WS13. According to the Horton’s law of stream lengths, Lsm of any given order is greater than that of lower order. This geometric relationship can be seen in Figure 3. A comparative analysis of Lsm and stream length ratio RL of all the subwatersheds is shown in Table 4.

3.3. Stream Frequency/Channel Frequency (Fs)

Stream frequency (Fs) is the total number of streams of all orders per unit area. The analysis of the results shown in Table 5 shows that Fs is maximum in subwatershed WS12 (8.19/km2), followed by WS8 (7.98/km2), WS9 (7.82/km2), WS13 (7.05/km2), WS7 (7.04/km2), WS5 (7.02/km2), WS10 (6.70/km2), WS6 (6.52/km2), WS16 (6.15/km2), WS1 (5.93/km2), WS4 (5.78/km2), WS15 (5.64/km2), WS3 (5.63/km2), WS2 (5.15/km2), WS11 (4.52/km2), WS14 (4.12/km2) and WS17 (3.47/km2). Fs has been related to permeability, infiltration capacity, and relief of watersheds [40, 41]. The values observed in West Lidder indicate that W12 is having rocky terrain and very low infiltration capacity out of all the 17 subwatersheds. Further, it is noted that Fs decreases as the stream number increases. Fs of WS17 reveals that this subwatershed is comparably covered with good amount of vegetation and has very good infiltration capacity. Overall, the results of Fs reflect early peak discharge for subwatersheds in order of their decreasing Fs resulting in flash floods, while the discharge from WS17 takes longer time to peak because of low runoff rates due to lesser number of streams [4244].

3.4. Bifurcation Ratio (Rb)/Mean Bifurcation Ratio (Rbm)

Analysis of the results shown in Table 5 shows that mean bifurcation ratio (Rbm) of the 17 subwatersheds is 3.6, 3.7, 4, 4.1, 2.7, 4.8, 4.7, 4.5, 6.4, 5.3, 3.6, 10.2, 9.1, 4.6, 4.2, 4.4, and 3.0, for WS1 to WS17, respectively. Rbm does not precisely remain constant from one order to the next, because of possibility of variations in watershed geometry and lithology, but it tends to be a constant throughout the series. High Rbm indicates early hydrograph peak with a potential for flash flooding during the storm events [45, 46]. Moreover, wherever in a basin, powerful geological control dominates; Rbm shows only a small variation for different regions on different environment [10]. Higher Rbm values are the characteristics of structurally more disturbed watersheds with a prominent distortion in drainage pattern and vice versa [14]. Maximum Rbm is seen for WS12 (Rbm = 10.2), and thus it will show early hydrograph peak (smaller basin lag time) which also indicates strong structural control on the drainage development for this watershed. The minimum Rbm is observed for WS5 (Rbm = 2.79), indicating delayed hydrograph peak.

3.5. Relief Ratio (Rh)

Quantitatively, it is the measurement of the overall steepness of a drainage basin. Also, it is an indicator of the intensity of erosion processes operating on the basin slopes [47]. Rh normally increases with decreasing drainage area and size of a given drainage basin [48]. In the present study as depicted in Table 7, Rh ranges from a minimum of 0.14 in WS9 to a maximum of 0.61 in WS17. Higher values of Rh indicate that intense erosion processes are taking place. This indicates that WS17 is more susceptible to erosion, and WS9 is the least among all the subwatersheds of the study area if this parameter alone is considered for erosion intensity analysis.

3.6. Drainage Density (D)

Factors which affect D are characteristically same as those affecting length of the stream, namely, resistance to weathering, permeability of rock formation, climate, vegetation, and so forth. The travel time by water within the basin is controlled by D [4951]. In humid regions, it varies between 0.55 and 2.09 km/km² [52]. Generally, D has a low value in regions underlain with highly resistant permeable material with vegetative cover and low relief. High drainage density is observed in the regions of weak and impermeable subsurface material and sparse vegetation and mountainous relief. Perusal of the D values from Table 5 indicates that WS17, WS14, and WS11 have low D values (below 2.0 km/km²), while WS1, WS15, WS16, WS4, WS3, WS2, WS5, and WS10 have moderate D values (2.0–2.5 km/km²). Subwatersheds WS6, WS9, WS12, WS7, WS8, and WS13 have high D values (between 2.5 and 3.0 km/km²). The low and moderate D watersheds reveal that they are composed of permeable subsurface material, good vegetation cover, and low relief which results in more infiltration capacity and comparably are good sites for ground water recharge as compared to high D watersheds [53, 54]. On the basis of D, WS17 will have the greatest basin lag time, while WS13 will demonstrate the shortest lag time.

3.7. Drainage Texture (Rt)

Drainage texture (Rt) is influenced by infiltration capacity [22]. There are five different texture classes: very coarse (<2), coarse (2–4), moderate (4–6), fine (6–8), and very fine (>8) [55]. According to this classification, WS17 has very coarse drainage texture, WS1, WS2, WS4, WS5, WS11, WS14, and WS16 have coarse Rt, WS3, WS6, WS12, WS13, and WS15 have moderate Rt, and WS7, WS8, WS9, and WS10 have very fine Rt (Table 6). Hydrologically very coarse texture watersheds have large basin lag time periods [56] followed by coarse, fine, and very fine texture classes. This indicates that WS17 (Rt = 1.62) shows longer duration to peak flow, while WS9 (Rt = 12.75) shows the shortest.

3.8. Form Factor (Rf)

The basins with high form factors have high peak flows of shorter duration, whereas elongated subwatersheds with low form factors have lower peak flow of longer duration [5759]. In this study, WS2, WS3, WS7, WS10, WS11, WS13, WS14, and WS15 show lower form factor (Table 7) indicating elongated shape and suggesting flat hydrograph peak for longer duration. Flood flows of such elongated basins are easier to manage than those of the circular basin. WS5, WS6, WS8, WS9, WS12, and WS16 have slightly circular shape as suggested by moderately higher Rf. Subwatersheds WS1, WS4, and WS17 have high Rf indicating that they have developed into quite circular to rectangular shape. Watershed morphology has profound impacts on the watershed hydrology [60].

3.9. Elongation Ratio (Re)

The value of elongation ratio (Re) generally varies from 0.6 to 1.0 associated with a wide variety of climate and geology. Values close to 1.0 are typical of regions of very low relief whereas, that from 0.6 to 0.8 are associated with high relief and steep ground slope [7]. These values can be grouped into three categories, namely, circular (>0.9), oval (0.9–0.7), and less elongated (<0.7). Re for the subwatersheds WS13, WS15, WS12, WS3, WS11, WS7, and WS10 is less than 0.70 which indicates that subwatersheds are elongated with high relief and steep slope. WS2, WS14, WS9, WS16, WS8, WS5, and WS6 sub-watersheds fall in oval category, while subwatersheds WS1, WS17 and WS4 fall in the circular category (Table 7). The circular watersheds have quick, though lower, hydrograph peak compared to the oval and elongated watersheds [61, 62]. This indicates that WS13 will show delayed time to peak flow, while WS4 will show shorter time to peak.

3.10. Circulatory Ratio (Rc)

It is that ratio of the area of the basin to the area of circle having the same circumference as the perimeter of the basin [25]. Many of the important basin characteristics influence Rc such as length and frequency of streams, geological structures, land use/land cover, climate, relief, and slope of the basin. In the present case, Rc for subwatersheds WS8, WS2, WS6, WS7, WS10, WS5, WS17, WS14, and WS16 is in the range from 0.53 to 0.62 indicating that the area is characterized by high relief, elongated, and permeable surface resulting in greater basin lag times, while WS3, WS9, WS1, WS13, WS12, WS15, WS4, and WS11 have lower circularity ratios indicating low relief and impermeable surface resulting in lower basin lag times (Table 7). This indicates that WS16 will show delayed time to peak flow, while WS3 will show shorter time to peak. Rf, Re, and Rc significantly influence the hydrological response of the watersheds as basin shape, and the arrangement of stream segments combine to influence the size and shape of flood peaks [63].

3.11. Length of Overland Flow (Lg)

Length of overland flow is one of the most important independent variables affecting both hydrologic and hydrographic development of drainage basins [24]. Table 7 reveals that Lg for the 17 subwatersheds in chronological order is 0.22, 0.16, 0.18, 0.18, 0.16, 0.14, 0.11, 0.11, 0.76, 0.16, 0.29, 0.11, 0.11, 0.32, 0.2, 0.18, and 0.71, respectively. The Lg for subwatersheds WS2-WS10, WS12, WS13, and WS16 indicates steep slopes and shorter flow paths, while as Lg for WS1, WS11, WS14, WS15, and WS17 indicates gentle slopes and longer flow paths.

3.12. Constant of Channel Maintenance (C)

In the present study, C varies from 0.33 for W8 and W9 to 0.84 for W17 as is shown in Table 6. The reciprocal of the drainage density (D) is constant of channel maintenance and signifies how much drainage area is required to maintain a unit length of channel. Low values of C in case of WS8 and WS9 indicates that among all the 17 subwatersheds, these two are associated with weakest or very low-resistance soils, sparse vegetation and mountainous terrain, while the watershed WS17 is associated with resistance soils, vegetation and comparably plain terrain [64].

3.13. Compactness Coefficient (Cc)

Compactness coefficient expresses the relationship of a basin with that of a circular basin having the same area. A circular basin yields the shortest time of concentration before peak flow occurs in the basin. indicates that the basin completely behaves as a circular basin. indicates more deviation from the circular nature of the basin. The values for all the subwatersheds range from 1.27 for WS16 to 2.20 for WS3, as seen in Table 6. Consequently, WS16 has the greatest deviation from the circular nature, and on the basis of this parameter alone, it will have the longest time of concentration before peak flow occurs in the basin out of all the subwatersheds [6062].

3.14. Shape Index (Sw)

Rate of water and sediment yield along the length and relief of the drainage basin is largely affected by the shape. The shape index values for subwatersheds of the study area range from 0.94 in WS4 to 4.02 in WS13 as shown in Table 6. In terms of Sw only, WS13 will have the shorter basin lag time, while WS4 will have the longer basin lag time.

3.15. Texture Ratio (T)

Texture ratio is one of the most important factors in the drainage morphometric analysis, which depends on the underlying lithology infiltration capacity, and relief aspect of the terrain is the texture ratio [65, 66]. As shown in Table 7, T varies between 1.21 in case of WS17 and 8.09 in case of WS9. Hydrologically, it can be said that WS17 in terms of T alone will have the longest basin lag times and WS9 will have the shortest. This reaffirms the results and conclusions drawn from the values of U, D, Fs, and Rt in case of W17.

3.16. Ruggedness Number (Rn)

The ruggedness number indicates the structural complexity of the terrain in association with the relief and drainage density. It also implies that the area is susceptible to soil erosion [67, 68]. In present study, the ruggedness number is minimum in case of WS2 ( ) and maximum in WS3 as seen in Table 6. This indicates that WS2 is least susceptible to erosion and WS3 is most susceptible among all the subwatersheds of the study area.

3.17. Shape Factor (Rs)

The shape factor observed in the study area is minimum in case of WS3 and maximum in case of WS16 , as shown in Table 6. This parameter is similar in interpretation to circularity ratio, elongation ratio, and form factor. It gives an idea about the circular character of the basin. The greater the circular character of the basin is, the greater is the rapid response of the watershed after a storm event [60]. Therefore, in terms of Rs only, WS3 has the longest basin lag time, and WS16 the shortest.

4. Conclusion

From the integrated analysis of the results on the morphometry of the 17 subwatersheds, it is concluded that the whole study area has uniform underlying lithology making the hydrological response in these subwatersheds a direct function of the geomorphology, the topography, and the existing vegetation conditions only. The spatial variability of the morphometric parameters analyzed in this study is quite significant. Since, the hydrology of a watershed changes significantly due to the spatial variations of the morphometric parameters, the subwatersheds will therefore also exhibit differential hydrological behaviour. Overall results indicated that comparatively WS16, WS13, WS12, and WS9 shall contribute dominantly to the stream runoff in the West Lidder catchment, and on the contrary, WS17, WS5, and WS3 shall contribute the least due to the cumulative and integrated effect of the morphometric parameters whose significance and influences on the hydrological behaviour have been discussed in detail in the preceding sections. Because of the close relationship between the morphometric parameters and the mean annual flood, it is concluded, on the basis of the range of the morphometric analysis, that WS16, WS13, WS12, and WS9, in that order, are more prone to flooding compared to the other 13 subwatersheds of the West Lidder catchment. The hydrological behaviour of these 17 subwatersheds shall have profound influences on the flood vulnerability and flood risk in the downstream portions of the Lidder River. It is hoped that information and knowledge generated under this study shall guide the informed planning and decision-making for flood disaster risk reduction in the Lidder watershed that forms an important tributary of the Jhelum basin.

Acknowledgments

This research work has been accomplished under a research grant provided by the Department of Science and Technology, Government of India for the project titled “Integrated Flood vulnerability Assessment for Flood Risk Management and Disaster Mitigation.” The authors express their gratitude to the funding agency for the financial assistance. The authors also express their sincere thanks to anonymous reviewers whose valuable comments have improved the quality of this paper.