Integrated Approach of Universal Soil Loss Equation (USLE) and Geographical Information System (GIS) for Soil Loss Risk Assessment in Upper South Koel Basin, Jharkhand ()
1. Introduction
Erosion is a natural geological phenomenon resulting from removal of topsoil by natural agencies like wind, water transporting them elsewhere while some human intervention can significantly increase erosion rates. It is a major agricultural problem and also one of the major global environmental issues. Erosion is triggered by a combination of factors such as slope steepness, climate (e.g. long dry periods followed by heavy rainfall), inappropriate land use, land cover patterns (e.g. sparse vegetation) and ecological disasters (e.g. forest fires) [1]. Moreover, some intrinsic features of a soil can make it more prone to erosion (e.g. a thin layer of topsoil, silty texture or low organic matter content [2]. Thus, some scholars even categorized soil erosion into two major types: geological erosion and man-induced erosion [3]. It has been estimated that about 113.3 m ha of land is subjected to soil erosion due to water and about 5334 million tonnes of soil is being detached annually due to various reasons in India [4]. The process of soil erosion involves detachment, transport and subsequent deposition [5]. Urbanization, deforestation leading to change in land use patterns is the major cause of soil erosion in recent times. Erosion process result in soil loss from a watershed and it is difficult to estimate soil loss as it arises from a complex interaction of various hydro-geological processes [6].
Estimating the soil loss risk and its spatial distribution are the one of the key factors for successful erosion assessment. Spatial and quantitative information on soil erosion on a regional scale contributes to conservation planning, erosion control and management of the environment. Identification of erosion prone areas and quantitative estimation of soil loss rates with sufficient accuracy are of extreme importance for designing and implementing appropriate erosion control or soil and water conservation practices [7]. Researchers have developed many predictive models that estimate soil loss and identify areas where conservation measures will have the greatest impact on reducing soil loss for soil erosion assessments [8]. These models can be classified into three main categories as empirical, conceptual and physical based models [9]. Despite development of a range of physical, conceptual based models, Universal Soil Loss Equation (USLE) and modified Universal soil loss equation (MUSLE) or revised Universal soil loss equation (RUSLE) are the most popular empirically based models used globally for erosion prediction and control and has been tested in many agricultural watersheds in the world. The main reason why empirical regression equations are still widely used for soil erosion and sediment yield predictions is their simplicity, which makes them applicable even if only a limited amount of input data is available.
Remote sensing and GIS techniques have become valuable tools specially when assessing erosion at larger scales due to the amount of data needed and the greater area coverage. For this reason use of these techniques have been widely adopted and currently there are several studies that show the potential of remote sensing techniques integrated with GIS in soil erosion mapping [10,11]. The main objective of the study involved estimation of soil loss using USLE model with the input parameters generated using remote sensing and GIS.
2. Geographical Setting of the Study Area
The present study area (Figure 1) is a part of South Koel basin with an area of 942.4 sq. km. , bounded by latitude 23˚17'16"N & 23˚32'16"N and longitude 84˚14'15"E & 85˚46'51"E, falling in Lohardaga and Ranchi districts of Jharkhand. South Koel is the main river flowing in the study area with Kandani and Saphi River as its major tributaries.
The climate of the area is subtropical. The watershed experiences precipitation with an annual average of 1380 mm. The mean maximum and mean minimum temperature in the watershed are 36˚C and 7˚C respectively. The study area is a part of Chotanagpur plateau with undulating topography. The elevation ranges from 640 m to 925 m above the mean sea level. The watershed is an agricultural watershed with 84% of the geographical area under agriculture.
3. Materials and Methodology
The soil loss of the watershed was calculated using
USLE [12] with some modification in its parameter estimation to suit Indian conditions. The general universal soil loss equation is as follows:
(1)
where, A is average annual soil loss (ton∙ha−1∙yr−1); R is the Rainfall and Runoff erosivity index (in MJ mm/ha/ hr/yr); K is the soil Erodibility factor (in ton/MJ/mm); LS is the Slope and Length of Slope Factor; C is the Cropping Management Factor; P is the is the supporting conservation practice factor.
3.1. Rainfall and Runoff Erosivity Factor (R)
R is the long term annual average of the product of event rainfall kinetic energy and the maximum rainfall intensity in 30 minutes in mm per hour [13]. Rainfall erosivity estimation using rainfall data with long-time intervals have been attempted by several workers for different regions of the world [14,15]. Using the data for storms from several rain gauge stations located in different zones, linear relationships were established between average annual rainfall and computed EI30 values for different zones of India and iso-erodent maps were drawn for annual and seasonal EI30 values [16]. The derived relationship is given below
(2)
where RN is the average annual rainfall in mm.
A 5-year average annual data has been used to calculate the average annual R-factor values. Since the rainfall data available for the study area is not homogenous, an interpolation of average annual rainfall data is applied to have a representative rainfall distribution map. This rainfall distribution map is used as input for calculation of R-factor.
3.2. Soil Erodibility Factor (K)
The soil erosivity factor, K, relates to the rate at which different soils erode. However, it is different than the actual soil loss because it depends upon other factors, such as rainfall, slope, crop cover, etc. K values reflect the rate of soil loss per rainfall-runoff erosivity (R) index.
Soil map of the study area was prepared from soil survey report prepared by the National Bureau of Soil and Landuse Planning. Details such as fraction of sand, silt, clay, organic matter and other related parameters information for different mapping units were taken from the same report.
In this study, Soil erodibility (K) of the study area can be defined using the relationship between soil texture class and organic matter content proposed by [17]. Table 1 presents the soil erodibility factor (K) based on the soil texture class.
Table 1. Soil erodibility factor K (by Schwab et al., 1981 [17]).
3.3. Topographic Factor (LS)
Topographic factor (LS) in USLE accounts for the effect of topography on sheet and rill erosion. The two parameters that constitute the topographic factor are slope gradient and slope length factor and can be estimated through field measurement or from a digital elevation model (DEM).
There are many relationships available for estimation of LS factor. Among these the best suited relation for integration with GIS is the theoretical relationship proposed by [18-20] based on unit stream power theory given below.
(3)
where A is the upslope contributing factor, B is the slope angle.
Slope shape, the interaction of angle and length of slope, has an effect on the magnitude of erosion. As a result of this interaction, the effect of slope length and degree of slope should always be considered together [21].
With the incorporation of Digital Elevation Models (DEM) into GIS, the slope gradient (S) and slope length (L) may be determined accurately The precision with which it can be estimated depends on the resolution of the digital elevation Model (DEM).Here in this study 30 m ASTER DEM has been used.
Raster calculator was used to derive LS map based on flow accumulation and slope steepness [22].The expression used is as follows
(3a)
3.4. Cropping Management Factor (C)
The C factors are related to the land-use and are the reduction factor to soil erosion vulnerability. It is an important factor in USLE, since they represent the conditions that can be easily changed to reduce erosion. C factor is basically the vegetation cover percentage and is defined as the ratio of soil loss from specific crops to the equivalent loss from tilled, bare test-plots. The value of C depends on vegetation type, stage of growth and cover percentage. Therefore, it is very important to have good knowledge concerning land-use pattern in the basin to generate reliable C factor values. The most widely used remote-sensing derived indicator of vegetation growth is the Normalized Difference Vegetation Index, which for Landsat-TM is given by the following equation:
(4)
where NIR: the reflection of the near infrared portion of the electromagnetic spectrum and IR: the reflection in the visible spectrum. NDVI values range between −1.0 and +1.0 where higher values are for green vegetation and low values for other common surface materials. Bare soil is represented with NDVI values which are closest to 0 and water bodies are represented with negative NDVI values [23-25].
Landsat TM data of the study area with spatial resolution of 30 m was used for generation of NDVI image. After the production of the NDVI image, the following formula was used to generate a C factor map from NDVI values [26-28]:
(5)
where α and β are unit less parameters that determine the shape of the curve relating to NDVI and the C factor. [26,27] found that this scaling approach gave better results than assuming a linear relationship. Finally, the values of 2 and 1 were selected for the parameters α and β, respectively which seems to give good results as mentioned in related literatures. As the C factor ranges between 0 and 1, a value of 0 was assigned to a few pixels with negative values and a value of 1 to pixels with value greater than 1.
3.5. Supporting Conservation Practice Factor (P)
The support practice factor P represents the effects of those practices such as contouring, strip cropping, terracing, etc. that help prevent soil from eroding by reducing the rate of water runoff. Table 2 shows the value of support practice factor according to the cultivating methods and slope [29]. The P value range 0 to 1 where 0 represents very good manmade erosion resistance facility and 1 represents no manmade erosion resistance facility.
Table 2. Support practice factor according to the types of cultivation and slope.
4. Result and Discussion
4.1. Rainfall and Runoff Erosivity Factor (R)
The average annual rainfall data of nine surrounding rain gauge stations was used to get the rainfall distribution map of the entire watershed. The average annual rainfall in the study area varied between 880 mm to 1480 mm.
The rainfall & runoff erosivity factor map (Figure 2) was generated in Arc-GIS from average annual rainfall map. The erosivity factor varied from 508 - 584 MJ ha/mm/hr/yr.
4.2. Soil Erodibility Factor (K)
Soil erodibility is a quantitative estimation of erodibility of particular soil type and the main factor affecting the capability of the soil to erode is its soil texture. However the other factors affecting K factor are soil structure, permeability and the organic matter content.
Soil erodibility factor map (Figure 3) was prepared from soil map of the study area based on different soil textures. The soil distribution in the study area and derived erodibility factor is given in Table 3.
4.3. Topographic factor (LS)
The topographic factors slope gradient and slope length significantly influence soil erosion. Slope and flow accumulation grid was prepared from ASTER DEM.
Since USLE is only applicable to rill and inter-rill erosion, there should be upper bound to slope length. To enforce such an upper bound the flow accumulation grid was modified to make sure the flow length did not exceed 180 m, which seems a reasonable flow length according to [30]. Since DEM has a resolution of 30 × 30 m; the flow accumulation cannot exceed 6 grid cells. Therefore using raster calculator all cells with flow accumulation greater than and equal to 6 were assigned the value 6. Slope was calculated using the available slope function in ARCGIS. Using ASTER DEM derived flow accumulation and slope grid, a LS-factor map was created, which is shown in Figure 4. LS factor values range from 0 - 20 but most of the area are in 0 - 2 range with some high values occurring in the hilly areas.
Figure 2. Rainfall & runoff erosivity factor map of study area.
Figure 3. Soil erodibility factor map of the study area.
Table 3. Soil types statistics and soil erodibility factor values of the study area.
Figure 4. Topographic factor (LS) map of study area.
4.4. Cropping Management Factor (C) & Conservation Practice Factor (P)
Using the NDVI image derived from Landsat TM satellite data of the study area C factor map was generated as shown in Figure 5. The C factor values in the study area varied between 0.058 to 1.The forest areas shows C values between 0.05 to 0.2 with waterbodies showing values approaching 1. Agricultural areas in the study area show c factor values varying between 0.3 to 0.6.
In the study area no major conservation practices are followed except that the agricultural plots under cultivation are bunded. On the basis of value suggested by [31], agricultural land were assigned P factor of 0.28 and other land uses were assigned P factor of 1.
4.5. Soil Loss Mapping
The main goal of this study was to test the USLE model in the study area. The main result of USLE is prediction of erosion risk. All the USLE parameters determined for the study area were either in spatial format and/or in numerical format. The spatial maps and other factors were integrated using USLE empirical formula Equation (1)
Figure 6. Average annual soil loss map of study area.
Table 4. Area under various soil loss zones of the study area.
and analyzed in Arc GIS module. The annual soil loss map obtained were grouped in the following scales of priority: Slight (0 to 5 ton/ha/yr), Moderate (5 to 10 ton/ha/yr), High (10 to 20 ton/ha/yr), Very High (20 to 40 ton/ha/yr), Severe (40 to 80 ton/ha/hr), Very Severe (>80 ton/ha/yr) annual erosion classes as per guidelines suggested by [32] for Indian conditions. Figure 6 shows annual erosion map of study area, helpful in identification of areas vulnerable to soil erosion. Table 4 shows the spatial distribution of different erosion classes in the study area. The mean annual soil loss estimated for the entire watershed is 12.20 ton /ha/yr. As expected hilltops and bare land witness the highest annual soil loss (>80 ton/ha/yr). There is an increase in soil loss with increased slope gradient. Major part of the watershed is under slight to moderate soil loss zone. Only 18.2% of the study area is under high or very high soil loss. We find that the agricultural lands with increase in slope gradient experience more soil loss. Therefore more emphasis should be provided to high slope gradient areas for soil conservation planning.
5. Conclusion
Creation of database through conventional methods is time consuming, tedious and is difficult to handle. Therefore various thematic layers representing different factors of USLE were generated and overlaid in GIS framework to compute the spatially distributed average annual soil erosion map for the Upper South Koel basin. A quantitative assessment of average annual soil loss was done for the study area. The study revealed that area covered under slight, moderate, high, very high, severe, very severe soil loss potential zones are 64.70%, 17.10%, 10.05%, 4.65%, 1.60%, 1.90% respectively. The average annual soil loss map will definitely be helpful in identification of priority areas for implementation of soil conservation measures and effective checking of soil loss.
NOTES