Next Article in Journal
Assessment of Soil–Structure Interaction Approaches in Mechanically Stabilized Earth Retaining Walls: A Review
Previous Article in Journal
Effect of Different Irrigation Managements on Infiltration Equations and Their Coefficients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Presentation of a Rainfall–Runoff Retention Model (3RM) Based on Antecedent Effective Retention for Estimating Runoff in Seven Basins in Iran

by
Shayan Shamohammadi
1,
Ahmad Reza Ghasemi
1,
Kaveh Ostad-Ali-Askari
2,3,* and
Saeedeh Izadi
1
1
Water Engineering Department, College of Agriculture, Shahrekord University, Shahrekord 8818634141, Iran
2
Department of Water Engineering, College of Agriculture, Isfahan University of Technology, Isfahan 8415683111, Iran
3
Department of Natural Sciences, Manchester Metropolitan University, Chester Street, John Dalton Building, Manchester M1 5GD, UK
*
Author to whom correspondence should be addressed.
CivilEng 2023, 4(3), 966-981; https://doi.org/10.3390/civileng4030052
Submission received: 6 May 2023 / Revised: 23 July 2023 / Accepted: 28 July 2023 / Published: 29 August 2023
(This article belongs to the Section Water Resources and Coastal Engineering)

Abstract

:
This study aims to correct and assess the SCS-CN model. In this research, the 3RM model (written by Shamohammadi) has been modified in such a way that the maximum primary retention (I), maximum secondary retention (Fmax), and basin potential retention (Smax) can be calculated using precipitation (Pa). The purpose of this study is to evaluate the total retention model ( S t = f ( F m a x , S m a x , p a ) ) and the runoff model ( Q = f ( S t , p a ) ) using the mountain basins of Iran, including Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin. The results showed that the primary retention, maximum secondary retention, and retention capacity are, respectively, 2.3, 30.4, and 32.7 mm in Imamah, 2.5, 48.6, and 51.1 mm in Kasilian, 2.4, 26.7, and 29.1 mm in Navrood, 3.2, 21.5, and 24.7 mm in Darjazin, 1.7, 15.0, and 16.7 mm in Kardeh, 2.5, 33.2, and 38.1 mm in Khanmirza, and 4.9, 44.5, and 50.6 mm in Mashine. Additionally, the λ (ratio of primary retention to potential retention) values for all basins are less than 0.2 (suggested by SCS) and vary between 0.05 in Kasilian and 0.1 in the Darjazin, Kardeh, and Mashine basins. The results of fitting the model to the rainfall-runoff data showed that the evaluation indices, including the coefficient of determination (R2), Nash–Sutcliffe (NS), and root mean square error (RMSE), for predicting the runoff in the basins varied between 0.78 to 0.96, 0.78 to 0.961, and 0.86 to 2.28, respectively. According to the obtained results, it can be concluded that the model has an acceptable ability to predict runoff for all the studied basins.

1. Introduction

With the progress of computer science, it is necessary to develop hydrological models along with the development of information technology, and instead of using descriptive methods and curve numbers, use data-driven models to predict floods and runoff [1]. One of the widely used methods in runoff estimation is the SCS-CN experimental method [2]. In this method, the classification of soil moisture conditions in each area is carried out only according to five-day antecedent rainfall and the growing and dormant season; this classification has a significant effect on determining the curve numbers and, finally, the runoff [3,4,5].
A lot of research has also tried to modify and optimize the SCS-CN method, especially in case of antecedent moisture conditions [6,7,8,9,10,11] indicated that, in the SCS-CN method, the use of a revised soil-moisture index instead of the use of antecedent moisture had better results. Verma et al. [12], developed a soil moisture accounting (SMA)-based SCS-CN method to avoid sudden jumps in estimated runoff. Ogarekpe et al. [13], used modified SCS models for predicting runoff in the Oraimiriukwa River Basin. They showed that the ratio of initial interception to maximum potential retention is not constant and is a function of precipitation.
Singh et al. [14], showed that one of the most important limitations of the SCS-CN method is that slope is not considered. In the modified SCS-CN method, the slope is considered to enhance the accuracy of the modeled results.
Many researchers also considered the SCS-CN model to be ineffective and showed that the SCS-CN method does not have the necessary efficiency for Karst landforms and forest areas [15,16,17]. Vaezi and Abbasi [18], investigated the efficiency of the SCS-CN curve number method for estimating the runoff in the Tehmchai Basin (in the northwest of Iran) and showed that the λ was <0.2. As a result, they suggested that the coefficient λ = 0.2 is not correct and the value of λ needs to be calibrated. These results were also repeated for the Kardeh Basin by Ebrahimian et al. [19].
Shamohammadi and Zomorodian [20], compared the performance of the SCS and Bennett Soil Moisture Accounting (SMA-B) models in estimating the flood in the Roodzard Basin in Iran. The results of this research indicated the inefficiency of the SCS model when compared to the SMA-B model. The results of the studies of Sadeghi et al. [21] also showed that the SCS-CN method is weak in estimating the runoff of the Emameh, Kasilian, Darjazin, and Khanmirza basins (in Iran). Similar research was also conducted by Asadi et al. [22], Arand and Torabi [23], Shammohammadi [24], and Orak and Farhadi [25] for the Mashin and Navrood basins, and it shows the weakness of the SCS-CN method for estimating runoff.
Shamohammadi [26], showed that the mathematical model of SCS-CN has a theoretical weakness and is basically unable to describe the conceptual model of SCS-CN. He stated that each basin has only one potential retention (Smax or Sp) and one secondary potential retention (Fmax). Both of them are obtained through their variable parameters (F and St, respectively).
Shamohammadi wrote his mathematical model based on the function variable F (secondary retention) and used the constant coefficient Ksh [26]. The most important weakness of the Shamohammadi model is that it was necessary to use the linear regression between rainfall-runoff data to calculate the initial retention (I). In the current research and in the Shamohammadi model (modified 3RM), instead of F, the total retention function (St) is used, and instead of the Ksh coefficient, the secondary maximum retention (Fmax) is used. As a result, while reducing the error in the calculation of initial retention (I), the values of I, Fmax, and Smax are directly extracted from the model (innovation of the present research), and the model also becomes simpler. In order to evaluate the new model, data from different areas of Iran (seven basins with different areas, different climates, and different landforms in the important mountain ranges of Iran) were used.

1.1. An Introduction to the Model

1.1.1. Definitions

I: Initial retention: this includes all the factors that prevent the creation of runoff (water runoff, impoundments, infiltration, evaporation, etc.). This parameter is constant for a basin and is equal to the sum of the initial retention (Ia) and the antecedent effective retention (IER) (I = Ia + IER). It is obvious that the value of I varies from each type of rainfall to another, but it is generally a constant value for a basin.
F: Secondary retention (in this study, secondary retention is used instead of actual retention): this includes all the factors that prevent runoff. The only difference between F and I is that F occurs at the beginning of the rainfall event, while I is related to retention before runoff. The value of F is variable in basins.
St: Total retention: this is equal to F + I, so St is also variable.
Smax: The maximum retention or potential retention (Sp) is the maximum retention after which rainfall is completely converted into runoff. Theoretically, when the rainfall tends to infinity, all the impoundments are filled, the soil is saturated, and even the air humidity reaches the saturation limit. As a result, after retention is equal to Smax, all rainfall converts to runoff. In this case, the slope of the rainfall-runoff curve will be equal to 1. It is obvious that the potential retention is an ideal and hypothetical condition for the dry state of the soil and can only be calculated through the model. In the SCS-CN method, the potential retention is estimated using the curve number, which, in the correct case, is a real number (it is part of the storage potential) and cannot be storage potential. This also applies to Fmax.
Fmax: It is defined as Smax and is linearly constant. Its relationship with Smax is as F m a x = S m a x I .
IER: Antecedent effective retention: this is a part of retention due to antecedent rainfall (PA), which is effective in the rainfall-runoff system and is calculated as per Equation (1) [27]:
I E R = P A Q A + n 1 n t E T 0 A
where PA, QA, and (ET0)A are the antecedent rainfall depth, antecedent runoff depth, and the daily antecedent potential evapotranspiration, respectively, and n1 and nt are the days of the beginning of the antecedent and the target, respectively.
Additionally, it can always be written as per Equation (2):
Q = P a S t
P a = P + I E R
where Q and Pa are the runoff depth (mm) and the modified target rainfall depth (mm), respectively.

1.1.2. Assumptions of the Shamohammadi Model

Equations (3) and (4) have been presented for the mathematical description of the conceptual curve in Figure 1.
d Q d p a = S t S m a x
d S t d p a = K s h K s h + p a
In Equation (4), the Ksh is the Shamohammadi constant. Additionally, d Q d p a and d S t d p a are the differential of runoff (compared to corrected rainfall) and the differential of total retention to rainfall, respectively. The other parameters are already defined.
On the other hand, based on the law of the survival of the fittest, it follows that
P a = S t + Q
Equation (5) is a general equation. Obviously, when the runoff is equal to 0, Pa = St = I, and when the total retention is completed and when rainfall is converted to runoff, then Pa = Q.
Using differentiation on both sides of Equation (5) with respect to rainfall, we obtain the following:
d P a d P a = d S t d P a + d Q d P a
Finally, Equation (6) can be written as Equation (7):
1 = d S t d P a + d Q d P a
Considering Equations (3), (4), and (7), we can write the following as Equation (8):
1 = k s h k s h + P a + S t S m a x
and then (Equation (9)):
S t = S m a x P a k s h + P a
According to Equation (9), when the rainfall tends to I (PaI), the value of Ksh will be equal to Fmax because when the rainfall is equal to I, the value of St will also be equal to I (F = 0) as a result:
{ I = S max l k sh + I S max = k sh + I k sh = S max I
However, based on the conceptual model in Figure 1 and the definition of Fmax, the value of ksh is equal to Fmax. Therefore, Equation (9) is modified as per Equation (10). Figure 2 shows the schematic of the maintenance of Equation (10).
S t = S m a x P a F m a x + P a
According to the relationship between runoff and total retention (Equation (2)), the amount of runoff depth is equal to what follows in Equation (11):
Q = P a S m a x P a F m a x + P a
The first advantage of model (11) compared to Shamohammadi’s model [27] is that, in the current model, the value of I is simply obtained (SmaxFmax). In this case, in addition to the error caused by linear regression being removed, the value of Fmax is also easily calculated by the model. As a result, model (11) is introduced as the modified model of Shamohammadi, and by applying boundary conditions, the correctness of the model (theoretically) is also proven because lim P a 0 ( Q ) = 0 , lim ( P a Q ) = and lim ( P a I Q ) = 0 .
Therefore, model (11) is also correct in limiting the conditions (based on the definitions and assumptions of the model) and is in harmony with the conceptual model (Figure 1) and the SCS-CN conceptual model [26].
This is despite the fact that, in the SCS-CN model, when the rainfall tends to 0 (Pa→0), the runoff depth (Q) does not become zero but becomes equal to 0.05S (Equation (12)). However, in reality, when the rainfall is zero, there is no runoff. As a result, it can be said that the SCS-CN model has a theoretical weakness in its hypotheses [26]. Even if it is claimed that the SCS-CN method is an experimental method, the influence of the theoretical weakness of the model on the experimental parameters cannot be ignored. Perhaps for this reason, the SCS-CN mathematical model is not defined for rainfall less than I (0.2S) following Equation (12):
Q = ( P 0.2 S ) 2 P + 0.8 S   P 0.2 S
Equations (10) and (11) are the basis equations for runoff estimation in this study. As mentioned before, the conceptual model of Figure 1 is the SCS-CN conceptual model. Only its hydrological parameters have been revised and redefined. Therefore, although Shamohammadi’s mathematical model is completely consistent with the SCS-CN conceptual model, it is completely different from the SCS-CN mathematical model (12) (both in terms of the relationships between the parameters and in terms of the method used for solving the problem). As a result, some of the SCS-CN limitations are also present in the current model, some of which are mentioned below:
  • The basin area should be so small that the whole basin is covered by rainfall (uniformly and at a constant intensity);
  • The intensity of different rainfalls should be equal to each other so that the runoff is only caused by the depth of the rainfall. Increasing the intensity of the rainfall causes runoff before the completion of primary retention; on the other hand, it causes an increase in mud in the runoff; thus, the amount of runoff at the hydrometric station is measured more than the actual runoff;
  • This model (SCS-CN) was developed for natural basins. In a basin (especially large basins), if there are man-made structures such as earthen dams, pools, large dry streams, or any unnatural structures, the data are not based on the theory of the conceptual model theory and, as a result, the models will not be able to correctly analyze the situation;
  • Accurate results can be expected from the model (11) when the rainfall-runoff curve follows the conceptual curve of Figure 1. Otherwise, the rainfall-runoff data must be checked carefully, and the incorrect data (for example, when the rainfall did not cover the entire basin) must be corrected or removed.
Of course, the limitations of the current model are much less than the SCS-CN model because it is based on historical data and operates as a black box model. In other words, in this model, in addition to the man-made structures, the changes in the parameters from one rainfall event to another also cause errors in the results.

2. Material and Methods

2.1. Studied Areas

In this study, the basins of Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin, which were equipped with rain gauges and hydrometer stations, were studied. The selection of these basins was based on the quality of rainfall and runoff data, the distribution of rainfall, climate diversity, vegetation diversity, area size, and geological conditions. Figure 3 shows the location of the studied basin in Iran. The general characteristics of the basins are also shown in Table 1.
All the meteorological and rainfall-runoff data were obtained from the Iranian Meteorological Organization and Iranian Water Resources Company, covering the period of 1985–2021. Due to the fact that an earthen dam was built in the Mashin area in 2001 and caused a change in the flood conditions of the basin, the data of the Mashin station from 1978 to 2001 were studied. In this research, Landsat 8 satellite images from 2013 were used to study land use and basin features.
In the present study, to simplify model (1), all the days between the beginning of the antecedent rainfall and the beginning of the target rainfall were added together, and 315 days were obtained. Additionally, the total potential evapotranspiration for all these days (315 days) was 159 mm. As a result, the average of potential evapotranspiration was approximately equal to 0.5 mm/day ( E T 0 = n 1 n t E T 0 A 315 0.5 m m / d a y ). Then, the average of potential evapotranspiration was used. On the other hand, no runoff has occurred in any of the antecedent rainfalls. As a result, based on [26], in this research, Equation (1) was used to simplify Equation (13).
I E R ( m m ) = P A ( m m ) 0.5 D ( d a y )
where D is the number of days of evapotranspiration (ntn1).
The maximum number of days for calculating evapotranspiration before the target precipitation was considered to be 5 days [2,28].

2.2. Model Evaluation

The following steps were taken to run the model:
  • The floods for which the associated rainfall covered the entire basin were selected. This depended on the number of rain gauge stations in the basin (between five rain gauges in the Emameh basin and 15 in the Mashin basin);
  • The rainfall data for which the amount of flood was as high as the target rainfall (P, Q) were selected. Then, based on the start date of the target rainfall (P), the rainfalls that occurred (at most) 5 days before (PA) were also determined. Because the antecedent rains did not have runoff, only the antecedent rain was recorded;
  • In the case of rainfall-runoff data that did not match the conceptual curve (mostly related to large basins, especially in the Mashin Basin), the data were corrected;
  • Model (11) was calibrated using half of the data. Then, using the calibrated model, the rest of the data were examined, and finally, the model was evaluated using the evaluation indices;
  • Three evaluation criteria, including the coefficient of determination (R2), root mean square error (RMSE), normalized root mean square error, and the Nash–Sutcliffe coefficient (NS), were used (Equations (14)–(16)).
R 2 = ( i = 1 N Q o b s Q o b s ¯ Q s i m Q s i m ¯ i = 1 N Q o b s Q o b s ¯ i 2 Q s i m Q s i m ¯ i 2 ) 2
R M S E = 1 N ( i = 1 N ( Q o b s Q s i m ) i 2
N S E = 1 i = 1 N Q o b s Q s i m i 2 i = 1 N Q o b s Q o b s ¯ i 2
where Qobs is the amount of observed runoff, Qsim is the amount of simulated runoff (model), and N is the amount of data.

3. Results and Discussion

The analysis of the rainfall-runoff data of the Emameh (35 data pairs), Kasilian (73 data pairs), Navrood (24 data pairs), Darjazin (66 data pairs), Kardeh (68 data pairs), Khanmirza (28 data pairs), and Mashin (data pair 54) basins showed that 34% of all data had antecedent rainfall. The antecedent rainfall values varied between 8.7 and 1.2 mm. Maximum and minimum values were observed in the Kasilian Basin (30 October 1999) and Khanmirza Basin (25 October 2005), respectively. Additionally, the highest and lowest IER values occurred in the Kasilian (8.4 mm) and Darjazin (0.5 mm) basins, respectively. The ET0 value ranged from 0.23 mm/day in the Navrood Basin (22 March 2001) to 0.65 mm/day in the Mashin Basin (28 September 1995).
The most important challenge of this research was to select precipitations that could meet the conditions of the model. This issue is more common in large basins. One of the most challenging basins was the Mashin Basin. Figure 4 shows the isohyetal map, hyetograph, and flood hydrograph for the Mashin Basin for 1–2 April 1996. As can be seen, the depth of precipitation increases from the outlet of the basin to the north and east (Figure 4a). Additionally, the rainfall intensity (Figure 4b) has severe fluctuations during the rainfall period (23 h): two peaks of precipitation at 4 to 8 h and 20 to 23 h can be observed. Therefore, both the non-uniform intensity and depth of precipitation over 23 h are the causes of error in the model. However, due to the fact that the rainfall covers the whole basin, this error did not have much effect on the performance of the model.
One of the factors that caused the flood lag time to decrease was antecedent rainfall. The average rainfall was 71.3 mm based on the isohyetal map. On the other hand, the antecedent rainfall with a depth of 7 mm happened 5 days before. The sum of potential evapotranspiration was calculated to be 2.5 mm (0.5D = 0.5 × 5 = 2.5 mm) and the IER = 4.5 mm. As a result, the total amount of corrected precipitation was equal to 75.8 mm. For this reason, the runoff coefficient in this event was equal to 0.57. It is noteworthy to mention that the Mashin Basin has a Karst landform [23,25].
On the other hand, due to the low slope of the basin (7.3%) and the large area (881 km2), the obtained runoff coefficient is very large. The most important reason for this issue is probably (1) the effect of the antecedent rainfall and (2) the high-intensity rainfall. This rainfall happened in March. It was expected that the melting of snow would be involved in the volume of runoff, but the investigations showed that the height of the snow border in April in the Mashin Basin was above 2700 m [24] as a result, the effect of snow melting on the volume of the flood was very low (see Figure 5). These results are consistent with [23,25,29].
Table 2 shows the results of the evaluation model using the rainfall-runoff data of the Imamah Basin. The basin has the smallest area among the studied basins (Table 1).
As can be seen in Table 2, out of the 35 rainfall-runoff data, 13 antecedent rainfall events occurred 5 days before the target rainfall (the rainfall events that have antecedent rainfall are marked with *). The results show that the potential retention and maximum secondary retention are estimated as 32.72 and 30.43 mm, respectively, and from the difference of these two values, the value of initial retention (I) is equal to 2.288 mm. Additionally, the value of λ (I/Smax) is calculated to equal 0.07. The value of λ for the Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins was 0.050, 0.08, 0.1, 0.1, 0.07, and 0.1, respectively. All the values obtained are less than the value recommended by the SCS-CN method. Additionally, based on the previous studies, the value of λ in the Emameh, Kasilian, Darjazin, Khanmirza, and Mashin basins have been reported to be about 0.09, 0.16, 0.2, 0.3 [30], and 0.2 [20], respectively. As can be seen, in all these studies, the values for λ are different, and it is not possible to consider a fixed number for it. Further investigations showed that there is no significant correlation between the values of λ and Smax [17,30]. Therefore, it is necessary to examine the values of I and Smax separately.
Figure 6 shows the initial storage value (I) and the potential storage, Smax (Fmax + I), in the basins. The highest I is related to Mashin Basin (4.90 mm), and the lowest I is related to the Kardeh Basin (1.74 mm). Additionally, the value of I for the Kasilian Basin is equal to 2.55 mm. Although the total area of soil hydrology groups A and B (Table 3), which have more permeability, is higher in the Kasilian Basin than all the other basins (93.7%), the value of I for the Mashin Basin is higher than the I for the Kasilian Basin. This is while the Kasilian Basin does not have soil hydrology group D. Additionally, the vegetation cover of the Mashin Basin (72.8% percent) is less than that of the Kasilian (90.6%). The reason for this difference is that the value of I is mainly related to the surface absorption of soil and plants, and a part of this is used to fill the impoundments in the basin.
As a result, the deep infiltration factor does not play a big role, but the evaporation and drainage factor plays a greater role in maintaining and preventing flooding. The area of the Mashin Basin is 13 times the area of the Kasilian Basin; on the other hand, the average air temperature in the warm areas of the Mashin Basin, which plays a large role in the model results [30], is much higher than the temperature of the Kasilian Basin. Therefore, in the area of the Mashin Basin, the water flow needs more time to reach the outlet. It is noteworthy that the slope of the Mashin Basin (7.3%) is much less than the slope of the Kasilian Basin (15.8%). Therefore, there is much more time for the evaporation of surface moisture in this basin. On the other hand, in large basins, the number of surface impoundments (i.e., pools, wells, streams, etc.) is more than in small basins.
Under these conditions, it is necessary to have more rainfall in the Mashin Basin so that the flood will flow. On the contrary, in small basins with a high slope, floods occur with a small amount of precipitation. For a better understanding of this issue, we compared the Smax values of these two basins. As seen in Figure 6, the Smax value in the Mashin Basin is 50.6 mm and is about 1% less than the Smax value in the Kasilian Basin. The value for Smax in the studies of [24] is 59.2 mm for the Mashin Basin.
Therefore, although a part of Smax is related to I, deep infiltration (Fmax) has a greater contribution to maximum retention (Smax). For this reason, the permeability of soil (Table 3) and vegetation cover (Table 4), which is higher in the Kasilian Basin than the other basins, has a greater contribution to Fmax and, consequently, to the potential retention. These results clearly show that there is no significant relationship between I and Smax. This also challenges the hypothesis of SCS-CN [17] in terms of the coefficient λ (I/Smax). In other words, I is more influenced by slopes, air temperature, and impoundments, and Fmax is more influenced by soil hydrological groups, vegetation cover, and impoundments. If the impoundments are very small, they will not play a role in the secondary retention value (Fmax), but if the impoundments are large, they play an important role in all three retentions, I, Fmax, and finally, Smax (which is the sum of all retentions).
The lowest retention values (I = 1.74 mm), Fmax (14.97 mm), and Smax (16.72 mm) are related to the Kardeh Basin (Figure 6); its area (547 square kilometers) is relatively large, and the slope of the basin (7.67%) is low. Therefore, it is expected that the value of I is high in the Kardeh Basin. The reason for the low retention in the Kardeh Basin is related to the soil hydrological groups and land use: when 55.2% of the basin is rocky (Group D) and the sum of groups C and D is equal to 86.8%, the impoundments are also greatly reduced, and floods will flow even with little rainfall. In other words, the value of I becomes much less than the other basins. On the other hand, the deep infiltration is also greatly reduced, and, as a result, Fmax and Smax are also reduced. Figure 7 also shows the land use of the Kardeh Basin.
According to the calibration of model (10) and total retention (St), Table 5 is adjusted based on the runoff-rainfall model (11).
Figure 8 shows the fitting of the total retention model (St) to the rainfall data for the Imamah, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins. As can be seen, with the increase in the amount of precipitation (Pa), the total retention slope line ( d S t d ( p a ) ) decreases so that, at the beginning of runoff, the retention slope line from one (slope line of 45 degrees) gradually starts to decrease until Pa tends to infinity. In this case, the retention slope line also approaches zero (the retention value has reached its maximum and remains constant, Equation (10)).
The comparison of rainfall-runoff curves in Figure 9 shows that the slope of the runoff in all basins is zero at first (Pa = I) and then gradually increases until it finally reaches 1. For a rainfall equal to Pa, the lowest value of I, the highest amount of runoff, and the highest slope line of runoff depth ( d Q d ( p a ) ) are related to the Kardeh Basin. As was said before, the reason for the lack of retention in the Kardeh Basin is related to the soil hydrological groups and land use [26]. Additionally, the lowest slope line is related to the Kasilian Basin. In the Kasilian Basin, due to land use (77% forest cover), secondary retention is higher than in the other basins and causes a decrease in the runoff slope line ( d Q d ( p a ) ).
In Figure 10, the total retention slope line ( d S t d ( p a ) ) and runoff slope line ( d Q d ( p a ) ), with respect to rainfall changes, are shown. At the beginning of a rainfall event, there is no runoff (all rainfall is kept by infiltration, surface absorption, impoundments, and evaporation), and the slope line of the total retention is one ( d S t d ( p a ) = 1 ). In fact, the runoff slope is equal to 0 ( d Q d ( p a ) = 0 ), while the retention slope is equal to 1. Further, with an increase in rainfall, the slope of the runoff increases, while, correspondingly, the slope of the retention decreases so that when the slope of the runoff tends to 1, the slope of the retention becomes 0. In other words, the sum of the runoff slope and retention slope is always equal to 1 (Equation (7)). These results are consistent with the results of [26].
The model evaluation results are shown in Table 6. As can be seen, the coefficient of determination (R2) for predicting runoff was between 0.78 and 0.96 in the Navrood and Darjazin basins, respectively, and are within an acceptable range. The root mean square errors (RMSEs) are also between 0.86 and 2.28, and the Nash–Sutcliffe efficiency (NSE) varied between 0.79 and 0.96, which are within the significant range.

4. Conclusions and Results

  • In this study, it was shown that the amount of primary retention (I) is more influenced by the slope of the basin, air temperature, and impoundment. While the secondary retention (Fmax) is more influenced by the soil hydrological group, vegetation cover, and impoundment. If the impoundments are very small, they will not play a role in the secondary retention value, but if the impoundments are large and deep, they play an important role in all three retentions (I, Fmax, and Smax);
  • Unlike the models that have fixed coefficients, in this model, all the parameters, including Pa, I, Fmax, and Smax, are defined. For this reason, the model has high sensitivity;
  • In the present model, the method of estimating the parameters (I, Fmax, and Smax) is completely different from the SCS-CN method and can be easily calculated;
  • The results of fitting the model to the rainfall-runoff data based on R2, RMSE, and NS showed that the model has an acceptable ability to predict runoff and retention in all studied basins, especially if the selection of rainfall-runoff data was carried out carefully. In any case, if the rainfall-runoff data are consistent with the conceptual rainfall–runoff curve, we can expect very good results;
  • The new model is completely consistent with the revised SCS-CN conceptual curve, but its mathematical model is different from the SCS-CN mathematical model (the SCS-CN model does not have boundary conditions ( lim P a 0 ( Q ) 0 )).

Author Contributions

S.S., A.R.G., K.O.-A.-A. and S.I. designed the study, collected the data, wrote the manuscript, and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Some of or all the data, models, or codes generated or used during the study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mosavi, A.; Ozturk, P.; Chau, K.W. Flood Prediction Using Machine Learning Models: Literature Review. Water 2018, 10, 1536. [Google Scholar] [CrossRef]
  2. Williams, J.R.; Kannan, N.; Wang, X.; Santhi, C.; Arnold, J.G. Evolution of the SCS Runoff Curve Number Method and its Application to Continuous Runoff Simulation. J. Hydrol. Eng. 2012, 17, 1221–1229. [Google Scholar] [CrossRef]
  3. Buszney, M. Improving the Efficency of SCS Runoff Curve Number. J. Irrig. Drain. Eng. 1989, 115, 798–805. [Google Scholar]
  4. Huang, M.; Gallichand, J.; Dong, C.; Wang, Z.; Shao, M. Use of Soil Moisture Data and Curve Number Method for Estimating Runoff in the Loess Plateau of China. Hydrol. Process. Int. J. 2007, 21, 1471–1481. [Google Scholar] [CrossRef]
  5. Nandhakumar, S.; Arsheya, S.; Kirtika Sri, V.K. Estimation of Precipitation Runoff Using SCS and GISApproach in Puzhal Watershed. Int. J. Civ. Eng. Technol. 2019, 10, 1978–1998. [Google Scholar]
  6. Lal, M.; Mishra, S.K.; Kumar, M. Reverification of Antecedent Moisture Condition Dependent Runoff Curve Number Formulae using Experimental Data of Indian Watersheds. Catena 2019, 173, 48–58. [Google Scholar] [CrossRef]
  7. Song, W.; Jiao, J.; Du, P.; Liu, H. Optimizing the Soil Conservation Service Curve Number Model by Accounting for Rainfall Characteristics: A Case Study of Surface Water Sources in Beijing. Environ. Monit. Assess. 2021, 193, 115. [Google Scholar] [CrossRef]
  8. Sharma, I.; Mishra, S.; Pandey, A. Estimation of Antecedent Soil Moisture Using Soil Conservation Service Curve Number (SCS-CN) Method Utilizing the Experimental Data of a Small Indian Watershed. In Proceedings of the American Geophysical Union (AGU) Fall Meeting, San Francesco, CA, USA, 9–13 December 2019. [Google Scholar] [CrossRef]
  9. Sharma, I.; Mishra, S.K.; Pandey, A.; Kumre, S.K. A Modified NRCS-CN Method for Eliminating Abrupt Runoff Changes induced by the Categorical Antecedent Moisture Conditions. J. Hydro-Environ. Res. 2022, 44, 35–52. [Google Scholar] [CrossRef]
  10. Sharma, N.K.; Mishra, S.K.; Pandey, A.; Verma, R.K.; Verma, S. Improved SCS-CN Model Incorporating Storm Intensity for Runoff Estimation. Eng. Rural. Dev. 2022, 25, 526–536. [Google Scholar]
  11. Upreti, P.; Ojha, C.S.P. Development and Performance Evaluation of SCS-CN Based Hybrid Model. Water Sci. Technol. 2022, 85, 2479–2502. [Google Scholar] [CrossRef]
  12. Verma, S.; Mishra, S.K.; Verma, R.K. Improved Runoff Curve Numbers for a Large Number of Watersheds of the USA. Hydrol. Sci. J. 2020, 65, 2658–2668. [Google Scholar] [CrossRef]
  13. Ogarekpe, N.M.; Nnaji, C.C.; Antigha, R.E.E. A Preliminary Case for Modification of the SCS-CN Hydrologic Model for Runoff Prediction in Imo River Sub-basin. Arab. J. Geosci. 2022, 15, 786. [Google Scholar] [CrossRef]
  14. Singh, N.M.; Winkangshu, T.; Devi, T.T. Comparison of Simple and Modified SCS-CN in Runoff Prediction in a Highly Flood Prone Zone. In Sustainable Water Resources Management; Springer: Singapore, 2023; pp. 133–143. [Google Scholar]
  15. Hawkins, R.H. A comparison of predicted and observed runoff curve numbers. In Proceeding of ASCE, Irrigation and Drainage Division, Special Conference; ASCE: Flagstaff AZ, USA, 1984; pp. 702–709. [Google Scholar]
  16. Sahu, R.K.; Mishra, S.K.; Eldho, T.I. An Improved AMC-coupled Runoff Curve Number Model. Hydrol. Process. 2010, 24, 2834–2839. [Google Scholar] [CrossRef]
  17. Caletka, M.; Šulc Michalková, M.; Karásek, P.; Fučík, P. Improvement of SCS-CN initial abstraction coefficient in the Czech Republic: A study of five catchments. Water 2020, 12, 1964. [Google Scholar] [CrossRef]
  18. Vaezi, A.; Abbasi, M. Efficiency of the SCS-CN Method in Estimating Runoff in Taham Cahi Watershed, North West of Zanjan. JWSS-Isfahan Univ. Technol. 2012, 16, 209–219. [Google Scholar]
  19. Ebrahimian, M.; Nuruddin, A.A.; Soom, M.A.B.M.; Sood, A.M. Application of NRCS-curve Number Method for Runoff Estimation in a Mountainous Watershed. Casp. J. Environ. Sci. 2012, 10, 103–114. [Google Scholar]
  20. Shamohammadi, S.; Zomorodian, M. Comparison of Soil Conservation Service Model SCS and the Bennett Soil Moisture Accounting Model (SMA-B) in the Flood Estimation Zard River Basin. Iran. J. Watershed Manag. Sci. Eng. 2013, 7, 9–17. [Google Scholar]
  21. Sadeghi, S.H.R.; Mahdavi, M.; Razavi, S.L. Importance of Calibration of Maximum Storage Index Coefficient and Curve Number in SCS Model in Amameh, Kasilian, Darjazin and Khanmirza Watersheds. J. Water Manage. Sci. Eng. 2008, 2, 12–24. (In Persian) [Google Scholar]
  22. Asadi, E.; Fakheri, A.; Gorbani, M. Development of Conceptual Rainfall-Runoff Model for Quick and Slow Runoff Simulation (Case Study: Navrood Basin, Iran). Water Soil Sci. 2012, 22, 61–75. (In Persian) [Google Scholar]
  23. Avarand, R.; Torabi, H. Estimation of Runoff Rainfall and Preparation of Land Use Maps and Agricultural Areas Using Remote Sensing Technology in the Yellow River Catchment in Khuzestan Province. In Proceedings of the First International Conference and the Third National Conference on Dams and Hydropower Plants, Tehran, Iran, 8 February 2011; pp. 1–12. (In Persian). [Google Scholar]
  24. Shamohammadi, S. Presenting the New Adsorption Isotherm Model. In Proceedings of the Second International Conference on Environmental Hazarde, Tehran, Iran, 29–30 October 2013. (In Persian). (In Persian). [Google Scholar]
  25. Orak, N.; Farhadi, N. Investigating the Factors Affecting Runoff in Rudzard Basin in Khuzestan Province. In Proceedings of the First International Conference on Geographical Sciences, Abadeh, Iran, 29 October 2015; pp. 1–9. (In Persian). [Google Scholar]
  26. Shamohammadi, S. Presenting the Mathematical Model to Determine Retention in the Watersheds. Eur. Water 2017, 57, 207–213. [Google Scholar]
  27. Izadi, S.; Shamohammadi, S. Evaluation of Rainfall-Runoff-Retention Model (3RM) in Kassilian and Darjazin Watersheds. Water Irrig. Manag. 2022, 12, 309–325. (In Persian) [Google Scholar]
  28. Mockus, V. Estimation of Total Surface Runoff for Individual Storms. Exhibit A, Appendix B, Interim Survey Rep.; (Neosho) River Watershed USDA: Washington, DC, USA, 1949; pp. 1–79. [Google Scholar]
  29. Jahanbakhsh, S.; Rezaee, M.; Haghighi, E.; Rosta, I. The Relationship between Large-scale Circulation Patterns of Sea Level with Snowfall in Northwestern Iran. Geogr. J. 2015, 12, 19–35. (In Persian) [Google Scholar]
  30. Gomis-Cebolla, J.; Garcia-Arias, A.; Perpinyà-Vallès, M.; Francés, F. Evaluation of Sentinel-1, SMAP and SMOS surface soil moisture products for distributed eco-hydrological modelling in Mediterranean forest basins. J. Hydrol. 2022, 608, 127569. [Google Scholar] [CrossRef]
Figure 1. Schematic of the conceptual model [26].
Figure 1. Schematic of the conceptual model [26].
Civileng 04 00052 g001
Figure 2. Schematic of the mathematical Shamohammadi model for retention.
Figure 2. Schematic of the mathematical Shamohammadi model for retention.
Civileng 04 00052 g002
Figure 3. Location of Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basin in Iran.
Figure 3. Location of Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basin in Iran.
Civileng 04 00052 g003
Figure 4. (a) Isohyetal map, (b) hyetograph, and (c) flood hydrograph of the Mashin basin.
Figure 4. (a) Isohyetal map, (b) hyetograph, and (c) flood hydrograph of the Mashin basin.
Civileng 04 00052 g004
Figure 5. DEM of the Mashin Basin.
Figure 5. DEM of the Mashin Basin.
Civileng 04 00052 g005
Figure 6. Comparison of initial retention (I), secondary retention (Fmax) and potential retention (Smax) in the Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza and Mashin basins.
Figure 6. Comparison of initial retention (I), secondary retention (Fmax) and potential retention (Smax) in the Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza and Mashin basins.
Civileng 04 00052 g006
Figure 7. The land use of the Kardeh Basin.
Figure 7. The land use of the Kardeh Basin.
Civileng 04 00052 g007
Figure 8. Shamohammadi model output, estimating the actual retention of the Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins.
Figure 8. Shamohammadi model output, estimating the actual retention of the Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins.
Civileng 04 00052 g008
Figure 9. The rainfall–runoff relationship in the Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins.
Figure 9. The rainfall–runoff relationship in the Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins.
Civileng 04 00052 g009
Figure 10. Changes in total retention and runoff compared to changes in precipitation in the studied basin.
Figure 10. Changes in total retention and runoff compared to changes in precipitation in the studied basin.
Civileng 04 00052 g010
Table 1. General characteristics of Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins in Iran.
Table 1. General characteristics of Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins in Iran.
Basin
Basin FeaturesMashinKhanmirzaKardehDarjazinNavroodKasilianEmameh
Longitude49°39′ 50°55′ 41°35′ 53°12′ 48°34′ 53°8′ 51°35′
50°10′51°18′47°39′53°29′48°54′53°15′57°35′
Latitude31°21′ 31°32′31°37′35°59′ 37°36′ 35°58′51°32′
31°41′31°37′35°51′36°59′37°47′36°07′51°38′
Area (km2)8823915473312606737
Mean Elevation (m)210123120872152157316202620
Slope (%)7.345.87.6714.67.6715.807.5
Mean annual precipitation (mm)7578141338385376570.6480
Mean annual temperature (°C)9.39.211.59.811.512.125/7
Climate (De Martonne) AridArid and coldArid and semi-aridSemi-aridSemi-humidVery humidHumid
Table 2. Calibration and evaluation of the Shamohammadi model for the Emameh Basin.
Table 2. Calibration and evaluation of the Shamohammadi model for the Emameh Basin.
IERCalibrationEvaluation
P a o b s (mm)
( P + I E R )
Q o b s   ( m m ) S t   m m
( P a Q )
Q (mm)
( P a S m a x P a F m a x + P a )
P a o b s (mm) Q o b s   ( m m ) Q m m
( P a 32.72 P a 30.43 + P a )
5.5033.70 *15.1018.6016.5128.1012.6012.39
-27.6013.2314.3712.0425.90 *11.3010.86
-26.5010.9015.6011.2722.50 *9.108.59
-21.809.4012.408.1421.606.708.02
3.7018.80 *8.5010.306.3018.50 *6.506.133
- 18.006.4411.565.8418.706.206.25
2.2017.40 *5.1212.285.5018.706.106.25
- 17.105.1012.005.3316.90 *5.705.22
- 16.50 *4.7811.725.0016.805.205.16
3.2015.50 *3.9811.524.4615.00 *4.304.20
- 14.803.2011.604.0914.804.004.09
- 13.402.7610.643.4014.50 *3.403.94
- 13.902.7311.173.6413.06 *3.923.49
1.9012.90 *2.6010.303.1613.203.003.30
- 10.402.507.902.0710.602.202.16
- 7.400.906.501.007.501.201.03
- 5.500.804.700.497.100.450.91
- 3.200.302.900.090.000.000.00
* rains that have been corrected.
Table 3. Soil hydrological groups for each basin (%).
Table 3. Soil hydrological groups for each basin (%).
Hydrological GroupsMashinKhanmirzaKardehDarjazinNavroodKasilianEmameh
A2.6016.733.843.270.0024.600.00
B10.3024.259.2732.9632.0069.105.00
C47.4028.4831.6034.3949.006.3044.00
D39.7030.5355.2029.3817.000.0051.00
Table 4. The land use in the basins (%).
Table 4. The land use in the basins (%).
Land UseMashinKhanmirzaKardehDarjazinNavroodKasilianEmameh
Rocky lands16.990.0068.5347.870.008.20.00
Pasture35.8420.4415.1645.2333.307.3062.65
Forest9.210.000.000.0066.7077.200.00
Agricultural12.565.673.104.450.006.102.36
Residential and River6.585.412.532.450.001.200.00
Rainfed14.638.4710.670.000.000.000.00
Wasteland4.250.000.000.000.000.000.00
Table 5. Results of Shamohammadi model fitting regarding the rainfall-runoff data of the Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins.
Table 5. Results of Shamohammadi model fitting regarding the rainfall-runoff data of the Emameh, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins.
Basin λ I (mm)Model (13)
Emameh0.072.3 Q = P a ( 32.7 P a 30.4 + P a )
Kasilian0.052.5 Q = P a ( 51.1 P a 48.6 + P a )
Navrood0.082.4 Q = P a ( 29.1 P a 26.7 + P a )
Darjazin0.103.2 Q = P a ( 24.7 P a 21.5 + P a )
Kardeh0.101.7 Q = P a ( 16.7 P a 14.97 + P a )
Khanmirza0.072.5 Q = P a ( 38.1 P a 33.2 + P a )
Mashin0.104.9 Q = P a ( 50.6 P a 44.5 + P a )
Table 6. Evaluation results of the modified 3RM model.
Table 6. Evaluation results of the modified 3RM model.
BasinNSRMSE (mm)R2
Emameh0.960.860.96
Kasilian0.940.930.86
Navrood0.781.10.78
Darjazin0.790.890.79
Kardeh0.960.890.87
Khanmirza0.902.280.84
Mashin0.911.420.88
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shamohammadi, S.; Ghasemi, A.R.; Ostad-Ali-Askari, K.; Izadi, S. Presentation of a Rainfall–Runoff Retention Model (3RM) Based on Antecedent Effective Retention for Estimating Runoff in Seven Basins in Iran. CivilEng 2023, 4, 966-981. https://doi.org/10.3390/civileng4030052

AMA Style

Shamohammadi S, Ghasemi AR, Ostad-Ali-Askari K, Izadi S. Presentation of a Rainfall–Runoff Retention Model (3RM) Based on Antecedent Effective Retention for Estimating Runoff in Seven Basins in Iran. CivilEng. 2023; 4(3):966-981. https://doi.org/10.3390/civileng4030052

Chicago/Turabian Style

Shamohammadi, Shayan, Ahmad Reza Ghasemi, Kaveh Ostad-Ali-Askari, and Saeedeh Izadi. 2023. "Presentation of a Rainfall–Runoff Retention Model (3RM) Based on Antecedent Effective Retention for Estimating Runoff in Seven Basins in Iran" CivilEng 4, no. 3: 966-981. https://doi.org/10.3390/civileng4030052

APA Style

Shamohammadi, S., Ghasemi, A. R., Ostad-Ali-Askari, K., & Izadi, S. (2023). Presentation of a Rainfall–Runoff Retention Model (3RM) Based on Antecedent Effective Retention for Estimating Runoff in Seven Basins in Iran. CivilEng, 4(3), 966-981. https://doi.org/10.3390/civileng4030052

Article Metrics

Back to TopTop