Sedimentation Yield Analysis of Kamthikhairy Dam Reservoir Using SWAT and SWAT-CUP Algorithm
Sedimentation Yield Analysis of Kamthikhairy Dam Reservoir Using SWAT and SWAT-CUP Algorithm
https://doi.org/10.22214/ijraset.2021.32920
International Journal for Research in Applied Science & Engineering Technology (IJRASET)
ISSN: 2321-9653; IC Value: 45.98; SJ Impact Factor: 7.429
Volume 9 Issue I Jan 2021- Available at www.ijraset.com
Abstract: Sedimentation is an immense problem that has threatened water resources development in the Kamthikhairy dam
reservoir, Nagpur, Maharashtra. Sediment yield is becoming a major problem now days which leads to siltation in reservoir and
not only reduces the capacity of dam but also can affect the power generation and safety of dam as well as aquatic life. In
current scenario all the dams in India are facing this problem. An insight into sedimentation mechanisms plays an important
role for managing the reservoir operations efficiently in the governing region. The objective of the study is to analyze the
sedimentation of the reservoir using Remote sensing and Soil and Water assessment Tool (SWAT). The SWAT can be the
powerful tool to analyze the sedimentation in the reservoir. The calibration using SWAT-CUP algorithm will help to make the
simulation more accurate to the observe results.
Keywords: ArcGIS, Remoe Sensing, SWAT, SWAT-CUP
I. INTRODUCTION
The sedimentation deposition of reservoirs in India has become a serious problem due to the increasing population, agricultural
activities, poor soil conservation, and deforestation. Sediment deposition creates many problems in the reservoirs. The
accumulation of sediments in reservoir has led to a decline in lake water storage capacity, navigation problems, and other water
quality and environmental issues.
To overcome the problems caused by the accumulation of reservoir sediments and to use these sediments, it is important to estimate
the sediment volume (quantity) and its distribution in these reservoirs. This is not only used to assess the capacity of reservoirs but
also to ensure the safest route of navigation through these reservoirs. The assessment of reservoir sediment deposits is important for
the operational planning and management of these reservoirs.
RS data has proven to be a useful tool for monitoring reservoir sediment deposits. Besides, GIS analysis tools are used to combine
satellite observations with field measurements to assess sediment deposits in reservoirs and to successfully apply them in different
parts of the world. Therefore, this study aims to estimate the total amount of accumulated sediment and its distribution on the
reservoir bed using RS and GIS techniques.
Typically, the reservoir is designed to accommodate the estimated settling rate. Remote sensing and GIS technology is a new
method for assessing sediment yields that provide high-precision results and spatial data over time.
Figure 1: Digital Elevation Map and Land Use and Land Cover Map
Digital Elevation Model is the geographic grid of an area where the contents of each grit cell describe the elevation of any point at a
given location and specific spatial resolution in the form of a digital form. It is one of the essential spatial inputs necessary for
delineation of the watershed into several sub-basin based on elevation in the SWAT model.
The SRTM DEM of 30m resolution was processed for the extraction of flow direction, flow accumulation, stream network
generation, and delineation of the watershed and sub-basins. The topographic parameters such as terrain slope, channel slope, or
reach length were derived from the DEM.
D. Soil Data
The soil map was obtained mainly from the FAO (Food and Agricultural Organization). The soils are compared with the watershed
area to use their properties to define HRUs. The soil associations are characterized in terms of the estimated proportions of the
component soil units, the presence of phases, and the estimated proportions of these soil units in one of 5 texture-slope classes such
as Shallow coarse soil: 8.60%, medium black soil: 5.65%, deep black soil: 2.87%, clayey medium black soil: 7.06%, clayey Lumi:
74.90%.
SWAT (Soil and Water Assessment Tool) is a model developed by the United States Department for Agriculture, Agricultural
Research Service (USDA-ARS) and it is a river basin scale model developed to quantify the impact of land management practices in
large, complex catchments (Neitsch et al., 2001). SWAT is a deterministic, physically-based, continuous simulation, watershed
model. The SWAT was developed principally to compute runoff and loadings from rural especially agriculture-dominated
watersheds and predicts the long term impacts of the land management practices on water, sediment in watersheds with varying
soils and land cover.
The semi-distributed model divides the watershed into multiple sub-basins up to ten thousand which are called hydrologic response
units (HRUs). The amount of water, sediment, nutrient, and pesticide loadings to the main channel from each sub-catchment is
controlled by the land phase of the hydrologic cycle. Aggregation of input data can be facilitated using a GIS interface (either
Geographic Resources Analysis Support System (GRASS) or Arc View) available with the model. The Arc View interface
developed for the SWAT model requires land use/land cover, soils, and DEM data as spatial inputs and can be used to prepare the
input files. SWAT computes the volume of surface runoff using the SCS curve number method or Green and Ampt infiltration
method; while the amended Rational Formula or the SCS TR-55 is used to count the peak runoff rate.
The SWAT model has become widely accepted as an effective tool to predict the effects of watershed management on runoff,
sediment, nutrients, etc. SWAT considers the base-flow contribution to total streamflow by routing a shallow aquifer storage
component to the stream.
F. Surface Runoff
The estimation of surface runoff can be performed by the model using the Soil Conservation Service (SCS) curve number method,
(S.M.Ghoraba) This method is widely used for the prediction of the approximate amount of runoff from a given rainfall event.
CN1 = CN2-( 20(100− 2) (100− 2+ (2.533−0.0636∗(100− 2)))
CN3 = CN2*(0.0063(100− ) (4)
Where;
CN1 is the moisture condition I curve number,
CN2 is the moisture condition II curve number, and
CN3 is the moisture condition III curve number.
G. Model Setup
1) Watershed Delineation by SWAT: Click on Watershed Delineator Automatic Watershed delineation, Once the watershed
delineation window is activated, click on the open file button in the DEM setup portion, and then click on select from map, and
then select the project DEM. Click OK. You will see that the name of the DEM will change to Source DEM. Then proceed to
create flow direction and flow accumulation grids. By leaving the default DEM-based option unchanged in the stream definition
portion, click on the flow direction and flow accumulation button. After flow direction and flow accumulation are complete, the
area box will be populated by some numbers. This is the critical stream area threshold that will be used to define the stream
network. Steps:
Click the create stream and outlets button to create a stream network, sub-basins, and outlets for each of the sub-basin. Next, in the
watershed outlets selection and definition, click on the select whole watershed outlet button, and then select the desired outlet
location using the cursor. If the outlet does not exist at the point you are interested in, use the add button in edit manually frame to
add the point, and then do the selection using the outlet button. After the outlet is selected, the Delineate watershed button will be
activated in green. After the watershed is delineated, the polygon feature class with sub-basins is added to the map document. Click
on the calculate sub-basin parameters button to populate the watershed feature class with necessary attributes for the SWAT model.
After all, attributes are calculated, click OK on the final message box. Now open the attribute table and you will see that it has more
attributes as shown below. The delineating watershed and creating sub-basins for a specific threshold table is done. Exit the
watershed delineation window by clicking on EXIT.
After processing land use data, click on the Soil data tab. After that click on the open soils grid button, land soils dataset from disk,
and browse to the clipped projected soil raster that was just created, and click Select. This step will input the raster and then clip that
raster to the watershed area. This process sometimes takes a long time depending on the size of the watershed. After that the soil
data are loaded, choose VALUE as the grid field and click OK. This will add all numerical values present in the land soil dataset to
the SWAT soil classification table. The next step is to link these values to a certain soil type. To use the lookup table, it is necessary
to define the attribute in the lookup table that corresponds to the soil type. Finally, click on reclassify which will finish the soil data
processing for the HRU analysis. The next step is to assign slope attributes to each HRU. Click on the slope tab. To keep the process
simple, use the single slope option, and click to reclassify. This will complete the processing of land use, soil, and slope data for
HRU analysis. Close the land use/ soil/ slope definition window. Leave the default options of multiple HRUs in the HRU definition.
Several hours will be added to the map document that would be described in the HRU report generated by this process.
Other input files after generating weather database, create all other input files by selecting write input table, and write SWAT input
tables. A new window will be displayed known as ‘Write SWAT Database Tables’. Then click on select all and click on create
tables. Also, choose YES in all the subsequent messages which will be showing up on the screen.
J. SWAT Simulation
To run SWAT Simulation click on SWAT Simulation and run SWAT output to set the period of simulation from 01/05/2003 to
31/12/2013 and change the printout settings to daily, and leave the options unchanged. After getting the finished SWAT setup
message, click on Run SWAT. While the SWAT is running, a DOS window is shown on the screen for the progress, and after the
simulation is done, the final message of the SWAT simulation is done is displayed on the screen. After the successful simulation
message, close the SWAT Simulation setup window and read and visualize the model output.
The Coefficient of determination R2 can be determined from the above equation and the value can range from 0 to 1. The value 1
shows complete correlation and a value near 1 can be considered as a good correlation. The graph below shows a correlation of
about 0.81. The value represents a good degree of correlation, but to improve that we can calibrate the output of the SWAT model
using the SUFI2 algorithm of SWAT-CUP.
Fig: The graph of observed VS SWAT simulated sediment with a coefficient of determination of 0.81
The sediment derived from the SWAT is calibrated using SWAT-CUP for different parameters which are from.SOL group. The.
SOL group of SWAT parameters represent the soil-related parameters like Soil Hydrologic group, Maximum rooting depth of soil
profile (mm), Soil CaCo3, Soil ph-value, Electrical conductivity, USLE equation soil erodibility (K) factor, Moist Soil Albedo, Clay
content(%s soil weight), Saturated hydraulic conductivity (mm/hr), Depth of soil surface to the bottom of layer (mm), Moist bulk
density (Mg/cu.M) for different maximum and minimum value to match the output values of the sediment derived from different
observed data and parameters. The coefficient of determination after calibration shows around 0.9469 shows a very good
correlation.
Values
Parameters Description
Min Max
r__CN2.mgt Initial SCS runoff curve number for moisture condition II -0.2 100
v__ALPHA_BF.gw Baseflow Alpha Factor (1/day) 0 1
v__GW_DELAY.gw The Delay time 30 450
Threshold depth of water in the shallow aquifer required for return
v__GWQMN.gw 0 2
flow to occur (mm H2O)
v__SNAM.sol Soil Name 0.1 100
r__HYDGRP.sol Soil Hydrologic group 0 100
r__SOL_ZMX.sol Maximum rooting depth of soil profile (mm) 0 100
r__SOL_CAL().sol Soil CaCo3 0 50
r__SOL_PH().sol Soil PH value 3 10
r__SOL_EC().sol Electrical conductivity 0 100
r__USLE_K().sol USLE equation soil erodibility (K) factor 0 1000
r__SOL_ALB().sol Moist Soil Albedo 0 1
r__SOL_CLAY().sol Clay content(%s soil weight) 0 100
r__SOL_K().sol Saturated hydraulic conductivity (mm/hr) 0 100
r__SOL_Z().sol Depth of soil surface to the bottom of the layer (mm) 0 100
r__SOL_BD().sol Moist bulk density (Mg/cu.M) 0 100
After getting the desired output, the calibration of the model is important for certain parameters like soil and land-use classes. The
sedimentation is majorly dependent upon the SOL parameters and RCH. For calibration of the model, SWAT-CUP is used. The
SWAT-CUP is equipped with the four algorithms for calibration i.e. SUFI2, GLUE, ParaSOL, and PSO. From the various kinds of
literature, the SUFI2 is suitable for analyzing and calibrating the sediment parameters. For calibrating the sediment value of Pench
dam the following parameters are used. The following table shows the minimum and maximum thresholds for all the factors along
with a description. Table: List of SWAT Model parameters for calibration along with a description and maximum and minimum
values for calibration in SUFI2 algorithm. The upper and the lower limit of each parameter are referred from the SWAT handbooks
for each parameter. The parameters which are responsible for the sediments are most of them from. Sol group i.e. soil group and it is
important to calibrate the soil parameters. The parameters like rooting depth of soil profile (SOL_ZMX) and soil surface depth
(SOL_Z) are the major parameters responsible and important to calibrate for assessing the soil and sediment in the reservoir. The
CN value of moisture zone II is also important to calibrate as the land-use data is responsible for soil fragment and looseness
The graph below shows the final sediment values including the data for calibration and validation. The sudden increase in the
sediment can be seen after the monsoon of 2011 which represents the full working condition of the Pench dam, as the year 2011 is
the year of attaining the fully functional level of Pench Dam. The year 2011 shows the actual sediment deposition and the data can
be seen increased after the year 2011. The study is important to analyze the change in the reservoir water capacity as the sediment
may reduce the capacity after saturating the dead storage.
IV. CONCLUSION
The SWAT model and SWAT-CUP algorithm model is developed and applied to case study of Kamthikhairy Dam in Pench river
basin in Maharashtra state, India. SWAT model calculated sediment yield data based on spatial datasets (DEM, Land cover
changes, and soil map) and weather data. This sediment yield data obtained from the SWAT analysis was used as target data in
SWAT-CUP model to validate the SWAT output data
The conclusions from SWAT analysis for the present study are given below:
A. It was observed that, the rate of sediment yield varied between 0.071 MCM/year in 1985 to 2.22 MCM/year in 2014
B. The Average rate of sediment yield was found to be 0.51 MCM / year from SWAT model.
C. For the duration between 1982 to 2014 almost 8.03% of the gross storage of dam was reduced due to problem of sediment yield.
The conclusions from SWAT-CUP for the present study are given below:
SWAT-CUP is a program for calibration of SWAT models. The program could be used to perform calibration, validation,
sensitivity analysis (one-at-a-time, and global) and uncertainty analysis. The program links SUFI2, GLUE, ParaSol, MCMC, and
PSO to SWAT. Any of the procedures could be used to perform calibration and uncertainty analysis of a SWAT model.
SWAT-CUP also has graphical modules to observe simulation results, uncertainty range, sensitivity graphs, watershed visualization
using Bing map, and statistical reports
REFERRENCES
[1] Bhere, S., & Patil, G. K. (2020). The Climate Change Analysis on Water Resources over Western Ghat by Deficit Calculation Method. 7(1), 16–35.
[2] Bonell, M., Purandara, B. K., Venkatesh, B., Krishnaswamy, J., Acharya, H. A. K., Singh, U. V., Jayakumar, R., & Chappell, N. (2010). The impact of forest
use and reforestation on soil hydraulic conductivity in the Western Ghats of India: Implications for surface and sub-surface hydrology. Journal of Hydrology,
391(1–2), 47–62. https://doi.org/10.1016/j.jhydrol.2010.07.004
[3] Butt, A., Shabbir, R., Ahmad, S. S., & Aziz, N. (2015). Land use change mapping and analysis using Remote Sensing and GIS: A case study of Simly
watershed, Islamabad, Pakistan. Egyptian Journal of Remote Sensing and Space Science, 18(2), 251–259. https://doi.org/10.1016/j.ejrs.2015.07.003
[4] Chakraborty, D., Reddy, M., Tiwari, S., & Umapathy, G. (2019). Land Use Change Increases Wildlife Parasite Diversity in Anamalai Hills, Western Ghats,
India. Scientific Reports, 9(1), 1–11. https://doi.org/10.1038/s41598-019-48325-8
[5] Chen, X. L., Zhao, H. M., Li, P. X., & Yin, Z. Y. (2006). Remote sensing image-based analysis of the relationship between urban heat island and land use/cover
changes. Remote Sensing of Environment, 104(2), 133–146. https://doi.org/10.1016/j.rse.2005.11.016
[6] Cohen, W. B. (1994). GIS Applications Perspective: Current Reseach on Remote Sensing of Forest Structure. Remote Sensing and GIS in Ecosystem
Management, 1990, 91–107.
[7] Freeman. (2013). 済無No Title No Title. Journal of Chemical Information and Modeling, 53(9), 1689–1699. https://doi.org/10.1017/CBO9781107415324.004
[8] Gambo, J., Mohd Shafri, H. Z., Shaharum, N. S. N., Abidin, F. A. Z., & Rahman, M. T. A. (2018). Monitoring and Predicting Land Use-Land Cover (Lulc)
Changes Within and Around Krau Wildlife Reserve (Kwr) Protected Area in Malaysia Using Multi-Temporal Landsat Data. Geoplanning: Journal of
Geomatics and Planning, 5(1), 17. https://doi.org/10.14710/geoplaning.5.1.17-34
[9] Giri, C., Zhu, Z., & Reed, B. (2005). A comparative analysis of the Global Land Cover 2000 and MODIS land cover data sets. Remote Sensing of
Environment, 94(1), 123–132. https://doi.org/10.1016/j.rse.2004.09.005
[10] Nayak, C. (2008). Comparing Various Fractal Models for Analysing Vegetation Cover Types at Different Resolutions with the change in Altitude and Season
Comparing Various Fractal Models for Analysing Vegetation Cover Types at Different Resolutions with the change in Altitud.