Terrain Characterization for Soil Resource Mapping Using IRS-P6 Data and GIS - A Case Study From Basaltic Terrain of Central India

In the present study, landforms and soils have been characterized in Borgaon Manju watershed of basaltic terrain located in Akola district, Maharashtra, Central India. Terrain characterization using Shuttle Radar Topography Mission (SRTM) elevation data (90 m) and IRS-P6 LISS IV data in conjunction with adequate field surveys shows nine distinct landforms. Soil resource inventory shows fourteen soil series in the study area. Soils formed on gently sloping (3–8 %) subdued plateau are very shallow (23 cm), moderately well drained, moderate (15–40 %) surface stoniness, severely eroded, clayey and slightly alkaline in reaction, whereas, the soils formed on level to nearly level (0–1 %) slope in the main valley are very deep (>150 cm), well drained, very slight (<3 %) surface stoniness, moderately eroded with clayey surface and moderately alkaline in reaction. Soils in the watershed are grouped into Lithic Ustorthents, Vertic Haplustepts, Calcic Haplustepts, Typic Haplustepts, Typic Haplusterts and Sodic Calciusterts. The study demonstrates that the analysis of SRTM elevation data and IRS P6–IV data in Geographic Information System (GIS) with adequate field surveys helps in characterization of landforms and soils in analysis of landscape-soil relationship.     

A method for groundwater prospect zonation in data poor areas using remote sensing and GIS: a case study in Kalikavu Panchayath of Malappuram district,Kerala, India

Abstract

The present study was an attempt to delineate potential groundwater zones in Kalikavu Panchayat of Malappuram district, Kerala, India. The geo-spatial database on geomorphology, landuse, geology, slope and drainage network was generated in a geographic information system (GIS) environment from satellite data, Survey of India topographic sheets and field observations. To understand the movement and occurrence of groundwater, the geology, geomorphology, structural set-up and recharging conditions have to be well understood. In the present study, the potential recharge areas are delineated in terms of geology, geomorphology, land use, slope, drainage pattern, etc. Various thematic data generated were integrated using a heuristic method in the GIS domain to generate maps showing potential groundwater zones. The composite output map scores were reclassified into different zones using a decision rule. The final output map shows different zones of groundwater prospect, viz., very good (15.57% of the area), good (43.74%), moderate (28.38%) and poor (12.31%). Geomorphic units such as valley plains, valley fills and alluvial terraces were identified as good to excellent prospect zones, while the gently sloping lateritic uplands were identified as good to moderate zones. Steeply sloping hilly terrains underlain by hard rocks were identified as poor groundwater prospect zones.     

Geoinformatics as a tool for the assessment of the impact of ground water quality for irrigation on soil health

Secondary salinisation is the most harmful and extended phenomenon of the unfavourable effects of irrigation on the soil and environment. An attempt was made to study the impact of poor quality ground water on soils in terms of secondary salinisation and availability of soil nutrients in Faridkot district of Punjab of northern India. Based on physiographic analysis of IRS 1C LISS-III data and semi-detailed soil survey, the soil map was finalized on a 1:50,000 scale and digitized using Arc Info GIS. Georeferenced surface soil samples (0–0.15 m) from 231 sites were collected and analyzed for available phosphorus (P) and potassium (K). Interpolation by kriging produced digital spatial maps of available P and K. Ground water quality map was generated in GIS domain on the basis of EC (electrical conductivity) and RSC (residual sodium carbonate) of ground water samples collected from 374 georeferenced tube wells. Integration of soil and ground water quality maps enabled generating a map showing degree (high, moderate and low) and type (salinity, sodicity and both) of vulnerability to secondary salinization. Fine-textured soils have been found to be highly sensitive to secondary salinisation, whereas medium-textured soils as moderately sensitive to secondary salinisation. The resultant map was integrated with available P and K maps to show the combined influence of soil texture and ground water quality on available soil nutrients. The results show that available P and K in the soils of different physiographic units were found in the order of Ap1 < Ap2 < Ap3. The soils of all physiographic units had sizeable area having high content of P (>22.5 kg / ha) and medium available K (135–335 kg ha⁻¹) in most of the test sites when irrigated with saline, sodic or poor quality water.     

Landscape analysis for pedo-geomorphological characterization in part of basaltic terrain,central India using remote sensing and GIS

In the present study, an attempt has been made to analyse IRS-ID LISS-III satellite data in conjunction with field observations for geomorphological mapping and pedo-geomorphological characterisation in Mohgaon area of Nagpur district, Maharashtra. Analysis of satellite data reveals distint geomorphological units viz., plateau top, isolated mounds, linear ridges, escarpments, plateau spurs, subdued plateau, rolling plains, pediments, narrow valleys and main valley floor. Soil profiles, studied on different identified landforms, showed variation in site and morphological charactaristics. Moderate soil erosion occurs on plateau top, isolated mouds, plateau spurs, rollinmg plains and pediments. Severe erosion was identified on escarpments and subdued plateau and narrow valleys suffer very slight erosional hazards. Moderately well drained soils were found on rolling plains, pediments, narrow valleys and main valley floor. Well drained soils were noticed on plateau top and plateau spurs. Very shallow soils were found on the plateau top and isolated mounds. Shallow soils are found in linear ridges, escarpments, plateau spurs and rolling plains. Moderately deep and deep soils are found on subdued plateau, pediments and main valley floor. The landform-soil relatioinship reveals that the soils on the plateau top and isolated mounds are very shallow, well drained, clay textured. The soils on the narrow valleys and main valley floor are deep, moderatly well drained, and clayey in texture. It also indicates that landform-soil processes are governed by physiographic position, drainage, slope and erosion conditions of the area. The present study reveals that the analysis of remotely sensed data in conjunction with field observations in GIS will be of immense help in geomorphology mapping, analysis of landform-soil relationships and generation of their geo-spatial database.     

Application of remote sensing data in characterization and mapping of soil resources for watershed planning in arid Western Rajasthan

Visual interpretation of IRS-L1SS-II (January, 1995) FCC (1:50,000 scale) of spectral bands 2, 3 and 4 was carried out for the identification and mapping of major physiographic units in an arid watershed of Jodhpur district (Rajasthan). Based on image characteristics and field traverses, seven major physiographic units identified are (1) hills (2) pediments, flat to undulating (3) buried pediments, moderately deep to deep, coarse textured (4) buried pediment, shallow to moderately deep and deep, medium to fine textured, saline (5) older alluvial plains, deep and very deep, coarse textured (6) younger alluvial plains, deep to very deep, very coarse textured and (7) dune complexes. Based on physiographicvariatton and soil or site characteristics such as texture, depth, slope, erosion and underneath substrata, 41 soil mapping units were identified and mapped. Final physiography, soil, slope, drainage and landuse maps were prepared on 1:5,000 scale. Taxonomically, the soils of the watershed were classified as Para-Lithic Torriorthents, coarse-loamy, Lithic/Typic Haplocambids, fine-loamy, Lithic/Typic Haplosalids and Typic Torrifluvents and Typic Torripsamments. Land suitability for various mapping units in the watershed have been assessed on the basis of soil physico-chemical characteristics.     

Assessment of groundwater potential zones using GIS technique

A case study was conducted to find out the groundwater potential zones in Kattakulathur block, Tamil Nadu, India with an aerial extent of 360.60 km². The thematic maps such as geology, geomorphology, soil hydrological group, land use / land cover and drainage map were prepared for the study area. The Digital Elevation Model (DEM) has been generated from the 10 m interval contour lines (which is derived from SOI, Toposheet 1:25000 scale) and obtained the slope (%) of the study area. The groundwater potential zones were obtained by overlaying all the thematic maps in terms of weighted overlay methods using the spatial analysis tool in ArcGIS 9.2. During weighted overlay analysis, the ranking has been given for each individual parameter of each thematic map and weights were assigned according to the influence such as soil −25%, geomorphology − 25%, land use/land cover −25%, slope − 15%, lineament − 5% and drainage / streams − 5% and find out the potential zones in terms of good, moderate and poor zones with the area of 49.70 km², 261.61 km² and 46.04 km² respectively. The potential zone wise study area was overlaid with village boundary map and the village wise groundwater potential zones with three categories such as good, moderate and poor zones were obtained. This GIS based output result was validated by conducting field survey by randomly selecting wells in different villages using GPS instruments. The coordinates of each well location were obtained by GPS and plotted in the GIS platform and it was clearly shown that the well coordinates were exactly seated with the classified zones.     

Land use changes in udumbanchola taluk, idukki district—Kerala: An analysis with the application of remote sensing data

Udumbanchola Taluk, located in the fragile zone of Western Ghats in the Idukki District of Kerala has undergone severe land use changes in the past century. The changing land use scenario of Udumbanchola Taluk was assessed using remote sensing technique. The historic land use map was derived from topographic maps of Survey of India, surveyed in 1910 and published during 1912–14. The land use map of 1997 was generated through the visual interpretation of IRS—1 C LISSIII images supported by ground truths and was observed that the original land use system was highly modified. The comparative study shows that the settlement areas that occupied hardly 0.73 % in 1910 were increased to 30.57 % in 1997 by transferring forest lands, grassland and cardamom plantations. Due to the depletion of natural vegetation and grassland system by the establishment of houses, cash crops and infrastructure facilities, combined with the absence of sufficient land management measures in the side slopes resulted severe land degradation and the present status of degraded lands is about 14.12 % of the total geographical area. It is also observed that the prevailing agro-climatic conditions favorable for the sustenance of sensitive crops like cardamom has undergone changes mainly due to the impact of large-scale land use modifications.     

Comparison of MODIS vegetation continuous field — based forest density maps with IRS-LISS III derived maps

The study compared forest cover maps derived using coarse resolution vegetation continuous fields (MODIS VCF; 500m resolution) with the maps derived from medium resolution (24m; IRS LISS-III) data. The comparison of VCF, per cent tree cover product, for the years 2000 to 2004 with LISS III forest density class maps of 2001 and 2003 was carried out for two sites representing hilly (Uttarakhand) and undulating terrains (Madhya Pradesh). Slicing VCF to corresponding forest crown cover, i.e., 0–10%, 10–40%, 40–70% and >70% produced considerable difference in forest area estimates when compared to original LISS III derived crown cover area. The corresponding value range in VCF for 0–10% of actual forest cover were 0–31% and 0–25% in 2 sites respectively, and the respective limit was consistent at 1–20% when VCF range were sliced with respect to upscaled LISS III at 500m resolution. Similarly, all other class limits were also found through iterative process. These limits were similar, within a site, across five years. Spatial Kappa match between these two data indicated higher match in 40–70% class, and also in undulating site. When compared at same resolution, similar forest area cover estimated with weighted area upscaling gave closest match. The study is useful in knowing the usability and limits of VCF product, and utility of spatial Kappa.     

Characterization of Landforms and Soils in Complex Geological Formations—A Remote Sensing and GIS Approach

The present investigation has been designed to analyze the landform and soil relationship in a geologically complex terrain of Tirora tahsil of Gondia district, Maharashtra using remotely sensed data and GIS technique. The geomorphologic units of the study area were delineated through visual interpretation of IRS–ID LISS-III data based on the spatial variation of the image characteristics. Thirteen landform units have been identified in the tahsil. The slope varied from level to nearly level with an area of about 63.76% of the tahsil. Rest of the area ranged from very gentle to moderately steep slopes. During soil survey, soil profiles were studied for morphological features. Horizon-wise soil samples were collected from the representative soil profiles on each landform unit. The depth of soil varied from 25 to 160 cm and colour from dark brown to very dark grayish brown. The texture ranged from clay loam to clayey in accordance with higher and lower topographic positions respectively. Higher available water holding capacity (AWC 285 mm) is found in low-lying area and low to medium AWC (140 mm) is noticed in the soils developed at higher elevation. The soils reaction (pH) is strongly acidic in nature (pH 5.2) on dissected hills, linear ridge and moderately weathered pediments, whereas, the soils are moderately to slightly acidic in nature (pH 5.5 to 6.5) on hills, shallow weathered pediments, moderately weathered pediments, deeply weathered pediments, narrow valleys, and broad valley floors. Slightly alkaline condition (pH 7.6) was observed on foot slopes and aggraded valley fills. The electrical conductivity of the soils is found almost same in all landforms. The cation exchange capacity of the area varies from 10.5 to 51.5 cmol(p⁺)kg^?1. The base saturation increases with decreasing elevation and slope. The four major soil orders viz, Entisols, Alfisols, Inceptisols and Vertisols are found in the study areas which are further classified into suborder and great group levels. The landform and soil relationship was analyzed to appraise the land resources in the tahsil. The study shows that the application of remotely sensed data and GIS are immensely helpful in land resources appraisal for their management on sustainable basis.     

Variations in soil spectral reflectance related to soil moisture,organic matter and particle size

Spectral reflectance measurements in the visible and near infrared wavelengths of alluvial, black cotton and lateritic soils under different conditions show that reflectance has negative association withsoil moisture and organic matter in all the three soils. In lateritic soils reflectance increases with decrease of particle size. Variations in reflectance due to changes in concentrations of parameters were generally restricted to certain concentration levels. The generally superior discriminant capability of band 4 (0.8 to 1.1 urn) is indicative of its utility in soil and soil characteristics mapping.     

Using InfoVal method and GIS techniques for the spatial modelling of landslide susceptibility in the upper catchment of river Meenachil in Kerala

A GIS-based statistical methodology for landslide susceptibility zonation is described and its application to a study area in the Western Ghats of Kerala (India) is presented. The study area was approximately 218.44 km² and 129 landslides were identified in this area. The environmental attributes used for the landslide susceptibility analysis include geomorphology, slope, aspect, slope length, plan curvature, profile curvature, elevation, drainage density, distance from drainages, lineament density, distance from lineaments and land use. The quantitative relationship between landslides and factors affecting landslides are established by the data driven-Information Value (InfoVal) — method. By applying and integrating the InfoVal weights using ArcGIS software, a continuous scale of numerical indices (susceptibility index) is obtained with which the study area is divided into five classes of landslide susceptibility. In order to validate the results of the susceptibility analysis, a success rate curve was prepared. The map obtained shows that a great majority of the landslides (74.42%) identified in the field were located in susceptible and highly susceptible zones (27.29%). The area ratio calculated by the area under curve (AUC) method shows a prediction accuracy of 80.45%. The area having a high scale of susceptibility lies on side slope plateaus and denudational hills with high slopes where drainage density is relatively low and terrain modification is relatively intense.     

Mapping of active tectonics intensity zones using remote sensing and GIS

The Bundi-Indergarh sector in southeast Rajasthan is characterized by folded and faulted Vindhyan rocks that are exposed as NE-SW trending long parallel ridges. The sector is separated from older rocks by the Great Boundary Thrust and is traversed by younger cross faults at several localities. The thematic maps of geomorphology, slope, vegetation index and morphotectonic parameters of Bundi-Indergarh sector have been prepared using IRS ID L1SS III and WiFS and, Landsat ETM digital data. These theme are integrated in GIS environment to assess the neotectonic potential in the area. The neotectonic potential map of the sector has been generated that indicates relative potential values as high (55–85), medium (35–55) and low (5–35) on 100-point scale. The observed four high potential zones in the area are located at the intersection of NE-SW and NW-SE lineaments. The study brings out methodology for assessing active tectonic potential of the area.     

Assessment of soil erosion using remote sensing and GIS in Nagpur district, Maharashtra for prioritisation and delineation of conservation units

In the present study, soil loss in Nagpur district of Maharashtra is predicted employing USLE method and adopting integrated analysis in GIS to prioritise the tahsils for soil conservation and for delineation of suitable conservation units. Remote sensing techniques are applied to delineate the land cover of the district and to arrive at annual cover factors. Results indicate that potential soil loss of very slight to slight (>5–10 tons/ha/year) exist in the valleys in north western, northern and in the plains of central and eastern parts of the district. Moderate to moderately severe erosion rates (10 to 20 tones/ha/year) is noticed in the southeastern and some central parts. Severe, very severe and extremely severe erosion types (20 to 80 tons/ha/year) are noticed in the northern, western, southwestern and southern parts of the district. The average soil loss is estimated to be 23.1 and 15.5 tons/ha/yr under potential and actual conditions respectively. Slight, moderate, moderately severe and extremely severe potential erosion covering about 41 per cent area of the district is reduced to negligible and very slight rates of actual erosion under the influence of present land cover leading to a reduction of 7421.2 tones of potential soil loss. Priority rating of the tahsils is evaluated from the area weighted mean quantum of soil loss. Multi-criteria overlay analysis with the parameters of soil erosion, slope, soil depth, land cover and surface texture with rating for the constituent classes has resulted in delineation of nine conservation units. Appropriate agronomic and mechanical practices are suggested in the identified units for minimizing the erosion hazard.     

Remote sensing and GIS in identification of soil constraints for sustainable development in Bhilwara district,Rajasthan

Soil is an integral part of ecosystem nurturing the biological system. Sustainable management of soil resources based on the consideration of constraints is the key to check land degradation and maintain productivity of biological system. To meet the objective remote sensing and GIS technology has been employed for identification of soil constraints in resource potential Bhilwara district. IRS LISS-III FCC images were interpreted for soil constraints using physiography soil approach, verified through field checking and laboratory analysis. On IRS LISS-III FCC images the salt affected soils of Kotri and Taswaria appeared in bright white to light grey tone, smooth texture with white mottles. These were also verified during ground truth and soil analysis for salinity (EC 2.90–3.32 dS m⁻¹) and sodicity (pH 9.50–9.86 and ESP 17.60–19.05). Similarly on the LISS III FCC, constraints due to water erosion near Bir, Sareri and Vijaypura soil series were apparent in light grey to whitish tone, intercepted by medium grey streaks indicating streams and exposed sub-soil. The constraints due to shallow depth associated with rock out crops and hilly areas of Balda and Delwara series appeared in greenish grey tone and coarse texture. There was close relationship between image characteristics, field observation and analytical data.     

Evaluation of groundwater potential zones in the sub-watersheds of north pennar river basin around Pavagada, Karnataka, India using remote sensing and GIS techniques

The study area is one of the watersheds of North Pennar basin, covering an area of 570 km² in Pavagada taluk of Tumkur district. The watershed has been subdivided into nine sub-watersheds namely Dalavayihalli, Maddalenahalli, Talamaradahalli, Puluvalli tank, Nagalamadike, Gowdatimmanahalli, Naliganahalli, Devadabetta and Byadanur. These nine sub-watersheds have been evaluated to delineate groundwater potential zones based on the characteristics of geomorphic units together with slope, geology, lineaments, borewell data using Remote Sensing and Geographic Information System (GIS) techniques. Slope varies from nearly level (0–1%) to very steep (>35%). The different geomorphic units in each sub-watershed consist of denudational hills, residual hills, inselbergs, pediment inselberg complex, pediments, shallow weathered pediplains, moderately weathered pediplains and valley fills. The lineament map for each sub-watershed has been prepared and the trends were analysed with rose diagrams. The analysis of borewell locations and their yield data in association with lineaments at subwatersheds level reveals that the lineaments are acting as a pathway for groundwater movement. The integrated map comprising groundwater potential zones prepared by “Union” function using GIS indicate that valley fills and moderately weathered pediplains are very good to good, shallow weathered pediplains are good to moderate, pediment inselberg complex and pediments are moderate to poor and denudational hills, residual hills and inselbergs are poor to very poor groundwater prospect zones.     

Selection of suitable sites for water harvesting structures in Soankhad watershed,Punjab using remote sensing and geographical information system (RS&;GIS) approach— A case study

A case study has been conducted to identify suitable sites for water harvesting structures in Soankhad watershed, Punjab using information technologies such as Remote Sensing and Geographical Information System (RS-GIS). The IRS-1C, P6 satellite imagery of the Soankhad watershed was used. The various Thematic maps such as land use map, hydrological soil group map, slope map and DEM map were prepared for selecting suitable site for construction of water harvesting structures. The suitable sites were not found for nala bunding and farm ponds due to steep slope, less soil thickness and high runoff velocity. Fourteen check dams and six percolation tanks were proposed for the construction as per Integrated Mission for Sustainable Development (IMSD) guidelines. The water balance study of the Soankhad watershed was also computed with monthly mean temperature and rainfall data using TM model. The average runoff for the wet season (July–September) 1996 was computed to be about 1543.82 mm and the total runoff volume from the Soankhad watershed was estimated to be about 143.52 Mm³.     

Land-Cover change analysis of garur Ganga watershed using GIS/Remote Sensing technique

Land-use change and Land-cover classes in Garur Ganga watershed of Bageshwar district in Uttranchal State during the periods 1963–1996 and 1986–1996 were analyzed through Survey of India Topographical Sheet and visual interpretation of LANDSAT 5 TM image bands 2, 3 and 4 using Geographical Information System (GIS). The detailed analysis have revealed that the area under agriculture and settlement increased from 34.98 to 42.34%. whereas the forest and barren land show a declining trend. Expansion of agriculture land and builtup areas have been found to be maximum in the 1200–1600 m elevation zone with 7–14° slope class. The loss of vegetation cover has been estimated to be 5.07% between 1963-1996 and 0.81% between 1986–1996.     

Estimation of Sediment Yield and Selection of Suitable Sites for Soil Conservation Measures in Ahar River Basin of Udaipur,Rajasthan using RS and GIS Techniques

The present study aimed at quantification of sediment yield for Ahar River basin of Udaipur district in Rajasthan, India by a regional empirical model using RS and GIS techniques. The land use/land cover (LULC) map of the study area was prepared by supervised classification using satellite imagery of IRS-P6 LISS III. Overall accuracy of the prepared LULC map was 90.78%. The major portion of the study area (49%) is covered with rangeland. Slope map for the study area was developed using digital elevation model. The slope in most of the study area (40% of the total area) ranges from 1% to 4%. In addition, drainage density map of the study area was generated on micro-watershed basis. The study area is covered by a dendritic pattern of drainage which shows that rocks in the area are homogeneous and uniformly resistant to water flow. The drainage density in the study area is 1.11 km km⁻². Annual sediment yield of the study area was quantified by Garde model. The mean annual runoff and sediment yield for the area was 37.58 million m³ and 8,760 m³/year, respectively. Finally, appropriate sites for construction of soil conservation measures are suggested using Boolean logic analysis method on combined slope and drainage maps.     

Hydromorphogeological investigations in a typical khondalitic terrain using remote sensing data

In order to examine the influence of tectonic and morphological characteristics on the occurrence and movement of ground water in Khondalitic (garnetiferrous sillimanite gneiss) suite of rocks, hydromorphogeological studies were carried out in a typical Khondalitic terrain situated in Viziangaram district of Andhra Pradesh, India. Different land forms have been identified with the aid of visual interpretation of Landsat imagery together with ground truth data in order to prepare hydromorphogeological and lineament maps. Drainage map and topographic slope map have been prepared using toposheets. These maps and other collateral data like well yields and geophysical data have been analysed to evaluate the ground water prospective geomorphic units. Ground water prospect areas are located in shallow buried pediplains and wash plains in such a way that they are identified on gently sloping uplands situated between the lineaments. Non potential areas are those, which are, low-lying areas near the streams and high slope areas near the residual hills. It is found near low lying areas i.e., beneath the streams that the khondalite must have transformed itself into kaoline and acting as barrier evidently preventing lateral movement of ground water forcing it to accumulate in flat upland areas between two streams or lineaments. From the lithologic cross sections it is found that there are four distinct subsurface layers namely (1) top soil, (2) highly weathered khondalite (kaolinised layer), (3) moderately weathered and fractured khondalite (aquifer layer) and (4) basement of granitic gneiss.     

Landscape ecological planning in a Basaltic terrain, central India, using Remote sensing and GIS techniques

In the present study, detailed field survey in conjunction with remotely sensed (IRS-1D, LISS-III) data is of immense help in terrain analysis and landscape ecological planning at watershed level. Geomorphologically summit crust, table top summits, isolated mounds. plateau spurs, narrow slopes, plateau side drainage floors, narrow valleys and main valley floor were delineated. The soil depth ranges from extremely shallow in isolated mounds to very deep soils in the lower sectors. Very good, good, moderate, poor and very poor groundwater prospect zones were delineated. By the integrated analysis of slope, geomorphology. soil depth, land use/land cover and groundwater prospect layers in GIS. 29 landscape ecological units were identified. Each landscape ecological unit refers to a natural geographic entity having distinctive properties of slope, geomorphology. soil depth, land use/ land cover and groundwater prospects. The landscape ecological stress zone mapping of the study area has been carried out based on the analysis and reclassification of tandscape ecological units. The units having minimum ecological impact in terms of slope, geomorphology, soil depth and land use/land cover were delineated under very low stress landscape ecological zones. The units having maximum ecological stress in the form of very high slopes, isolated mounds, table top summits and summit crust, extremely shallow soils, waste lands and very poor groundwater prospects were delineated into very high stress landscape ecological zones. The integrated analysis of remotely sensed data and collateral data in GIS environment is of immense help in evaluation of landscape ecological units and landscape ecological stress zones. The delineated landscape ecological stress zones in the watershed have been recommended for landscape ecological planning for better utilization of natural resources without harming the natural geo-ecosystem of the area.