An Improved SCS-CN Method Incorporating Slope, Soil Moisture, and Storm Duration Factors for Runoff Prediction

: Soil Conservation Service Curve Number (SCS-CN) is a popular surface runoff prediction method because it is simple in principle, convenient in application, and easy to accept. However, the method still has several limitations, such as lack of a land slope factor, discounting the storm duration, and the absence of guidance on antecedent moisture conditions. In this study, an equation was developed to improve the SCS-CN method by combining the CN value with the tabulated CN 2 value and three introduced factors (slope gradient, soil moisture, and storm duration). The proposed method was tested for calibration and validation with a dataset from three runoff plots in a watershed of the Loess Plateau. The results showed the model efficiencies of the proposed method were improved to 80.58% and 80.44% during the calibration and validation period, respectively, which was better than the standard SCS-CN and the other two modified SCS-CN methods where only a single factor of soil moisture or slope gradient was considered, respectively. Using the parameters calibrated and validated by dataset of the initial three runoff plots, the proposed method was then applied to runoff estimation of the remaining three runoff plots in another watershed. The proposed method reduced the root-mean-square error between the observed and estimated runoff values from 5.53 to 2.01 mm. Furthermore, the parameters of soil moisture ( b 1 and b 2 ) is the most sensitive, followed by parameters in storm duration ( c ) and slope equations ( a 1 and a 2 ), and the least sensitive parameter is the initial abstraction ratio λ on the basis of the proposed method sensitivity analysis. Conclusions can be drawn from the above results that the proposed method incorporating the three factors in the SCS method may estimate runoff more accurately in the Loess Plateau of China.


Introduction
The ability to make surface runoff estimations has become an essential part of development strategies in water resources management, flood control, and water and soil conservation [1]. A multitude of hydrologic models have been developed to predict direct runoff. Among these methods, the Soil Conservation Service Curve Number (SCS-CN) method [2] is one of the most widely used. The method was originally developed for surface runoff prediction, but has been applied to several other areas such as the infiltration, sediment yield, and transport of pollutants [3,4]. Moreover, it has been extensively integrated into many hydrological and ecological models [5], including AnnAGNPS (Annualized Agricultural Nonpoint Source Pollution Model) [6], CREAMS (Chemicals, Runoff, and Erosion from Agricultural Management Systems) [7], EPIC (Environmental Policy Integrated Climate) [8], SWAT (Soil and Water Assessment Tool) [9], and EBA4SUB (Event-Based Approach for Small and Ungauged Basins) [10,11].
The popularity of the SCS-CN method is due to its simplicity, convenience, widespread acceptance, applicability to ungauged watersheds, and requirement of only one parameter (CN), which is determined by four readily grasped watershed characteristics (soil group, land cover, surface condition, and antecedent moisture condition (AMC)) [12,13]. However, the SCS-CN also has some limitations such as lack of a land slope factor, discounting the storm duration, and the absence of guidance on antecedent moisture conditions [14][15][16][17][18][19]. These advantages and disadvantages made the SCS-CN method a perennial topic of discussion over the last four decades [20][21][22][23][24][25].
Storm duration is a vital part of the rainfall-runoff process, which could greatly influence runoff estimations [14]. However, the storm duration is not considered in the SCS-CN method, which leads to uncertainty in runoff prediction because of the spatiotemporal variability of rainfall [12]. Several enhanced SCS-CN methods have been proposed to overcome the problem of discounting the storm duration. Jain et al. [17] suggested a new model formulation with storm duration-based rainfall adjustment; Mishra et al. [18] presented a rainfall duration-dependent procedure by developing an equation for the new CN value with an introduced minimum CN and rainfall duration; Petroselli et al. [26] presented a model for continuous simulation and the SCS-CN method was modified to calculate the net rainfall by resolving the problem of insufficient including infiltration; Sahu et al. [27] incorporated the storm duration factor into the SCS-CN method by dividing the cumulative infiltration amount F into dynamic infiltration and static infiltration (Fc), which is a product of the minimum infiltration rate and rainfall duration; Shi et al. [25] introduced the static infiltration into the soil moisture accounting (SMA)-based SCS-CN method to improve runoff predictions on the Loess Plateau. However, all these methods have no contact with the CN value of the original SCS-CN method, which has limited the application of the models.
Antecedent moisture condition (AMC) is a significant factor determining the initial abstraction of the SCS-CN method. The AMC prior to a storm event is divided into three levels: dry (AMC 1), average (AMC 2), and wet (AMC 3) based on the five-day antecedent rainfall amount in the traditional SCS-CN method, which will cause the CN value to suddenly jump from one level to another [28]. Soil moisture seems to be a better choice for selecting an appropriate AMC value as compared with the five-day antecedent rainfall depth, and improved knowledge of the relationship between antecedent moisture and the CN value would greatly improve the runoff estimation of the SCS-CN method [29,30]. Saxton [31] redefined the antecedent moisture condition with a step function of soil water content. Koelliker [32] suggested that CN changes linearly with the soil moisture. The EPIC and SWAT models also adopted a continuous equation between soil moisture and CN for runoff prediction. Jacobs et al. [33] used remote sensing of soil moisture for CN adjustment to predict runoff in five basins of Oklahoma. Based on the climatic conditions of the Loess Plateau in China, Huang et al. [34] established an equation between soil water content and CN value for better runoff estimation.
The land slope is another main factor affecting flow movement within the landscape. The CN values in the USDA-NRCS (US Department of Agriculture-Natural Resources and Conservation Service) handbook table are based on 5% slope, and should be adjusted according to the actual slope. Several efforts have also been made on consideration of a slope factor in the CN method. Sharpley and Williams [8] introduced a slope-adjusted CN2 in the calculation of the runoff volume, but it has not been tested in the field. Based on the data of experimental plots with slopes varying from 14% to 140% on the Loess Plateau of China, Huang et al. [15] proposed an equation for considering the influence of slope on CN value. However, there is no study that couples these factors (storm duration, soil moisture, and slope) in the SCS-CN method for predicting runoff, which would be the focus of this study.
Therefore, the objectives of this study are (1) to develop an equation of the CN value with original CN2 value, slope gradient, soil moisture, and storm duration in the conventional SCS-CN method; (2) to compare the performance of the standard SCS-CN method, those from Huang et al. [15] and Huang et al. [34], and the proposed method by observing three experimental plots in the Loess Plateau region; (3) to apply the proposed method to predict runoff from three runoff plots in the other watershed.

The Original SCS-CN Method
The SCS-CN method [2] was originally developed based on two fundamental hypotheses and one simple water balance equation, which can be written as follows: where P and Q are the depth of observed rainfall and direct runoff, respectively (mm); Ia and  are the initial abstraction (mm) and coefficient of initial abstraction (dimensionless), respectively; F is the cumulative amount of infiltration (mm); and S is the maximum potential retention (mm), which can be calculated by where CN is in the range of 0-100 (dimensionless). This CN value is determined by the CN value (CN2) of the average moisture condition (AMC 2), which depends on land cover, soil group, and hydrologic conditions using a table from the SCS handbook [2], and is then converted to AMC 1 or AMC 3 based on the five-day-prior precipitation.

The Proposed Method
In the proposed method, slope gradient, soil moisture, and storm duration factors are considered. The CN value for each rainfall-runoff event can be improved by multiplying the reference value of SCS handbook CN2, and the function of slope, soil moisture, and storm duration: where α is the slope gradient (m m −1 ), θ is soil moisture (cm 3 cm −3 ), and t is storm duration (h). f(α), f(θ), and f(t) are functions of α, θ, and t, respectively.
The slope equation proposed by Huang et al. [15] was used to consider the influence of slope on the CN2 value. They used the modified SCS-CN method to evaluate runoff prediction in runoff plots with an 11-year observation experiment and a slope range of 14%-140% in Xifeng City of the Loess Plateau. It is adopted in this study: A nonlinear equation between CN and soil moisture was developed by Huang et al. [34], who found that the V0 value prior to each storm event can be expressed by the water storage in the upper 15 cm soil. It is also adopted and can be expressed as follows: For the storm duration function, we proposed a one-parameter (c) equation for f(t) after several trials: where a1, a2, b1, b2, and c are the empirical coefficients. Therefore, Equations (5)-(8) combined with Equations (1)-(4) make up the proposed method for runoff prediction. The flowchart of the whole methodology applied in this study is shown in Figure 1.

Performance of the Methods
The following two statistical indices were used to evaluate the performance of the methods:  (10) where NSE is the Nash-Sutcliffe efficiency [35,36] and RMSE is the root-mean-square error; i Q and i Q  are the i th observed and estimated runoff, respectively; Q is the average measured runoff of events; and N is the total events. Higher NSE and lower RMSE values indicate that the model has better agreement with the observations.

Study Area
This study was conducted in six test plots, three each in the Xindiangou (XDG) and Chabagou (CBG) watersheds, which are both in the hilly region of the Loess Plateau ( Figure 2).

Method 1
The proposed method with λ = 0.2 and slope parameters fixed

Method 2
The proposed method with slope parameters fixed

Method 3
The proposed method with all parameters optimized The XDG watershed (latitude: 37°29′ N, longitude: 110°08′ E, elevation: 990-1010 MASL, area: 30 km 2 ) is located in Suide County. The average annual temperature is 8 °C and the mean annual precipitation is 470 mm, mostly between June and September. The soil type of the XDG watershed is mainly silty loam. The characteristics of soil physical and particle size distribution are homogeneous in the top 30 cm soil [37].
The CBG watershed (latitude: 37°31′ N, longitude: 109°47′ E; elevation: 900-1100 MASL; area: 205 km 2 ) is located in Zizhou County. The mean temperature is also 8 °C and the average annual precipitation is 450 mm. The precipitation from July to August accounts for more than 70%, most of which are heavy and short-term rainstorms. The main soil in CBG watershed is Malan loess soil, in which the content of clay particles is less than 40%, leading to the large porosity of soil and the vulnerability to erosion [38].

Data Collection
Rainfall characteristics in the six plots, including depth and duration of each rainfall, were monitored by self-recording rain gauges. For all the experimental sites in the two watersheds, the rain gauges were located less than 500 m from each experimental plot. At the downstream end of each plot, a funnel-type collector was adopted to collect the surface runoff ( Figure 3). Soil samples in each plot were collected at an interval of 10 cm for 0-50 cm soil layer at 10:00 a.m. during the monitoring period, and the gravimetric water content of the collected soil samples was then determined by the oven-drying method. The water content of the top 15 cm soil could be calculated from the mean soil water measurements at 10 and 20 cm soil depth. These datasets have been compiled based on rainfall-runoff events and printed for internal use by the Yellow River Administration Committee of the Ministry of Water Resources and the Shaanxi Institute of Soil and Water Conservation.
The standard SCS-CN and proposed methods were tested using data derived from the three plots in the XDG watershed for calibration and validation, and the tested methods were then applied to the remaining three runoff plots in the CBG watershed.

Parameter Estimation
The rainfall-runoff events of the three plots in the XDG watersheds were divided into two groups: the datasets of cropland and one grassland plot were used for calibration, and the remaining datasets of another grassland plot were used for validation. The Marquardt [39] algorithm for solving constrained least squares problems was adopted to optimize all the parameters of the test methods. In order to examine the applicability of the models, λ = 0.2, assigned in its original development, is adopted for the traditional SCS-CN, Huang et al. [15], Huang et al. [34], and the proposed method (Method 1). The slope parameters (a1 and a2) obtained from Huang et al. [15] were also used for Huang et al. [15] and the proposed method (Method 1). Moreover, the effect of initial abstraction and slope on the proposed method was tested with the optimized λ (Method 2) and both optimized λ and slope parameters (Method 3).

Model Calibration and Validation
The number of storm events observed for calibration and validation in watershed XDG was 47 and 29, respectively. Table 2 presents the parameters of the tested methods optimized with the calibration dataset. The overall performance of all tested models based on statistical indicators is compared in Table 3.   Figure 4 presents the runoff estimated by the original method plotted against the corresponding observed values for the full datasets in XDG watershed. The standard SCS-CN method underpredicted the large and overpredicted as well as underpredicted some small storm-runoff events (Figure 4a). This is also confirmed by the statistical indexes in Table 3. There is a similar pattern in the calibration and validation sub-datasets, where the slopes of regression are 0.83 and 0.77, and the NSE values is −118.64% and −182.60%, respectively (Table 3), which illustrates the poor performance of the conventional SCS-CN method. The Huang et al. [15] and Huang et al. [34] methods incorporated the slope gradient and soil moisture condition factors into the standard SCS-CN method. The prediction of large storm-runoff events was improved with a higher slope of 0.90 using the Huang et al. [15] method as compared with the conventional SCS-CN method (0.80) (Figure 4b). However, it also overpredicted some small storm-runoff events when the slope factor is considered in the original SCS-CN method. Therefore, it can be seen that the Huang et al. [15] method did not improve the prediction accuracy of the full dataset much, as verified by the NSE value, which decreased from −148.17% (original SCS-CN method) to −177.63% (Table 3) Table 3 also shows that the NSE value (−9.8%) of the Huang et al. [34] method during the calibration period was significantly improved as compared to the SCS-CN (−118.64%) and Huang et al. [15] method (−133.00%). The overprediction of the latter two methods incurred for some small storm-runoff events were mitigated. However, the regression slope of the Huang et al. [34] method deviates more from the 1:1 line more than the traditional SCS-CN and Huang et al. [15] method, which suggests that the former still underestimated large and overestimated as well as underestimated small runoff events (Figure 4c).

The Proposed Method (Methods 1-3)
The runoff estimated by Method 1 is closer to the observed values, with the coefficient of determination R 2 increasing to 0.79 from 0.21 (the SCS-CN method), and the regression line of which is closer to the perfect line, where the regression line slope and intercept is 0.90 and 0.30, respectively (Table 3). Moreover, the prediction results of Method 1 are satisfactory because most of the values are closer to the perfect prediction line, and in good agreement with observed runoff even for large runoff events (Figure 5a). Method 1 yielded a larger NSE value of 75.84% and a smaller RMSE value of 3.26 mm in the calibration case, whereas the NSE value of Method 1 increases from −182.60% (the traditional SCS-CN method) to 77.45%, and the RMSE value decreased from 13.18 to 3.72 mm in the validation case (Table 3). Thus, we can conclude that Method 1 performed the best of the four methods for both the calibration and validation cases. However, the points of Method 1 around the 1:1 line are still scattered, which might be inherent in the CN method where the initial abstraction ratio λ is set to a standard value of 0.2. The assumption that λ = 0.2 is often questioned by many scholars [12,41]. Therefore, in order to improve the SCS-CN method, we introduced Method 2 with an optimized λ based on Method 1 to test the effect of initial abstraction on the proposed method. Moreover, we also optimized the parameters in Equation (6) (a1 and a2) (Method 3) based on Method 2 to test the slope factor on the proposed method as compared with the fixed values obtained from Huang et al. [15]. Figure 5b,c presents the comparison of runoff predicted by Methods 2 and 3 and the corresponding observed values, respectively. It can be seen that the runoff estimation is a great improvement over Method 1, with observation points close to the perfect line, especially when the depths of runoff are between 10 and 20 mm. The statistics also demonstrate the predictive capacity of Methods 2 and 3 with a higher model efficiency of 80.50% and 80.95%, respectively, in Table 3. However, there is a slight improvement with Method 3, which uses an optimized slope parameter (a1 = 213.99 and a2 = 25.38), as compared with Method 2, with NSE values of 80.73% vs. 80.44% in calibration and 81.21% vs. 80.58% in validation. The results indicated that the proposed methods with optimized λ (Methods 2 and 3) could further improve the SCS-CN method for runoff prediction in this study area.

Model Application
The tested standard SCS-CN and the proposed method (Methods 1-3) using the parameters calibrated and validated by the dataset of the three plots in XDG watershed were applied to estimate runoff of the remaining three plots in the CBG watershed. Figure 6 shows the measured versus estimated runoff values for the conventional SCS-CN and proposed method with the dataset in plots of the CBG watershed. The original SCS-CN method continued to underpredict most of the large as well as small storm-runoff events, whereas Method 1 performed better than the original SCS-CN method, with the former yielding a higher NSE (75.33%) and a smaller RMSE value (2.38 mm) ( Figure 6). However, Method 1 with λ = 0.2 still underestimated some small storm-runoff events as compared with Method 2 (λ = 0.001), which achieved a higher NSE value of 82.29% and had most of its data points visibly closer to the perfect line. Moreover, Method 2, using the parameters obtained from Huang et al. [15], with data from slopes ranging from 14% to 140%, performed better than Method 3 (NSE = 79.22%), which used optimized parameters with a narrow range of slope (26.7% and 70%). Therefore, from the results of the experimental sites, Method 2 is the most suitable SCS-CN method for surface runoff estimation in the plots of the two watersheds, which may also be applicable in the Loess Plateau.

Sensitivity Analyses
The above results indicate that Method 2 predicts runoff with greater accuracy than the other SCS-CN methods. A sensitivity analysis can identify the primary parameters that affect the performance of the model. In this study, sensitivity analyses were conducted by observing the effect of a variation in the calibrated parameter on the runoff prediction in terms of NSE, using the datasets of the plots in XDG watershed. Figure 7 presents the sensitivity of the runoff estimations to each parameter of Method 2. A sensitive variable is a parameter for which NSE changes dramatically as the parameter increases or decreases from the calibrated value. The most sensitive variable is apparently parameter b1, where the value of NSE sharply decreases from 80.58% to 47.83% as the b1 value changes in a small range of 130%-80% of the calibrated value. However, the initial ratio λ appears to be the least sensitive. In general, the parameters of soil moisture (b1 and b2) are the most sensitive, followed by parameters in storm duration (c) and slope equations (a1 and a2), and the least sensitive parameter is the initial abstraction ratio λ.

Discussion
From the above results, the original SCS-CN method had no good performance for the three plots of the XDG watersheds in both calibration and validation cases; it still underpredicted numerous rainfall-runoff events, and it also performed poorly when applied to the plots of the CBG watershed. The poor quality of the SCS-CN method might be due to the CN value that was taken from the USDA tables [42], which need improvement either by accounting for different factors that influence runoff or rebuilding the CN values based on different climate, soil, land use, and slope conditions from a large set of monitored rainfall-runoff data. In this study, several factors that influence runoff are considered to improve the CN value. Although the Huang et al. [15] and Huang et al. [34] methods, in which single factors of slope gradient and soil moisture were considered, respectively, could enhance the performance to some extent, the improvement is limited. Overall, the above three methods did not perform as well as the proposed method. The better performance of the proposed method shows it can accurately estimate runoff in the six plots with different vegetation types in the XDG and CBG watersheds using the optimized parameters. The results are also given some understanding of the physical processes represented by incorporating the slope, soil moisture, and storm duration factors into the proposed method, which eliminated the main drawbacks of the conventional SCS-CN model.

The Effect of Rainfall Duration
The characteristics of rainfall are the main driving force of runoff process, and different rainfall regimes will lead to different capacities of runoff production for each plot.
The 148 rainfall events of the six plots in the two watersheds (XDG and CBG) were divided into three categories using K-means clustering [43] based on rainfall depth and duration to test the performance of the SCS-CN method for different types of rainfall regimes ( Table 4). The classification conformed to the ANOVA criterion of a significant level (***p < 0.001). Rainfall regime 1 is a group of storm events with heavy intensity, short rainfall duration, and lower precipitation, and its occurrence frequency is the highest, accounting for 79.73% of the total events. Rainfall Regime 3 consists of storm events with low intensity, long-term duration, and low frequency (3.38%), while Rainfall Regime 2 is composed of storm events with moderate storm characteristics, i.e., shorter rainfall duration and lower precipitation than Rainfall Regime 3, but longer rainfall duration and higher precipitation than Rainfall Regime 1. In general, the descending order of average P and D is as follows: Rainfall Regime 3 > Rainfall Regime 2 > Rainfall Regime 1. This classification of rainfall regimes is consistent with Wei et al. [44] and Fang et al. [45], who also classified the rainfall into three rainfall regimes with 14 and 9 years of field measurements on the Loess Plateau, respectively.   (Figure 8a). This is because the SCS-CN method accounts for the rainfall amount but ignores the storm duration. The runoff predicted by the conventional SCS-CN method increased as the rainfall amount increased from Rainfall Regimes 1 to 3, while the storm duration also increased from Rainfall Regimes 1 to 3. Therefore, the measured runoff did not increase with the rainfall amount, which ultimately meant that the SCS-CN method performed poorly. However, when incorporating the storm duration in the proposed method with Equation (8), the underprediction of Rainfall Regime 1 and overprediction of Rainfall Regime 3 were removed and the performance of predicting Rainfall Regime 2 was improved, with most of its data points lying close to 1:1 as compared with the traditional SCS-CN method (Figure 8b). The better performance indicated that storm duration plays a vital role in rainfall-runoff generation and estimation, and the proposed method can account for different rainfall regime types of varying duration. Others have also confirmed that storm duration plays a significant role in both runoff production and prediction [45][46][47].

The Effect of Slope
As the slope increased, the initial abstraction, infiltration, and recession time of overland flow decreased, which resulted in less of an opportunity for infiltration and more surface runoff as compared with a horizontal surface. These drawbacks have been demonstrated by many researchers [48][49][50][51][52]. However, the CN2 of the standard SCS-CN method obtained from the handbook table [2] are based on 5% slope. Therefore, a slope equation (Equation (6)) was developed that incorporated the slope in the standard SCS-CN method for steep slope conditions.
In this study, the effect of slope is not reflected as the Huang et al. [15] method performed worse than the original SCS-CN method when accounting for the slope factor in the method. The reason is that the traditional SCS-CN has already overpredicted the storm events of Rainfall Regime 3, and the overprediction intensified when the slope factor was taken account in the Huang et al. [15] method because the runoff increased with slope gradient (Figure 4). However, when the overprediction was removed by the proposed method, the effect of slope could be found easily and as the slope increased the runoff increased.

The Effect of Soil Moisture
The antecedent moisture condition (AMC) plays a significant role in the progress of runoff. The conventional SCS-CN method divides the AMC into three levels based on the five-day antecedent rainfall amount. The measured values of CN should be continuous values for runoff prediction and are not limited to the three discrete AMC levels determined by the five-day AMC, which causes a sudden jump in runoff prediction [25]. Moreover, Koelliker [53] found that the probability of AMC2 occurring before a rainfall-runoff event was only 7%, rather than 50% [54], which indicated that the AMC value was underestimated and thereby resulted in runoff underprediction by the original SCS method. Huang et al. [34] indicated that there was no correlation between the CN value and the fiveday AMC, and the tabulated CN values were less than the measured ones in most cases, which further confirmed the runoff underprediction. Therefore, the five-day AMC used in the original CN method is unreasonable.
An alternative to soil moisture for determining the AMC would be to use direct measured or predicted soil moisture prior to a storm event, instead of antecedent rainfall, which is a vital factor in determining the initial abstraction of the SCS method. The rainfall-runoff prediction will be greatly enhanced by improving knowledge of antecedent moisture [29,30]. The results of the Huang et al. [34] method performed better than Huang et al. [15] and the standard SCS-CN method (Figure 4), which indicated that the equation of soil moisture (Equation (7)) is more suitable for AMC than the antecedent five-day rainfall.
In the new CN value, three factors of slope, soil moisture, and storm event were introduced, which represented the land condition, antecedent moisture, and characteristics of the storm event, respectively; all three factors have a great influence on runoff. However, as compared with the Huang et al. [15] and Huang et al. [34] methods, the performance of Method 1 was significantly improved after accounting for the storm duration factor (Table 3). Thus, storm duration has a greater effect on the runoff prediction of the SCS-CN method, which has frequently been ignored in other modified SCS-CN methods, followed by soil moisture and slope. Moreover, λ = 0.001 is more appropriate for the initial abstraction coefficient in this study area according to the comparison between Method 1 and 2; the same conclusion was also obtained by Huang et al. [34] for four plots in the Loess Plateau. The slope parameters in Method 2, which is obtained from Huang et al. [15] with data from slopes ranging from 14% to 140%, are also more suitable for application as compared with Method 3. This inference can be further confirmed by the sensitivity analysis, which also suggested that the initial abstraction ratio λ and the slope parameters (a1 and a2) are less sensitive, so Method 2 can be applied to similar sub-humid, semi-arid, and arid regions with the optimized parameters, but may need adjustment for humid regions because the soil moisture and rainfall may differ from the tested results in this study.

Conclusions
In this study, an equation was developed to improve the SCS-CN method by combining the CN value with the tabulated CN2 value and three introduced factors (slope gradient, soil moisture, and storm duration). Six models including the original SCS-CN method, Huang et al. [15], Huang et al. [34], and the proposed method with different optimized parameters (Methods 1-3) were used to test the reliability of the data from three runoff plots in the XDG watersheds of Loess Plateau. Subsequently, using the parameters calibrated and validated by dataset of the initial three runoff plots, the proposed method was then applied to runoff estimation of the remaining three runoff plots in CBG watershed. High NSE and low RMSE values during the calibration, validation, and application period indicated the proposed method can estimate runoff accurately for six plots in two watersheds and had greater reliability than the standard SCS-CN, Huang et al. [15], and Huang et al. [34] methods. Moreover, Method 2 with an initial abstraction ratio of λ = 0.001 rather than the standard value of 0.2 and slope parameters obtained from Huang et al. [15] with data from slopes ranging from 14% to 140% seems to be the most promising SCS method for runoff estimation in the experimental plots of the Loess Plateau. Furthermore, the parameters of soil moisture (b1 and b2) are the most sensitive, followed by parameters in storm duration (c) and slope equations (a1 and a2), and the least parameter is the initial abstraction ratio λ on the basis of the proposed method sensitivity analysis.