Next Article in Journal
Built Environment Typologies Prone to Risk: A Cluster Analysis of Open Spaces in Italian Cities
Next Article in Special Issue
Old Wine in New Bottles: Exploiting Data from the EU’s Farm Accountancy Data Network for Pan-EU Sustainability Assessments of Agricultural Production Systems
Previous Article in Journal
Development of a Novel Helical-Ribbon Mixer Dryer for Conversion of Rural Slaughterhouse Wastes to an Organic Fertilizer and Implications in the Rural Circular Economy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prioritization of Watershed Using Remote Sensing and Geographic Information System

1
Haryana Space Applications Center, Hisar 125004, Haryana, India
2
Department of Soil and Water Engineering, Punjab Agricultural University, Ludhiana 141004, Punjab, India
3
Department of Agrometeorology, G.B. Pant University of Agriculture and Technology, Pantnagar 263153, Uttarakhand, India
4
Department of Soil and Water Engineering, CCS Haryana Agricultural University, Hisar 125004, Haryana, India
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(16), 9456; https://0-doi-org.brum.beds.ac.uk/10.3390/su13169456
Submission received: 5 January 2021 / Revised: 22 April 2021 / Accepted: 30 April 2021 / Published: 23 August 2021

Abstract

:
Soil erosion is becoming a major concern at the watershed scale for the environment, natural resources, and sustainable resource management. Therefore, the estimation of soil loss through this phenomenon and the identification of critical soil erosion-prone areas are considered to be key tasks in the soil conservation programme for the design and implementation of best management practices for specific regions or areas. In the present study, revised universal soil loss equation (RUSLE) modelling is combined with remote sensing (RS) and geographical information system (GIS) techniques and used to predict soil erosion and the prioritization of watersheds in Nainital district Uttarakhand, India. For the estimation of soil loss, different factors, namely, rainfall-runoff erosivity (R) factor, soil erodability (K) factor, slope length steepness (LS) factor, cover management (C) factor, and the erosion control practices (P) factor were computed. The data on various other aspects such as land use/land cover (LU/LC), the digital elevation model (DEM), slope, contours, drainage network, soil texture, organic matter, and rainfall were integrated to prepare a database for the RUSLE equation by employing ENVI & QGIS software. The results showed that a major portion (70.26%) of Nainital district is covered with forest, followed by area under fallow and agricultural land. Annual average soil loss ranged between 20 to 80 t ha−1 yr−1 in the study area. Out of 50 watersheds in the study area, 7 watersheds were given top priority for conserving natural resources, while 11 watersheds, mostly in the east-central part of Nainital, were kept under the next priority category. Only 4 watersheds of the total were given lowest priority. Moreover, it was concluded that major portions of Nainital district were in a severely prone category of soil erosion, and therefore required immediate action plans to check soil erosion and evade the possibility of landslides.

1. Introduction

India supports 16 percent of the world’s population on 2.42 percent of the world’s land area, and global population growth has led us to make effective use of diminishing natural resources, particularly for livelihood needs and agricultural practices. Natural resource development programs are generally applied at a water level, because catchments and sub-catchments are the basic units for land and water resource management [1]. An estimated 147 million hectares (M ha) of land comprising approximately 45% of India’s total geographical area suffers from soil erosion and land degradation problems [2]. Water and wind alone account for the degradation of 150 M ha of land in the form of active erosion, resulting in losses of approximately 5300 million tons of top soil; in addition, 25 M ha of land has been eroded due to ravines and gullies, shifting crops, salinity/alkalinity, water logging, etc. [3,4].
Nainital, a popular hill station in the Uttarakhand state of India, has mountains with a height of 2615 m and a sub-tropical to sub-humid climate with annual average rainfall of 1500 mm with high intensity. For the prevention of soil erosion, better practices and prioritization will be necessary for natural resources management.
Watershed management planning highlights the management techniques to control erosion in the catchment/watershed area [3,5]. Land and water resources are scarce and their widespread use is imperative, especially in developing countries such as India, where population pressure is steadily increasing while natural resources are dwindling day after day [6,7]. In addition to industrial expansion and the consequent need for infrastructure services, this trend will eventually lead to increased pressure on land for grain, fibre, and fodder while land is being converted for urbanization/settlement [8]. The watershed is an ideal unit for the planning and management of land and water resources [9,10]; however, realistic assessments of the hydrological behaviour of a watershed [6,7] and recognition of the inter-relationships among the linkages between soil, slope, uplands, lowlands, land use, and geomorphology [11] are important to the development of an effective management plan.
Remote sensing technology deals with the requirements of reliability and speed, and is an ideal tool for generating spatial information that is considered to be a prerequisite for planned and balanced development at the watershed level [12]. The geographical information systems (GIS) technology provides suitable alternatives for efficient management of large databases whereas the drainage characteristics of basins and sub-basins have been studied using conventional methods [13,14,15]. Morphometric analysis using remote sensing and GIS techniques has been attempted by a number of researchers [16,17]; the results of their studies concluded that remote sensing and geographical information systems are powerful tools for the study of basin morphometry [18,19,20]. Additionally, the integration of remote sensing data and GIS technologies has proved to be an efficient tool for water resources development and management projects as well as for watershed characterization and prioritization [21,22,23,24].
The integration of the RUSLE model with GIS enables detailed analysis of soil erosion because the mechanism has a spatially distributed character. GIS and remote sensing (RS) provide the model with spatial input data, whereas the RUSLE can be used to estimate watershed soil erosion [25,26].
Nevertheless, management interventions from the highest priority micro-watershed available should always be initiated [27]. Watershed prioritization is the ranking of different micro-watersheds of a watershed according to the order in which they have to be addressed with treatment and soil conservation measures [28]. The lowest-ranking micro-basin has the highest erosion risk and is urgently recommended for conservation measures [29]. Keeping in view the above concerns, the present study examined the utility of RS and GIS in the prioritization of natural resources in a hilly catchment using the RUSLE model inputs related to land use/land cover, soil etc.

2. Materials and Methods

2.1. Study Area

Nainital is a popular hill station in the Indian state of Uttarakhand in, and a sparkling jewel in the Himalayan necklace, blessed with scenic natural splendour and varied natural resources. The study area extended from 29°05′ N to 29°05′ N and 80°14′ E to 78°80′ E with an altitude of 2084 m above mean sea level. Figure 1 shows the location map of the study area. It is surrounded by mountains, the highest of which are Naina on the north (2615 m), Deopatha on the west (2438 m), and Ayarpatha on the south. The district’s total geographic area is 4064.33 km2. The study area has sub-tropical to sub-humid climate. The maximum temperature in the plain areas ranges from 42 °C to 46 °C and the minimum between 1 °C and 9 °C. In the hilly areas, the minimum temperature falls below the freezing point during winter. The rainfall rate usually increases from north to south. Quite steep to steep hills and glaciofluvial valleys are dominated by very shallow to moderately shallow, heavily drained, sandy-skeletal to loamy skeletal, neutral to slightly acidic soil with low water-holding capacity. Major crops grown in the study area are rice, maize, mandua, wheat, soybean, lentil, peas, and tomato.

2.2. Data Use

The 14-year rainfall data (2001–2014) were obtained from the Department of Meteorology, GBPUA&T, Pantnagar with location 29°24′00.0″ N 79°28′00.1″ E. The data for soil such as silt, sand and clay in soil, percentage of organic matter, soil structure code, and permeability code were collected from the National Bureau of Soil Survey and Land Use Planning (NBSS&LUP) Nagpur. Satellite imagery from LANDSAT 8 (30 m Resolution) was taken on 3 October 2015; the path and row of the satellite were 145 and 40, respectively, from https://landsat.usgs.gov/ site. ENVI and QGIS were used to create specific thematic layers such as contour, drainage, LU/LC, and other study area data sets.

2.3. Methodology

The RUSLE model quantifies soil erosion using five factors [30]
A = R K LS C P
where A is the average annual soil loss per unit area (t/ha/y)
R is the rainfall and runoff factor (MJ mm/ha/h/y)
R = i = 1 12 1.735 × 10 ( 1.5   log 10 Pi P 0.08188 )
where Pi is the monthly rainfall (mm), P is the annual rainfall (mm).
K is the soil erodibility factor (t.ha.h/ha/MJ/mm)
K = [2.766 × M1.14 × 10−7 × 12 − OM) + 4.28 × 10−3 × (S − 2) + 3.28 × 10−3(P − 3)] × 0.1317
where M = particle size parameter (% silt + % very fine sand) × (100—% clay)
OM = organic matter (%), S = soil structure code
P = soil permeability class and 0.1317 is used for the conversion to SI units.
LS is the combined slope length and slope steepness factor (dimensionless)
LS = Flow   accumulation × Cell   size 22.13 0.4 × Sin slope 0.0896 1.3
where cell size = size of grid cell (for this study 30 m)
sin slope = Z slope degree value in sin
C is the cover and management factor (dimensionless, ranging between 0 and 1.5)
C = exp × NDVI β NDVI
where α and β are parameters that determine the shape of the NDVI–C curve.
Values of α = 2 and β = 1
P is the support practice factor (dimensionless, ranging between 0 and 1)
P = 0.2394 × (Slope) × 0.4061
where slope is in percent, slope more than 30% will have the value 1. The range of P is 0–1.
The details of the above-mentioned factors and their effect on soil erosion processes are discussed below. The methodology flowchart is shown in Figure 2. For the estimation of soil erosion, the rainfall data for a 14-year period (2001–2014) were obtained from the Meteorological Department, GBPUA&T, Pantnagar, and used to estimate the rainfall runoff erosivity factor (R-factor). The soil map was procured from the National Bureau of Soil Survey & Land Use Planning and converted to digital format on 1:50,000 scale. The soil erodibility factor (K-factor) map was then digitized and prepared in QGIS using the soil map procured from NBSS&LUP; the values of K-factor were then assigned to the different soil types in the region according to U.S. customary units.
The ASTER data were used to draw up the slope map and the study field flow accumulation map, which were then used to prepare the LS factor map in QGIS. The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) instrument on-board the National Aeronautics and Space Administration’s Terra spacecraft has an along-track stereoscopic capability using a near-infrared spectral band to acquire the stereo data. With a base-to-height ratio of 0.6, ASTER has two telescopes, one for nadir viewing and the other for backward viewing. In the horizontal plane, the spatial resolution is 15 m. During the initial operation time, parameters such as the line-of-sight vectors and the pointing axis were adjusted to generate Level-1 data products with high-quality stereo system performance. Both Japanese and American science teams independently evaluated the digital elevation model (DEM) data, using different DEM generation tools and reference databases. Without ground control point (GCP) correction, the vertical accuracy of the DEM data produced from Level-1A data is 20 m with 95 percent confidence for individual scenes. The essential geolocation accuracy for DEM datasets is better than 50 m. This tends to be constrained by the precision of the spacecraft’s location. Furthermore, using GCPs to produce the stereo data results in a small improvement in accuracy. The DEM was utilized to deduce the flow direction to prepare the flow accumulation map. The LS factor map was then prepared using the slope and flow accumulation map in the raster calculator in QGIS. The drainage of the study area was delineated from the Survey of India topographical map in 1:50,000 scale and was updated by using higher resolution Landsat 8 satellite images. Radiometric and geometric corrections were made to remove errors and distortions during image pre-processing. The land-use/land-cover (LU/LC) map was developed from Landsat 8 satellite data, which, due to their high spatial resolution, provided accurate mapping of various LU/LC categories. The accuracy was assessed by computing the Kappa coefficient using ground truth data. The LU/LC map was used to prepare the land cover and management factor (C-factor) map. The values of C-factor were assigned to the different land-use/land-cover classes in the study area. The P-factor map was also prepared using the land-use/land-cover map of the watershed, and the values of P-factor were assigned to the different features based on the soil conservation practices addressed in the study area with reference to previous studies.
All five parameter maps (having the same coordinate system) were used to predict soil erosion and the spatial distribution of different soil erosion zones in the Nainital district. Raster Calculator was used to build the expression (R × K × LS × C × P), which, when applied to all cells in a raster coverage of the watershed, produced a map of average annual soil erosion.

2.4. Prioritization of Resources

The prioritization of natural resources is limited and their wide utilization is imperative and essential. The resource-based approach is considered to be practical for watershed prioritizations as it requires an integrated approach. In the present study, the knowledge-based weighting method was adopted for sub-watershed prioritization on the basis of study variables and after careful observation of field conditions. The basis for assigning weightage to different themes was according to the relative importance assigned to each parameter in the study area as depicted in Table 1.
The weightage system adopted here depends entirely on the local terrain and can vary from place to place. The water-holding capacity or porosity and permeability of the formation in the area were also considered while assigning weightage.

3. Results and Discussion

The prioritization of the watersheds of the Nainital district was carried out using remote sensing and GIS technique by assessing potential soil losses arising due to the rainfall, slope, and soil factors and reduced by land cover and management practices. This technique is quite robust in evaluating the influence of different land covers and soil management factors on quantitative estimations of soil loss. The remotely sensed data plays an important role in delineating the land cover with greater accuracy of form and extent and in assessing the related annual cover factors. Implementing the revised universal soil loss equation using the GIS integration procedures enabled soil loss predictions at the sub-water level. The estimation of spatially varying soil losses requires the generation of the thematic layers such as land use and land cover, land slope, soil characteristics, land cover, and management practices.

3.1. Land Use and Land Cover

The supervised classification technique was used to produce maps of Nainital district land use and land cover (LU/LC). The minimum distance algorithm was employed to achieve the task. The ground truth was performed for verification of objects at different locations. A perusal of data presented in Table 2 indicated that a major portion of the Nainital district is under forest cover (70.26%), while agriculture is over a small region (9.32%). The other major land uses include fallow land (14.72%), water bodies (0.12%), rivers (3.92%), rocks (0.67%), and built-up (1.01%), respectively. Data clearly indicated that a major portion (70.26%) of Nainital is covered with forest, followed by agriculture at only 9.32% of total area, while the smallest percentage of area was occupied by water bodies accounting for only 0.12% of the total area. The value of the Kappa coefficient was found to be 0.9926, or 99.26%, which indicated that the accuracy of the classification is very high. The classified images of the Nainital district are presented in Figure 3. The classified results are listed in Table 2. Nainital has its highest area under forest cover, much than others features, and the probable reason behind this is that Nainital has totally undulating topography; therefore, land cannot be used for other purposes or particularly for agricultural or settlement purposes. Secondly, the weather conditions and especially the high rainfall also favour higher vegetation growth, which results in thick forest cover. Moreover, the southern region of Nainital has relatively more area in plain; therefore, most agricultural and built-up lands are found in this region. High concentrations of built-up land are found in the southern region as it is considered a gateway to the Kumaon region. Most of the facilities and installations whether industrial, educational or health related, have been developed in this region. As the region is topographically undulated, a significant amount of water is drained off, creating many seasonal and perennial water streams. Therefore, approximately 3.92 percent of Nainital’s total area is occupied by water channels, although no reservoirs have been constructed.

3.2. NDVI Map of Study Region

The NDVI image was generated using a combination of the NIR band and Red band. NDVI is a direct indication of vegetation cover in the region. The value of NDVI increased with the increase in vegetation cover. The NDVI image of Nainital district is shown in Figure 4. The value of NDVI was at its maximum in the lower and middle regions of Nainital, but gradually decreased with increasing elevations toward the northern part of the district. The region of higher vegetation in the lower part of Nainital district has comparatively less slope, whereas the lower density of vegetation in the higher regions of Nainital is due to the steep slopes. Water flow down the steep slopes is fast, which contributes to heavy soil loss. Therefore, the topsoil depth is very shallow, which does not allow root penetration into the deeper soil layers. Other major factors behind the higher density of vegetation in the lower regions of Nainital include the higher soil water-retaining capacity, higher nutrient accumulations, etc.

3.3. Land Slope Layer

The slope features of the land are one of the most important parameters that strongly influences soil losses. Therefore, we generated a slope layer map of the Nainital district using ASTER data. A digital elevation map of the region was created, and the slope was then derived using the hydrological module available in Quantum GIS. The slope layer map of Nainital is shown in Figure 5. The figure indicated that the Nainital district has varying degrees of slope, from 1° to more than 40°. The lower areas of Nainital exhibit gradual sloping land with slopes of less than 10°; however, upper areas of the district (northern Nainital) exhibit steep slopes mostly of more than 40°. The elevation of the Nainital district gradually increases toward the northern region and the change in slope exhibits a similar trend.

3.4. Average Annual Rainfall

The rainfall layer was generated by interpolating point weather data from various stations, primarily in Nainital, which served as the interpolation centre with location 29°24′00.0″ N 79°28′00.1″ E. the inverse distance weightage (IDW) method embedded in QGIS software was employed to transform point data into the spatial layer of rainfall. The rainfall layer of Nainital district is shown in Figure 6. Average annual rainfall in Nainital is very high, almost 1500 mm in most parts of the district, and around 2000 mm in some areas of the central part of the district. The higher rainfall in the north-central region of Nainital is due to the steep slopes in this region [31]. The steep slopes induce orographic effects that lead to higher rainfall [20].

4. Soil Layers

4.1. Organic Matter

An organic matter map of Nainital district was extracted from the database created by NBSSLUP, Nagpur. The hardcopy maps created by NBSSLUP were digitized in the QGIS environment, and the attributes were assigned as given in the legend. The map originally contained values for organic carbon, which were converted into organic matter. Organic matter is higher throughout the district of Nainital, except in some pockets of southern and western Nainital. The fairly high values for organic matter in the district are probably due to the higher amounts of forest cover. The forest trees shed leaves that ultimately are mixed into the soil and maintain the higher values for organic matter [32]. The organic matter layer of Nainital district is presented in Figure 7.

4.2. Soil Structure Class

The soil structure codes were extracted from the soil database of NBSSLUP, Nagpur. The legend of the map shows the information on soil texture. The soil textural classes were derived with the help of soil texture triangles, and the structure codes were accordingly assigned. The values of the soil textural codes were their lowest in the case of sandy soil, whereas the maximum was observed for fine silty and fine loamy soils. The sandy soil exhibits less adhesive properties and is more prone to soil erosion, while silty and loamy soils exhibit a higher soil structural class values and are therefore more effective in reducing soil erosion [32]. The soil structural classes for Nainital are shown in Figure 8. A major portion of the Nainital district shows soil structure codes of 2 or less, which suggests that the soils of Nainital district are more prone to soil erosion.

4.3. Particle Size Parameter

The particle size parameter is a function of the texture of the soil and mostly depends on the silt, fine sand, and sand content. Higher values for the particle size parameter lead to heavier soil losses. The spatial distribution of the particle size parameter over the district of Nainital is shown in Figure 9. Closer scrutiny of the map revealed that the values for the particle size parameter over most of the district of Nainital are quite high. Some areas of Nainital district, especially western Nainital, have lower parameter values attributable to thicker forest cover. The distribution of particle size parameters over two districts showed that Nainital district is more prone to soil losses.

5. Hydrological Modelling

5.1. Flow Accumulation

The water from precipitation flows over land masses from the upper to lower regions in a gradual manner. It accumulates as it flows from one region to another. The water flow from regional points of origin is of low velocity, but gradually increases velocity as it passes through lower regions. The flow gradually increases not only in velocity but also in volume. Therefore, the lower regions of the topographical space are more prone to soil losses than relatively higher regions [33]. The flow accumulation over the district of Nainital is shown in Figure 10 and exhibits fairly extensive spatial distribution.

5.2. Delineation of Watershed Boundaries and Identification of Water Channel

The watershed boundaries in the district of Nainital were delineated using a digital elevation map (DEM) prepared with ASTER data as shown in Figure 11. ASTER-DEM is available at 30 m resolution and is widely acclaimed for its three-dimensional accuracy. The Grass module embedded in Quantum GIS was used to delineate the watershed boundaries. Grass computes the direction of water flow from each cell and thereafter summarizes the accumulated water flow. Based on the accumulated water flow, major water streams were also identified. Grass produces its results in a raster data model, which was transformed into a vector data model using the raster-to-vector conversion module available in Grass.

5.3. Estimation of Soil Loss Using RUSLE Model

Soil loss assessment was carried out using the revised soil loss equation (RUSLE) for Nainital district. Various factors as discussed below are required to compute soil losses using the RUSLE model.

5.4. Rainfall-Runoff Erosivity Factor (R)

The quantitative expression of the rainfall-runoff parameter causing soil erosion is the R factor. It is greatly influenced by the intensity and duration of precipitation events and by the amount of resulting runoff [34]. The R value incorporates the total precipitation intensity and duration pattern of rainfall. The annual value of the R factor related to annual precipitation was computed. The value of the R factor for Nainital varies from 1000 to 5000 MJ mm/ha/h/year, as shown in Figure 12.

5.5. Soil Erodibility Factor (K)

Factors such as texture, structure, organic matter content, and permeability are highly significant determinants of soil erodibility. Soil erodibility is regulated by a complex set of physical and chemical properties and is usually determined empirically [32]. Soil analysis data were available for all soil types found in Nainital. Using K factor values, a K factor map was prepared in QGIS, and is shown in Figure 13. The value of the K factor was found to range between 0.002 and 0.012 t ha h/ha/MJ/mm for Nainital district.

5.6. Slope Length Steepness Factor (LS)

The slope length created by DEM is based on the assumption that each slope plane consists of a homogeneous type of slope and vegetation cover, which may not be the case in practice. GIS techniques tend to predict very long slope lengths on flat to very gentle slopes while deriving topographic factors, which can result in overestimations of soil loss. As a result, the LS factor does not fully account for the hydrological processes that influence runoff and erosion, compromising its significance as a measure of the landscape’s sediment carrying capacity [20]. The slope map shown in Figure 5 was derived from ASTER DEM in percentages and degrees. The LS factor map was prepared using QGIS and is shown in Figure 14. The value of the LS factor for the study area ranged from 0 to 60 for Nainital district.

5.7. Cover Management Factor (C)

Remote sensing and GIS technologies are capable of generating a region’s LU/LC thematic layer [35]. The study area was classified into seven land use classes shown in Figure 3. The factor of crop management was assigned to various patterns of land use using the values given in Table 2. The C factor map was prepared in QGIS using an LU/LC map and C factor values and is shown in Figure 15. It was found that the crop management factor ranges between 0.2 and 1.00. The conservation practice factor for micro watersheds was assigned on the basis of percent slope. Its P factor map was prepared in QGIS and is shown in Figure 16.

5.8. Erosion Control Practices Factor (P)

The P factor value ranges from 0.3 to 1.2 and it is clear from Figure 16 that the higher the value of slope, the higher the value of the P factor. It was found to be more than 1 when the slope is above 30%. The contribution of the P factor is highly significant in Nainital district.

5.9. Soil Losses

Annual soil loss for micro-watersheds was estimated using the annual average R (based on data from daily rainfall), K, LS, C, and P variables. All the different layers, namely, R, K, LS, C, and P, were generated in GIS and overlays to obtain the product that produces the Nainital annual soil erosion map (Figure 17). The soil erosion rate (t/ha/yr) for a given watershed was measured as total soil loss. Micro watersheds in each basin were prioritized on the basis of soil erosion risk.

5.10. Watershed Level Soil Losses

The watershed erosion maps show the ranges of erosion in different watersheds and thus illustrate the nature of the watersheds based on the severity of soil erosion within them. According to the soil erosion-prone rankings, various watersheds in Nainital fall into the severely prone range. Watershed boundaries were overlaid on raster-based soil loss erosion images for computation of average soil loss from each watershed, as depicted in Figure 18. Watershed-level soil erosion in Nainital ranges from 0 to 80 tonnes per year (Figure 17). Soil erosion is strongly influenced by the slope of the land and soil organic matter. Higher levels of soil erosion are found in regions with steeper slopes. Steeper slopes not only cause higher runoff but also add more momentum to flowing water, which erodes the soil almost every year [36]. The continuous erosion of soil does not allow accumulation of organic matter; the bonding power of the soil therefore deteriorates, which again leads to higher soil erosion. The northern hills of the Nainital district exhibit steep slopes and therefore high soil erosion, i.e., 60–80 ton/ha/year. High soil erosion is also responsible for landslides, thereby making this area highly prone to natural disasters. Moving toward the southern part of Nainital, the slope gradually decreases, and soil erosion is therefore at its minimum. In most parts of southern Nainital, soil erosion varies from 0 to 20 ton/ha/year. Human settlements around the southern fringes of the district have further lowered the risk of soil erosion. Higher forest cover and agricultural lands intersecting the soil slopes in the southern region of the district also reduce the impact of rainwater.
The cropping intensity of the district is close to 200. That means farmers usually harvest two crops per year. Therefore, agricultural fields remain occupied by crops throughout the year, which helps reduce the risk of erosion. Thick forest cover, human settlements, and waterbodies further contribute to reduced soil erosion.

5.11. Prioritization of Watersheds for Sustainable Development

The watersheds identified in Figure 11 were prioritized on the basis of their potential for soil losses. Soil is assumed to be one of the most important natural resources due to its importance in providing support to vegetation, nutrients to plants, water retention capacity, and many other functions. Water is another important natural resource that should be conserved and judiciously used for sustainable development. However, water in the study region is sufficiently available in terms of precipitation; therefore, prioritization based on soil losses assumed prime importance [37]. The geospatial database was created as mentioned above in different sections. Soil erosion risk was calculated using the revised universal soil loss equation (RUSLE). Pixel-level soil erosion was averaged at the watershed boundaries so that watersheds could be prioritized. The critical level of the sub-watersheds was then determined using the criteria suggested by the authors of [38]. Each sub-watershed was ranked according to its soil erosion risk. The Nainital district was divided into 50 watersheds for prioritization. The watersheds were prioritized as listed in Table 3 and Table 4 and as shown in Figure 19. Six different categories for potential soil losses—of less than 5 t/ha/yr (slight), 5–10 (moderate), 10–20 (high), 20–40 (very high), 40–80 (severe), and more than 80 (very severe)—were defined for watershed prioritization and were assigned priority rankings of sixth, fifth, fourth, third, second, and first, respectively. That means watersheds with higher potential soil losses were assigned higher priority for the conservation of precious soil resources, while those with lower soil loss potential received less priority, in descending order. The map of prioritized watersheds in Figure 19 for Nainital district shows that the top-priority watersheds (those ranked first or second)) were in the northern part of Nainital district. Extensive soil losses triggered by steep slopes, high rainfall, and low vegetation cover are the reason for the top priority given to the watersheds in the northern region. The watersheds in the southern part of Nainital district take low priority due to lower soil losses. Seven watersheds out of a total of 50 took top priority, while 11 mostly in the east-central part of Nainital received second priority.
These watersheds should be supported with the development of physical structures to reduce soil losses. The slope gradient should be pieced into relatively low sloping zones to reduce water-induced soil losses. Massive plantations of fast-growing tree species should be extended to enhance vegetation cover and root-based soil retention capacity. Four watersheds, three in the eastern part and one in the western part of Nainital, take third priority. Fourteen watersheds were given fourth priority. In the third and fourth priority categories, contouring, terracing, and bund formation should be adopted for agricultural lands to check large-scale soil erosion. Plantations should be developed on non-agricultural lands. Ten watersheds were assigned fifth priority, while the remaining four watersheds were assigned last priority. These watersheds should be targeted for contour farming, proper bunding, etc. Additional applications of organic matter in agricultural fields will not only increase the water retention capacity of the soil, but also the percolation rate and particle adhesiveness. These practices together with persistent plant cover that helps control soil losses will be crucial for sustainable development [39]. Therefore, controlling large amounts of runoff and enhancing organic matter buildup in the soil will suffice for the purpose [32,38].

6. Conclusions

The creation of a geo-database on land use and land cover, slope, soil properties, climatic variables, and hydrological characterization was found to be very effective and useful for estimating soil losses. Landsat 8 (OLI) optical data of the 145/40 path and row were used to create a land use and land cover map of the Nainital district. According to the findings of the prioritization analyses, 32 watersheds had very high priority and 28 had high priority. These watersheds should be thoroughly investigated for soil and water conservation, water resource conservation, and proper land-use planning. This will aid in the development of a comprehensive database for each natural resource theme, which is essential for successful and efficient watershed management and the protection of deteriorating watersheds.
The proposed action plans are based on the prioritization outcomes. We recommend the construction of major check dams in high-priority watersheds and minor check dams in low-priority watersheds. Additionally, percolation tanks and recharge wells should be built and forest and grassland expansion and improvement projects should also be implemented in these watersheds. Water harvesting techniques are being implemented to prevent rainwater from being wasted in the watersheds. In addition to providing supplementary irrigation, this will increase groundwater recharge. Farmers should be encouraged to use soil conservation techniques and build farm ponds. Contour bunding is a good example of a soil and water conservation measure. Fodder, silvo pasture, and social forestry can be used to cultivate lands with low productivity and severe soil erosion. These growth activities aid in soil erosion reduction, moisture conservation, and enhanced soil productivity. The study area mostly comprised hard rock terrain. Water harvesting and artificial recharge sites would raise the water table level and meet water demand for crop irrigation during the summer season, thus meeting the needs of both urban and rural areas.
Based on the results of our study, it may be concluded that a major portion of Nainital district comes under the severely prone category of soil erosion, and therefore demands immediate attention as well as appropriate measures by government agencies and civil society to control soil erosion and avert the possibility of landslides.

7. Future Scope of Research

In the future, a segments and sub-watershed study will be carried out to further investigate the area’s infiltration potential in extensive detail.

Author Contributions

D.K., A.D. and A.S.N. were involved in the conceptualization, literature review, methodology design, investigation, data analysis, and review. A.S.N. and M.P.S. supervised this research. A.S. was involved in writing this manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of Institute and University and approved by the examiner committee.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gumma, M.K.; Birhanu, B.Z.; Mohammed, I.A.; Tabo, R.; Whitbread, A.M. Prioritization of watersheds across Mali using Remote sensing data and GIS techniques for agricultural development planning. Water 2016, 8, 260. [Google Scholar] [CrossRef] [Green Version]
  2. Mythili, G.; Goedecke, J. Economics of Land Degradation in India. In Economics of Land Degradation and Improvement—A Global Assessment for Sustainable Development; Nkonya, E., Mirzabaev, A., von Braun, J., Eds.; Springer: Cham, Switzerland, 2016. [Google Scholar]
  3. Gajbhiye, S.; Mishra, S.K.; Pandey, A. Simplified sediment yield index model incorporating parameter curve number. Arab. J. Geosci. 2015, 8, 1993–2004. [Google Scholar] [CrossRef]
  4. Meshram, S.G.; Sharma, S.K. Prioritization of watershed through morphometric parameters: A PCA-based approach. Appl. Water Sci. 2017, 7, 1505–1519. [Google Scholar] [CrossRef] [Green Version]
  5. Gajbhiye, S.; Sharma, S.K.; Awasthi, M.K. Application of principal components analysis for interpretation and grouping of water quality parameters. Int. J. Hybrid Inf. Technol. 2015, 8, 89–96. [Google Scholar] [CrossRef]
  6. Sharma, S.K.; Gajbhiye, S.; Tignath, S. Application of principal component analysis in grouping geomorphic parameters of a watershed for hydrological modelling. Appl. Water Sci. 2014, 5, 89–96. [Google Scholar] [CrossRef] [Green Version]
  7. Sharma, S.K.; Gajbhiye, S.; Nema, R.K.; Tignath, S. Assessing vulnerability to soil erosion of a watershed of tons river basin in Madhya Pradesh using remote sensing and GIS. Int. J. Environ. Res. Dev. 2014, 4, 153–164. [Google Scholar]
  8. Biswas, S.; Sudhakar, S.; Desai, V.R. Prioritisation of subwatersheds based on morphometric analysis of drainage basin: A remote sensing and GIS approach. J. Indian Soc. Remote Sens. 1999, 27, 155. [Google Scholar] [CrossRef]
  9. Gajbhiye, S.; Mishra, S.K.; Pandey, A. Prioritization of shakkar river catchment through morphometric analysis using remote sensing and gis techniques. J. Emerg. Technol. Mech. Sci. Eng. 2013, 4, 129–142. [Google Scholar]
  10. Gajbhiye, S.; Mishra, S.K.; Pandey, A. Prioritizing erosion-prone area through morphometric analysis: An RS and GIS perspective. Appl. Water Sci. 2013, 4, 51–61. [Google Scholar] [CrossRef] [Green Version]
  11. Sebastian, M.; Jeyaraman, V.; Chandrasekhar, M.G. Space Technology Applications for Sustainable Development of Watersheds; Technical Report, ISROHQTR10495; ISRO: Bengaluru, India, 1995. [Google Scholar]
  12. Ravindran, K.V.; Jeyaram, A. Groundwater prospects of Shahbad tehsil, Baran district, Eastern Rajasthan: A remote sensing approach. J. Indian Soc. Remote Sens. 1997, 25, 239–246. [Google Scholar] [CrossRef]
  13. Horton, R.E. Erosional development of streams and their drainage basins: Hydrophysical approach to quantitative morphology. Geol. Soc. Am. Bull. 1945, 5, 275–370. [Google Scholar] [CrossRef] [Green Version]
  14. Miller, V.C. A quantitative geomorphic study of drainage basin characteristics in the Clinch Mountain area, Varginia and Tennessee. J. Geol. 1957, 65, 1. [Google Scholar]
  15. Strahler, A.N. Handbook of Applied Hydrology. In Quantitave Geomorphology of Drainage Basins and Channel Networks; Section 4-II; McGraw Hill Book Company: New York, NY, USA, 1964. [Google Scholar]
  16. Reddy, G.P.O.; Maji, A.K.; Gajbhiye, K.S. Drainage morphometry and its influence on landform characteristics in a basaltic terrain, Central India–a remote sensing and GIS approach. Int. J. Appl. Earth Obs. Geoinf. 2004, 6, 1–16. [Google Scholar] [CrossRef]
  17. Singh, S.; Singh, M.C. Morphometric analysis of Kanhar river basin. Natl. Geogr. J. India 1997, 43, 31–43. [Google Scholar]
  18. Chalam, B.N.S.; Krishnaveni, M.; Karmegam, M. Correlation of runoff with geomorphic parameters. J. Appl. Hydrol. 1996, 9, 24–31. [Google Scholar]
  19. Chaudhary, R.S.; Sharma, P.D. Erosion hazard assessment and treatment prioritization of Giri river catchment, North western Himalayas. Indian J. Soil Conserv. 1998, 26, 6–11. [Google Scholar]
  20. Kumar, R.; Lohani, A.K.; Kumar, S.; Chatterjee, C.; Nema, R.K. GIS based Morphometric Analysis of Ajay river basin uptoSarathh gauging site of South Bihar. J. Appl. Hydrol. 2000, 14, 45–54. [Google Scholar]
  21. Ali, S.; Singh, R. Morphological and hydrological investigation in Hirakud catchment for watershed management planning. J. Soil Water Conserv. 2002, 1, 246–256. [Google Scholar]
  22. Pandey, A.; Behra, S.; Pandey, R.P.; Singh, R.P. Application of GIS for Watershed Prioritization and Management: A Case Study. Int. J. Environ. Sci. Dev. Monit. 2011, 2, 25–42. [Google Scholar]
  23. Pandey, A.; Chowdary, V.M.; Mal, B.C.; Dabral, P.P. Application of remote sensing and GIS for identification of suitable sites for soil and water conservation structures. Land Degrad. Dev. 2010, 22, 359–372. [Google Scholar] [CrossRef]
  24. Pandey, A.; Mathur, A.; Mishra, S.K.; Mal, B.C. Soil erosion modeling of a Himalayan watershed using RS and GIS. Environ. Earth Sci. 2009, 59, 399–410. [Google Scholar] [CrossRef]
  25. Chopra, R.; Dhiman, R.D.; Sharma, P.K. Morphometric analysis of sub-watersheds in Gurdaspur district, Punjab using remote sensing and GIS techniques. J. Indian Soc. Remote Sens. 2005, 33, 531. [Google Scholar]
  26. Tirkey, A.S.; Pandey, A.C.; Nathawat, M.S. Use of satellite data, GIS and RUSLE for estimation of average annual soil loss in Daltonganj watershed of Jharkhand (India). J. Remote Sens. Technol. 2013, 1, 20–30. [Google Scholar] [CrossRef]
  27. Naik, S.D.; Das, S.N. Watershed Prioritization of Dhansiri Catchment, Brahmaputra Basin, using Remote Sensing Techniques. Available online: http://library.cept.ac.in/cgi-bin/koha/opac-detail.pl?biblionumber=37241 (accessed on 30 April 2021).
  28. Kumar, B.; Kumar, U. Micro watershed characterization and prioritization using Geomatics technology for natural resources management. Int. J. Geomat. Geosci. 2011, 1, 789–802. [Google Scholar]
  29. Khare, D.; Mondal, A.; Mishra, P.K.; Kundu, S.; Meena, P.K. Morphometric analysis for prioritization using remote sensing and GIS techniques in a hilly catchment in the state of Uttarakhand, India. Indian J. Sci. Technol. 2014, 7, 1650–1662. [Google Scholar] [CrossRef]
  30. Benavidez, R.; Jackson, B.; Maxwell, D.; Norton, K. A review of the (Revised) Universal Soil Loss Equation ((R) USLE): With a view to increasing its global applicability and improving soil loss estimates. Hydrol. Earth Syst. Sci. 2018, 22, 6059–6086. [Google Scholar] [CrossRef] [Green Version]
  31. Rawat, J.S.; Joshi, R.C.; Mesia, M. Estimation of erosivity index and soil loss under different land uses in the tropical foothills of Eastern Himalaya (India). Trop. Ecol. 2013, 54, 47–58. [Google Scholar]
  32. Shinde, V.; Tiwari, K.N.; Singh, M. Prioritization of micro watersheds on the basis of soil erosion hazard using remote sensing and geographic information system. Int. J. Water Resour. Environ. Eng. 2010, 5, 130–136. [Google Scholar]
  33. Zhao, J.; Oost, K.V.; Chen, L.; Govers, G. Moderate topsoil erosion rates constrain the magnitude of the erosion-induced carbon sink and agricultural productivity losses on the Chinese Loess Plateau. Biogeosciences 2016, 13, 4735–4750. [Google Scholar] [CrossRef] [Green Version]
  34. Khadse, G.K.; Vijay, R.; Labhasetwar, P.K. Prioritization of catchments based on soil erosion using remote sensing and GIS. Environ. Monit. Assess. 2015, 187, 333. [Google Scholar]
  35. Tamene, L.; Adimassu, Z.; Aynekulu, E.; Yaekob, T. Estimating landscape susceptibility to soil erosion using a GIS-based approach in Northern Ethiopia. Int. Soil Water Conserv. Res. 2017, 5, 221–230. [Google Scholar] [CrossRef]
  36. Abdullah, A.F.; Aimrun, W.; Nasidi, N.M.; SAF, K.H.; Sidek, L.M.; Selamat, Z. Modelling erosion and landslides induced by farming activities at Hilly Areas, Cameron Highlands, Malaysia. J. Teknol. 2019, 81, 195–204. [Google Scholar]
  37. Gupta, P.; Uniyal, S. Watershed Prioritization of Himalayan Terrain, Using SYI Model. Int. J. Adv. Remote Sens. GIS Geogr. 2013, 1, 42–48. [Google Scholar]
  38. Vemu, S.; Pinnamaneni, U.B. Sediment yield estimation and prioritization of watershed using remote sensing and GIS. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2012, 39, 529–533. [Google Scholar] [CrossRef] [Green Version]
  39. Uddin, K.; Murthy, M.S.R.; Wahid, S.M.; Matin, M.A. Estimation of soil erosion dynamics in the Koshi basin using GIS and remote sensing to assess priority areas for conservation. PLoS ONE 2016, 11, e0150494. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location map of study area.
Figure 1. Location map of study area.
Sustainability 13 09456 g001
Figure 2. Methodology flow chart.
Figure 2. Methodology flow chart.
Sustainability 13 09456 g002
Figure 3. LU/LC map of Nainital district.
Figure 3. LU/LC map of Nainital district.
Sustainability 13 09456 g003
Figure 4. NDVI map of Nainital district.
Figure 4. NDVI map of Nainital district.
Sustainability 13 09456 g004
Figure 5. Slope map of Nainital district.
Figure 5. Slope map of Nainital district.
Sustainability 13 09456 g005
Figure 6. Average annual rainfall map of Nainital district.
Figure 6. Average annual rainfall map of Nainital district.
Sustainability 13 09456 g006
Figure 7. Organic matter map of Nainital district.
Figure 7. Organic matter map of Nainital district.
Sustainability 13 09456 g007
Figure 8. Soil structure map of Nainital district.
Figure 8. Soil structure map of Nainital district.
Sustainability 13 09456 g008
Figure 9. Soil particle size parameter map of Nainital district.
Figure 9. Soil particle size parameter map of Nainital district.
Sustainability 13 09456 g009
Figure 10. Flow accumulation map of Nainital district.
Figure 10. Flow accumulation map of Nainital district.
Sustainability 13 09456 g010
Figure 11. Watershed boundaries and water channel map of Nainital district.
Figure 11. Watershed boundaries and water channel map of Nainital district.
Sustainability 13 09456 g011
Figure 12. R factor map of Nainital district.
Figure 12. R factor map of Nainital district.
Sustainability 13 09456 g012
Figure 13. K factor map of Nainital district.
Figure 13. K factor map of Nainital district.
Sustainability 13 09456 g013
Figure 14. LS factor map of Nainital district.
Figure 14. LS factor map of Nainital district.
Sustainability 13 09456 g014
Figure 15. C factor map of Nainital district.
Figure 15. C factor map of Nainital district.
Sustainability 13 09456 g015
Figure 16. P factor map of Nainital district.
Figure 16. P factor map of Nainital district.
Sustainability 13 09456 g016
Figure 17. Soil erosion potential map of Nainital district.
Figure 17. Soil erosion potential map of Nainital district.
Sustainability 13 09456 g017
Figure 18. Watershed-level soil loss map of Nainital district.
Figure 18. Watershed-level soil loss map of Nainital district.
Sustainability 13 09456 g018
Figure 19. Watershed prioritization map of Nainital district.
Figure 19. Watershed prioritization map of Nainital district.
Sustainability 13 09456 g019
Table 1. Priority classes and the criteria adopted for watershed prioritization.
Table 1. Priority classes and the criteria adopted for watershed prioritization.
Sr. NoPriority ClassSoil Loss (t ha−1yr−1)Class
16<5Slight
255–10Moderate
3410–20High
4320–40Very High
5240–80Severe
61>80Very Severe
Table 2. Area of different land use and land cover features during Kharif season.
Table 2. Area of different land use and land cover features during Kharif season.
Sr. No.FeatureArea in (ha) of NainitalArea in (%) of Nainital
1Water bodies452.070.12
2Forest285,472.1770.26
3Rocks2714.040.67
4Builtup Land4084.471.01
5Fallow Land59,795.8214.72
6Agricultural Land37,849.509.32
7River15,918.123.92
8Total406,286.2100
Table 3. Quantitative analysis of watersheds with priority value.
Table 3. Quantitative analysis of watersheds with priority value.
Watershed No.Soil Loss
(t/h/year)
Perimeter (m)Area (m2)Slope (Degree)Priority
129.2239,887639617.993
287.40132,49737,37623.701
360.0552,08410,61324.152
442.6656,320921720.122
513.8850,141272710.534
613.4065,385732410.044
77.6786,07018,4547.615
88.00107,05122,4047.375
916.0240,417622616.964
1046.7889,20914,72919.872
119.1921,17446722.885
1261.7921,752187426.602
1313.1659,45110,03215.584
1474.7952,217861426.002
155.6763,66910,7639.815
1684.7689,95424,12526.651
17155.2041,906725427.531
1819.4568,55911,94422.244
1981.7375,33016,59825.821
2096.56102,65022,25216.351
2118.9212,67269219.634
2267.3252,268916924.812
2368.07105,17827,68613.332
240.8714,8469055.266
2512.9872,64817,3396.864
2667.5459,52010,87923.932
273.0541,91327996.336
283.2615,4056805.806
296.7058,07267656.615
3088.4822,27392421.891
3114.1894,55296757.644
3244.53127,01434,21521.872
3354.9787,04817,14014.322
34100.8817,030117029.521
3514.9892,62888875.394
3611.6338,79713735.574
3737.4738,575560328.773
3814.6726,49916875.794
395.1069,26683223.415
407.5543,43349255.375
415.9371,24179813.575
4210.3162,28572574.564
435.3658,43573864.285
4411.5244,86161025.134
4537.6039,633540511.933
4642.8651,585806515.852
4734.6956,825774314.053
4815.6367,22061548.674
498.4816,00312014.305
504.5019,3276534.696
Table 4. Priority classes, criteria adopted for watershed prioritization, and the total number of watersheds in each category.
Table 4. Priority classes, criteria adopted for watershed prioritization, and the total number of watersheds in each category.
Sr. NoPriority ClassSoil Loss (t/ha/yr)ClassNo. of Watersheds
16<5Slight4
255–10Moderate10
3410–20High14
4320–40Very High4
5240–80Severe11
61>80Very Severe7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kumar, D.; Dhaloiya, A.; Nain, A.S.; Sharma, M.P.; Singh, A. Prioritization of Watershed Using Remote Sensing and Geographic Information System. Sustainability 2021, 13, 9456. https://0-doi-org.brum.beds.ac.uk/10.3390/su13169456

AMA Style

Kumar D, Dhaloiya A, Nain AS, Sharma MP, Singh A. Prioritization of Watershed Using Remote Sensing and Geographic Information System. Sustainability. 2021; 13(16):9456. https://0-doi-org.brum.beds.ac.uk/10.3390/su13169456

Chicago/Turabian Style

Kumar, Devendra, Arvind Dhaloiya, Ajeet Singh Nain, Mahendra Paal Sharma, and Amandeep Singh. 2021. "Prioritization of Watershed Using Remote Sensing and Geographic Information System" Sustainability 13, no. 16: 9456. https://0-doi-org.brum.beds.ac.uk/10.3390/su13169456

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop