1
1
com
ScienceDirect
Aquatic Procedia 4 (2015) 1322 – 1330
Abstract
Watershed prioritization has gained importance in natural resources management, especially in the context of watershed
management. Delineation of watersheds within a large drainage basin and their prioritization is required for proper planning and
management of natural resources for sustainable development. Delineation of potential zones for implementation of conservation
measures above the entire watershed at similar occurrence is inaccessible as well as uneconomical; therefore it is a prerequisite to
apply viable technique for prioritization of sub watersheds (SW). The present research attempted to study various morphological
characteristics and to implement Geographical Information System (GIS) and Multi Criteria Decision Making (MCDM) through
Fuzzy Analytical Hierarchy Process (FAHP) techniques for identification of critical sub watersheds situated in transaction zone
between mountainous and water scarcity region of kallar watershed, Tamil Nadu. The morphometric characterization was
obtained through the measurement of three distinct linear, areal and relief aspects over the eleven sub-watersheds. The
morphometric characterization showed vital role in distinguishing the topographical and hydrological behavior of the watershed.
Each morphometric parameters were ranked with respect to the value and weightings obtained by deriving the relationships
between the morphometric parameters obtained through classification of the SW by associating the strength of Fuzzy Analytical
Hierarchy Processes (FAHP). Based on FAHP approach, sub watersheds were evaluated and divided into five prioritization
zones: very less, less, medium, high and very high classes. The FAHP techniques is a practical approach for identification of the
2214-241X © 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license
(http://creativecommons.org/licenses/by-nc-nd/4.0/).
Peer-review under responsibility of organizing committee of ICWRCOE 2015
doi:10.1016/j.aqpro.2015.02.172
S. Abdul Rahaman et al. / Aquatic Procedia 4 (2015) 1322 – 1330 1323
sensitive priority zones and is useful for better management practices such as implementation of land and water resource
management, conservation and sustainable agricultural development.
© 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license
© 2015 The Authors. Published by Elsevier B.V.
(http://creativecommons.org/licenses/by-nc-nd/4.0/).
Peer-review under
Peer-review under responsibility
responsibility of organizing
of organizing committee
committee of ICWRCOE
of ICWRCOE 2015 2015.
Keywords: Prioritization, Morphometric, Fuzzy Analytical Hierarchy Process, Geographic Information System
1. Introduction
A drainage basin or watershed is an extent or an area of land where surface water from rain, melting snow, or ice
converges to a single point at a lower elevation, usually the exit of the basin, where the waters join another water
body, such as a river, lake, reservoir, estuary, wetland, sea, or ocean. The water shed plays a dominant role in the
development of landforms and therefore, the study of drainage basin has a great significance in geomorphic studies.
A watershed is an ideal unit for management of Natural resources like land and water and for mitigation of the
impact of natural disasters for achieving sustainable development. The watershed management concept recognizes
the interrelationships among the linkages between uplands, low lands, land use, geomorphology, slope and soil.
Morphometry is the measurement and mathematical analysis of the configuration of the earth's surface, shape and
dimension of its landforms (Agarwal, 1998; Obi Reddy et al., 2002). A major emphasis in geomorphology over the
past several decades has been on the development of quantitative physiographic methods to describe the evolution
and behaviour of surface drainage networks (Horton, 1945; Leopold & Maddock, 1953). Morphometric analysis of a
watershed provides a quantitative description of the drainage system, which is an important aspect of the
characterization of watersheds (Strahler, 1964).
Pioneering work on the drainage basin morphometry has been carried out by Horton (1932, 1945), Miller (1953),
Smith (1950), Strahler (1964) and others. In India, some of the recent studies on morphometric analysis using
remote sensing technique were carried out by Natuiay (1994), Srivastava (1997), Nag (1998) and Srinivasa et al
(2004). In the earlier researches, prioritization of watershed was accomplished through different approaches to
instance soil erosion or sediment yield indexing (SYI), morphological characterization, socio-economic aspects, etc.
some other studies focused on soil erosion and SYI modelling aspects by classifying the erosion affected priority
areas (Suresh et al., 2004; Ratnam et al., 2005; Kalin and Hantush, 2009; Pandey et al., 2009; Niraula et al., 2011;
Pai et al., 2011). Land deterioration as well as land use change impacts were also measured for evaluation of
prospective zones of watersheds (Adinarayana, 2003; Deb and Talukdar, 2010; Kanth and Hassan, 2010; Javed et
al., 2011; Sarma and Saikia, 2011).
The remote sensing and GIS technique is a convenient method for morphometric analysis as the satellite images
provides a synoptic view of a large area and is very useful in the analysis of drainage basin morphometry. GIS and
remote sensing (RS) techniques are proved to be proficient tools for morphometric characterization of sub-
watersheds (Singh, 1994; Grohmann, 2004; Sreedevi et al., 2009; Aher et al., 2010; Rao et al., 2011). Mishra et al.
(2007) carried out prioritization of sub-watersheds through morphological characteristics by using Soil and Water
Assessment Tool (SWAT) model in the small multi-vegetated watershed of a sub-humid subtropical region in India.
In recent research on prioritization of watersheds using multi-criteria evaluation through fuzzy analytical
hierarchy process (Aher, P. D., 2013). Assessment of different vulnerability producing factors is the decision
making process associated with formation of system knowledge database which involves multiple criteria and
alternatives, which results in great degree of complexity. Therefore, in this research an attempt has been made for
prioritization of sub-watersheds through analysis of the natural drainage system that implements a novel approach
by investigating the fuzzy analytical hierarchy process (FAHP) to circumvent the complex information associated
with various morphological characteristics by accomplishing better accuracy in identification and prioritization of
sub-watersheds with earlier approaches.
1324 S. Abdul Rahaman et al. / Aquatic Procedia 4 (2015) 1322 – 1330
The Kallar watershed situated in Eastern slope of Western Ghats stretching from West to the East. It is located
between 11˚17’0” to 11˚31’0” N Latitude and 76˚ 39’ 0” to 77˚ 8’ 45” E Longitude and covers an area of 1281.24
Sq.Km. It comprises of three districts The Nilgiris, Coimbatore and Erode, covering 6 taluks (Coonoor, Kothagiri,
Udhagamandalam, Mettupalayam, Coimbatore north, and Sathyamangalam); it consists of 89 Revenue villages (Fig.
1. (a). The maximum and minimum elevation encountered in the watershed about 177m and 2615m above MSL
(Fig. 1. (b). About 50% of areas are mountains covered with diverse plant communities that form various types of
forest and other agricultural activities, especially Tea, coffee plantation, vegetables and orchards, which are
normally cultivated in the upper and the lower area. The climate of this area is temperate and salubrious for more
than half of the year. The average day temperature of the sub watershed is 16.15° C and the average rainfall is about
901.65mm. The winter is relatively cool. The maximum rainfall is received during the month of October and
November. The Kallar streams flow from Southwest to Northeast and it connects the Bhavani River, which finally
joins with Bhavanisagar Dam, it was built in the Northeastern part of the watershed. Which is primarily serves as
source of irrigation. The area represent by clay soil, loam soil and rock out creep on steep to narrow sloping
landform. Geologically, the area represented by charnockite and fissile hornblende-biotite gneiss. The present study
area falls under the Survey of India toposheets (1:50,000 scale) 58 A/11, 15, 16 and 58E/3 & 4. The watershed and
sub watershed boundary has been demarcated using this drainage network.
3. Methodology
Watershed is a natural hydrological unit, topographically delineated area drained by a stream system, from which
runoff resulting from precipitation flow past from a point into single stream. The drainage network was derived from
Survey of India toposheet, and the sub watershed boundaries were demarcated on the basis of contour and drainage
flow direction. The Kallar watershed was divided into ten sub watersheds named as SW1 to SW10 (Fig. 2. (a) and
(b)). The smallest (SW9) and the largest (SW1) sub watersheds measures 40.9 and 210 Sq.km, respectively. The
morphometric parameters broadly classified into three aspects, viz... linear, areal, and relief, the sub parameters such
as stream order, number of streams, stream length, bifurcation ratio, drainage density, stream frequency, drainage
texture, relief ratio, basin shape, form factor, circularity ratio, elongation ratio, and length of overland flow.
The naming of stream order is the first step in morphometric analysis of drainage basin, based on the hierarchic
of stream proposed by straheler (1964). In the study area SW9 and SW10 are of forth order, SW4, SW5, SW6, and
SW8 are of fifth order, SW1, SW2, and SW3 are of sixth order, whereas , SW7th have seventh order. The entire
Kallar watershed is falls under seventh order. The morphometric characterization in the form of linear, areal and
relif aspects for the delineated sub watershed was calculated based on given formula (Table 1 & 2).
Constant of Channel Maintenance(C) 0.37 0.32 0.39 0.29 0.38 0.47 0.41 0.38 0.4 0.47 0.38
Lemniscate’s (k) 3.56 3.47 3.15 3.31 3.3 3.27 3.32 3.3 2.85 3.53 4.56
Circularity Ratio (Rc) 0.47 0.36 0.42 0.41 0.39 0.57 0.54 0.39 0.48 0.27 0.36
Form Factor Ratio (Rf) 0.28 0.29 0.32 0.3 0.3 0.31 0.3 0.3 0.35 0.28 0.22
Elongation Ratio (Re) 0.29 0.32 0.48 0.39 0.4 0.41 0.39 0.4 0.74 0.30 0.10
Shape Factor Ratio (Rs) 3.56 3.47 3.15 3.31 3.3 3.27 3.32 3.3 2.85 3.53 4.56
Min Elevation (h) (M) 293 291 302 327 310 273 247 258 251 177 177
Max Elevation (H) (M) 2614 2242 812 1578 1499 450 1369 2116 581 1887 2615
Total Basin Relief (R) (m ) 2321 1951 510 1251 1189 177 1122 1858 330 1710.00 2438
Relief Ratio (Rr) 84.8 79.9 31.2 62.2 59.8 9.26 54.8 93.4 30.8 65.29 31.75
Relative Relief (RhP) 31.1 25 10.1 20.4 19.2 3.58 20.7 30 10.1 18.01 11.47
Ruggedness number (Rn) 6274 6159 1323 4329 3115 378 2747 4868 834 3644.19 6483.40
Gradient Ratio (Rg) 84.8 79.9 31.2 62.2 59.8 9.26 54.8 93.4 30.8 65.29 31.75
Melton Ruggedness Number (MRn) 160 149 55.3 113 109 16.7 99.9 170 52 122.59 67.80
Texture ratio (T) 7.65 7.42 4.13 7.05 4.12 3.74 4.52 4.12 3.14 3.83 15.25
In order to prepare the prioritize watershed, various methods such as quantitative, fuzzy logic, statistic methods
and AHP where used by several researches. The present study of Fuzzy Analytical Hierarchical Process (FAHP)
with extent analysis method was used to prioritize the watersheds. The Analytical Hierarchy Process (AHP), which
is a multiple criteria decision making tools, especially in the problems with spatial nature or GIS-based. In addition
AHP, its manner to apply and its weaknesses and strengths and ultimately the fuzzy modified Analytical Hierarchy
Process (FAHP) which is proposed after that the concepts like of fuzziness, uncertainty and vagueness was broadly
posed in expert's decision making. Fuzzy AHP is described to obtain a crisp priority vector from a triangular fuzzy
comparison matrix. In spite of popularity of AHP, this method is often criticized for its inability to adequately
handle the inherent uncertainty and imprecision associated with the mapping of the decision maker’s perception to
exact numbers (Deng, 1999).
To apply the process depending on this hierarchy, according to the method of extent analysis, each criterion is
taken and extent analysis for each criterion, gi; is performed on, respectively. Therefore, m extent analysis values for
each criterion can be obtained by using following notation.
M 1gi , M gi2 , M gi3 , M gi4 , M gi5 ,}}, M gin Where gi is the goal set (i = 1, 2, 3, 4, 5 … n) and all the M gij (j = 1, 2, 3, 4,
5…m) are Triangular Fuzzy Numbers (TFNs). The steps of Chang’s analysis can be given as in the following:
Step 1: The fuzzy synthetic extent value (Si) with respect to the ith criterion as follows:
S. Abdul Rahaman et al. / Aquatic Procedia 4 (2015) 1322 – 1330 1327
1
m ªn m º
Si ¦M j
gi «¦ ¦ M » j
gi (1)
j 1 ¬i 1 j 1 ¼
This involves computation of
m
¦M
j 1
j
gi (2)
Perform the “fuzzy addition operation” of m extent analysis values for a particular matrix given in (3) below, at
the end step of calculation, new (l, m, u) set is obtained and used for the next:
m m m m
¦M
j 1
j
gi (¦ l j , ¦ m j , ¦ u j )
j 1 j 1 j 1
(3)
Where “l” is the lower limit value, “m” is the most promising value and “u” is the upper limit value and to obtain
(4); this involves computation of
1
ª n m º
«¦¦ M gi »
j
(4)
¬i 1 j 1 ¼
Perform the “fuzzy addition operation” of M gij (j = 1, 2, 3, 4, 5…m) values given
n m m m m
¦¦ M
i 1 j 1
j
gi (¦ li , ¦ mi , ¦ ui )
j 1 j 1 j 1
(5)
and then compute the inverse of the vector in (5) . (6) Is then obtained such that:
ª º
1
ª n m º « 1 1 1 »
«¦¦ M gi » « n , n , n »
j
(6)
« u »
«¬ ¦ i ¦ mi ¦ li »
¬i 1 j 1 ¼
i 1 i 1 i 1 ¼
Step 2: The degree of possibility of M 2 (l 2 , m2 , n2 , ) t M 1 (l1 , m1 , n1 , ) is defined as equation (7):
and x and y are the values on the axis of membership function of each criterion. This expression can be equivalently
written as given in equation (8) below:
° 1, if
° m2 t m1 ,
° (8)
V (M 2 t M 1 ) ® 0, if m2 t m1,
°
° l1 u 2
° (
¯ 2 m u 2) ( m1 l1)
Step 3. The degree possibility for a convex fuzzy number to be greater than k convex fuzzy numbers Mi (i = 1, 2,
3, 4, 5 … k) can be defined by
V ( M t M 1 , M 2 , M 3 ,}} M k )
V [( M t M 1 ) and ( M t M 2 ) and ( M t M 3 ) and }} and ( M t M k ) (9)
min V ( M t M i ), i 1,2,3,4,} k
Assume that equation (9) is:
1328 S. Abdul Rahaman et al. / Aquatic Procedia 4 (2015) 1322 – 1330
d ( Ai ) minV ( S i t S k )
For k = 1, 2, 3, 4, 5…n; k≠ i. Then the weight vector is given by equation (10):
W1 (d ( A1 ), d ( A2 ), d ( A3 ), d ( A4 ), d ( A5 )}d ( An ))T (10)
Where Ai (i = 1, 2, 3, 4, 5, 6… n) are n elements.
Step 4. Via normalization, the normalized weight vectors are given in (11):
W (d ( A1 ), d ( A2 ), d ( A3 ), d ( A4 ), d ( A5 )}d ( An ))T (11)
Where W is non-fuzzy numbers.
The morphometric parameters i.e., bifurcation ratio (Rb), drainage density (Dd), stream frequency (Fs), drainage
texture (Dt), circularity ratio (Rc), relief ratio (Rr), ruggedness number (Rn) basin relief (Rh) and basin shape (Bs)
are also termed as erosion risk assessment parameters and have been used for prioritizing watersheds. The linear
parameters such as drainage density, stream frequency, bifurcation ratio, drainage texture, have a direct relationship
with erodibility, higher the value, more is the erodibility. Hence the prioritization of sub watershed, the highest
value was rated as score or rank based on Triangular Fuzzy Numbers (TFN) as very strong importance (VSI -
2,5/2,3), second highest value was rated as strongly importance (SI – 3/2,2,5/2), and so on, the least value was rated
as weakly more importance (WMI-1,3/2,2), and both parameters meet an same influence was rated as just equal (JE
– 1,1,1). Shape parameter such as circularity ratio, basin shape and compactness coefficient have an inverse
relationship with erodibility, lower the value of shape parameters was rated as very strong importance (VSI), nest
lower value was rated as ranked strongly importance (SI – 3/2,2,5/2), the highest value was rated as weakly more
importance (WMI-1,3/2,2), and both parameters meet an same influence was rated as just equal (JE – 1,1,1), sub
watershed wise morphometric variables used for pair-wise comparison matrix in fuzzy analytical hierarchy process
to generate the morphometric criterion weight. Same way each morphometric parameters compare with every sub
watershed (Table.3 & 4).
Finally we obtained criterion weight and alternative weight, for prioritization of sub watershed criterion weights
were multiply with alternative weight, and overlaid all sub watershed criterion in GIS environment to prepare a
prioritization of sub watershed. The final score ranges in Fuzzy Analytical Process is 0.074 to 0.167, where highest
rank assigned to the sub watershed having the highest analysis value, the same way ranks were assigned to each sub
watershed (Fig. 3. (a)) In this study SW10 was allotted for highest priority having FAHP value of 0.167 followed by
SW4, SW1, SW7, SW9, SW6, SW3, SW5, SW8, SW2 is the least priority in this watershed (Table. 5).
The prioritizes FAHP score was divided in to five priority classes from very less to very high based on the overall
weights(Table.6). The Kallar watershed categorized priority classes, very less priority is SW2, it covers an area of
13.15%, and very high priority class is SW4, SW10, it covers an area of 24.50%. (Fig. 3. (b)). The SW10 falls under
very high priority category, which is located near to the reservoir, it contain high amount of erosion. Hence the
SW10 need immediate concern to protect and control the erosion. The SW1 is falls under the high priority category,
which is located on high altitude area with huge amount of population and this sub watershed falls under the highly
landslide hazard zone, this sub watershed also need immediate action to sustain the environmental condition.
6. Conclusion
The present study demonstrate the utility of Geographical Information System (GIS) and Multi Criteria Decision
Making (MCDM) techniques in prioritizing sub watershed based on morphometric as well as integration of Fuzzy
Analytical Hieratical Process (FAHP) with extent analysis. In this study, a new and consistent approach of MCDM
processes i.e. FAHP analysis based prioritization was formulated successfully which plays a very important role in
illustrating the problem through integration of risk assessment factors causing natural resources. This may be one of
the feasible and efficient techniques, particularly in conventional watershed prioritization approaches for designing
and developing the efficient sustainable development and management practices. In this study, it has found that
SW4 and SW10 are falls in the high priority category based on morphometric analysis; hence, these sub watersheds
may be taken for conservation measurement by decision maker for planning and development.
References
Agarwal CS 1998. Study of drainage pattern through aerial data in Naugarh area of Varanasi district, U.P. Jour Indian Soc Remote Sensing 26:
169-175.
Strahler AN 1964. Quantitative geomorphology of drainage basins and channel networks: Handbook of applied hydrology. McGraw Hill Book
Company, New York, Section 4-76.
Horton, R. E. 1932. Drainage basin characteristics. Am.Geophys. Union, Trans. 13, 348–352.
Biswas S, Sudhakar S, Desai VR 1999. Prioritization of sub-watersheds based on Morphometric Analysis of Drainage Basin, District Midnapore,
West Bengal. J Indian Soc Remote Sensing 27: 155-166.
Nag SK 1998. Morphometric analysis using remote sensing techniques in the Chaka sub-basin, Purulia district, West Bengal. J Indian Soc
Remote Sensing 1: 69-76.
Singh, S. and Singh, M.C. 1997 Morphometric analysis of Kanhar river basin. National Geographical J. of lndia, (43), 1: 31-43.
M. Imran Malik, M. Sultan Bhat and Nissar A. Kuchay., 2011. Watershed based drainage Morphometric Analysis of Lidder catchment in
Kashmir Valley using Geographical Information System. Recent Research in Science and Technology 2011, 3(4): 118-126.
Adinarayana, J. 2003. Spatial decision support system for identifying priority sites for watershed management schemes. In Proc. First Interagency
Conference on Research in the Watersheds (ICRW), 405-408. Benson, Arizona, USA: U. S. Department of Agriculture, Agricultural
Research Service, 27-30 October.
Aher, P. D., J. Adinarayana, and S. D. Gorantiwar. 2013. Prioritization of watersheds using multi-criteria evaluation through fuzzy analytical
hierarchy process. Agric Eng Int: CIGR Journal, 15(1): 1118.
Chang, D. Y. 1996. Applications of the extent analysis method on fuzzy AHP. European Journal of Operational Research, 95(3): 649-655.
Saaty, T. L. 1980. The analytic hierarchy process. New York, U.S.: Mc Graw Hill.
Grohmann, C. H. 2004. Morphometric analysis in Geographic Information Systems: applications of free software GRASS and R star. Computers
and Geosciences, 30(10): 1055-1067.
Hlaing, K., S. Haruyama, and A. Maung. 2008. Using GIS-based distributed soil loss modeling and morphometric analysis to prioritize watershed
for soil conservation in Bago river basin of Lower Myanmar. Frontiers of Earth Science in China, 2(4): 465-478.
Javed, A., M. Y. Khanday, and R. Ahmed. 2011. Watershed prioritization using morphometric and land use/land cover parameters: a remote
sensing and GIS based approach. Journal Geological Society of India, 78(1): 63-75.
Mishra, A., S. Kar, and V. P. Singh. 2007. Prioritizing structural management by quantifying the effect of land use and land cover on watershed
runoff and sediment yield. Water Resources Management, 21(11): 1899-1913.
Kanth, K. A., and Z. U. Hassan. 2010. Prioritization of watersheds in Wular catchment for sustainable development and management of natural
resources. Recent Research in Science and Technology, 2(4): 12-16.
Pandey, A., V. M. Chowdary, B. C. Mal, and M. Billib. 2009. Application of the WEPP model for prioritization and evaluation of best
management practices in an Indian watershed. Hydrological Processes, 23(21): 2997-3005.
Sarma, S., and T. Saikia. 2011. Prioritization of Sub-watersheds in Khanapara–Bornihat area of Assam–Meghalaya (India) based on land use and
slope analysis using remote sensing and GIS. Journal of the Indian Society of Remote Sensing, 40(3): 435-446.
Srivastava, R. K., H. C. Sharma, and A. K. Raina. 2010. Suitability of soil and water conservation measures for watershed management using
geographical information system. Journal of Soil and Water Conservation, 9(3): 148-153.
Ratnam, N. K., Y. K. Srivastava, V. V. Rao, E. Amminedu, and K.S. R. Murthy. 2005. Check dam positioning by prioritization micro-watersheds
using SYI model and morphometric analysis – remote sensing and GIS perspective. Journal of the Indian Society of Remote Sensing, 33(1):
25-38.
Thakkar, A. K., and S. D. Dhiman. 2007. Morphometric analysis and prioritization of mini-watersheds in a Mohr watershed, Gujarat using
remote sensing and GIS techniques. Journal of the Indian Society of Remote Sensing, 35(4): 313-321.
Singh, N. 1994. Remote sensing in the evaluation of morphohydrological characteristics of the drainage basin of Jojri catchment. Annals of Arid
Zone, 33(4): 273-278.
Vittala, S., Govindaiah, S. and Honne, G. H. 2004. Morphometric analysis of sub-watersheds in the Pavagada area of Tumkur district, South
India using remote sensing and GIS techniques Jour. Indian Soc. Remote Sensing (32), 4: 351-362.
Schumn, S. A. 1956. Evolution of drainage systems and slopes in badlands at Perth Amboy, New Jersey. Geol. Soc. Am. Bull. 67: 597-646.
Oguchi, T. 1997. Drainage density and relative relief in humid steep mountains with frequent slope failure. Earth Surf. Processes Landforms,
22:107–120.