1Krish Vigyan Kendra, Tamil Nadu Agricultural University, Tiruchirappalli– 639115, Tamil Nadu, India

2Tamil Nadu Agricultural University, Coimbatore – 641 003, Tamil Nadu, India

3AEC&RI, Tamil Nadu Agricultural University, Coimbatore – 641 003, Tamil Nadu, India

4Institute of Forest Genetics and Tree Breeding, Coimbatore – 641 002, Tamil Nadu, India

Corresponding Author Email: guhanthiran@gmail.com

DOI : https://doi.org

Keywords

Calibration, Sensitivity, Streamflow, SUFI-2, SWAT-CUP, Uncertainty, Validation

Download this article as:

Abstract

For efficient agriculture and flood risk reduction, the Palar sub-basin requires strong water asset management. For hydrological demonstration, A semi-circulated based model, SWAT, was chosen and created in the Palar sub-basin. As part of the validation and calibration process, the Sequential Uncertainty Fitting (SUFI-2) approach was applied to the SWAT CUP model (SWAT-Calibration and Uncertainty Programs).The seven-year observed daily series of inflow data was used to calibrate and validate the SWAT model (2003-2012). The initial three years (2003-2005) were used to set the warm-up period in the SWAT model. The subsequent four years (2006-2009) were used for calibration and the remaining three years (2010-2012) for validation. Two variables, the p-factor, and the r-factor, were used to assess the competence of model validation and calibration. The calculated P-factor and r-factor for validation were 0.84 and 0.24, respectively. The P-factor with high value (0.84) expressed that 84 percent of observed inflow fell within the model uncertainty range (95 PPU). As the r-factor describes the thickness of the uncertainty band, less r-factor of 0.24 indicated the lower predictive uncertainty of the model during the validation stage. The further predictive capability of the model was quantified by PBIAS, RSR, and R2. The presence of high agreement between the observed and simulated stream flow was shown by R2 values of more than 0.8 for calibration (0.81) and validation (0.85). The NS value above 0.7 during calibration (0.75) and validation (0.8) also stated the good performance of the model. PBIAS demonstrated a smaller deviation of simulated values from the observed values by 13.3 and 11.5 % for the calibration and validation phase respectively. The discoveries would be valuable to the hydrological local area, water resource administrators working in Agriculture and soil water protection, policymakers and managing catastrophes like dry spells and floods.

Introduction

Due to the peculiarity of the research elements, the Soil and Water Assessment Tool (SWAT) Model was utilized in this study to understand the hydrological process of the Palar sub-basin. The goal of this study is to test and assess SWAT model use and performance, as well as to model hydrological parameters with judicious accuracy. Hydrological understanding could open the gates of agricultural productivity estimates using Crop Simulation Models as key components. This is extremely helpful in determining the best time for deficit irrigation to save the crop [3]. Thinking about this, the essential objective of this work was to reproduce streamflow in the Palar sub-basin utilizing the SWAT model, which was joined with model alignment and vulnerability investigation utilizing the SUFI-2 methodology, and to assess its reasonableness for the Palar sub-basin. This modeling work also helps water resource managers plan and manage agricultural water resources more effectively.

Materials and Methods

2.1 Study area

Palar sub-basin is one among the sub-basins of the Parambikulam Aliyar Project  (PAP) basin, which is an inter-state water distribution project located in the southwestern part of the peninsular India as a collaboration between Kerala and Tamil Nadu States. This project diverts the water from the eight rivers. The PAP basin area lies in coordinates of 10° 10’ 00” N to 10°57’20” N latitude, 76°43’00” E to 77° 12’30” E longitudes and covers 2388.72 sq. km area (Fig.1).  .In this study Palar sub-basin which was considered to be the major sub-basin of PAP basin was taken for SWAT simulations.

Figure 1. The study area of ParambikulamAliyar basin

Climate data

Data from ECMWF reanalysis Interim (ERA-I) with 0.75° horizontal resolution was used to create the starting and boundary conditions for the Weather Research and Forecasting model (WRF) [2]. To adjust for bias, the observed data was compared to the data from the simulation. Bias-corrected data was used for this study.

SWAT model

The Schematic diagram of the Input and Output files of the SWAT model is described in Fig 2.

Description of datasets used to set up the SWAT model

It was ArcSWAT 2012 that we employed in this research [6]. Soil, land use, and climate are all factors that must be taken into account while constructing the SWAT model. For the SWAT model, the Digital Elevation Model (DEM) is essential because it offers digital topographic information that helps with stream flow and watershed delineation (Fig 3.3). Shuttle Radar Topographic Mission’s Digital Elevation Model of 1 Arc Second (30m) was utilized in the investigation[6]. The Parambikulam-Aliyar basin’s SRTM-DEM elevations ranged from 140 to 2498 meters.

Soil data

The soils of the Tamil Nadu section of the basin were defined using a digitized soil map of Tamil Nadu at a scale of 1:50,000 acquired from the Department of Remote Sensing and Geographical Information System, Tamil Nadu Agricultural University (TNAU). 

Data on land usage and land cover

The National Remote Sensing Centre provided data on land use/land cover for the year 2011. (NRSC). The NRSC created this using multi-temporal Resourcesat-1 AWiFS (Advanced Wide Field Sensor) data.

SWAT model setup

Through the GIS interface ArcSWAT2012, the streams of the Palar sub-basins catchment area were delineated using SRTM DEM and a digitised stream network from Google Earth. The Palar sub-basin was split into three sub-basins and three hydrological response units.

Figure 2.  Schematic diagram of the Input and Output files of the SWAT model

2.10 Calibration and validation of SWAT model

SWAT model calibration and validation were performed using the seven years observed daily series of inflow data (2003-2012). The initial three years (2003-2005) were used to set a warm-up period in the SWAT model. The subsequent four years (2006-2009) were used for calibration (Fig 2.) and the remaining three years (2010-2012) for validation (Fig 3.).SWAT-CUP is an automated model calibration tool for the SWAT model. It stands for Calibration and Uncertainty Programs. This tool was created by the Eawag Aquatic Research Institute in Switzerland [1]. SWAT-CUP is a freeware application with a generic interface. Within SWAT-CUP, several sensitivity analyses, calibrations, validations, and uncertainty analyses are possible. SWAT-CUP includes five distinct uncertainty algorithms (SUFI-2, PSO, MCMC, ParaSol, and GLUE) [1]. The unknown model parameters are systematically modified in SWAT-CUP, and the model is run. The appropriate outputs from the model output files are then extracted and compared to the observed data.SWAT-CUP was calibrated and validated using the Sequential Uncertainty Fitting (SUFI-2) method (Calibration and Uncertainty Programs). SWATCUP was calibrated by using p-factor and r-factor, which reflect model prediction uncertainty (95 PPU) as well as r-factor to determine uncertainty band thickness  (PBIAS), Nash-Sutcliffe efficiency (NSE) [5], and Coefficient of determination (R2) were also calculated using the equation from 1,2,3.

The percent BIAS was computed using the following formula — (1)

The Nash-Sutcliffe Efficiency (NSE) was computed using the following formula —(2)

The coefficient of determination (R2) was computed using the following formula — (3)

3. Result and Discussion

3.1 Sensitive analysis of parameters in SWAT cup in Palar sub-basin of PAP basin

For SWAT model calibration 15 parameters were considered based on the sensitivity levels tested in SWAT-CUP.  The initial minimum and maximum values, methods used in parameter change and fitted value of individual parameters during calibration for Palarsub-basin are furnished in Table 1.

Table 1. Sensitive hydrological parameters included for calibrating the SWAT model in Palar sub-basin

Parameters nameDescriptionMethod of changeMinimumMaximumFitted value
CN2.mgtThe soil water state is represented by a curve number.Relative-0.120.410.09
ALPHA_BF.gwThe alpha factor for base flow (days)Replace0.0360.530.19
GW_DELAY.gw   Delay in groundwater (days) Replace45465409
GWQMN.gw   Return flow requires a minimum depth of water in the shallow aquifer (mm).Replace0.080.310.16
ESCO.bsnCompensation factor for soil evaporationReplace0.150.810.64
OV_N.hruFor overland flow, Manning’s n valueReplace-0.09-0.18-0.06
EPCO.hruCompensation factor for plant uptakeReplace0.350.810.45
REVAPMN.gwThe minimum depth of water in the shallow aquifer required for “recap” to take place (mm).Replace1811450308
CH_K2.rteAlluvium in the main channel has a high effective hydraulic conductivity.Replace18.5108.123.1
CH_N2.rteFor the primary channel, Manning’s “n” value.Replace0.150.630.53
MSK_CO1.bsnThe calibration coefficient is used to manage the storage time’s influence. constant for normal flow.Replace0.856.13.56
SURLAG.bsnlag time for surface runoffReplace0.8315.83.9
GW_REVAP.gwCoefficient of “revap” in groundwaterReplace0.0150.280.23
RCHRG_DP.gwPercolation fraction in deep aquifersReplace0.050.450.33
SOL_AWC(1).solThe initial soil layer’s available water capacity (mm/mm)Replace-0.150.180.15

Table 2.  SWAT model statistics for calibration and validation of streamflow in Palar sub-basin

ParametersP-factorr-factorR2NSPBIASRSR
Calibration0.780.260.810.7513.30.6
Validation0.840.240.850.811.50.45

The calculated P-factor and r-factor for validation were 0.84 and 0.24, respectively. The strong P-factor (0.84) meant that 84% of the measured inflow was within the model’s uncertainty range (95 PPU). Because the r-factor defines the thickness of the uncertainty band, a lower r-factor of 0.24 indicates that the model’s prediction uncertainty during the validation stage is smaller.

The objective function used for parameterization was the Nash Sutcliffe coefficient (NS) [5]. The further predictive capability of the model was quantified by PBIAS, RSR and R2. The R2 value was more than 0.8 for calibration (0.81) and validation (0.85) indicating the existence of good agreement between the observed and simulated stream flow (Table 2). The model’s high performance was further demonstrated by NS values over 0.7 during calibration (0.75) and validation (0.8).PBIAS demonstrated a smaller deviation of simulated values from the observed values by 13.3 and 11.5 % for the calibration and validation phase respectively.

Fig . 2. Calibration of simulated stream flow with the observed stream flow (2006-2009) in Thirumoorthy reservoir in the Palar sub-basin

Fig . 3. Validation of simulated stream flow with the observed stream flow (2010-2012) in Thirumoorthy reservoir in the Palar sub-basin

Conclusion

Rigorous model calibration is essential for hydrological prediction, such as discharge, to achieve an efficient outcome. It is necessary to report the uncertainty in the model forecast along with the outcomes for appropriate modeling practice. As part of this work, the SWAT model was used to simulate streamflow in the Palar sub-basin of the PAP basin, after a thorough calibration and validation method. When it comes to forecasting water flow and quantifying uncertainties and related assumptions, SWAT CUP has been demonstrated to be a significant tool. As per the consequences of the alignment and approval, the model precisely duplicated the noticed streamflow. The discoveries of this study might be helpful in water management in agriculture, soil and water preservation, and lessening catastrophic events like dry spells and floods, among different applications.

References

  1. Abbaspour, Karim C, ElhamRouholahnejad, Srinivasan Vaghefi, Raghavan Srinivasan, Hong Yang, and BjørnKløve. 2015. A continental-scale hydrology and water quality model for Europe: Calibration and uncertainty of a high-resolution large-scale SWAT model.  Journal of Hydrology 524:733-752.
  2. Dee, Dick P, SM Uppala, AJ Simmons, Paul Berrisford, P Poli, S Kobayashi, U Andrae, MA Balmaseda, G Balsamo, and d P Bauer. 2011. The ERA‐Interim reanalysis: Configuration and performance of the data assimilation system.  Quarterly Journal of the royal meteorological society 137 (656):553-597.
  3. Farahani, Hamid J, Gabriella Izzi, and Theib Y Oweis. 2009. Parameterization and evaluation of the AquaCrop model for full and deficit irrigated cotton.  Agronomy Journal 101 (3):469-476.
  4. https://earthexplorer.usgs.gov/
  5. Nash and Sutcliffe, 1970Gurtz, Joachim, AndriBaltensweiler, and Herbert Lang. 1999.Spatially distributed hydrotope‐based modeling of evapotranspiration and runoff in mountainous basins.Hydrological Processes 13 (17):2751-2768.
  6. Winchell, M, R Srinivasan, M Diluzio, and J Arnold. 2013. Arcswat Interface for Swat 2012: User Guide.  Blackland Research Center, Texas AgriLife Research.

Similar Posts