A CN-Based Ensembled Hydrological Model for Enhanced Watershed Runoff Prediction

A major structural inconsistency of the traditional curve number (CN) model is its dependence on an unstable fixed initial abstraction, which normally results in sudden jumps in runoff estimation. Likewise, the lack of pre-storm soil moisture accounting (PSMA) procedure is another inherent limitation of the model. To circumvent those problems, we used a variable initial abstraction after ensembling the traditional CN model and a French four-parameter (GR4J) model to better quantify direct runoff from ungauged watersheds. To mimic the natural rainfall-runoff transformation at the watershed scale, our new parameterization designates intrinsic parameters and uses a simple structure. It exhibited more accurate and consistent results than earlier methods in evaluating data from 39 forest-dominated watersheds, both for small and large watersheds. In addition, based on different performance evaluation indicators, the runoff reproduction results show that the proposed model produced more consistent results for dry, normal, and wet watershed conditions than the other models used in this study.


Introduction
The one-parameter traditional curve number (CN) model (CN model) developed by the U.S. Soil Conservation Service (SCS) now known as Natural Resources Conservation Service (NRCS), has enjoyed a long history of application as a lumped hydrological model.Its simplicity, versatility, and the availability of the necessary data have made it popular worldwide, as reported in [1], and an essential component in various hydrologic models, including water balance and storm routing models.Obviously, the CN model is reputable in the realm of applied hydrology, and CN is a mature concept that will remain in the forefront of engineering design.Its different versions, despite the complicated forms that result from introducing new parameters, cannot apply directly in approximating real situations because of the problem of model closure.Therefore, the CN model has been one of the most appealing and popular models for watershed runoff estimation for more than six decades [2].The unavailability of any other simple contender has allowed this model to enjoy a long application history [3].The widespread application of the CN model has led to its inclusion in hydrological software for surface runoff computations, such as CREAMS, CELTHYM, EPIC, HELP, L-THIA, PRZM, SWAT, SWIM [4], AGNPS, EPA-SWMM, GLEAMS, HEC-HMS, NLEAP, WinTR20, and WinTR55 [5].Furthermore, the absence of hydrologic gauging stations and the high cost of gauging station installations lead hydrologists to estimate the surface direct runoff for ungauged watersheds using various techniques, and the CN model plays the leading role in such a situation [6].Nonetheless, the CN model is not without its shortcomings.Despite its wide applications, the CN model has some limitations and misapprehensions caused mainly by its basic empirical assumptions [5].For example, it considers a three-step CN variation based on three antecedent moisture conditions (AMCs), dry, normal, and wet, which permit unreasonable sudden jumps in runoff estimation [7].Using tabulated CNs also underestimates most of the storm surface direct runoff in steep slope watersheds [8].Furthermore, application of the CN model for forested watersheds can result in an inaccurate estimate of runoff from a given storm rainfall [9][10][11].Thus, despite the long history, the CN model requires further study to develop and introduce a more robust model to better estimate runoff.Among other researchers, [1] and [11] modified the CN model to improve its runoff prediction capability.The modification proposed by [1] introduced storm-to-storm variation by incorporating a new parameter (M c ) to account for soil moisture before rainfall occurrence.However, this parameter sometimes results in negative values, which limits its applications [1].The simplified one-parameter modified model suggested by [11] claimed improved runoff prediction.Nevertheless, that model has been criticized for not preventing the sudden jumps in runoff estimation [4,5,12].To avoid sudden jumps in runoff estimation, simulated soil moisture (SM) has been employed using the soil water balance equation [13].Further the parameter S (maximum potential retention) of the CN model has been determined by means of an experimentally derived relationship between S and SM as S = a(1-θ e ), where θ e is the simulated relative SM at the beginning of the rainfall storm event and a is a parameter to be optimized.However, due to limited data in ungauged watersheds, neither θ e nor a can be estimated and hence this reduces the application of the CN model only for gauged watersheds.Recognizing the limitations of the traditional CN model and its modified versions, we set out to develop a new ensembled lumped model that accounts for a continuous initial abstraction and prevents sudden jumps in runoff estimation.In addition, by using a small number of required parameters, we intended for our new conceptualization to be less sensitive to CN variation and structurally more consistent than the traditional CN model and its modified versions.

Study Area and Data
To analyze the runoff reproduction, we selected 39 South Korean forest-dominated mountainous watersheds.The selected watersheds represent the overall hydro-meteorological setting of South Korea and their corresponding characteristics are described in [14].The selected watersheds vary in size from 42.32 to 888.01 km 2 and are characterized by low to high elevations (26 to 911 m above mean sea level) with average slopes from 7.50% to 53.53%.In Table 1, we differentiated the watersheds as small (Area ď 250 km 2 ) and large (Area > 250 km 2 ) following [9].The land cover is mainly forests, followed by agricultural and urbanized land as illustrated in Table 1.Loam and sandy loam are the two major soil textures with some fractions of silt loam.
We collected 30 min time-step rainfall data from the Korea Meteorological Administration (KMA) whereas observed discharge data at the same time step came from the Hydrological Survey Center (HSC) of Korea.We used land cover information from the Ministry of Land, Infrastructure, and Transport (MOLIT).The collected measurements constitute 1804 rainfall-runoff events from 39 watersheds between 2005 and 2012.Figure 1 gives the locations of the large watersheds in the study area; the small watersheds can be seen in [14].
To separate baseflow and direct runoff from discharge, we used the straight-line hydrograph method [15].We first screened the measured data to exclude small storm rainfall events (P < 25.4 mm) to prevent bias in estimating runoff using the traditional CN model [16].We applied the P 5 (the accumulated prior five days rainfall) criterion to determine the watershed antecedent condition to adjust the CN and its corresponding S values from normal to dry and wet conditions [12].We took the composite CN by applying the procedure documented in the NRCS [17] based on the watersheds' land cover characteristics, which we considered to be the same for all models.According to [18], Water 2016, 8, 20 of 17 hydrological models are more reliable when they reproduce satisfactory measurements, especially from watersheds with limited data.To make the models applicability more realistic in ungauged watersheds, we did not calibrate any of the parameters used in this study.Some statistics regarding the measured and predicted data are shown in Table 2.  Other parameters are defined in the text, T is the storm duration.
where Q, P, Ia, λ and S are surface direct runoff (mm), total rainfall (mm), initial abstraction (mm), initial abstraction coefficient (dimensionless), and potential maximum retention (mm), respectively.The standard Ia = 0.20S [19] and its one modified version for comparable improved runoff prediction is Ia = 0.05S [20].Different researchers have found structural inconsistencies in the traditional CN model, such as sudden jumps in S with corresponding runoff values for dry, normal, and wet conditions [11,12].Similarly the standard Ia = 0.20S has been found to be unreasonably high; and researcher across the globe suggested Ia < 0.2S (e.g., [1,14,[20][21][22][23]).For the CNM, we estimated the weighted composite CN for watershed characteristics and calculated the transformed S using the following Equation:

Mishra et al. Model (MRM)
To prevent unreasonable sudden jumps in runoff estimation, the relationship for Ia was modified by incorporating a new parameter (Mc) to account for soil moisture content before the rainfall occurrence [1].The modified Ia is expressed as:

The Traditional CN Model (CNM)
For a given P and its corresponding initial abstraction, I a , the traditional CN model in general is expressed as: where Q, P, I a , λ and S are surface direct runoff (mm), total rainfall (mm), initial abstraction (mm), initial abstraction coefficient (dimensionless), and potential maximum retention (mm), respectively.The standard I a = 0.20S [19] and its one modified version for comparable improved runoff prediction is I a = 0.05S [20].Different researchers have found structural inconsistencies in the traditional CN model, such as sudden jumps in S with corresponding runoff values for dry, normal, and wet conditions [11,12].
Similarly the standard I a = 0.20S has been found to be unreasonably high; and researcher across the globe suggested I a < 0.2S (e.g., [1,14,[20][21][22][23]).For the CNM, we estimated the weighted composite CN for watershed characteristics and calculated the transformed S using the following Equation:

Mishra et al. Model (MRM)
To prevent unreasonable sudden jumps in runoff estimation, the relationship for I a was modified by incorporating a new parameter (M c ) to account for soil moisture content before the rainfall occurrence [1].The modified I a is expressed as: Water 2016, 8, 20

of 17
The new parameter (M c ) can be found from the prior five days accumulated rainfall (P 5 ), λ and where β = 0.72 is the optimized coefficient.Despite some improvements in runoff prediction, the application of M c in [1] sometimes causes negative Q values, which is the major drawback of that formulation.Therefore, the model should be revisited to eliminate this complication when estimating runoff.For runoff estimation, the MRM can be expressed as: In Equation ( 4), M C " 0.72 ?P 5 S which is one of the best versions presented in [1] and I a is calculated using Equation (3) (λ = 0.2).

Michel et al. Model (MCM)
For λ = 0.2, the simplified forms of the MCM for dry (AMC-I), normal (AMC-II), and wet conditions (AMC-III), respectively, are as follows: Q " P ˆ0.48S `0.72P S `0.72P ˙( 6) Similarly, a simplification of the model parameters reported in [11] also considers a three-step S parameter and cannot eradicate the problem of sudden jumps in runoff prediction.

The Proposed Model (AJM)
To circumvent the sudden jumps and other inconsistencies, we here derived a new model ensembling the traditional CN and GR4J models [23] with the pre-storm moisture content concept [11] and further introduced a continuous S formulation.The GR4J (which stands for modèle du génie rural à 4 paramé tres journalier) is a French daily lumped four-parameter rainfall-runoff model that belongs to the family of soil moisture accounting models.It determines four parameters for daily runoff prediction: maximum storage capacity, groundwater exchange coefficient, one day prior maximum routing storage capacity, and time base of unit hydrograph.More details on the GR4J model are available in [23].
We replaced the M c in Equation ( 3) with V o adopted from [11], to make the following relationship: where V o is the soil moisture content before a storm event.Different values of V o are found for different AMCs [11] which are differentiated based on P 5 [14].Using the available CN adjustment formulae from normal (AMC-II) to dry (AMC-I) and wet (AMC-III) conditions, the expression for V o becomes [4]: where γ is the fraction to link V o and S that can be optimized from measured data.After testing the measured rainfall-runoff data from 140 French watersheds with different surface and variable moisture conditions, V o as some fraction of S was simplified for dry, normal and wet conditions as; V o = 0.33S, V o = 0.61S, V o = 0.87S, respectively.Taking λ = 0.20, as generally assumed for practical applications in the traditional CN model, and substituting V o from Equations ( 12)-( 14), respectively into Equation (8) gives: I a " 0.2S 2 S `0.61S " 0.124S for AMC-II ( 13) To incorporate the new initial abstraction formulation, we conceptualized a new model by combining the CN model, the PSMA procedure [11], and the event-based empirical GR4J model.To do this, we initialized V o as the soil moisture store level at the onset of an event and V as the moisture level corresponding to the total rainfall P [11].
Substituting Q from Equation (1) into Equation ( 15) yields Now the simplified form of the GR4J runoff model described in [11] can be expressed in its cumulative form as: Here PE is the daily potential evapotranspiration and is assumed negligible because the runoff from rainfall usually lasts for an event of sufficiently limited duration [24].Hence, the above equation becomes: Equation ( 18) yields Q = P for V = S + S a as a maximum capacity of V, where S a is an intrinsic parameter equal to; S a = V o + I a [11].Substituting the expression for V from Equation ( 16) into Equation (18) and simplifying yields Now, substituting I a = 0.15S from Equation ( 12) and V o = 0.33S (assuming the watershed dry condition), after simplification Equation (19) becomes: It is interesting that Equation ( 20) is a simple, one-parameter model like the traditional CN model.Using the traditional CN model, as in Equation (1), after substituting I a from Equation (13) could Water 2016, 8, 20 7 of 17 improve runoff prediction if S is taken as the transformed value from the dimensionless CN assessed as a function of land use/cover, soil types, and AMCs (differentiated from P 5 ).However, many researchers of various Mediterranean watersheds reported limitations and drawbacks to that approach, and it inconsistently displays pre-storm soil moisture temporal variations (e.g., [19,[25][26][27]).Some researchers (e.g., [26,28]) have specifically mentioned that the cumulative rainfall in the previous five days is not a good proxy of the wetness conditions of the catchment before a flood event.Similarly, the use of a fixed I a in the traditional CN model causes unusual sudden jumps in runoff estimation [1,4,11,12] For P 5 Ñ 8 then Q Ñ P which means the watershed is fully saturated by pre-storm soil moisture content and it has no more capacity to retain water.Hence all rainfall will be converted to runoff.In the case of S Ñ 0, then Q = P, whereas, if S Ñ P and P 5 = 0 (dry condition), then Q = 0.40P.Using this new model (Equation ( 22)) does not require any adjustment for S values from normal to dry or wet conditions, because the P 5 variation accounts for the I a adjustment associated to different pre-storm soil moisture responses in a watershed.
From this point, we investigated four models in this study: the traditional CN model (CNM), the model proposed in [11] as Michel model (MCM) and in [1] Mishra model (MRM), and our proposed model (AJM).Details of the different models are given in Table 3.  2), ( 5), ( 6) and ( 7) 2) and ( 22) Proposed Model -do-means as above.

Models' Goodness-of-Fit Evaluation
A vital step in evaluating and inferring model results is creating graphs to display them and determining goodness-of-fit indicator values to quantify model performance in terms of prediction accuracy [29].No formal standards evaluate goodness-of-fit to data visually or numerically through quantitative assessment.This leads to considerable variability in the techniques used, which are frequently decided based on applications in previous similar studies.According to [30], model evaluation (validation) and sometimes performance are judged by comparing calculated values with measured benchmark data.Commonly, evaluations of model performance are based on several statistics and techniques.To assess models' prediction accuracy quantitatively, the most commonly used methods are the root mean squares error (RMSE) [15,[31][32][33], Nash-Sutcliffe efficiency (NSE) [32][33][34], and Percent Relative Error (PRE) [35] expressed as: PRE " where Q oi and Q ci are the observed and estimated runoff values for storm events i to N and SD is the standard deviation of observed data.Smaller RMSE values (optimum = 0.0) depict better model runoff estimation.Here, Q o represents the mean of observed runoff values for storm events i to N. A model is deemed satisfactory if its NSE > 0.50 [35].Another recent study considered a hydrologic model satisfactory if NSE > 0.65 [36].
We further tested the appropriateness of each model based on the best-fit distribution and its corresponding statistic for the cumulative observed and modeled runoff from the 39 watersheds.To do so, we analyzed the observed and modeled data using three different tests: Kolmogorov-Smirnov test, Anderson Darling test, and chi-squared test [37].We used those tests to identify that the data belong to a specific distribution.We identified the best-to-worst fit distribution for those three tests from rank 1 to rank 18 after analyzing for 18 different distribution fittings.The model depicting the same distribution fit as exhibited by the measured runoff would be the best alternative for runoff prediction in the study area.Significant agreement between the measured and modeled runoff can be further validated in comparing models based on closeness of statistical values from any of the aforementioned three tests.
Beside the above evaluation indicators, to examine the model performance, the time series data for measured and modeled runoff were compared for two example watersheds selected from the study area.In addition, percent errors on peak discharge (E Q P ) [13] were also calculated to identify the physical behaviors of the models in estimating runoff.
where max (Q o ) and max (Q c ) are the maximum peak discharge values for the observed and modeled runoff.

Results and Discussion
To select an appropriate physical relationship between P 5 and I a we assessed the watersheds' runoff-producing response for the common discrete (fixed) I a = 0.2S and the proposed variable (continuous) formulation.This assessment was accomplished once the pre-storm soil moisture was adjusted in I a based on the recommendation from [11] and was subsequently adopted by other researchers in their studies (e.g., [4,5,12]).In our analysis, we kept the CN as a major parameter obtained from watershed characteristics under normal conditions fixed for all models.The measured rainfall-runoff events we used in the model assessments were characterized by a complex pattern to represent a diverse response.Considering P and CN as the fundamental and common parameters and other physiographic characteristics as shown in Table 1, all the models compared in this study both structurally differ from each other and lie within the parameter space of the traditional CN model.
Figure 2 illustrates the variation of fixed and variable I a used in the CNM and AJM respectively.For the CNM, the evident high I a and sudden jumps in S were the two major reasons for underestimating runoff from steep slope watersheds like in South Korea.As illustrated [9], the traditional CN model should be employed cautiously in watersheds larger than 250 km 2 .Keeping in mind this criteria, the model's analysis was carried out for small (Area ď 250 km 2 ) and large (Area > 250 km 2 ) watersheds as well as based on pre-storm watershed hydrological conditions (dry, normal and wet).As shown in Figure 3a, for small watersheds, the MCM showed no evident improvement in terms of mean (median) RMSE (mm) values of 22.66 (23.25) compared to the CNM values of 22.07 (20.15).However, the spread of RMSE values in those watersheds between the 25th and 75th percentiles depicted improvement followed by the MRM with mean (median) RMSE of 20.79 (18.78).Contrarily, the AJM indicated more accurate mean (median) RMSE of 13.13 (12.34).The improved results by the AJM indicate the reliability of the proposed variable I a in estimating runoff from mountainous watersheds compared to the CNM and the models presented in [1] and [11].
Water 2016, 8, 20 9 of 17 mind this criteria, the model's analysis was carried out for small (Area ≤ 250 km 2 ) and large (Area > 250 km 2 ) watersheds as well as based on pre-storm watershed hydrological conditions (dry, normal and wet).As shown in Figure 3a, for small watersheds, the MCM showed no evident improvement in terms of mean (median) RMSE (mm) values of 22.66 (23.25) compared to the CNM values of 22.07 (20.15).However, the spread of RMSE values in those watersheds between the 25th and 75th percentiles depicted improvement followed by the MRM with mean (median) RMSE of 20.79 (18.78).Contrarily, the AJM indicated more accurate mean (median) RMSE of 13.13 (12.34).The improved results by the AJM indicate the reliability of the proposed variable Ia in estimating runoff from mountainous watersheds compared to the CNM and the models presented in [1] and [11].To assess the model's efficacy in large watersheds, as shown in Figure 3b, it is similar to that for small watersheds, as depicted in Figure 3a.The mean (median) RMSE for the CNM, MCM, MRM, and AJM were 27.94 (24.85), 27.00 (25.86), 25.50 (23.48) and 17.40 (16.16), respectively.The MCM showed almost no improvement over the CNM and the MRM exhibited modest improvement, whereas the AJM indicated significantly lower RMSE values in all watersheds.The cumulative RMSE values for all types of events, shown in Figure 3c, support that the AJM has more reliable results than the other contenders.The consistently lower RMSE values evidenced by the AJM for all kinds of events (combined dry, normal, and wet conditions from 39 watersheds) indicate its superiority.To assess the model's efficacy in large watersheds, as shown in Figure 3b, it is similar to that for small watersheds, as depicted in Figure 3a AJM indicated significantly lower RMSE values in all watersheds.The cumulative RMSE values for all types of events, shown in Figure 3c, support that the AJM has more reliable results than the other contenders.The consistently lower RMSE values evidenced by the AJM for all kinds of events (combined dry, normal, and wet conditions from 39 watersheds) indicate its superiority.
Analyzing the runoff prediction efficacy illustrated by different models using NSE as the performance indicator demonstrated modest improvement by the MCM followed by the MRM, as shown in Figure 4a.The mean NSE (0.57) for the MCM showed almost no improvement compared to the CNM mean NSE (0.57), but the spread of the NSE from this figure indicates good model performance by the MCM.Moreover, the MRM mean NSE (0.63) indicated some increases in the model's overall efficacy compared to the CNM and MCM.Nevertheless, the AJM had the highest NSE values (mean 0.85) indicating its reliability for accurate runoff estimation in small watersheds.Figure 4b evidences a similar pattern of performance for large watersheds.Here, the mean NSE for the CNM, MCM, MRM, and AJM were 0.64, 0.65, 0.71, and 0.86, respectively.The higher NSE values in both small and large watersheds confirmed the AJM enhanced runoff prediction.According to the NSE criteria in [35], the CNM, MCM, MRM, and AJM had respectively, good to very good results 46.67%, 33.33%, 60.00%, and 100% of small watersheds.Similarly, the models respectively fulfilled the same performance rating in 41.67%, 54.17%, 58.33%, and 100% of large watersheds.Considering the performance rating criteria in [36], the good to very good performance by the CNM, MCM, MRM, and AJM respectively, was in 26.67%, 6.67%, 26.67%, and 66.67% of small watersheds.Likewise, those ratings respectively were met in 25.00%, 16.67%, 37.50%, and 87.50% of large watersheds.
Water 2016, 8, 20 10 of 17 Analyzing the runoff prediction efficacy illustrated by different models using NSE as the performance indicator demonstrated modest improvement by the MCM followed by the MRM, as shown in Figure 4a.The mean NSE (0.57) for the MCM showed almost no improvement compared to the CNM mean NSE (0.57), but the spread of the NSE from this figure indicates good model performance by the MCM.Moreover, the MRM mean NSE (0.63) indicated some increases in the model's overall efficacy compared to the CNM and MCM.Nevertheless, the AJM had the highest NSE values (mean 0.85) indicating its reliability for accurate runoff estimation in small watersheds.Figure 4b evidences a similar pattern of performance for large watersheds.Here, the mean NSE for the CNM, MCM, MRM, and AJM were 0.64, 0.65, 0.71, and 0.86, respectively.The higher NSE values in both small and large watersheds confirmed the AJM enhanced runoff prediction.According to the NSE criteria in [35], the CNM, MCM, MRM, and AJM had respectively, good to very good results 46.67%, 33.33%, 60.00%, and 100% of small watersheds.Similarly, the models respectively fulfilled the same performance rating in 41.67%, 54.17%, 58.33%, and 100% of large watersheds.Considering the performance rating criteria in [36], the good to very good performance by the CNM, MCM, MRM, and AJM respectively, was in 26.67%, 6.67%, 26.67%, and 66.67% of small watersheds.Likewise, those ratings respectively were met in 25.00%, 16.67%, 37.50%, and 87.50% of large watersheds.For cumulative events, all AMC-I, AMC-II and AMC-III events from 39 watersheds, the CNM showed improvement from AMC-I to AMC-II and AMC-III.However, the MCM exhibited improvement in AMC-II events compared to AMC-I but worse results in AMC-III.The MRM increased in efficiency from AMC-I to AMC-II and AMC-III events.Nonetheless, the higher and consistent NSE values exhibited by the AJM indicated superior performance in all watersheds for all types of events.This performance can be attributed to the incorporation of the variable Ia in the CN model in the AJM.The RMSE and NSE cannot indicate over-or underestimation of the model runoff.Therefore, we assessed the models' performances using the PRE statistic.Figure 5a,b show that the CNM most underestimated all events, followed by the MRM.Contrarily, the MCM overestimated runoff for the majority of the events.Both the CNM and MCM were inferior because their PRE values are in the unsatisfactory range (PRE ≥ +25 or PRE ≤ −25) [14,35].The mean PRE (%) in small (large) watersheds for CNM, MCM, MRM and AJM respectively were 36.55 (34.64), −23.01 (−20.83),13.04 (16.59), and 5.66 (8.26).Evidently, both in small and large watersheds, the mean PRE values for the AJM were in the very good performance range [14,35].Using the PRE statistics, all models show better performance in large watersheds than in small watersheds in the mountainous parts of the study area.For not only the cumulative events from all 39 watersheds, but also the cumulative AMC-I, II, and III events, the AJM more accurately and consistently predicted the runoff.On the contrary, the other models were inconsistent in their runoff predictions.The performance ratings of different models for runoff prediction are given in Table 4. Using different statistical indicators, the proposed model performed significantly better than the other models investigated.One reason for the CNM's inferior results might be the watersheds slope difference.The CNM was derived from rainfall-runoff data from watersheds with up to 5% slope [15], whereas the watersheds in our study area range between 7.50% and 53.53% slope (refer to Table 1).In addition, no statistical indicator showed any significant overall performance difference in predicting runoff from small or large watersheds, which For not only the cumulative events from all 39 watersheds, but also the cumulative AMC-I, II, and III events, the AJM more accurately and consistently predicted the runoff.On the contrary, the other models were inconsistent in their runoff predictions.The performance ratings of different models for runoff prediction are given in Table 4. Using different statistical indicators, the proposed model significant results for runoff estimation.The close agreement between the measured runoff values and those modeled by the AJM can be seen from the best-fit line approaching the 1:1 line with a high coefficient of determination (R 2 = 0.88) compared to that for the CNM (R 2 = 0.75), MRM (0.76), and MCM (R 2 = 0.82).
A time series plot (Figure 7a,b) compares the measured and modeled runoff values for the selected events from Cheoncheon and Donghyang watersheds highlighting the performance of different models used in this study.It was found that the AJM is comparatively more consistent than the other models.
Figure 8 shows the model performance in individual watersheds and represents their cumulative frequency distribution for the total 39 watersheds in this study.The maximum peak discharge values were calculated after estimating the peak discharge for the individual storm events in the respective watersheds.The evaluation based on the percentage error in peak discharge, E Q P (Figure 8a) as well as the percentage absolute error in peak discharge, ˇˇE Q P ˇˇ(Figure 8b), indicated that the runoff estimated using the proposed model is in better agreement with the measured data when compared to other models.Similarly, the CNM, MCM, MRM, and AJM depicted the overall mean (median) percentage absolute peak discharge ˇˇE Q P ˇˇfor the 39 watersheds as 16.82 (19.19), 26.15 (25.66), 18.30 (15.99), and 14.33 (12.81), respectively.From these statistical scores the proposed model can be encouraged for application in the study area as well as other hydrological similar areas.

Conclusions
We tested the rainfall and runoff data collected for 1804 selected events from 39 mountainous watersheds for runoff reproduction using our proposed model and three earlier contenders.Because it underestimated the runoff for the majority of events, the CNM was found inferior to the other models.This might be attributable to the fixed initial abstraction and the absence of a credible parameter for pre-storm soil moisture, which results in sudden jumps in runoff estimation.After testing data from 140 French watersheds [11], researchers conceptualized an intrinsic parameter and incorporated it into the CN model to enhance its runoff prediction capability.However, their simplified model was unable to circumvent sudden jumps in runoff prediction and hence overestimated the majority of events.The model's amendment presented in [1] was also unable to handle this inconsistency and provided only modest improvement on the CNM.
In the proposed model, we conceptualized a CN-based ensembled approach by amending the previously suggested formulation [1] for initial abstraction after incorporating the pre-storm soil moisture adjustment [11].To circumvent sudden jumps in runoff estimation, we used a continuous and variable initial abstraction that depends not only on S, but also takes into account the storm magnitude and the prior rainfall effects irrespective of the three AMCs.Our modification resulted in a structurally more consistent model.Our proposed model has shown a high degree of reliability in predicting runoff from a majority of watersheds in the study area irrespective of the watershed size.Because of its simplicity and statistically reliable performance, the proposed model can be incorporated fruitfully in continuous hydrologic modeling.Nevertheless, a thorough investigation is required for watersheds of other biomes.In addition, the current study is limited to the prediction of runoff from large storm rainfall events (P ≥ 25.4 mm).The application of the proposed model can be tested in future studies by estimating runoff from small storm rainfall events (P < 25.4 mm).It could be more interesting to compare runoff estimation based on watershed scale and the partial source areas.

2. 2 .
Development of a New Hydrological Model2.2.1.The Traditional CN Model (CNM)For a given P and its corresponding initial abstraction, Ia, the traditional CN model in general is expressed as:

Figure 2 .
Figure 2. Variation of fixed and variable initial abstraction (Ia) for rainfall events in (a) dry condition; (b) normal condition; and (c) wet condition.

Figure 2 .
Figure 2. Variation of fixed and variable initial abstraction (I a ) for rainfall events in (a) dry condition; (b) normal condition; and (c) wet condition.

Figure 3 .
Figure 3. Model performance based on root mean squares error (RMSE) (a) in small watersheds (Area < 250 km 2 ); (b) in larger watersheds (Area > 250 km 2 ) and (c) for different types of events (dry, normal, and wet).

Figure 3 .
Figure 3. Model performance based on root mean squares error (RMSE) (a) in small watersheds (Area < 250 km 2 ); (b) in larger watersheds (Area > 250 km 2 ) and (c) for different types of events (dry, normal, and wet).

Figure 4 .
Figure 4. Model performance based on Nash-Sutcliffe efficiency (NSE) (a) in small watersheds (Area ≤ 250 km 2 ); (b) in larger watersheds (Area > 250 km 2 ) and (c) for different types of events (dry, normal and wet).The circled NSE value for the CNM shows the second overall value (the minimum was −0.05, which is not shown here).

Figure 4 .
Figure 4. Model performance based on Nash-Sutcliffe efficiency (NSE) (a) in small watersheds (Area ď 250 km 2 ); (b) in larger watersheds (Area > 250 km 2 ) and (c) for different types of events (dry, normal and wet).The circled NSE value for the CNM shows the second overall value (the minimum was ´0.05, which is not shown here).

Figure 6 .
Figure 6.Comparison of events from all watersheds for (a) observed runoff (Qo) vs. predicted with CNM; (b) Qo vs. predicted with MCM; (c) Qo vs. predicted with MRM; and (d) Qo vs. predicted with AJM.

Figure 7 .
Figure 7.Comparison between the measured runoff time series values and those obtained from application of four models for (a) Cheoncheon watershed; and (b) Donghyang watershed.

Figure 6 .
Figure 6.Comparison of events from all watersheds for (a) observed runoff (Q o ) vs. predicted with CNM; (b) Q o vs. predicted with MCM; (c) Q o vs. predicted with MRM; and (d) Q o vs. predicted with AJM.

Figure 6 .
Figure 6.Comparison of events from all watersheds for (a) observed runoff (Qo) vs. predicted with CNM; (b) Qo vs. predicted with MCM; (c) Qo vs. predicted with MRM; and (d) Qo vs. predicted with AJM.

Figure 7 .
Figure 7.Comparison between the measured runoff time series values and those obtained from application of four models for (a) Cheoncheon watershed; and (b) Donghyang watershed.

Figure 7 .
Figure 7.Comparison between the measured runoff time series values and those obtained from application of four models for (a) Cheoncheon watershed; and (b) Donghyang watershed.

Figure 8 .
Figure 8. Model evaluation for 39 watersheds based on: (a) cumulative distribution of percentage error in peak discharge ( P Q E ), and (b) cumulative distribution of absolute percentage error in peak

Table 1 .
Watersheds and data description.

Table 2 .
Statistics of model outputs from the combined data set from 39 watersheds.
. To circumvent this issue, we present a continuous spatio-temporal expression for the initial abstraction by replacing the S value in the proposed model with a new expression; S

Table 3 .
Models and their corresponding parameters identification.

Table 4 .
Model performance ratings in 39 watersheds based on different statistical indicators.

Table 5 .
Comparison of observed and modeled runoff based on best-fit distributions (statistic).