In terms of the land and climate characteristics of the southern Uttar Pradesh Basin, tendencies toward sudden floods and soil erosion present some of the most extreme issues. Because of soil erosion, the soil of this watershed is continuously degrading and losing nutrients, thus affecting agriculture in the state. The severity of the issue is intensified in arid and semi-arid lands, where short-duration intense rainfall and unsustainable land use have quickened soil losses by erosion. Detailed and accurate information regarding soil erosion and surface water discharge is helpful for watershed managers to better manage and conserve natural resources such    as soil and water, and to promote sustainable development. At present, various procedures are used for hydrological modeling, utilizing precipitation, land use, and soil characteristics information. Among these are the Water Evaluation and Planning (WEAP) system, Agricultural Non-Point Source Pollution (AGNPS) system, Areal Nonpoint Source Watershed Environment Response Simulation (ANSWERS), Soil and Water Assessment Tool (SWAT), and Water Erosion Prediction Project (WEPP). Among these strategies, the SWAT is a process-based hydrological model developed by the United States Department of Agriculture (USDA) Agricultural Research Service (ARS) [1]. Additional algorithms have been developed to best estimate the parameters used in hydrological modeling.

 The SWAT is widely used for hydrological and sediment yield modeling. It is a time-consistent and spatially appropriate test model that was developed to assist watershed managers in anticipating the effects of land use management activities on runoff, soil erosion, and agricultural chemical yields [2]. Specialists have effectively used the SWAT for runoff assessments [3,4], water quality modeling [5,6], hydrological and river basin modeling [7–9], sediment yield modeling [10–13], and the management of erosion-prone areas [1,14]. Psomas et al. [15] utilized the SWAT and WEAP to develop water efficiency measures in Greece.

To modify the elements that influence SWAT yields, calibration using observed values and evaluated estimations of runoff, evapotranspiration, and  other SWAT outputs is necessary. Validation is  a way that compares SWAT results with the observed data, without modifying the values of the influencing factors. Calibration of various parameters is necessary for proper hydrological modeling. There are 26 parameters distinguished for runoff, more than 30   for soil disintegration, and 41 for water quality, all of which can     be utilized for calibration. The calibration of such a circulated parameterized watershed model entails some significant issues  that require careful attention on the part of investigators, particularly in regard to uncertainty [9]. Qiu and Wang [5] performed validation and calibration for the periods 1997–2002 and 2003–2008, respectively, for the hydrological and water quality evaluation of central New Jersey. Vigiak et al. [16] performed calibration and validation for hillslope erosion modeling of the Danube Basin.

Earlier researchers applied deterministic approaches such as the trial and error method for uncertainty  analysis,  calibration, and validation. When using such methods, it is necessary to continue to alter the parameters until a sensible match is achieved between the simulation and observations. However, these approaches are now outdated; scientists have developed many stochastic algorithms for calibration, validation, and uncertainty analysis. Some of the widely used algorithms for calibration are sequential uncertainty fitting version 2 (SUFI-2), parallel solution (ParaSol), particle swarm optimization (PSO), generalized likelihood uncertainty estimation (GLUE), artificial neural network (ANN), and Markov chain Monte Carlo (MCMC).

Vilaysane et al. [17] used the SUFI-2 technique to calibrate the values for hydrological stream flow modeling of the Xedone River Basin. Yesuf et al. [12] utilized SUFI-2 and SWAT-CUP for sediment yield modeling in Ethiopia; they reported that SUFI-2 utilized a computerized advancement system that was equipped to perform affectability, alignment, approval, and vulnerability examinations with enhanced model run time proficiency. Ercan et al. [18] used cloud techniques for SWAT model calibration. Talebizadeh et al. [19] used ANN and the SWAT for sediment load modeling and uncertainty analysis. Calibration and validation have also been done using genetic algorithms and Bayesian  model  averaging  [20]. Tuo et al. [21] applied the SUFI-2 algorithm for uncertainty analysis when evaluating precipitation input for the SWAT model. Noori and Kalin [4] coupled ANN with the SWAT and SWAT-CUP for daily streamflow prediction. Similarly, runoff [3] and sediment fluxes modeling [13] have been done using the SUFI-2 algorithm. An approximation of the SWAT model can also be constructed using ANN and a support vector machine (SVM) [22]. Out of all these algorithms, it is difficult to  determine  the  best  algorithm  for the calibration of SWAT outputs. Therefore, in the present study, we consider the three calibration-uncertainty analysis algorithms of GLUE, SUFI-2, and ParaSol.

In the present study, SWAT 2012 was used for hydrological modeling. For calibration, the three widely used algorithms of GLUE, SUFI-2, and ParaSol were applied and compared in order   to determine the best calibration technique. In this investigation, the ArcSWAT Interface was associated with the Ganga River Basin for depicting and organizing sub-watersheds in order to gauge and model runoff, sediment yield, and evapotranspiration. The model required a digital elevation model (DEM), a land use/land cover (LULC) map, and a soil map as inputs for watershed delineation and hydrological response unit (HRU) analysis. Daily meteorological data from the four stations of Balia, Varanasi, Babatpur, and Gazipur were utilized from 1996 to 2015. The major target of this examination was to delineate the watersheds, divide  the  study  area into sub-watersheds, determine the number of streams, and estimate their order, and then to establish information regarding the outlet points and reservoirs available in the watershed in order to accurately divide it into HRUs having similar but unique land types, soil types, and elevation properties. These units permit simpler modeling. Finally, after providing the meteorological data, the SWAT was run to estimate runoff, sediment yield, and evapotranspiration for each HRU and sub-watershed. The model was then calibrated for the observed data from the years 2004–2009 at the sub-watershed and watershed level using the SUFI-2, GLUE, and ParaSol algorithms. By analyzing the results, assumptions, and limitations of these calibration techniques, the best technique among the three was determined based on five categories, as described in Section 4.2. The model calibrated using the best technique  was then validated for the period of 2009–2015. Finally, the estimated and corrected results of the hydrological parameters were calculated, and modeling was conducted.

2. Study area

The study area for which modeling was performed in the present study is part of the Ganga River Basin that occupies the southern of Uttar Pradesh State of India. The length of the river in this watershed is about 50 km, and the total area of the watershed is 15 621.612 km2.  The  area  lies  between  latitude  82°1052.43900 E  to 83°55010.63 E   and   longitude   26°2'7.842'' N   to   24°22'53.034'' N. The major stations covered by the study area are the Varanasi, Balia, Babatpur, Gazipur, Mirzapur, and Chandauli Districts. The mainstream flow of the study area lies in  the  Varanasi  District. The average rainfall in this area is 941.2 mm, with the maximum rainfall occurring in July, August, and September and the minimum rainfall occurring in June and October. From November to May, the study area receives negligible rainfall. This part of the basin is barren, and urban lands cover more than 50% of the basin.

3. Data

Elevation data, LULC data, soil data, and daily meteorological data were the prerequisites for the present modeling. The Shuttle Radar  Topology  Mission  (SRTM)  DEM,  with  a  resolution  of   90 m 90 m, was used for elevation data. For the LULC map, image classification  was performed using  satellite imagery from  Landsat 8. For meteorological input, data for the daily rainfall, temperature, solar radiation, and pressure for 20 years were used. Detailed information on the data used and on their locations, periods, and the organizations from which they were procured is given in Table 1.

Table 1 Details of the raw data used for modeling.

IMD: Indian Meteorological Department; USGS: United States Geological Survey; NBSSLUP: the National Bureau of Soil Survey and Land Use Planning.

4. Methods and procedures

4.1. Soil and Water Assessment Tool (SWAT)

In this study, SWAT 2012 was utilized for hydrological modeling. ArcGIS 10 software was used along with its extension, ArcSWAT. The catchment was split into various sub-catchments, which were then further categorized into HRUs. Each HRU comprises a unique combination of soil characteristics, elevation, and LULC. The water budget was the primary impetus behind all forms, and the HRU was the unit used to estimate the hydrologic parameters. These procedures were separated into two stages: ① A land stage, wherein the SWAT mimics the catchment loadings of the stream, residue, and supplements from each HRU, which are then locale-weighted to the sub-basin level; and ② an in-stream stage, wherein the model tracks the course of the catchment loadings from each sub-basin throughout the channel organization.

A solitary plant-improvement model was utilized as a part of  the SWAT in order to reenact extraordinary land-cover types including expansive information yield. This development model was utilized to evaluate the expulsion of water and supplements from the root zone, transpiration, and biomass generation. Planting, reaping, culturing passes, and supplement and pesticide applications could be recreated for each  trimming  framework  with particular dates [23]. Once the data for each HRU were determined at the sub-basin level, the runoff, sediments, nutrients, and pesticides were routed through channels, ponds, reservoirs, and wetlands to the watershed outlet.

The SWAT theoretical documentation version 2009 by Neitsch et al. [24] can be used as a reference to learn more about SWAT modeling.

4.2. Sequential uncertainty fitting version 2 (SUFI-2)

The SUFI-2 technique is a stochastic algorithmic approach that is most frequently used by scientists to evaluate uncertainty. In this algorithm, the extent to which all uncertainties are represented is assessed by a parameter referred to as the P-factor. This parameter is the percentage of estimated information sectioned by the 95% prediction uncertainty, which is also called the 95PPU. Another parameter evaluating the quality of the calibration is the R-factor, which is the standard thickness of the 95PPU band separated by the standard deviation of the predetermined information. Hypothetically, the estimation of the P-factor ranges from 0 to 100%, while that of the R-factor ranges from 0 without limit. A P-factor of 1 and an R-factor of 0 signify a recreation that precisely compares with the predetermined information. In this study, SWAT-CUP was the software used to consolidate the SWAT 2012 simulated values with the observed values. Using these values and the SUFI-2 algorithm, the uncertainty analysis and calibration were performed.

 A short well-ordered depiction of the SUFI-2 algorithm is as follows:

Step 1: The objective function (g_i) is characterized. Next, the minimum and maximum absolute ranges (θ_j) of the physically important parameters being optimized are identified.

Step 2: A sensitivity analysis for each of the parameters is completed; afterward, the initial uncertainty ranges are relegated to the parameters for the first round of Latin hypercube testing.

Step 3: Latin hypercube sampling is performed, and the corresponding objective functions are assessed. The sensitivity matrix J_ij and the parameter covariance grid C are calculated by the following:

where C_m is the number of rows in the sensitivity matrix, and p is the number of parameters to be estimated.

where S_g² is the variance of the objective function values resulting from the m model runs.

Step 4: The 95PPU is then calculated, followed by the P-factor and R-factor.

4.3. Generalized likelihood uncertainty estimation (GLUE)

The GLUE was applied in order to consider the non-uniqueness of the parameter sets during the estimation of model parameters in over-parameterized models. This method is simple; for vastly overparameterized models, GLUE accepts that there is no unusual arrangement of parameters that optimizes the decency-of-fit criteria. In GLUE, all sources of parameter uncertainties are considered in the parameter uncertainty. The probability esteem, which is related to a parameter set, mirrors all sources of error and the impacts of covariation of parameter esteem on the model execution. A GLUE examination comprises the following three stages:

(1) After defining the "generalized likelihood measure,”  , many parameter sets are randomly sampled from the earlier dissemination; each parameter set is assessed as either "behavioral" or "non-behavioral” through a comparison of the "likelihood measure” with a selected threshold value.

(2) Each behavioral parameter set is given a "likelihood weight”, W_i, according to the following:

where n is the number of behavioral parameter sets.

(3) Finally, the prediction uncertainty is depicted as a forecast quintile from the aggregate  dissemination  acknowledged  from  the weighted behavioral parameter sets. The most frequently utilized likelihood measure for GLUE is the Nash-Sutcliffe coefficient of efficiency (NSE)

4.4 Parallel solution (ParaSol)

ParaSol depends on an adjustment to the worldwide optimization calculation, SCE-UA. The motivation is to utilize the simulation performed during optimization to derive prediction uncertainty. The simulations accumulated by SCE-UA are extremely important, as the calculation tests over the whole parameter space with an emphasis on arrangements close to the ideal/optima.

The process of ParaSol is as follows:

(1) After streamlining the application of the changed SCE-UA (the arbitrariness of the SCE-UA calculation is expanded to enhance  the scope of the parameter space), the reproductions performed are partitioned into  "great” recreations and "not great” recreations using a limit estimation of the target work, as in GLUE. This yields "great” parameter sets and "not great” parameter sets.

(2) A prediction of uncertainty is built by similarly weighing every "great” recreation. The target work utilized as a part of ParaSol is the sum of the squares of the residuals (SSQ):

Brief descriptions of the SWAT-CUP, SUFI-2, GLUE, and ParaSol model components are found in the SWAT-CUP user manual by Abbaspour [25] and in the SWAT model use, calibration, and validation document by Arnold et al. [26].

4.5. Data pre-processing

The input data required for SWAT modeling are DEM, a LULC map, a soil map, and meteorological data. First, the raw data procured from the different organizations described in Table 1 were pre-processed. The different DEMs downloaded from the United States Geological Survey (USGS) were mosaicked using the mosaicking tool in ArcGIS 10. In the mosaicking tool, multiple images are joined into one unit. To obtain a LULC map, the satellite imagery was processed. Image classification was performed using the ERDAS Imagine tool. Any satellite image comprises multiple bands in the electromagnetic spectrum. All the bands together form an image; however, this provides only spectral information. Converting the spectral information into LULC information classes or into any other particular information class requires image classification. Information classes are categorical classes, in which a pixel may represent a water area, forested area, or urban area. Based on a spectral curve or the spectrum of the pixel, classification is possible. There are two types of automated image classification techniques: supervised and unsupervised digital image classification. In supervised digital image classification, the software is guided by researchers or image interpreters with expertise in specifying the land cover classes of interest as a signature dataset, which is then automatically used by the software to create spectral classes. In unsupervised digital image classification, an interpreter only specifies the number of classes required to classify the image; the system then classifies the image without researchers’ expertise. The LULC map of this region (Fig. 1) was extracted by performing supervised classification on the optical satellite imagery. A signature dataset was created for six classes—water (WATR), forested areas (FRSD), urban areas (URBN), range (RNGE), agriculture (AGRL), and barren land (BARR)—using the visual interpretation keys; approximately 150 training datasets were used for a single class. The classification algorithm that was used was maximum likelihood. Most of the watershed comprises urban or barren land.

《Fig.1》

Fig. 1. LULC map of the study area.

To construct the soil map shown in Fig. 2, soil data were collected using National Bureau of Soil Survey and Land Utilisation Planning (NBSS & LUP) data and laboratory test results. The soil was divided into five layers. A user soil database describing each soil class was constructed for use in HRU analysis, as shown in Table 2. A climate database was created using Indian Meteorological Department (IMD) data; the latitude and longitude were given along with the rainfall, temperature, and solar radiation data in separate text files.

《Fig.2》

   Fig. 2. Soil map of the study area.

Table 2 Details of user soil table.

^a Hydrological soil groups are of four types: Type A, having an infiltration rate of 7.6–11.4 mm·h^-1 ; Type B, having an infiltration rate of 3.8–7.6 mm·h^-1 ; Type C, having an infiltration rate of 1.3–3.8 mm·h^-1 ; and Type D, having infiltration rate of 0–1.3 mm·h^-1 .

^b Texture L: loamy; CL: clay loamy; SCL: silty clay loamy; SL: silty loamy.

^c Soil ALB is the soil albedo: ratio of amount of solar radiation reflected by a body to the amount incident upon it.

4.6. SWAT modeling

The watershed is defined to include both the catchment and the drainage channel within a single morphometric divide. It is a nat-urally occurring hydrologic unit that is covered by natural boundaries and characterized by similar conditions such as physical characteristics, the topology of the land surface, and climatic conditions. Watershed delineation implies the drawing of lines on a map to indicate a watershed’s limits; these  are commonly drawn  on maps utilizing data from DEM or contour maps.

Watershed delineation is the first step of SWAT modeling. DEM is used as the input in this step. Along with the slope of the watershed, the streams and outlet points are generated. According to the outlet points selected by the researcher, the watershed is delineated. The user can also provide data related to the reservoir and predefined streams as input. In this project, the watershed was delineated by dividing the  watershed  into  46  sub-watersheds.  Fig. 3 shows the delineated watershed, streams, and monitoring points of the study area. The next step is  the HRU analysis. In   HRU analysis, the watershed is divided into units having similar but unique land types, soil types, and elevation properties. In this step, using the soil map, user soil table, LULC  map,  and  slope map (Fig. 4) of the study area, the  watershed  was  divided  into 160 HRUs. In the next step, the climate database was provided to the model. Finally, the SWAT model was run to yield the runoff, evapotranspiration, and sediment yields of each sub-watershed  and HRU for 20 years.

《Fig.3》

Fig. 3. Watershed delineation of the study area.

《Fig.4》

Fig. 4. Slope map of the study area.

4.7. Calibration and validation

As noted earlier, calibration is the modification of the elements influencing the SWAT yields, and is assisted by observed values  and evaluated estimations of runoff, evapotranspiration, and other SWAT outputs. Validation compares the SWAT results with the observed data without modifying the values of the influencing factors. Calibration is necessary for proper hydrological modeling. There are 34 parameters distinguished for runoff and soil erosion that can be utilized for calibration. Table 3 describes the parameters used in the present study for calibration, along with their minimum and maximum absolute SWAT values. In this step, SWAT-CUP was used as an interface between the SWAT and calibration algorithms to perform the uncertainty analysis, calibration, and validation of the hydrological modeling outputs. The TxtInOut folder of the SWAT model was imported into the SWAT-CUP software for input. The observed values of surface water discharge of the Varanasi sub-watershed were calculated using the acoustic Doppler current profiler (ACDP) for the years 1996–2015. These values were used for comparison and calibration. The input files of the SWAT-CUP model were updated as necessary. Fig. 5 shows the detailed procedure as a flowchart.

《Fig.5》

Fig. 5.  Calibration flowchart.

《Table 3》

Table 3 Parameters used for calibration, with their absolute ranges.

SCS: soil conservation service; USLE: universal soil loss equation.

4.8. Methodology and criteria for comparison

Different challenges appear when comparing calibration techniques used for soil and water modeling. The most important concerns that we addressed while comparing these techniques are as follows:

(1) Most algorithms are diverse in their theories. Subjective decisions must be made in their formulation concerning prior parameter distributions and objective functions. In this study, we addressed this issue by choosing objective functions for each algorithm as they would be utilized in hydrological applications. Fundamentally, this leads to various objective functions for various algorithms. In discussing the outcomes, we indicate whether an issue is caused by the theoretical definition of a specific algorithm or by the choice of the function.

(2) Every algorithm has its own fundamental concept and objective functions, which impairs comparison. To solve this problem, the value of each function behind every algorithm was calculated to permit fair comparison. Similarly, we utilized  the  measures of computational proficiency and an appraisal of the conceptual basis criteria for better comparison.

(3) It is fairly obvious that each algorithm yields different results. To address this problem, we considered the results of all  the algorithms for all possible criteria; we then outlined them all   so that the reader can draw his or her own particular conclusions.

(4) The outcomes of the comparison innately depend upon the applications. To solve this issue, the specific results of each application of the algorithm are separated from the generic ones.

The following five categories were used for the comparison of the three (SUFI-2, GLUE, and ParaSol) calibration algorithms: The calibration techniques use several parameters for calibration. For each technique, the best estimate and uncertainty range of these parameters differ, so the first comparison is based on the best estimate, the minimum and maximum uncertainty ranges of each algorithm, and the parameter correlation. The various algorithms use various objective functions, so the second comparison is based on the values of NSE, R 2, and other objective functions. The third comparison was based on the values of the R-factor, which is the average width of the band divided by the standard deviation of the corresponding measured variable, and on the P-factor, which is the percentage of data bracketed by the 95PPU band of each algorithm. The fourth comparison was based on theoretical concepts, testability, and the fulfillment of statistical assumptions. The final  comparison was  based on the difficulty of implementation.

5. Results

5.1. SWAT output

In the present study, the SWAT model was run for a watershed located in the northern part of India in the state of Uttar Pradesh. The SWAT divided the watershed into  46  sub-watersheds  and  760 HRUs for easy and accurate modeling. It was estimated that the average annual precipitation of the basin is 941.24 mm, snowfall is 0 mm, snowmelt is 0 m, surface runoff (SUR Q) is 358.56 mm, lateral discharge is 1.39 mm, and groundwater discharge for shallow and deep aquifers is 138.69 and 8.53 mm, respectively. The average value for total aquifer recharge is 180.65 mm, the total water yield of the basin is 507.17 mm, and evapotranspiration is 411.7 mm. A pictorial representation of the SWAT output regarding runoff and evapotranspiration is shown in Fig. 6. The results indicated that, on average, more than 45% of the total precipitated water is lost in runoff and evapotranspiration. The  average monthly values of all the parameters of the basin—that is, rainfall, snowfall, surface runoff, lateral runoff, water yield, and evapotranspiration—are given in Table 4. It was also estimated that there are an average of 52.11 annual water-stress days and 10.54 temperature-stress days.

《Fig.6》

Fig. 6.  Pictorial representation of the SWAT output.

《Table 4》

Table 4 Average monthly values of the parameters for the watershed.

Q: discharge or runoff; ET: evapotranspiration; PET: potential evapotranspiration; SED: sediment yield.

During watershed image classification, the watershed was classified into six classes according to the land use. It was found that most of the watershed is urban or barren land. Table 5 gives the annual average values of the parameters according to each land use type. It was also concluded from the results that for urban land, surface runoff may be excessive, and less than 22% of the water yield is base flow. For range and agricultural land, surface runoff may be excessive. For barren land, the sediment yield may be  overly high, more than half of the precipitation is lost to runoff,   the surface runoff is the highest, and less than 22% of the water yield is base flow.

《Fig.7》

Fig. 7.  SWAT output for sediment  yield.

5.2. Comparison output

Table 6 summarizes the comparison of the calibration algorithms. The comparison was performed over the five categories described earlier, with the following results:

Category 1: According to this category, GLUE is better than SUFI-2 and ParaSol, as GLUE has the widest range of uncertainties. Thus, most of the intervals considered by SUFI-2 and ParaSol are covered by GLUE. In this category, SUFI-2 is the second best.

Category 2: According to this category, ParaSol is the best for its objective function, that is, NSE, as it is based on a global optimization algorithm. As shown by the table values, it can be concluded that SUFI-2 and GLUE have similar values, so they hold the same position.

Category 3: According to this category, ParaSol is worse than  the other two algorithms because of its narrow prediction uncertainty bands. From the other values, it can be concluded that  SUFI-2 works better than GLUE.

Category 4: According to this category, all sources of uncertainties are considered by SUFI-2 and GLUE, whereas ParaSol only considers parameter uncertainty and ignores all other uncertainties. On a conceptual basis, SUFI-2 is better than GLUE and ParaSol because it gives the most efficient results.

Category 5: According to this category, GLUE is the best, as it is the easiest of the three algorithms.

From the comparison and description given in Section 5.2, we can conclude that the performances of SUFI-2 and GLUE are very similar. However, by looking at the 95PPU plot (Figs. 8 and 9),     we concluded that SUFI-2 performs better than GLUE, as the results are more efficient and the algorithm considers uncertainties better. SUFI-2 can be run with the smallest number of parameters. GLUE works with a large number of simulations, but cannot provide better results than SUFI-2. The main drawback of GLUE is the length of the computational techniques because of its random sampling strategy. Thus, SUFI-2 is considered to be the best for this study, and further analysis of sediment yield calibration and validation was done using only the SUFI-2 algorithm.

Table 6 Comparison of calibration algorithms, according to the criteria given in Section 4.8.

^a Best estimate values of parameters and their minimum and maximum range.

^b Values of objective functions.

^c Uncertainty described by parameter uncertainty.

《Fig.8》

Fig. 8. Plot of estimated and observed values after calibration using GLUE.

《Fig.9》

Fig. 9.  Plot of estimated and observed values after calibration using SUFI-2.

5.3. SUFI-2 output

In the present analysis, the model was calibrated using the sensitive parameters given in Table 3. As indicated by the execution evaluation criteria of the model, which is suggested for a monthto-month time cycle, five parameters were determined to be the most sensitive parameters for modeling the present watershed. These parameters are the effective channel hydraulic conductivity (CH_K2), universal soil-loss equation (USLE), support parameter (USLE_P), Manning’s n value for the main channel (CH_N2), surface runoff lag time (SURLAG), and available water capacity of the soil layer (SOL_AWC). The behavioral threshold for calibration was set to 0.500000. The period from 2000 to 2004 was considered to be a warm-up period. Data from the period 2004–2009 were used for calibration, and the data period 2009–2015 was used for validation. The results of the calibration and validation performed by the SUFI-2 algorithm are shown in Figs. 9 and 10. Fig. 9 shows a plot of the estimated and observed values of discharge and Fig. 10 shows a plot of the estimated and observed values of sediment data after calibration. Estimated outputs were compared at the same outlet point in sub-basin #27.

《Fig.10》

Fig. 10. Plot of estimated and observed values of sediment yield (mm).

Table 7 lists the evaluation coefficients of the simulated  monthly sediment yield of different objective functions for the study watershed. The assessment coefficients of the simulated month-to-month sediment yield of various  target  functions  for  the study area are recorded in Table 6. The P-factor, or the percentage of observations bracketed by the 95PPU, was 0.69. The R-factor was 0.763. The outcomes demonstrated that SUFI-2 sectioned a lower  percentage  of  the   observed   sediment   yield.   Similarly,  R² = 0.78, NSE = 0.76, percent bias (PBIAS) = 2.4×10², and the observation standard deviation ratio (RSR) = 0.49 are not strong indicators of the goodness of fit; however, they are within adequate evaluation ratings. The outcomes demonstrate that  the SWAT can simulate the hydrological attributes of the Indian watershed extremely well. Hence, this model can be used for further hydrological studies in the basin.