Morphometric characteristics are the dimension and mathematical evaluation of the configuration of the earth’s surface [1]. Moreover, the morphometric characteristics may control imperative information concerning its formation and growths due to all hydrologic and geomorphic processes occur within the watershed. Morphometric techniques have been applied to a variety of earth surface features and proven to be very effective in planning and management of land surface processes [2,3]. There have been few quantitative investigations of the Lateritic upland region of West Bengal. While, more attention has been given to the morphometric and spatial distribution parameters of its surroundings. Furthermore, recognition of drainage networks within the watershed region can be obtained using conventional methods like field observations and topographic maps or instead with remote sensing and GIS technology [3,4]. However, due to absence of topographic maps in adequate scale, or the availability of data, it is difficult to observe all drainage networks from field observations by conventional methods in rough terrain and inaccessible areas.

The Remote sensing and GIS technology are a convenient method to study the morphometric characteristics as the satellite images provides detailed information of earth surface features with its synoptic coverage, high receptivity, cost effectiveness [5]. Furthermore, Digital elevation models (DEMs) are progressively more used for visual and mathematical analysis of topography, landscapes and landforms, as well as modeling the surface processes [6-8]. Actually, the automatic generation of a DEM from remotely sensed data with sub-pixel precision is promising [9] Application of GIS to spatial data has proven to be instrumental in the analysis of complex problems in the earth and environmental sciences.

The varying nature of surface characteristics due to natural hazard (e.g. agricultural drought) is a persistent phenomenon in the western part of West Bengal during the last two-three [10,11]. Due to uneven distribution of the surface properties, soil moisture and inadequate rainfall played an important role for the growth of healthy crop and cause extreme crop stress and wilt. The study of morphometric characteristics of the lateritic upland area will reveal an apparent idea about the spatial distribution of relief and hydrological characteristics and will provide the knowledge to devise and accomplish an appropriate plan to progress agriculture and others allied activities. Hence, in the present study we investigated the morphometric characteristics of Kangshabati-Darkeswar Interfluves Area using remote sensing and GIS technology.

The Kangshabati-Darkeswar Interfluves Area is extended between 24°04’02.76” N- 21°44’45.33” N latitude and 85°34’27.02” E - 87°59’46.05” E longitude (Figure 1). The Interfluves area is bounded in the Hazaribagh in North, Baripada in the South, Bardwan in the East and Ranchi in the West. It is surrounded by lateritic soil, older and younger alluvial soil, red gravel soil and red sandy soil. Geological properties of the study area showed fine and medium sands, unconsolidated sands, silt and clay, fragment of pebbles and boulder, granite gneiss, quartzite and mica schist. Geomorphologically, the study area is characterized by floodplains, upland plains, badlands, duricrust, paradeltaic fan surfaces, pediments, pediplains, ridges and hills [12].

Figure 1: Location map of the Kangsawati-Darkeswar Interfluve area in West Bengal, India.

The morphometric analysis of the study area was carried out using the Advance Space Thermal Emission and Radiometer (ASTER, Spatial resolution-30 m). ASTER GDEM data was downloaded from the ASTER GDEM Project of Japan-US ASTER Science Team (http://gdem.ersdac.jspacesystems.or.jp/). For the present study a total of 16 ASTER GDEM scenes (i.e., ASTGTM2 - N21E086, N21E087, N22E085, N22E086, N22E087, N22E088, N22E089, N22E090, N23E085, N23E088, N23E089, N23E090, N24E087, N24E088 and N24E089). 10 km² grids were generated to delineate the relief characteristics. A total of 363 grids were generated using GIS to extract the morphometric characteristics of the study sites. For each grid maximum, minimum and mean elevation was obtained separately. According to [13], the maximum rate of change in value from each cell to its neighbors is identified as slope in grid areas. The digital Elevation Model (DEM) was computed to deduce the morphometric measurement like, absolute relief, relative relief [14], dissection index [15], ruggedness index [16,17], drainage characteristics (e.g., drainage density, drainage frequency) [18], infiltration number [19], overland flow [20]. Drainage networks were automatically extracted from digital elevation models coupled with the constant stream threshold value based on the method proposed by [21].

Relative relief was derived as the difference in height between the highest and the lowest points (height) in a unit area. The study area has been categorized into ten relative relief characteristics based on the geometric interval methods. Dissection index is calculated, as the ratio between relative relief and absolute relief. Ruggedness index was computed to evaluate the complexity of topography and the roughness of the terrain. Frequency of drainage can be defined as the number of drainage segments per 10 km² grid area. Drainage density (DD) refers as the ratio of total length of streams of all orders with in the basin to the basin area or per unit area, is expressed in terms of km/ sq. km [18]. Drainage texture ratio (RT) refers to the relative spacing of drainage lines. Stream ordering is defined as a measure of the position of a stream in the hierarchy of tributaries [22]. The length of overland flow approximately equals to half of the reciprocal of drainage density [20]. According to Smith’s classification of drainage texture, the texture value below 4 is designated as coarse; 4-10 as intermediate; above 10 as fine and above 15 as ultra-fine texture. Infiltration number is the product of drainage density and stream frequency which provides an idea of the infiltration characteristics of the basin. Drainage intensity refers to the ratio of the stream frequency to the drainage density [19].

Second order local polynomial (LP) interpolation technique was used to estimate the surface characteristics of the study area using ArcGIS v 9.3. The advantage of this technique is that it is a quite quick deterministic interpolator that is smooth and more flexible than the global polynomial method [23]. As such, there is no measurement of prediction errors and this method does not let you to examine the autocorrelation of the data, making it less flexible and more automatic than Kriging techniques.

Absolute relief

The descriptive characteristics of each morphometric analysis of the study area are represented in Table 1. The absolute elevation of the study region is between 4.0-949.29 m with an average elevation of 484.50 m (S.D. ± 277.85). The Maximum elevation is recorded in the extreme North- eastern corner and Eastern part of the study site and the minimum elevation is recorded in the Western and Southwestern part of the study site (Figure 2). The absolute mountainous relief of the study area range between 8.61-589.6 m (mean ± standard deviation 141.59 m ± 112.17). The average elevation value is recorded in the central and Northern part of the study site.

Figure 2: Absolute relief features of the Kangsawati-Darkeswar Interfluves Area.

Relative relief

Relative relief is an important morphometric variable used for the overall assessment of morphological characteristics of terrain. The highest relative relief is calculated as 833.69 m, while the lowest value is recorded as 20 m (Figure 3). The low relative relied demonstrated that the region is almost flat and appear like mature stage of geomorphic evolution. Therefore, the region could be used for agricultural purpose for its flat nature depending on water accessibility.

Figure 3: Relative relief characteristics of the Kangshabati-Darkeswar Interfluves Area.

Dissection index

The highest computed dissection index value was 0.99, while the lowest value is 0.20 (Figure 4). Moreover, the higher dissection index in the study area highlight that both a low river erosion and a transition of the area toward a mature stage of development.

Figure 4: Dissection index characteristics of the Kangshabati-Darkeswar Interfluves Area.

Slope

The slope of a terrain refers to the amount of inclination of physical feature, topographic landform to the horizontal surface. Slope analysis is an important parameter in morphometric studies. The slope elements, in turn are controlled by climato-morphogenic processes in areas having rock of varying resistance. The degree of slope exhibited by Kangshabati-Darkeswar interfluves area varies from 0º to 65º. The average slope is 26.53 º (Table 1). The slope map of the Kangshabati- Darkeswar interfluves area (Figure 5). The higher slope gradient in the study area belongs to Chotanagpur hilly Eastern part. Higher slope gradient will result in rapid runoff with potential soil loss.

Figure 5: Slope map of the Kangshabati-Darkeswar Interfluves Area.

Stream order

The Kangshabati-Darkeswar stream interfluves area have been demarcated according to the Strahler’s system of stream ordering has been designated as a fifth-order basin (Table 2, Figure 6) the present investigation, maximum frequency is observed in the firstorder streams. The first-order streams are those that do not have any tributary and these channels normally flow during the wet weather [24]. More number of first-order streams can be observed in the hilly region of the study area, which points towards terrain complexity and compact nature of the bedrock lithology. Moreover, the presence of large number of streams in the basin indicates that the topography is still undergoing erosion and at the same time, less number of streams indicates mature topography. The calculated result matched with [16,17], which described that the total number of streams gradually decreases as the stream order increases.

Figure 6: Stream Ordering of the Kangshabati-Darkeswar Interfluves Area.

Drainage frequency

Maximum drainage frequency of the study area is recorded as 2.78/10 km² grid area. Very small pockets of high drainage frequency areas were found in the southern side of the study area (Figure 7). A higher stream frequency points to a larger surface runoff and steeper ground surface. According to [20], stream frequency (Fs) is defined as the ratio of the total number of stream segments of all the orders in the basin to the total area of the basin. ‘Fs’ is an index of the various stages of landscape evolution. The occurrence of stream segments depends on the nature and structure of rocks, vegetation cover, nature and amount of rainfall and soil permeability. The stream frequency for the Kangshabati-Darkeswar interfluves area is 0.62 (Table 1). The stream frequency is dependant more or less on the rainfall and the physiography of the region.

Figure 7: Drainage frequency/10km2 grid area of the Kangshabati-Darkeswar Interfluves Area.

Table 1: Descriptive characteristics of morphometric parameters in the Kangsawati-Darkeswar Interfluve area in West Bengal, India.

Stream length ratio (RI)

It is the ratio between the lengths of streams in a given order to the total length of streams in the next order [18]. The Rl values for the Kangshabati-Darkeswar interfluves area vary widely from 0.5 to 0.9 (Table 2) and are strongly dependant on the topography and the slope. It shows an important relationship with the surface flow discharge and the erosional stage of the basin or interfluves area.

Table 2: Results of morphometric analysis of Kangshabati-Darkeswar interfluves area.

Bifurcation ratio (Rb)

Rb is related to the branching pattern of a drainage network and is defined as the ratio of the number of streams of any given order to the number of streams in the next higher order in a drainage basin [14]. It is a dimensionless property and shows the degree of integration prevailing between streams of various orders in a drainage basin. Rb shows a small range of variation for different regions or for different environments except those where the powerful geological control dominates. The Rb for the Kangshabati-Darkeswar interfluves area varies from 1.1 to 2.1 (Table 2) with a mean Rb of 1.67. Low Rb value indicates poor structural disturbance and the drainage patterns have not been distorted [16,17], whereas the high Rb value indicates high structural complexity and low permeability of the terrain. A low mean Rb value (0.64) indicates less structural disturbances in Kangshabati- Darkeswar interfluves area.

Drainage density

Drainage density (Dd) is one of the important indicators of the landform element and provides a numerical measurement of landscape dissection and runoff potential [25]. Dd is related to various features of landscape dissection such as valley density, channel head source area, relief [26], climate and vegetation [27], soil and rock properties [28,29] and landscape evolution processes. A low drainage density indicates permeable sub-surface strata and has a characteristic feature of coarse drainage, which generally shows values less than 5.0. [16,17] noted that low drainage density is favored where basin relief is low and vice versa. In this study highest drainage density is recorded as 4.55/sq. km and the average drainage density is 0.73/sq. km with a standard deviation of ± 0.58 which indicates that the study area has a weak or permeable subsurface material with intermediate drainage and low relief. Highest DD is recorded in the extreme northeastern and southern part of the study area. Some small pockets of highest drainage density are recorded in the central part also. Low drainage densities are often associated with widely spaced streams due to the presence of less resistant materials. Moreover, the lowest DD is recorded in the eastern and northern part of the study site (Figure 8). The result of our analysis showed that the area may have significantly moderate to good amount of surface water resource. However, it may be due to the origin of some first order streams from the highlands and it does not add significantly to the accessibility of water in the region.

Figure 8: Drainage density/10km2 grid area of the Kangshabati-Darkeswar Interfluves Area.

Ruggedness index

The average ruggedness index of the study area is recorded as 0.09 per 10 km² area. Maximum portion of the study area is estimated as lower ruggedness index, comprised with the central and southern part (Figure 9). The results of the study indicated that the minimum effort is needed to manage the surface water and irrigational facilities which can make the good agricultural region.

Figure 9: Ruggedness index characteristics of the Kangshabati-Darkeswar Interfluves Area.

Drainage texture

According to [28], drainage texture (T) is a product of stream frequency and drainage density (Table 1). The ‘T’ index depends on underlying lithology, infiltration capacity and relief aspect of the terrain. The drainage texture value of the Kangshabati-Darkeswar interfluves area basin is 0.48, and indicates coarser drainage texture and due to the presence of high relief in the western part of the study area (Table 2). Most of the investigated area showed coarser drainage texture per 10 km² grid area (Figure 10).

Figure 10: Drainage texture/10 km2 grid area of the Kangshabati-Darkeswar Interfluves Area.

Length of overland flow

This depends on the rock type, permeability, climatic regime, vegetation cover and relief as well as duration of erosion [14]. The high overland flow value indicates that the rainwater had to travel relatively longer distance before getting concentrated into stream channels [30,31]. The average length of overland flow of the study area is shown as 2.56/10 km², it may be due to under the influence of high structural disturbance, low permeability, steep to very steep slopes and high surface runoff. The maximum value of overland flow is recorded in northern and some small pockets of eastern corner of the study site (Figure 11). However, the low value of overland flow concentration is recorded from the southwestern parts of the study site.

Figure 11: Length of overland flow/10km2 grid area of the Kangshabati-Darkeswar Interfluves Area.

Infiltration number

The higher value indicates low infiltration and high runoff. The infiltration number of the study area is the highest infiltration number 3.64/10 km² grid area in the study site (Figure 12). The average infiltration number of the study site is estimated at 0.43 with a standard deviation of ± 0.44.

Figure 12: Infiltration number/10km² grid area of the Kangshabati-Darkeswar Interfluves Area.

Remote sensing and GIS application employed in the demarcation and modernize of the relief and drainage characteristics in the present study. The variation in drainage characteristics might be due to changes in slope and topography. It is contingent from the study that the streams have created a mature stage which is a sign of good geomorphic development. Furthermore, by the complete analysis of drainage basin characteristics of Kangsabati-Darkeswar interfluves area can be considered as an area with moderate to poor groundwater prospect as the area has less permeable subsurface and the conditions somewhat favorable for infiltration of surface water. However, it is also required to observe the other factors like, land use, climate, soil type, etc. which are identified to influence hydrologic processes to disentangle the multivariate nature of the predicament with the view to finding a holistic elucidation.

We are very much thankful to the Ministry of Economy, Trade and Industry of Japan (METI) and the National Aeronautics and Space Administration (NASA) for freely providing the ASTER GDEM data. We also thanks to the Department of Geography, Vidyasagar University (Midnapur) for their laboratory facilities to carry out the research work.

1. Pareta K, Pareta U (2011) Quantitative Morphometric Analysis of a Watershed of Yamuna Basin, India using ASTER (DEM) Data and GIS. International Journal of Geomatics and Geosciences 2: 248 – 269.
2. Bhunia GS, Bramhya S, Mishra AK (2013) Spatial analysis of relief characteristics in Kondagaon District (Chattishgarh, India) – A study based on remote sensing and Geographic Information system. International Journal of Physical and Social Sciences 3: 145 – 157.
3. Bhunia GS, Samanta S, Pal B. Quantitative analysis of relief characteristics using space technology. International Journal of Physical and Social Sciences 2: 350-365.
4. Macka, Z (2001) Determination of texture of topography from large scale contour maps. Geografski Vestnik 73: 53–62.
5. Sreedevi PD, Owais S, Khan HH, Ahmed S (2009) Morphometric Analysis of a Watershed of South India Using SRTM Data and GIS. J Geol Soc 73: 543-552.
6. Kamp U, Bolch T, Olsenholler J (2005) Geomorphometry of Cerro Sillajhuay (Andes, Chile/Bolivia): Comparison of Digital Elevation Models (DEMs) from ASTER Remote Sensing Data and Contour Maps. Geocarto International 20: 23-33.
7. Bishop MP, Shroder JF (2000) Remote sensing and geomorphometric assessment of topographic complexity and erosion dynamics in the Nanga Parbat massif. USA.
8. Treloar PJ, Searle MP, Khan MA, Jan MQ Tectonics of the Nanga Parbat Syntaxis and the Western Himalaya: an introduction. Geological Society London, Special Publications, London 170: 1-6.
9. Tucker GE, Catani F, Rinaldo A, Bras RL (2001) Statistical analysis of drainage density from digital terrain data, Geomorphology 36: 187-202.
10. Krzystek, P., (1995). New investigations into the practical performance of automatic DEM generation, Proceedings, ACSM/ASPRS Annual Convention, Charlotte, North Carolina, American Society for Photogrammetry and Remote Sensing 2: 488-500.
11. Bhunia GS, Samanta S, Pal DK, Pal B (2012) Assessment of Groundwater Potential Zone in Paschim Medinipur District, West Bengal – A Meso-scale study using GIS and Remote Sensing Approach. Journal of Environment and Earth Science 2: 41-59.
12. Maity SK, Maiti R (2013) Hydrodynamics at the Junction of Silabati, Dwarakeswar and Rupnarayan Rivers at Bandar, Paschim Medinipur, West Bengal, India. Earth Science India 6: 77-89.
13. Singh LP, Parkash B, Singhvi AK (1998) Evolution of the Lower Gangetic Plain landforms and soils in West Bengal, India. Catena 33: 75–104.
14. Burrough PA (1986) Principles of geographical information systems for land resources assessment. Geocarto International 1: 54.
15. Schumm SA (1956) Evolution of Drainage Systems and Slopes in Badlands at Perth Amboy, New Jersey. Geological Society of America Bulletin 67: 597-646.
16. Singh S, Dubey A (1994 ) Geo -environmental Planning of Watershed in India, Chugh Publications, Allahabad, India.
17. Strahler AN (1964). Quantitative geomorphology of drainage basins and channel networks.
18. Chow VT (2010) Handbook of applied hydrology. International Association of Scientific Hydrology. Bulletin 10: 82-83.
19. Horton RE (1932) Drainage Basin Characteristics. Transactions, American Geophysical Union, 13: 350-361.
20. Faniran A (1968) The Index of Drainage Intensity - A Provisional New Drainage Factor. Australian Journal of Science 31: 328-330.
21. Horton RE (1945) Erosional Development of Streams and their Drainage Basins; Hydrophysical Approach to Quantitative Morphology. Geological Society of America Bulletin 56: 275-370.
22. Lin WT, Chou WC, Lin CY, Huang PH, Tsai JS (2005) Automated suitable drainage network extraction from digital elevation models in Taiwan’s upstream watersheds. Hydrological Processes, 20: 289-306.
23. Leopold LB, Wolman MG, Miller JP (1964). Fluvial processes in geomorphology. WH Freeman and Company, San Francisco and London.
24. Schaum A (2008) Principles of local polynomial interpolation. Proceeding AIPR '08 Proceedings of the 2008 37th IEEE Applied Imagery Pattern Recognition Workshop. IEEE Computer Society Washington, DC, USA.
25. Chow VT, David RM, Larry WM (1988) Handbook of applied hydrology. TATA McGraw Hill Inc., New York.
26. Chorley RJ (1969) Introduction to physical hydrology. Methuen and Co. Ltd.
27. Montgomery DR, Dietrich WE (1989) Source areas, drainage density and channel initiation. Water Resour Res 25:1907–1918.
28. Moglen GE, Eltahir EAE, Bras RL (1998) On the sensitivity of drainage density to climate change. Water Resour Res 34: 855–862.
29. Smith KG (1958) Erosional processes and landforms in badlands national monument, South Dakota. Geol Soc Am Bull 69: 975–1008.
30. Kelson KI, Wells SG (1989) Geologic influences on fluvial hydrology and bedload transport in small mountainous watersheds, northern New Mexico, USA. Earth Surf Processes 14: 671–690.
31. Chitra C, Alaguraja P, Ganeshkumari K, Yuvaraj D, Manivel M (2011) Watershed characteristics of Kundah sub basin using remote sensing and GIS techniques. Int J Geomatics Geosci 2: 311–335.