Abstract

The soil and water assessment tool (SWAT) is a physically based model that is used extensively to simulate hydrologic processes in a wide range of climates around the world. SWAT uses spatial hydrometeorological data to simulate runoff through the computation of a retention curve number. The objective of the present study was to compare the performance of two approaches used for the calculation of curve numbers in SWAT, that is, the Revised Soil Moisture Index (SMI), which is based on previous meteorological conditions, and the Soil Moisture Condition II (SMCII), which is based on soil features for the prediction of flow. The results showed that the sensitive parameters for the SMI method are land-use and land-cover features. However, for the SMCII method, the soil and the channel are the sensitive parameters. The performances of the SMI and SMCII methods were analyzed using various indices. We concluded that the fair performance of the SMI method in an arid region may be due to the inherent characteristics of the method since it relies mostly on previous meteorological conditions and does not account for the soil features of the catchment.

1. Introduction

The soil and water assessment tool (SWAT) is a semiphysical and semiempirical modeling tool developed by the United States Department of Agriculture (USDA) [1]. SWAT can be developed for rainfall-runoff simulation with different short-term and long-term objectives, for example, prediction of stream flow and sedimentation, identification of best management practices (BMP), and planning water resources of large and complex watersheds [24]. In recent years, the application of SWAT has increased rapidly as a modeling tool for various hydrological processes under different climate conditions around the world. It is anticipated that SWAT will continue to evolve in response to the need for improved understanding of key hydrological processes [5].

SWAT uses the soil conservation services’ curve number (SCS-CN) to estimate the volume of stream flow. The SCS-CN method uses soil classes and land-use/land-cover information for the prediction of runoff [6]. The method separates effective rainfall from total rainfall by using the water balance concept [7]. The SCS-CN method is used extensively for consecutive flow simulation for long-term runoff analysis [8]. It can also be used for optimal estimation of surface runoff in ungauged watersheds [9]. Due to its feasibility and uncomplicated characteristics, the method has been integrated into several hydrological models [7]. SWAT is one such hydrological modeling tool that uses the SCS-CN method to predict water flows in large and ungauged watersheds.

SWAT uses two methods for the calculation of the retention parameter and estimation of the runoff curve number. One method is related to soil profile water content, and it is known as the Soil Moisture Condition II (SMCII) method [10]. The other method is related to accumulated plant evapotranspiration. Usually, SMCII method predicts much more runoff in shallow and/or low storage soil areas [11]. In other words, this method simulates stream flows in rivers or channels by overestimation, especially for soil types and land use that involve high-permeability mixtures of sand and silt. The Revised Soil Moisture Index (SMI) method was developed for SWAT to account for the daily curve number as a function of plant evapotranspiration. This method is less dependent on soil storage and more dependent on the previous climate [12].

In recent years, extensive research has been conducted to examine the capability of CN methods and/or to develop modified CN methods according to local conditions. A group of studies focused on exploring the behavior of the SMCII and SMI methods in the SWAT model for flow prediction and the contribution of the regional climate on stream flow [1217]. The results showed that the use of SMI method in SWAT model yielded promising performance in calculating discharges from plant evapotranspiration and the previous climate. Kim et al. [7], White et al. [14], and Easton et al. [15] proposed new methods for high slopes and monsoon climates to explore the optimal estimation of flood events. Their results indicated that modification of the SMCII method was required for reliable simulation of stream flow in monsoon climates in order to avoid overprediction of extreme events, such as floods. Therefore, Kim and Lee [16] used weighted average curve number to improve the poor performance of SWAT in predicting extreme rainfall events. The performance of SMI method in SWAT has also been assessed for simulating stream flow in a low-gradient, forested climate, and it has been reported that the SMI method can overcome the complexities of calculating stream flow [17]. However, the performance of the SMI method is still unknown in a thin forest climate, so it is not known whether it is viable for practical use in hydrological sciences in such a climate [17]. This emphasizes the need for additional research on exploring the behavior of the SMCII and SMI methods in the SWAT model in various climates, such as arid and semiarid zones [12, 13].

During the past few decades, SWAT has been used in different catchments of Iran for various purposes. Most of the studies addressed the impact of land-use changes, water resource management, pollution tracking, impact of effective precipitation on stream flow, and the impact of global climate change on river discharge [1825]. Generally, all the studies used the SMCII method (the conventional method) for calibration of the SWAT model, and the impact of the SMI method has not been reported for discharge simulation. Moreover, most of the studies mentioned above used manual calibration.

This research is a supplementary concept for developing a SWAT model in Roodan Plain by Jajarmizdeh et al. [26, 27]. The overall objectives of the study are to evaluate the performance of SMI and SMCII methods in SWAT for stream flow in a large arid catchment in order to have a better understanding of the sensitive parameters and the behavior of SWAT in such climate. Furthermore, the present study includes the application of a semiautomatic calibration algorithm and uncertainty analysis for the SWAT model to compare the performances of the SMCII and SMI methods in approximating discharges.

2. Methodology

SWAT consists of several different modules, including simulation of weather, crop growth, evapotranspiration, runoff, soil erosion, and transport of nutrients and pesticides. SWAT divides a watershed into subcatchments, and the subcatchments are divided into hydrological response units, which are individual combinations of soil, land cover, and slope. In SWAT, surface runoff and pollutants from hydrological response units (HRUs) within a subcatchment are routed to the catchment reach by using a lumped procedure for simulating the runoff of the catchment. Two methods are available in SWAT for estimating runoff, that is, the modified NRCS-CN method (or Direct Soil Moisture (SMCII)) method and the Revised Soil Moisture Index (SMI) method.

The SMCII method is based on the runoff curve number (CN) method [6], which uses watershed land-cover and soil characteristics to estimate direct runoff from a rainfall event [13]. CN values in this model are calculated on a daily basis using soil retention parameter values as a proxy for soil water status characterized by field capacity, wilting point, and saturation water content. When water content varies with the soil profile, the following equation is used to calculate the retention parameter: where is the retention parameter for a given day (mm), is the maximum value the retention parameter can achieve on any given day (mm), SW is the soil water content of the entire profile, excluding the amount of water held in the profile at wilting point (mm), and and are shape coefficients. The shape coefficients can be obtained by solving (1) with the following assumptions: (a) the retention parameter for moisture condition I curve number corresponds to wilting point of soil profile water content; (b) the retention parameter curve number for moisture condition III corresponds to field capacity soil profile water content; and (c) the soil has a curve number of 99 when completely saturated. This means that the retention parameter is 2.54 as the soil is fully saturated [1]. Based on the above assumptions, (2) can be obtained by solving (1) to derive shape coefficients that reflect the nonlinearity of water content in the soil layers as follows: where is the first shape coefficient, is the second shape coefficient, is the amount of water in the soil profile at field capacity (mm), is the retention parameter for the moisture condition III curve number, is the retention parameter for the moisture condition I curve number, is the amount of water in the soil profile when completely saturated (mm), and 2.54 is the retention parameter value for curve number 99.

The Revised Soil Moisture Index (SMI) method is related to the continuous CN approach. The curve number is modified by replacing the five-day antecedent rainfall with a Soil Moisture Index as part of a daily water-yield model. The SMI is driven by rainfall, runoff, and lake evaporation [13]. By calculating the daily CN as a function of plant evapotranspiration, the value of CN becomes less dependent on soil storage and more dependent on antecedent climate. In the SMI, the retention parameter is updated at the end of every day according to the following equation: where is the retention at the present time step, is the retention at the previous time step, is the potential evapotranspiration for the day , is the depletion coefficient, is the rainfall depth on the day, is the runoff depth, and is the maximum value of the retention. The concept of SMI method is developed to overcome the disadvantages of SMCII method which overestimate runoff in low storage/or shallow soils. The SMI method allows the retention parameter to vary with accumulated plant evapotranspiration and was found to overcome the disadvantage of SMCII method. Full description for runoff estimation methods by using curve number approaches in SWAT can be found in [1, 7, 1113].

In the present study, inverse modeling (IM) was used for the calibration and validation of SWAT model. Sequential uncertainty fitting algorithm (SUFI-2) was used for uncertainty analysis by considering all sources of uncertainties, namely, uncertainty in driving variables, conceptual model, and parameters and measured data. Two indices were used to measure uncertainty, namely, -factor and -factor. The -factor is the mean thickness of the 95% prediction uncertainty (95PPU) band and the -factor is the percentage of measured data within 95% prediction uncertainty. Theoretically, the value for -factor lies between 0 and infinity, while -factor ranges between 0% and 100%. If the -factor is equal to 100% and the -value is 0, it is considered that the simulation exactly matches the observed data [2830]. In addition to -factor and -factor, accuracy of SWAT model was evaluated by using Nash-Sutcliffe (NS) coefficient and Root Mean Square Error as given in (4) and (5), respectively, as follows: where is the number of time steps, and are the simulated and the observed stream flows, respectively, at time step , and is the average observed stream flow over the simulation period. Usually, if the value of the Nash-Sutcliffe coefficient is more than 50%, the SWAT model can be accepted [3].

The study area known as the Roodan watershed is located in the south of Iran between Hormozgan and Kerman Provinces. The area of catchment is 10570 km2, and it lies between the northern geographical latitude of 26°57′ to 28°31′ and the eastern longitude of 56°47′ to 57°54′. The location of the study area in Iran is shown in Figure 1. The average annual precipitation of the study area is 215 mm. Most of the precipitation occurs between October and March, and no substantial amount of precipitation occurs during the warm months. The predominant soil type is a heterogeneous mix of clay, silt, and sand in the northern and center parts of the watershed, whereas the soils at the southern and eastern part are mostly silt and clay. According to the reports of an agricultural organization of the Hormozgan Province of Iran, approximately 30% area of the Roodan watershed is covered with low-storage soil, which is commonly a mixture of sand and silt [31]. The mechanical compositions of different types of soil and their coverage in Roodan watershed are given in Table 1. Some important and dominating land covers of the Roodan watershed are shrub land (range brush) and a mixture of grassland with shrub land.

SWAT requires the following data for runoff modeling: (1) digital elevation map (DEM); (2) land-use map; (3) soil map; and (4) meteorological data. In the present study, DEM was prepared from 1 : 25000 topographic maps with a resolution of 90 m. A technique known as digital burning of river network [32] was applied on the DEM for the delineation of subbasins. Smallest subbasins cover around 3 to 4% of the total watershed area. The soil map was obtained from the global soil map archives of the Food and Agriculture Organization (FAO) of the United Nations. Land-use information of Roodan was prepared from a Landsat-7 image in conjunction with an available land-use map (1 : 25000) and statistical records provided by the agricultural organization of Hormozgan Province. Daily precipitation and temperature data for the period 1988–2008 were collected from a local meteorological station. Potential evapotranspiration was calculated by using the Hargreaves method.

In the present study, 26 parameters of SWAT having important roles in the estimation of daily stream flow were used for sensitivity analysis [33]. First, the sensitivity analysis was performed by using the Latin hypercube-one-factor-at-time (LH-OAT) method before the calibration of the model. Details of LH-OAT method can be found in [34]. Next, the SUFI-2 algorithm was used for uncertainty analysis. Data from 1996 to 2002 were used for calibration, and data from 2003 to 2007 were used for validation. SWAT is a continuous stream flow simulator and, therefore, it requires continuous data. Uncertainty in model simulation is supposed to increase if one or more years are excluded during model calibration or validation. Therefore, continuous daily time series data over the time period 1996–2007 was used for the study.

3. Result and Discussion

Table 2 shows the sensitive parameters for the SMI and SMCII methods in Roodan watershed obtained by using SUFI-2 algorithm. The research outcomes showed that the effective hydraulic conductivity of the main channel was the most sensitive parameter for both SMI and SMCII methods. This might be due to arid climate of the catchment. The aridity in climate results in many transmission losses, which are related to the effective hydraulic conductivity.

It can be seen from Table 2 that the parameters related to soil, such as soil evaporation compensation factor and soil conductivity, do not have any notable significance in the SMI method because the method is based on plant evapotranspiration and antecedent climate. In contrast, the plant evaporation compensation factor does not have any substantial impact on the SMCII method. This is due to the fact that SMCII method calculates runoff based on land-use/land-cover and soil-texture characteristic. Previously, it was reported by Kannan et al. [35] that the most sensible parameters are minimum depth of water in soil for base flow to occur, the groundwater reevaporation coefficient, and the threshold depth of water in the shallow aquifer for percolation to the deep aquifer. Amatya and Jha [17] and Williams et al. [13] found the most sensitive parameters to be soil evaporation compensation factor, plant evaporation compensation factor, curve number, manning coefficient for channel, surface runoff lag coefficient, and maximum canopy index. These findings collaborate to the results obtained in the present study.

Table 3 shows the values of the RMSE and NS coefficients obtained during calibration and validation of the model using SMCII and SMI methods. Comparison between the SMCII and SMI methods in simulating daily stream flow showed that the SMCII performed better than SMI. The RMSE was 50.2 during calibration and 36.2 during validation for SMI, and it was 36.9 during calibration and 27.6 during validation for SMCII. According to the obtained values of Nash-Sutcliff coefficient (NS) [36], the performance of SMI methods in simulating daily stream flow can be classified as fair during calibration and good during validation. On the other hand, the performance of SMCII method can be classified as good during both calibration and validation. Amatya and Jha [17] developed a daily flow simulation model by using SWAT with the SMI method for estimating the retention parameter in a large forested watershed in Turkey and reported NS coefficients of 0.59 and 0.7 during model calibration and validation, respectively. In comparison to their results, the NS coefficients obtained for Roodan watershed can be considered promising for an arid region. Since the -factor should be near to 100% and the -factor should be close to zero [30], it can be concluded that the performances of both the SMCII and SMI methods are acceptable in accordance with model uncertainty in prediction. It appears that SMI had slightly better values for the uncertainty factors ( and ). The reason might be better prediction of low flows (near zero) by this method. The Roodan watershed includes mostly intermittent or ephemeral flows. Therefore, better adjustment for low flows resulted in better simulation and consequently less uncertainty for SMI.

Figures 2 and 3 show observed and predicted stream flow obtained by using the SMI and SMCII methods in the Roodan watershed during model calibration and validation, respectively. The results show that SMI method overestimated and SMCII method underestimated peak stream flow during calibration. The largest event during calibration was on March 1, 1996 (observed = 1462 m3/s), which was estimated by SMI as 2018 m3/s and SMCII as 1193 m3/s. On the other hand, both methods underestimated the largest peak flow during the validation period (in February 2005, observed = 1248.5 m3/s, SMI = 664.7 m3/s, and SMCII = 746.5 m3/s). It can be speculated that the SMI method is sensitive to the amount of precipitation and the area of low-storage or shallow soils. It can be noted that about 30% of the area of the Roodan watershed has low-storage soil, and the annual rainfall is less than 250 mm. However, SMCII still was found to provide better results. This may be due to the lower percentage of low-storage soil in the Roodan watershed, or it may be possible that the low-storage soil only has an impact on runoff for a brief period of time.

Figures 4 and 5 show the cumulative daily stream flow obtained by using the SMI and SMCII methods during model calibration and validation, respectively. Generally, the cumulative flow graph helps in estimating runoff for watershed planning or management. Figure 4 shows that the SMI method estimated closer cumulative discharge compared to SMCII method at the end of calibration period (1996–2002). The SMCII method showed that the trend of the fluctuation of the cumulative discharge was in better agreement with the observed data (Figure 4). On the other hand, SMCII showed overestimation for cumulative discharge.

Figure 5 shows the underestimation of cumulative flow by SMI method and overestimation by SMCII method during model validation. The SMCII was found to provide excellent prediction of cumulative stream flow for the period of 2003–2006. During the rest of the validation period (2006-2007), SMCII was found to yield slight overestimations. On the other hand, SMI method underestimated the cumulative discharge seriously during model validation, especially from 2005 to 2007.

Kannan et al. [12] assessed the performance of SMCII and SMI methods in the catchments of the United States where the average annual rainfall was 292 mm and found that the SMI method underestimated the stream flow. They concluded that the SMI method is not as capable as the SMCII method in simulating stream flow in low precipitation catchments; this is in agreement with the finding of present study. Figure 5 also shows that the SMI method predicted stabilization of cumulative stream flow at the end of validation period. This may be due to more sensitivity of SMI model to climate and the amount of area that is covered by low-storage soil in a catchment.

Table 4 illustrates the percentile analysis of stream flow during model validation. It can be seen from the table that the SMI method yields closer values for percentile of stream flow during model validation. For 25 or less percentile data, both SMI and SMCII were found to predict similar stream flow values. Amatya and Jha [17] also reported that SWAT could simulate low stream flows (near zero) in a forested catchment. This means that the forested area in the catchment improves the predictions for the low flows by SWAT. For 50 percentile (value of median), SMI had fair agreement with the observed data in comparison with the SMCII method.

Finally, the behavior of the SMI and SMCII methods is assessed by analyzing the relative errors (%) for stream flow exceeding 500 m3/s during model validation (Table 5). The table shows that SMCII yielded lower relative errors for high stream flows during model validation. On the other hand, both methods were found to underestimate the highest discharge during model validation. Both methods also displayed an underestimation trend for flow that exceeded 500 m3/s. Table 6 presents the general statistics of stream flows during model validation. The SMI method underestimated the mean flow (observed = 5.4 m3/s, SMI = 2.5 m3/s, and SMCII = 5.7 m3/s). This may be due to the inherent characteristics of the SMI method that offset high predictions of stream flow. In general, the SMCII method was more promising than those of the SMI method. However, SMI showed better estimation for low flows.

4. Conclusions

Two approaches of curve number generation embedded in the SWAT model were discussed for hydrologic prediction. The first was the Soil Moisture index (SMI), which is based on previous meteorological conditions, and the second was the Soil Moisture Condition II (SMCII), which is based on soil features. The study showed that the SMI method relies mostly on previous climate conditions. It does not account for soil features in the catchment in accordance with the sensitivity of related parameters. However, the SMCII method is highly sensitive to soil features. The present study showed that both SMI and SMCII methods can predict daily steam flow with reasonable accuracy; however, SMCII showed better total performance for flow prediction. However, low flows were modeled better via the SMI method in SWAT. The results suggested that the results of the SMI method can be affected by the amount of area in a catchment that has low-storage soils. Regional climate also might be an important factor that determines the performance of the SMI method. Therefore, moderately accurate simulation of stream flows was by the SMI method for flows greater than 50 percentile (value of median). Statistical evaluation showed that the SMCII yielded better values in simulating average flow during calibration and validation. It can be suggested that the performance of SMI and SMCII methods should be analyzed in other arid catchments with low-storage soil coverage for better assessment of these methods in an arid region. Also, there is a need to conduct additional evaluations to explore the behavior of the SMI and SMCII methods in various regional climates.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors sincerely appreciate the Research Management Center (RMC) of Universiti Teknologi Malaysia for funding them to carry out this research under postdoctoral fellowship scheme as well as through university research Grant (vot no. Q.J130000.2522.06H36). The authors are thankful to all members of Consultant Engineers of Ab Rah Saz Shargh Corporation in Iran and the Regional Water Organization, Agricultural Organization, and Natural Resources Organization of Hormozgan Province, Iran.