Estimation of Sediment Yield of Govindsagar Catchment, Lalitpur District, (U.P.), India, Using Remote Sensing and GIS Techniques ()
1. Introduction
Soil erosion is universally recognized as a serious threat to human being [1] . It is a widespread problem threating the sustainability of agriculture productivity and causing deterioration of both land and water resources. The severity of this problem is more pronounced in arid and semi-arid regions especially on grazing land, where high rainfall intensities of short duration, susceptibility of soil to erosion and human mismanagement of land have accelerated soil losses through erosion leading to downstream sedimentation. Catchments and watersheds have been identified as planning units for administrative purpose to conserve land and water resources [2] - [5] . The concept of watershed management recognizes the inter-relationships between land use, soil and water and the linkage between uplands and downstream areas [6] . Soil and water conservation are key issues in watershed management in India for demarcating the priority watersheds. The average soil loss rate has been computed as 16.75 t・ha−1・year−1 which translates into approximately 1 mm per year [4] . Besides degrading the land, erosion has direct bearing on the quantity and quality of water as well, since water carries with it loose soil that later gets deposited in the reservoirs and dams reducing both their storage capacity and life span. Compared to humid regions, the problem of soil and water erosion is altogether different in arid regions where flash floods are a common phenomenon [7] - [9] . Sediment yield from any area is a resultant of the interaction of meteorological factors with the land surface. Rainfall is the most important single meteorological factor that through in arid areas wind power assumes equal significance. The land factors include physiography, slope, soil, land use, vegetation and present status of erosion. Soil factors comprise the broad nature of the soil, effective depth, texture of the surfaces oil, soil reaction and stone content. It is the combined and reciprocal influence of these factors that determine the magnitude of the sediment yield from any locality.
Numerous attempts have been made to explain global and regional patterns of sediment yield in terms of climate and topography [10] - [15] . It has been estimated that in India about 113.3 million hectares of land are subjected to soil erosion due to water and about 5334 million tons of soil are being detached annually due to various reasons [16] . Since a host of variables are involved in erosion, it becomes difficult to measure or predict erosion accurately. The latest advances in remote sensing technology and geographical information sciences (GIS) provide real-time information on various aspects of the watershed such as land use, physiography, soil, relief, drainage characteristics, etc. It also assists in identification of existing or potential erosion-prone areas and provides data inputs to many of the soil erosion and runoff models. The rate of soil loss is judged by the Sediment Yield Index (SYI), which can be derived through various empirical formulae. In India, SYI model is developed by the All India Soil and Land Use Survey (AISLUS), and is commonly employed in providing criteria for priority delineation in river valley projects and flood-prone rivers. Several attempts have been made to estimate sediment yield in the recent past in various basins of India, using SYI and other conventional methods [17] - [23] .
The present study makes an attempt to assess the sediment yield rate of the Govindsagar Catchment by using standard method with the aid of remote sensing data and drainage parameters.
2. Study Area
Govindsagar catchment falls in Lallitpur district which is the southernmost part of Uttar Pradesh (India) and lies between 24˚24'N to 24˚43'N latitude and 78˚17'E to 78˚28'E longitude, covering an area of about 384.35 Km2 (Figure 1). The area falls in Survey of India (SOI) Toposheet numbers 54L/6 and 54L/7 on 1:50,000 scale. The maximum and minimum elevations in the catchment are 531 m and 336 m above mean sea level (amsl) respectively. The main Shahzad river and its tributaries form dendritic to sub dendritic drainage pattern, representing more or less gentle slope and plain topography. Geologically, the study area comprises of granites, gnisses, schist, sandstone, limestone, shale and alluvium. The soil in the catchment comprises of sand, silt, clay, red soil and laterites. Govindsagar reservoir is formed due to the dam constructed on
Shahzad river near Lalitpur town. The catchment is predominated by agricultural land, however, forest cover is present in the southern part. The main crops grown are wheat, barley, maize, chana, mustard etc. which are irrigated through canals, tube wells and surface streams. Though the catchment receives average annual rainfall of more than 1044 mm, but major part of the catchment faces acute shortage of water for drinking and irrigation purpose during summer. Socio-economically, this region known as “Bundelkhand”, has remained backward educationally, infrastructurally and growth in human development. Though the Union government had announced special economic package for the region few years back, however, the fruits of development are still awaited. The area is typically known for recurring droughts and scarcity of water during the last couple of decades, and in dire need of overall development.
3. Data Used and Methodology
Survey of India (SOI) Toposheet numbers 54L/6 and 54L/7 on 1:50,000 scale, were obtained from Survey of India, Dehradun, and used for the base map preparation. The toposheets were scanned and georeferenced in Erdas Imagine 14 software with the help of Ground Control points (GCPs) and were projected into Universal Transverse Mercator (UTM) projection, taking World Geodetic System (WGS84) as the datum. The study area lies in UTM zone 44 N. For geo-referencing, the latitudes and longitudes of a reference map were converted into X-Y co-ordinates, and the corresponding X-Y co- ordinates were transferred on the map. After georeferencing, these toposheets were mosaiced and watershed boundary was delineated from Advanced Space borne Thermal Emission and Reflection Radiometer (ASTER) data of 30 meter resolution, which is freely available and downloaded from http://www.earthexplorer.usgs.gov, using SAGA software, and this data was subsequently utilized for the preparation of Digital Elevation Model (DEM).
IRS-P6 LISS III data, having a spatial resolution of 23.5 meter, (Path/Row: 98/54) of 27th February, 2005 was procured from National Remote Sensing Centre (NRSC), Hyderabad (Figure 2).
The other relevant information was obtained from various government sources in form of maps, reports and ancillary data set.
The drainage network was created manually by digitizing drainage lines and overlaid on DEM in ArcGIS 10 (Figure 3). The drainage network of the basin and the stream ordering and morphometric parameters were computed using standard methods as adopted by Horton and Strahler [24] [25] . Standard visual image interpretation method based on photo-recognition elements such as tone, texture, size, shape, pattern, association and field knowledge was followed. Land use/land cover categories such as agriculture land, dense forest, open forest, open scrub, settlement, stone quarry, exposed rock, waste land and water body, were delineated on the basis of image interpretation, which was supplemented by ground truth verification. Land use/cover map of 2005 was imported to ArcGIS for digitization, editing and topology building. Polygon ids were assigned for each Land use/cover category and area both in km2 as well as in percentage
![]()
Figure 2. IRS-P6 LISS III standard FCC, 2005.
was computed in GIS.
In order to understand slope characteristics of the watershed, slope map was derived from DEM. Employing the standard procedure for calculating the slope in degrees, three classes as gentle (0˚ - 13˚), moderate (13˚ - 25˚), and steep (25˚ - 38˚) were assigned (Figure 4).
Sediment Yield Index (SYI) was computed taking inputs from drainage parameters, land use/cover and slope of area to assess sediment yield rate, in Govindsagar catchment.
4. Results and Discussion
4.1. Morphometric Analysis
Morphometric analysis provides quantitative description of basin geometry to understand initial slope or inequalities in the rock hardness, structural controls, recent diastrophism, geological and geomorphic history of a drainage basin [25] . The drainage in a
![]()
Figure 3. Drainage network superimposed on DEM.
basin is largely governed by the rocks, soil and rock structures, and can be used to draw inferences about lithology, structure and soil types. The higher order streams are often controlled by the rock structures, while the behavior of the lower order streams and their inter-relationships provide information about nature of rocks and soils [26] . The parameters utilized for the assessment of morphometric analyses are categorized into Linear, Areal and Relief aspects. The linear aspects include stream order, stream length, mean stream length, stream length ratio, and bifurcation ratio, while areal aspects include drainage density, steam frequency, drainage texture, basin shape, from factor, circularity ratio and elongation ratio, whereas relief aspects consists of relief ratio.The designation of stream order is based on a hierarchic ranking of streams, and is a useful indicator of stream size, discharge and drainage area [27] . Measurement of these linear, areal and relief parameters envisages to understand basin morphology [28] .
Stream ordering is the first step of quantitative analysis of the drainage basin as proposed by Horton [24] , and later modified by Strahler [25] . The study area has 720
various orders of streams and the main Shahzad river is of the 6th order stream. These streams contribute their discharge to Govingsagar reservoir. The main Shahzad river is of the 6th order (Table 1). The stream length ratio of different order in the Govindsagar catchment area reveals variation in slope and topographic conditions. Accordingly, stream length ratio has an important relationship with the surface flow discharge and erosional stage of the basin [29] . The bifurcation ratio values of Govindsagar catchment area range from 2.00 to 5.23, which indicate less structural control on the drainage development i.e. the drainage pattern of the study area has not been distorted because of structural disturbances.Drainage density is the total length of streams of all orders divided by the area of drainage basin and indicates closeness of spacing of channels [24] . The catchment has a drainage density of 4.87 Km−1 falling in high drainage density category which may be attributed to impermeable subsurface materials. The drainage texture value is 4.53, which indicates moderate texture. Basin shape is 2.54, which indicates peaked flood discharge and high stream flow carrying sediment load. Form
![]()
Table 1. Result of the morphometric analysis of Govindsagar catchment.
factor value of 0.39, indicates that watershed is less elongated in shape with flat peak flow for longer duration. Circularity ratio (0.19) suggests catchment area is less circular with low to moderate relief (6.26) [30] .
4.2. Land Use/Cover Analysis
Land use/cover analysis was carried out using IRS-P6 LISS III FCC data. Visual interpretation of FCC data led to the identification and delineation of various land use/cover categories such as agriculture land, dense forest, open forest, open scrub, settlement, stone quarry, exposed rock, waste land and water body (Figure 5). The area statistics of Land use/cover compited in GIS suggests that the catchment occupies an area of about 192.32 km2 under agriculture land, which is 50.08% of the total watershed area. Dense forest occupies an area of about 68.45 km2, which is largely confined to the southern periphery of the watershed. The area under open forest is about 15.89 km2, whereas open scruboccupies an area of about 91.21 km2. The other land use/cover categories such as settlement, exposed rock, stone quarry, waste land and water body occupy areas of 2.86 km2, 0.42 km2, 4.75 km2, 0.71 km2 and 7.71 km2 respectively (Table 2).
4.3. Sediment Yield Index (SYI) Analysis
A number of sediment yield models, both empirical and conceptual, are in practice to address wide ranging soil and water management problems. Most conservation planning for erosion control, however, uses empirical model(s) to estimate average annual soil loss. However such empirical models require input parameter in terms of spatial
![]()
Figure 5. Land use/cover map based on IRS-P6 LISS III.
![]()
Table 2. Land use/cover statistics.
information of land use, vegetation cover, soils, slope, drainage density, besides runoff and rainfall intensity. Soil and water conservation over large areas is difficult, expensive and unmanageable, hence requires a selective approach to demarcate smaller hydrological units watershed/sub-watersheds, for more efficient and targeted resource management programs is best suited for conservation planning [22] . The quantitative estimation of Sediment Yield Index (SYI) were carried out by Bali and Karale [17] and updated by National Remote Sensing Agency [31] . A common empirical model for SYI employed under Indian conditions, was carried out following Kumar [32] besides Rao and Mahabaleswara [21] , and is as follows:
(1)
Vs = Sediment Yield, (Mm3/year × 103)
P = Annual precipitation, (cm)
A = Catchment area, (km2)
D = Drainage density, (km/km2)
S = Average slope of watershed, (degrees)
F = Vegetative cover factor, (km2)
where, F can be determined as:
(2)
F1 = Agricultural land
F2 = Dense forest
F3 = Open forest
F4 = Open scrub
F5 = Exposed rock
F6 = Stone query
F7 = Wasteland
F8 = Settlement
F9 = Waterbody
The parameters like A, D, S and F are essential units, which can be derived from stream drainage map, slope map and land use/cover map. However, the parameter F needs to be redefined on land use/cover information that has been extracted from satellite image. Sediment Yield Index (SYI) of the Govindsagar catchment has been derived using Equations (1) and (2). The value of every parameter in Equations (1) and (2) have been put and the resultant value of Vs has come out to be 2.74 M・m3/year × 103 (Table 3). The sediment yield rate in Govindsagar reservoir is found to be 0.07 ha・m/year.
5. Conclusions
Soil erosion study and measurements of sediment are quite cumbersome since many factors influence the erosion rates. Precise field observations such as sediment measurement, soil properties, rainfall duration and intensities, vegetation densities and agricultural practices are required. The sediment yield rate of Govindsagar reservoir is found to be 0.07 ha・m/year which falls under very low category [33] . Thus, the catchment
![]()
Table 3. Sediment yield and its parameters in Govindsagar reservoir.
is subjected to least soil erosion, and has little effect on water storage capacity. This suggests that the catchment possesses gentle and uniform topography where soil erosion is not a serious threat. Moreover, land use/cover analysis also suggests that the catchment is predominantly agriculture based besides having sufficient forest cover, which restricts the soil erosion.
The result of the study clearly demonstrates that the catchment is less prone to soil erosion but requires regular monitoring from the planners and decision makers for making it better in terms of land and water conservation. Remote sensing data can be employed to derive surface features and can be used as an important input in such studies.
Acknowledgements
The authors are thankful to the Chairman Department of Geology, A.M.U, Aligarh for extending infrastructure facilities. Senior author wishes to thank Uttar Pradesh Council of Science and Technology, Lucknow for financial support (CST/SERPD/D-1225).