Surface Runoff Estimation using Remote Sensing & GIS based Curve Number Method Ishtiyaq Ahmad, Dr. M. K. Verma Department of Civil Engineering, National Institute of Technology, Raipur, India Abstract One of the important components of hydrologic cycle is runoff and influenced by various factors including precipitation and watershed characteristics. Numbers of mathematical models are available to quantify runoff. National Resources Conservation Service (NRSC) has developed a Geographic Information System (GIS) based method known as Soil Conservation Services-Curve Number method (SCS-CN) for computing the runoff depth based on the rainfall depth. This method has universal acceptance as it is simple, predictable and stable method for computing runoff depth. This method is based on one parameter i.e. curve number, which is basically a coefficient that reduces the rainfall to runoff. In the present study SCS-CN method has been applied to estimate the runoff depth in Kharun River basin, a subbasin of Sheonath river in Chhattisgarh state. Various layers has been prepared namely base map, soil map, land use map and other map of the study area using GIS and remote sensing data. Based on the rainfall data of 21 raingauge stations in and around the study area, daily runoff depth has been estimated. Keywords Curve Number, Hydrologic Soil Group, Geographic Information System, Remote Sensing, Runoff I. INTRODUCTION This study aims to compute the runoff depth using Soil Conservation Service-Curve Number (SCS-CN) method using Remote Sensing and Geographic Information System (GIS). The SCS-CN is a quantitative description of land use / land cover / soil complex characteristics of a watershed. This model is a widely used hydrological model for estimating runoff using runoff and curve number (CN). The CN is an index that represents the watershed runoff potential. In the present study GIS based SCS-CN method is used for estimating the runoff depth in the Kharun River Sub-Basin of Sheonath river Basin of Chhattisgarh State of India. The present study reveals that the remote sensing and GIS based SCS-CN can be effectively used to estimate the runoff from the river basins of similar geo-hydrological characteristics. II. KHARUN RIVER SUB-BASIN Kharun River sub-basin a major tributary to Sheonath river in Chhattisgarh State was considered for this study. The study area extends between latitudes 20º32'9.6'' N and 21º39'25''N, and longitudes 81º12'54'' E and 81º58'26'' E. As per GIS total area of Kharun river subbasin is about 4178.33 sq.km. It comprises of Balod (area= 500.19 sq.km.), Dhamtari (area= 593.70898 sq.km.), Raipur (area= 1709.9 sq.km.), Bemetara (area= 327.936 sq.km.)and Durg (area= 1046.6 sq.km.) Districts of Chhattisgarh State. There are about 21 rain gauge stations recording the rainfall data in the study area. Some of the data needed for the study were available from various sources and some of them were procured. The Study area map is shown in Figure 1. The following paragraph gives brief information on the data sources. The Indian Remote Sensing satellite with Linear Imaging Self Scanning sensors (IRS LISS III) satellite data of scale 1:50000 were collected from Bhuvan portal of Indian Space and Research Organization (ISRO), to use land use/ land cover of the study area. Daily rainfall data for all the 21 rain gauge stations from Water Resources Department Chhattisgarh were used. The soil data from National Bureau of Soil Survey & Land Use Planning (NBSS & LUP). III. SCS-CN METHOD Curve Number is basically a coefficient which reduces the rainfall amount. Soil Conservation Services (SCS) CN method is based on two concepts. www.ijaers.com 73
Fig 1. Location Map of Study Area The first concept is that the ratio of actual amount of runoff to maximum potential runoff is equal to the ratio of actual infiltration to the potential maximum retention. This proportionality concept is expressed as ( )/ =/( ) (1) Where P = precipitation in millimeters (when P Q); Q = runoff in millimeters; S = potential maximum retention in millimeters; I a = Initial Abstraction The second concept is that the amount of initial abstraction is some fraction of the potential maximum retention and thus expressed as: Solving equation (1) and (2) we have =( )^2/( + ) (3) For Indian condition I a =0.3S Above equation is used in the estimation of daily runoff from the storm rainfall. Hydrologic Soil Group As per National Engineering Handbook (NEH) developed by USDA, soils are classified in four groups A, B, C and D based upon the infiltration and other characteristics. = (2) Where S = 25400/CN- 254 Group A: Soils in this group have low runoff potential and high infiltration rate when thoroughly wet. Water is transmitted freely through the soil; Group B: Soils in this group have moderately low runoff potential and moderate www.ijaers.com 74
infiltration rate when thoroughly wet. Water transmission through the soil is moderate; Group C: Soils in this group have moderately high runoff potential and low infiltration rate, when thoroughly wet. Water transmission is somewhat restricted through the soil; Group C: Soils in this group have high runoff potential and low very low infiltration rate, when thoroughly wet. Water transmission is restricted through the soil. Antecedent Moisture Condition AMC indicates the moisture content of soil at the beginning of the rainfall event. The AMC is an attempt to account for the variation in curve number in an area under consideration from time to time. Three levels of AMC were documented by SCS AMC I, AMC II & AMC III. The limits of these three AMC classes are based on rainfall magnitude of previous five days and season (dormant season and growing season). AMC for determination of curve number is given in Table 1. Fig 2.Flowchat for computing runoff AMC Table 1.AMC for determination of CN value Total Rain in Previous 5 days Dormant Season Growing Season I Less than 13 mm Less than 36 mm II 13 to 28 mm 36 to 53 mm III More than 28 mm More than 53 mm A conclusion section must be included and should indicate clearly the advantages, limitations, and possible applications of the paper. Although a conclusion may review the main points of the paper, do not replicate the abstract as the conclusion. A conclusion might elaborate on the importance of the work or suggest applications and extensions. IV. METHODOLOGY The methodology adopted in assessing the runoff potential of the study area is explained in the following steps. The same is shown with the help of flowchart given in Figure 2. 1. Preparation of Land use/land cover information of the study area using the satellite imageries in GIS. Land use / Land cover map of the study area is shown in Figure 3. 2. Soil information of the study area obtained is used for making appropriate hydrological soil classification A, B, C & D as shown in Figure 4. 3. Theissen polygons are established for each identified rain gauge station. The weightage of each rain gauge stations are given in Table2. Table 2.Raingauge Stations Weight S.No. Raingauge Station Area_sqkm Weightage 1 Balod 19.5329 0.004675 2 Banbarod 144.056 0.034477 3 Bhalukona 0.194043 0.000046 4 Chandi 141.60899 0.033891 5 Dhamtari 452.12701 0.108208 6 Durg 14.6084 0.003496 7 Gangrel 103.672 0.024812 8 Gunderdehi 187.96899 0.044987 9 Gurur 429.397 0.102768 10 Khapri 129.34 0.030955 11 Kharun_Amdi 289.685 0.069330 12 Kondapar 194.424 0.046532 13 Newara 87.061501 0.020836 14 Oteband 178.08 0.042620 15 Patan 248.50999 0.059476 16 Patharidih 303.95801 0.072746 17 Pindrawan 449.845 0.107662 18 Raipur 382.04901 0.091436 19 Simga 124.366 0.029765 20 Sond 40.871799 0.009782 21 Thanod 256.96899 0.061500 Total 4178.324633 1 www.ijaers.com 75
For each Theissen cell, area weighted CN (AMC II) and also CN (AMC I) and CN (AMC III) were determined. CN for AMC II is given in Table 3. Table 3.Curve Number for HSG under AMC II Conditions Land Use Hydrologic Soil Group A B C D Agriculture Land 76 86 90 93 Buid Up 49 69 79 84 Tree cover 41 55 69 73 Forest 26 40 58 61 Wasteland 71 80 85 88 Water bodies 97 97 97 97 CN for AMC I iscalculated as: = /(2.281 0.01281 ) CN for AMC I iscalculated as: = /(0.427+0.00573 ) Fig 4. Hydrologic Soil Group Map SCS runoff CN for hydrologic soil cover complex under AMC II condition for the study area is given in Table 2. Area weighted composite curve number for various conditions of land use and hydrologic soil conditions are computed as follows: CN=(CN A )+(CN A )+ +(CN " A " )/A Where A 1, A 2, A 3,..., A n represent areas of polygon having CN values CN 1, CN 2, CN 3,..,CNn respectively and A is the total area. Composite curve number for different AMC conditions computed is tabulated in Table 4. Table 4.Raingauge Stations Weight AMC Condition Composite CN AMC I 74.66 AMC II 87.05 AMC III 94.03 4. Using equation (3) with rainfall data, corresponding runoff series is derived Fig 3. Land Use map of Study Area www.ijaers.com 76
Table 5.Sample of Daily Rainfall Runoff Computation of Study Area Day Month Year Daily rainfall (mm) 5-day cumulative rainfall Season AMC Condition Composite Curve number (CN) Surface retention (S) Daily runoff (mm) 15 7 2013 1.571988 67.02837 D AMC III 94.03 16.12656 0 16 7 2013 0 61.58728 D AMC III 94.03 16.12656 0 17 7 2013 0 27.89094 D AMC III 94.03 16.12656 0 18 7 2013 47.38861 112.0837 D AMC III 94.03 16.12656 30.86 19 7 2013 17.38403 112.847 D AMC III 94.03 16.12656 5.49 20 7 2013 1.020052 94.30126 D AMC III 94.03 16.12656 0 21 7 2013 20.8479 100.4585 D AMC III 94.03 16.12656 7.98 22 7 2013 3.975076 90.61567 D AMC III 94.03 16.12656 0 23 7 2013 34.87773 75.02617 D AMC III 94.03 16.12656 19.55 24 7 2013 19.31186 92.73276 D AMC III 94.03 16.12656 6.85 25 7 2013 16.26383 94.818 D AMC III 94.03 16.12656 4.74 26 7 2013 27.0258 109.9985 D AMC III 94.03 16.12656 12.85 27 7 2013 32.02194 129.5012 D AMC III 94.03 16.12656 17.06 28 7 2013 2.11032 96.73375 D AMC III 94.03 16.12656 0 29 7 2013 2.737779 80.15967 D AMC III 94.03 16.12656 0 30 7 2013 58.42269 122.3185 D AMC III 94.03 16.12656 41.19 31 7 2013 121.7108 217.0035 D AMC III 94.03 16.12656 102.70 1 8 2013 1.685128 12.7042 G AMC I 74.66 86.20895 0 2 8 2013 16.96575 29.64545 G AMC I 74.66 86.20895 0 3 8 2013 11.32686 32.57644 G AMC I 74.66 86.20895 0 4 8 2013 1.436419 31.63369 G AMC I 74.66 86.20895 0 5 8 2013 25.82693 57.24109 G AMC I 74.66 86.20895 0 6 8 2013 0.10277 36.05113 G AMC I 74.66 86.20895 0 7 8 2013 0 9.642675 G AMC I 74.66 86.20895 0 8 8 2013 56.87893 56.9817 G AMC I 74.66 86.20895 8.21 9 8 2013 42.53611 99.51781 G AMC I 74.66 86.20895 2.70 10 8 2013 0.283746 99.80155 G AMC I 74.66 86.20895 0 11 8 2013 0.11904 99.81782 G AMC I 74.66 86.20895 0 12 8 2013 0 61.6229 G AMC II 87.05 37.78633 0 13 8 2013 0 48.46627 G AMC II 87.05 37.78633 0 14 8 2013 27.92665 212.9082 G AMC III 94.03 16.12656 13.59 15 8 2013 16.56669 227.3646 G AMC III 94.03 16.12656 4.94 www.ijaers.com 77
V. CONCLUSION GIS based curve number method along with daily rainfall data were used for computing the daily runoff. Antecedent moisture condition plays an important role in the estimation of runoff as it provides the information on moisture content of the land surface for previous the five days rainfall data. Land use layer and soil layer of the study area were prepared and merged in GIS to identify the suitable curve number. Weighted or composite curve number for the study area were calculated and found to be 74.66, 87.05 and 94.03 for AMC I, AMC II and AMC III conditions respectively. Now using equation (3) the daily runoff depth were computed for the year 2013. The sample of runoff computation is shown in Table 5. From this daily runoff event, monthly and yearly runoff can be computed. For those rainfall events whose intensity is less than 0.3 times the surface retention, runoff is taken as zero.monthly runoff generated in Kharun river sub-basin is shown in Table 6. Table 6.Monthly runoff depth for the year 2013 Month Rainfall (mm) Runoff (mm) Jan 0.086 0 Feb 11.59 0 Mar 0.02 0 Apr 9.27 0 May 6.83 0 Jun 251.55 61.5 Jul 589.69 281.86 Aug 458.14 149.26 Sep 167.60 20.66 Oct 187.71 37.76 Nov 0 0 Dec 0 0 Runoff being the important component in planning and management of watershed, its proper quantification is necessary. With the availability of remote sensing data in public domain and GIS, its precise quantification is possible. In The present study reveals that the GIS based SCS-CN method proves to be suitable tool for runoff computation, which helps in proper planning of watershed and its management ACKNOWLEDGEMENTS The authors acknowledge the support provides by Chhattisgarh Council of Science & Technology Raipur, Chhattisgarh & State Data Centre, Water Resources Department Chhattisgarh Raipur, NBLSS & LUP Nagpur. REFERENCES [1] Ahmad I, Verma V, VermaM. K Application of Curve Number Method for Estimation of Runoff. Potential in GIS Environment International Proceedings of Chemical, Biological and Environmental Engineering, Vol. 80: 16-20. [2] Chapter 7, Hydrologic Soils Groups, National Engineering Handbook, National Resources Conservation Services, USDA, May 2007. [3] Jena SK, Tiwari KN, Pandey Ashish, Mishra SK RS and Geographical Information System-Based Evaluation of Distributed and Composite Curve Number Techniques Journal of Hydrologic Engineering, ASCE, Vol. 17, No. 11, November 1, 2012: 1278-1286. [4] Jena SK, Tiwari KN, Pandey Ashish, Mishra SK Runoff Estimation by Distributed Curve Number Technique using Remot Sesnsing and GIS Journal of Indian Water Resources Society, Vol. 30, No. 1, January, 2010: 31-38. [5] Ministry of Agriculture, Govt. of India, Handbook of Hydrology, New Delhi 1972. [6] Nagarajan N, Poongothai Spatial Mapping of Runoff from a Watershed Using SCS-CN Method with Remote Sensing and GIS. Journal of Hydrologic Engineering, ASCE, Vol. 17, No. 11, November 1, 2012: 1268-1277. [7] Reza Kabiri, 2014, Simulation of Runoff using SCS-Cn Method using GIS System, Case Study: Klang Watershed in Malaysia. Research Journal of Environmental Science, 8: 178-192. [8] Seth SM, Kumar Bhism, Thomas T, Jaiswal RK Rainfall- Runoff Modelling for Water Availability Study in Ken River Basin Using SCS-CN Model and Remote Sensing Approach Technical Reports, National Institute of Hydrology, Roorkee, No. CS/AR-12/97-98. [9] SherifMM, Mohamed MM, Sheety Amapr, Almulla M Rainfall-Runoff Modeling of Three Wadis Journal of Hydrologic Engineering, ASCE, Vol. 16, No. 1, January 1, 2011: 10-20. [10] Subramanya K, Engineering Hydrology. Fourth Edition McGraw-Hill Education (India) Private Limited, New Delhi. [11] Xiao Bo, Wang Qing Hai, Fan Jun, Han Feng Peng, Quan Hou Application of the SCS-CN Model to Runoff Estimation in a Small Watershed with High Spatial Heterogeneity Pedosphere, 21 (6), 2011: 738-749. www.ijaers.com 78