Assessing the vulnerability to soil erosion of the Ukai Dam catchments using remote sensing and GIS

Transcription

1 Hydrologicat Sciences Journal-des Sciences Hydrologiques, 7( ) February 00 Assessing the vulnerability to soil erosion of the Ukai Dam catchments using remote sensing and GIS SANJAY K. JAIN & M. K. GOEL National Institute of Hydrology. Roorkee 7667 (UP). India sanjavfc/nih.emet.in Abstract The investigation of basins for planning soil conservation requires a selective approach to identify smaller hydrological units, which would be suitable for more efficient and targeted conservation management programmes. One criterion, generally used to determine the vulnerability of catchments to erosion, is the sediment yield of a basin. In India, sediment yield data are generally not collected for smaller sub-catchments and it becomes difficult to identify the most vulnerable areas for erosion that can be treated on a priority basis. An index-based approach, based on the surface factors mainly responsible for soil erosion, is suggested in this study. These factors include soil type, vegetation, slope and various catchment properties such as drainage density, form factor, etc. The method is illustrated with a case study of subcatchments immediately upstream of the Ukai Reservoir located on the River Tapi in Gujarat State, India. The area is divided into 6 watersheds and different soil, vegetation, topography and morphology-related parameters are estimated separately for each watershed. Satellite data are used to evaluate the soil and vegetation indices, while a GS system is used to evaluate the topography and morphology-related indices. The integrated effect of all the parameters is evaluated to find different areas vulnerable to soil erosion. Two watersheds were identified as being most susceptible to soil erosion. Based on the integrated index, a priority rating of the watersheds for soil conservation planning is recommended. Key words soil erosion; soil conservation: sediment yield; GS; remote sensing; NDVI; slope; India Estimation de la vulnérabilité à l'érosion des sols des bassins du Barrage Ukai à l'aide de la télédétection et d'un SIG Résumé La planification de la conservation des sols nécessite une approche sélective pour identifier des unités hydrologiques plus petites que les bassins versants, dont l'échelle serait plus pertinente pour des programmes plus efficaces et ciblés de conservation. Un critère, généralement utilisé pour déterminer la vulnérabilité d'un bassin versant à l'érosion, est son apport solide. En Inde, les données sur l'apport solide ne sont généralement pas disponibles par sous-bassins si bien qu'il est difficile d'identifier les zones les plus vulnérables à l'érosion, à traiter en priorité. Nous proposons une approche basée sur un indice exprimé en fonction de facteurs caractérisant la surface du sol. Ces facteurs englobent le type de sol, la végétation, la pente et différentes propriétés du bassin comme la densité de drainage, un facteur de forme, etc. La méthode est illustrée avec un ensemble de sous-bassins à l'amont immédiat du Barrage Ukai, sur la Rivière Tapi, dans l'état de Gujarat en Inde. La zone est divisée en 6 sous-bassins, et différents paramètres caractérisant le sol, la végétation, la topographie et la morphologie sont estimés pour chacun d'entre eux. Des données satellitaires sont utilisées pour évaluer les indices liés au sol et à la végétation, tandis qu'un système d'information géographique est utilisé pour évaluer ceux qui sont liés à la topographie et à la morphologie. L'effet intégré de tous les paramètres est estimé afin d'identifier les différentes zones vulnérables à l'érosion. Deux sous-bassins particulièrement concernés ont ainsi été identifiés. Sur la base de l'indice intégré, un ordre de priorité pour la conservation des sols est recommandé. Mots clefs érosion des sols; conservation des sols; apport solide; SIG; télédétection; pente: Inde Open for diseussion until I August 00

2 Sanjay K. Jain & M. K. Goel INTRODUCTION The development of land and water resources on a sustainable basis without deterioration and with a constant increase in productivity is the mainstay of mankind. Soil erosion is a complex dynamic process of land denudation by which productive surface soils are detached, transported and accumulated at a distant place. The detachment of soil particles occurs either by hydrological (fluvial) processes of sheet, rill or gully erosion, or through the action of wind. Soil erosion results in loss of precious soil resources for cultivation and causes siltation of reservoirs and natural streams (Kothyari, 996; Biswas et al., 999; Jain & Dolezal, 000). In India, about % of the total land area is prone to erosion (Narayan Dhurva & Rambabu, 98). The formulation of proper basin management programmes for sustainable development requires an inventory of the quantitative soil loss erosion and the priority classification of basins. A basin with a higher rate of erosion needs to be given higher priority for soil conservation measures to be adopted. The All India Soil and Land Use Survey (AISLUS), established in 988, has been assigned the task of priority delineation. Initially, the AISLUS conducted soil surveys in the upper parts of catchments using Survey of India (SOI) topographic maps and village cadastral maps. Erosion control treatments were started in the upper parts of the catchments with a view that treatments taken up in the downstream catchments at later stages would not be adversely affected by unprotected upper reaches. During the early and mid 960s, aerial photographs were used to identify severely eroded areas in different catchments. Subsequently, soil conservation works were initiated in close vicinity of the reservoirs and sediment carrier streams (Bhan, 997). Sediment yield from a catchment is one of the main criteria for assessing the vulnerability of a watershed to soil erosion. However, this criterion requires continuous monitoring of sediment samples at the catchment outlet. Such data are hardly available in India for small watersheds. Although the sediment yield from large basins can be obtained from such observations, it is not possible to ascertain the vulnerability to soil erosion of small watersheds within a basin. A soil conservation programme is an expensive and cumbersome process, carried out in steps starting from the most vulnerable (highest sediment producing) region. Therefore, there is a need to assign relative priorities to different regions within a catchment. The major factors responsible for soil erosion include rainfall, soil type, and vegetation, topographic and morphological characteristics of the basin. (Kothyari, & Jain, 997). Where there is a lack of data on rainfall and sediment yield, the relative vulnerability of watersheds can be assessed with respect to time-independent factors (soil type, topography and morphology). With the advancement of remote sensing techniques and data acquisition, it is now possible to generate and revise vegetation resource maps at the scale of even a few metres. The effect of vegetation can also be incorporated in such analysis. A geographic information system (GIS) is a computerbased system designed to store, process and analyse geo-referenced spatial data and their attributes. Using a GIS, the topographic and morphometric analysis may be carried out efficiently and different layers of information can be integrated, geographic information systems (using traditional and remotely sensed data) have already proved to play a very important role in analysing soil erosion and sediment yield, as evident from recent studies in the Indian Peninsula and southeast Asia (Jain & Kothyari, 000; Baban & Yusof, 000).

3 Assessing the vulnerability to soil erosion of the Ukai Dam catchments In the present study, digital analysis of remotely sensed data has been carried out to assess the vegetation and soil related indices. Topographic and morphometric indices have been generated in a GIS using the topographic information (contour and drainage) from the SOI toposheets. Finally, all the indices have been combined to prioritize different watersheds in the area. METHODOLOGY In assessing the relative vulnerability of different watersheds to soil erosion, the major factors responsible for soil erosion were considered using the Watershed Erosion Response Model (WERM). The determination of the different factors is briefly discussed below. Rainfall The amount and intensity of rainfall affect the sediment yield from a basin. Rainfall is a random meteorological phenomenon. If a dense network of raingauge stations and long-term rainfall data are available, then the effect of rainfall characteristics on soil erosion may be taken into consideration. However, within small regions, rainfall characteristics do not vary to a large extent and can be assumed to be similar over a larger time span. Vegetation Land use may have a major influence on erosion. Vegetation reduces the raindrop's capability to detach soil particles and significantly affects the erosion process. The effectiveness of vegetation depends on the height and continuity of canopy, density of ground cover and the root density. Generally, forests are most effective in reducing erosion because of their large canopies (Jain & Dolezal, 000). The classification and mapping of vegetation are fundamental tools for obtaining knowledge about vegetation cover and its relationship to the environment. A number of methods have been used to identify different phenological stages of vegetation, including the application of the normalized difference vegetation index (NDVI), which is used as an indicator of vegetation condition (Ramamoorthi et al., 99) and is expressed as: ND VI = (MR - VIS)I{NIR + VIS) ( ) where, Mi? and VIS are reflectance in near-infrared and visible red bands, respectively. The normalization minimizes the effect of illumination geometry as well as surface topography. Soil type The physical properties of soil affect its infiltration capacity and the extent to which the soil can be detached, dispersed and transported. The properties which most influence

4 San]ay K. Jain & M. K. Goel erosion include soil structure and texture, organic matter content, moisture content, density (compactness), shear strength, as well as chemical and biological characteristics (Sharma et al., 990). To study the effect of soil conditions in the watersheds, the soil brightness index (SBI) for each watershed was estimated. The variability in agriculture scene data can be captured in two dimensions by soil brightness index, as shown by Sharma et al. (990), who gave coefficients for the calculation of brightness using the data of IRS-B LISS-II. To compute these coefficients, Sharma et al. (990) took about 60 soil samples with a wide range of physical properties. A principal components analysis was performed on the combined soil data to determine the overall distribution of the soil spectral data in four-dimensional space. A positive linear relationship was found among all the four bands (B-B), with maximum correlation between bands and. The formula used for the calculation of soil brightness index (SBI) is as follows: SBI= () Topography One topographic feature that mostly influences the erosion process is the degree of slope: a higher slope resulting in higher erosion. With the advent of GS techniques, it is now possible to prepare the digital elevation model (DEM) of an area. Using the DEM, the slope on all the grids in the area can be calculated and this information can be utilized for assessing the relative vulnerability to soil erosion. Morphology The relationships between the morphology of streams and sediment yield have been considered important for many decades, especially when changes in morphology might somehow be linked to changes in sediment yield from the landscape. Important morphological characteristics of a watershed include drainage density, form factor, elongation ratio and circulatory ratio. Drainage density is defined as the quotient of the cumulative length of the streams to the total drainage area and is expressed in length per unit area (Choubey & Jain, 99). A higher drainage density represents a relatively higher number of streams per unit area and thus a rapid storm response. It also represents conditions favourable for higher erosion from the catchment. The form factor is defined as the ratio of the basin area to the square of basin length. The circulatory ratio is defined as the ratio between the area of the basin and the area of the circle having the same perimeter as that of the basin (Choubey & Jain, 99). A higher circulatory ratio induces lesser erosion. The elongation ratio is the ratio between the diameter of a circle having the same area (as that of the basin) and the maximum length of the basin. A higher elongation ratio induces lesser erosion. The vast size of areas, time constraints and yearly variations demand a fast inventory to be made of the situations. Land-use changes with time need continuous monitoring. In all these circumstances, remote sensing can be looked upon as an aid in land-use mapping. Land-use features can be identified, mapped and studied on the

5 Assessing the vulnerability to soil erosion of the Ukai Dam catchments basis of their spectral characteristics: healthy green vegetation has considerably different characteristics in visible and near-infrared regions of the spectrum, whereas dry bare soil has a relatively stable reflectance in both regions of the spectrum. Thus, by using multispectral data appropriately, different ground features can be differentiated from each other and a thematic map depicting land use can be prepared. Using a GIS database, topographical and morphological characteristics of the watersheds can be estimated quite easily and all such information can be integrated for assessing the vulnerability of different watersheds to soil erosion. Using this methodology, the relative vulnerability of different watersheds has been assessed in the immediately contributing catchment area of Ukai Reservoir. CASE STUDY OF UKAI CATCHMENT The study area The Ukai is the largest multipurpose project so far completed in Gujarat State, India. The Ukai Dam is located across the River Tapi near Ukai village in Surat District, Gujarat State. It is located between longitudes 7 '"-78 6'0"E and latitudes 0 '0"~ '0"N. The total catchment area of the Ukai Reservoir (6 km ) lies in the Deccan Plateau. For the present study, the catchment area immediately upstream of the Ukai Reservoir was investigated to assess the relative vulnerability to soil erosion. The study area is shown in Fig.. The AISLUS is engaged in conducting rapid reconnaissance surveys for prioritization of smaller hydrological units within the catchment areas of river valley projects and flood prone rivers. It has developed a system for delineating and codifying the catchment areas into smaller hydrological units. The > JF ) () w. Indus drainage. Ganges drainage. Brahmaputra drainage. All drainageflowinginto Bay of Bengal except those of and. All drainage flowing into Arabian Sea except that at 6. WesternRajasthan mostly ephemeral drainage (i) Fig. The location of the study area.

6 6 Sanjay K. Jain & M. K. Goel V ""Y Scale : Fig. The sub-basin map of the study area, entire country has been divided into six major water resources regions. There are basins, catchments, 00 sub-catchments and 7 watersheds, following a fivestage delineation approach. The study area considered herein falls under region in the Tapti basin. The area has 86 watersheds and covers 6.9 million ha (ASLUS, 988). The total area is covered in two catchments, viz. left bank Ukai Dam to Puma confluence and right bank Ukai Dam to Purna confluence, with five and seven watersheds respectively in these two catchments. In the original codification scheme, the study area is covered in watersheds, but, to handle the data, four watersheds were further subdivided into CA, CA, CA and CA6 making the total number of watersheds 6 (Fig. ). The total study area up to the dam site is km", Data used for the study Various watersheds in the area of interest were marked using the : scale Watershed Atlas of India maps (AISLUS, 988). For preparation of vegetation and soil indexes, the remote sensing data of the LISS-II sensor of the RS-IB satellite, at spatial resolution of 6 m, were used. The data of 7 December 99 were used in the present study, since the vegetation is more developed after the monsoon season and appears clearly in the remote sensing image. The remote sensing data of all the watersheds studied are covered in path/row 0- in quadrants Al and Bl of the satellite. For the preparation of drainage and contour maps at higher scales, Survey of India toposheets at scales of :0 000 and :0 000 were used. The area was covered by 9 toposheets. Interpretation and analysis The Watershed Erosion Response Model (WERM) used in the study considers four parameters, namely, morphological and topographical features, vegetation index and

7 Assessing the vulnerability to soil erosion of the Ukai Dam catchments 7 soil brightness index (CWC, 999). The image processing was carried out on the ERDAS IMAGINE 8.. system. The remote sensing image was geometrically rectified with respect to the Survey of India (SOI) :0 000 topographical maps. The geo-referenced images of different quadrants (Al and Bl) were stitched together using the MOSAIC function. A base map of the area of interest was prepared using the :0 000 toposheets. This base map contained the overall catchment boundary and the boundaries of various watersheds. This base map was used to extract the remote sensing data of the area of interest from the full scene of the satellite. Different watersheds were also separated from the full catchment area. The generation and calculation of all these parameters is explained below. Preparation of vegetation maps Vegetation is important in assessing the relative vulnerability of different watersheds to soil erosion. The presence of vegetation reduces both the detachment of sediments and their transportation. Remote sensing is a fundamental tool for the classification and mapping of vegetation. The normalized difference vegetation index (NDVI) as referenced to IRS-B LISS-II data, is expressed as: ND VI = ( - B)I{B + B) () where B and B are the reflectance in band ( \im) and band ( jam), respectively. From the geo-referenced raw image, the NDVI image was derived using the above equation. Using the base map, the NDVI images were separated for each basin, and then organized into five different classes using the technique of unsupervised classification. These five classes represented the extent of NDVI at each grid in the watershed. Since each basin was to be considered as a single identity, the area-weighted-vegetation {AWV) was calculated for each watershed as follows:, A\-wV\ + A-wV + A-wV + A-wV + A-wV AWV = () AÏ + A + A + A + A where AI...A are the areas under each vegetation class, and wvl...wv are the weights for each vegetation class. The range of AWV was classified into five classes with weights varying from to. High AWV h given the lowest weight, based on the reasoning that the watershed with a higher vegetation amount must suffer lesser erosion. Thus, a watershed with higher vegetation content is given low vulnerability. Preparation of soil brightness index maps Similar to the case of NDVI, the SBI images obtained for each watershed were divided into five different classes using unsupervised classification. These five classes represented the extent of SBI. In a similar method to that used for vegetation, area weighted soil (A WSo) values were calculated for all the watersheds. Then the range of A WSo values was divided

8 8 Sanjay K. Jain & M. K. Goel into five classes with weights varying from to. A high AWSo was given higher weight. This signifies that the higher the soil content (barren land), the higher is the vulnerability to erosion. Estimation of morphological/topographical parameters For the estimation of morphological and topographical parameters, Survey of India toposheets were used. The drainage network was derived from the :0 000 Survey of India toposheets: the area is covered in part or in full by 9 topographical maps. All the maps were digitized and different thematic layers were generated for the contours and the drainage pattern. For drainage networking and creation of DEMs, the GIS software ILWIS. (Integrated Land and Water Information System) was used. The Strahler system (Strahler, 96) was used for stream ordering. All the watersheds were found to be in the range of sixth or seventh order. After networking, the lengths of streams of each order were evaluated separately for each watershed. Using this GIS database, the physical characteristics of the watersheds, such as the drainage density, form factor, circulatory ratio and elongation ratio were estimated. Based on the length of different order streams, drainage density was calculated for each watershed separately and divided into five classes (weights varying from to ). The higher the drainage density, the higher will be the vulnerability to erosion and, hence, the greater the weight. Form factor, circulatory ratio and elongation ratio were calculated similarly for each watershed and divided into five classes. Higher values of the form factor, circulatory ratio and elongation ratio induce lesser erosion and higher values in these cases was assigned less weight. The different weights obtained for each morphological parameter were averaged out and again divided into five different classes. Thus, a single weight was assigned for all the morphological parameters taken together. Weights were assigned to each different range of average morphological weight, assuming that higher morphological weight induces higher erosion. Slope layer In addition to the basin characteristics, slope is another prominent factor for soil erosion. The higher the slope, the greater will be the erosion. In the present case, the ILWIS system was used for evaluation of slope maps. Using the contour information for each watershed, a DEM was generated, and used for the estimation of slope factor for all the watersheds. Since a watershed may contain many slope categories, area weighted slope (AWS) was calculated for each watershed in the same way as for soil and vegetation. Relative vulnerability of watersheds To account for the integrated effect of all the four parameters considered in this study, the individual weights of all the parameters were added together. This sum was further sub-divided into four different categories for the purpose of assessing the relative

9 Assessing the vulnerability to soil erosion of the Ukai Dam catchments 9 vulnerability, and codes were assigned to each category. Watersheds with higher final weights were considered to be most vulnerable to soil erosion. Thus, a watershed with a code of is highly vulnerable to soil erosion and must be given the highest priority for the purpose of basin treatment and for the adoption of soil conservation measures. Table presents the relative prioritization status of the watersheds, and the range of accumulated weights and the corresponding priority code is given in Table. Table Priority code derived for all watersheds, using 7 December 99 data. Watershed number CA() CA() CA CA() CA() CA CA CA CA() CA() CA CA CA CA6() CA6() CA7 Weight: Morphology Soil Vegetation Slope Sum of weights Priority code ->. Tabic Range of accumulated weights and the corresponding priority code. Range of accumulated weights > - - < Vulnerability Very High High Medium Low Code CONCLUSION Investigation of basins with a view to planning soil conservation measures is an important aspect of basin management. Sediment yield data are generally not collected for smaller sub-catchments and it becomes difficult to identify the areas most vulnerable to erosion that should be treated on the basis of priority. An index-based approach is suggested using the surface factors mainly responsible for soil erosion soil type, vegetation, slope and various catchment properties such as drainage density, form factor, etc. With the advent of remote sensing data and GS capabilities, it has become very convenient to calculate indices for different basins such that their relative vulnerability to soil erosion can be assessed. A case study is presented of vulnerability assessment for the catchment of the Ukai Reservoir, in Gujarat, India. Watersheds with different soil erosion potentials have

10 0 Sanjay K. Jain & M. K. Goel been assessed with a view to adopting soil conservation measures. Watersheds such as CA(), CA, CA and CA, were classified as Very High priority from the point of view of soil conservation. Similarly, watersheds CA(), CA and CA6() were classified as High priority. The major factors that contribute to soil erosion potential in these watersheds were analysed. For the Very High priority watersheds, it can be seen that all factors apart from soil/vegetation in watershed CA(), slope in watersheds CA and CA and morphology in watershed CA, are the dominant factors of soil erosion. Of the High priority watersheds, CA() is affected by all factors except soil, while soil plays dominant role in watershed CA and morphology in watershed CA6(). Therefore, if some conservation measures have been adopted in the catchment, or land-use practices have undergone some change, their impact on soil erosion may be analysed using the remote sensing technique. REFERENCES AISLUS (All India Soil and Land Use Survey) (988) Watershed Alias of India Dept. of Agriculture and Co-operation. Government of India ( : M scale). Baban, S. M... & Yusof, K. W. (00) Modelling soil erosion in tropical environments using remote sensing and geographical information systems. Hydro/. Sci. J. 6() Bhan S. (997) Strategy for watershed development and management../. Soil Wat. Conserv. S (&), 0-. Biswas S., Sudhakar, S. & Desai, V. R. (999) Prioritisation of sub-watersheds based on morphometrie analysis of drainage basin: a remote sensing and GIS approach../ Indian Soe. Remote Sens. 7(), -66. CWC (Central Water Commission) (999) Watershed prioritization using remote sensing and GIS: Hirakud Dam. Unpublished report of CWC, New Delhi, India. Choubev V. K. & Jain, S. K, (99) Morphometrie analysis of Sabarmati basin using remotely sensed data. Hydro! J. Indian Assoc. Hydrologists XV(&), -0. Jain, S. K. & Dolezal, F. (000) Modelling soil erosion using EPIC supported by GIS, Bohemia, Czech Republic../. Environ. Hydro!. 8 Paper, -. Jain, M. K. & Kothvari. U. C. (000) Estimation of soil erosion and sediment vield usint» GIS. Hydro! Sei. J. () Kothvari, U. C. {996) Erosion and sedimentation problems in India. In: Erosion and Sediment Yield: Global and Regional 'Perspectives (ed. by D. E. Wailing & B. W. Webb) (Proc. Exeter Symp.. July 996), -0. IAIIS Publ. no. 6. Kothvari, U. C. & Jain, S. K. (997) Sediment yield estimation using GS. Hydro! Sci. J. 6() Naravan Dhurva, V. V. & Rambabu (98) Estimation of soil erosion in India. J. Irrig Drain. Ençnç ASCE 09(), '9-. Ramamoorthi, A. S., Thiruvengadaehari, S. & Kulkarni, A. V. (99) IRS-A applications in hydrology and water resources. J. Current Sci. 6(&), Sharma, S. A., Bhatt, II. P. & Ajai (990) Generation of brightness and greenness transformations for RS LISS II data. J. Indian Soc. Remote Sens. 8(), -. Strahler, A. N. (96) Quantitative ueomorphologv of drainaue basins and channel networks. In: Handbook of Applied Hydrology (ed. by V. T. Chow), vol., McGraw-Hill, New York, USA. Received 8.Inly 000; accepted 7 August 00