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

: This study aims to correct and assess the SCS-CN model. In this research, the 3RM model (written by Shamohammadi) has been modiﬁed in such a way that the maximum primary retention ( I ), maximum secondary retention ( F max ), and basin potential retention ( S max ) can be calculated using precipitation ( P a ). The purpose of this study is to evaluate the total retention model ( S t = f ( F max , S max , 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 ﬁtting the model to the rainfall-runoff data showed that the evaluation indices, including the coefﬁcient of determination ( R 2 ), 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.


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 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 (I a ) and the antecedent effective retention (I ER ) (I = I a + I ER ).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.
S t : Total retention: this is equal to F + I, so S t is also variable.S max : The maximum retention or potential retention (S p ) 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 S max , 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 F max .
F max : It is defined as S max and is linearly constant.Its relationship with S max is as F max = S max − I.
I ER : Antecedent effective retention: this is a part of retention due to antecedent rainfall (P A ), which is effective in the rainfall-runoff system and is calculated as per Equation (1) [27]: where P A , Q A , and (ET 0 ) A are the antecedent rainfall depth, antecedent runoff depth, and the daily antecedent potential evapotranspiration, respectively, and n 1 and n t are the days of the beginning of the antecedent and the target, respectively.Additionally, it can always be written as per Equation (2): where Q and P a are the runoff depth (mm) and the modified target rainfall depth (mm), respectively.

Assumptions of the Shamohammadi Model
Equations ( 3) and ( 4) have been presented for the mathematical description of the conceptual curve in Figure 1.
dQ dp a = S t S max (3) CivilEng 2023, 4, FOR PEER REVIEW 3 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  =  − .
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]: 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): where Q and Pa are the runoff depth (mm) and the modified target rainfall depth (mm), respectively.

Assumptions of the Shamohammadi Model
Equations ( 3) and ( 4) have been presented for the mathematical description of the conceptual curve in Figure 1.In Equation (4), the K sh is the Shamohammadi constant.Additionally, dQ dp a and dS t dp 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 Equation ( 5) is a general equation.Obviously, when the runoff is equal to 0, P a = S t = I, and when the total retention is completed and when rainfall is converted to runoff, then P a = Q.
Using differentiation on both sides of Equation ( 5) with respect to rainfall, we obtain the following: Finally, Equation ( 6) can be written as Equation ( 7): Considering Equations ( 3), (4), and ( 7), we can write the following as Equation ( 8): and then (Equation ( 9)): According to Equation ( 9), when the rainfall tends to I (P a →I), the value of K sh will be equal to F max because when the rainfall is equal to I, the value of S t 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 F max , the value of k sh is equal to F max .Therefore, Equation ( 9) is modified as per Equation (10). Figure 2 shows the schematic of the maintenance of Equation (10).According to the relationship between runoff and total retention (Equation (2)), the amount of runoff depth is equal to what follows in Equation (11): 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 (S max − F max ).In this case, in addition to the error caused by linear regression being removed, the value of F max 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 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 (P a →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): 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: 1.
The basin area should be so small that the whole basin is covered by rainfall (uniformly and at a constant intensity); 2.
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; 3.
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; 4.
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.

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.
because it is based on historical data and operates as a black box model.In in this model, in addition to the man-made structures, the changes in the pa one rainfall event to another also cause errors in the results.

Studied Areas
In this study, the basins of Emameh, Kasilian, Navrood, Darjazin, K mirza, and Mashin, which were equipped with rain gauges and hydrometer studied.The selection of these basins was based on the quality of rainfall an the distribution of rainfall, climate diversity, vegetation diversity, area size, conditions.Figure 3 shows the location of the studied basin in Iran.The g teristics 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 (ET 0 = (∑ n t n1 ET 0 ) A 315 ≈ 0.5 mm/day).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).
where D is the number of days of evapotranspiration (n t − n 1 ).The maximum number of days for calculating evapotranspiration before the target precipitation was considered to be 5 days [2,28].

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); 2.
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 (P A ) were also determined.Because the antecedent rains did not have runoff, only the antecedent rain was recorded; 3.
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; 4.
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; 5.
where Q obs is the amount of observed runoff, Q sim is the amount of simulated runoff (model), and N is the amount of data.

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.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 km 2 ), 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].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 I ER = 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 km 2 ), 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).
rainfall happened in March.It was expected that the melting of snow would in the volume of runoff, but the investigations showed that the height of the in April in the Mashin Basin was above 2700 m [24] as a result, the effect of s on the volume of the flood was very low (see Figure 5).These results are con [23,25,29].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/S max ) 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 S max [17,30].Therefore, it is necessary to examine the values of I and S max separately.
Figure 6 shows the initial storage value (I) and the potential storage, S max (F max + 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  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 S max values of these two basins.As seen in Figure 6, the S max value in the Mashin Basin is 50.6 mm and is about 1% less than the S max value in the Kasilian Basin.The value for S max in the studies of [24] is 59.2 mm for the Mashin Basin.
Therefore, although a part of S max is related to I, deep infiltration (F max ) has a greater contribution to maximum retention (S max ).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 F max and, consequently, to the potential retention.These results clearly show that there is no significant relationship between I and S max .This also challenges the hypothesis of SCS-CN [17] in terms of the coefficient λ (I/S max ).In other words, I is more influenced by slopes, air temperature, and impoundments, and F max 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 (F max ), but if the impoundments are large, they play an important role in all three retentions, I, F max , and finally, S max (which is the sum of all retentions).The lowest retention values (I = 1.74 mm), F max (14.97 mm), and S max (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, F max and S max 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).
the sum of groups C and D is equal to 86.8%, the impoundments are a and floods will flow even with little rainfall.In other words, the value less than the other basins.On the other hand, the deep infiltration is al and, as a result, Fmax and Smax are also reduced.Figure 7 also shows Kardeh Basin.According to the calibration of model (10) and total retention (St) based on the runoff-rainfall model (11).Figure 8 shows the fitting of the total retention model (S t ) 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 (P a ), the total retention slope line ( dS 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 P a 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)).
Imamah, Kasilian, Navrood, Darjazin, Kardeh, Khanmirza, and Mashin basins.As c seen, with the increase in the amount of precipitation (Pa), the total retention slop ( ( ) ) decreases so that, at the beginning of runoff, the retention slope line from one ( line of 45 degrees) gradually starts to decrease until Pa tends to infinity.In this cas retention slope line also approaches zero (the retention value has reached its maxi and remains constant, Equation ( 10)).The comparison of rainfall-runoff curves in Figure 9 shows that the slope of th off in all basins is zero at first (Pa = I) and then gradually increases until it finally re 1.For a rainfall equal to Pa, the lowest value of I, the highest amount of runoff, an highest slope line of runoff depth ( ( ) ) are related to the Kardeh Basin.As was sa fore, the reason for the lack of retention in the Kardeh Basin is related to the soil h logical groups and land use [26].Additionally, the lowest slope line is related to th silian Basin.In the Kasilian Basin, due to land use (77% forest cover), secondary rete is higher than in the other basins and causes a decrease in the runoff slope line ( ( The comparison of rainfall-runoff curves in Figure 9 shows that the slope of the runoff in all basins is zero at first (P a = I) and then gradually increases until it finally reaches 1.For a rainfall equal to P a , the lowest value of I, the highest amount of runoff, and the highest slope line of runoff depth ( dQ 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 ( dQ d(p a ) ).
CivilEng 2023, 4, FOR PEER REVIEW  In Figure 10, the total retention slope line ( dS t d(p a ) ) and runoff slope line ( dQ 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 ( dS t d(p a ) = 1).In fact, the runoff slope is equal to 0 ( dQ 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].In Figure 10, the total retention slope line ( ( ) ) and runoff slope line ( ( respect to rainfall changes, are shown.At the beginning of a rainfall event, the runoff (all rainfall is kept by infiltration, surface absorption, impoundments, and ration), and the slope line of the total retention is one ( ( ) = 1).In fact, the runo is equal to 0 ( ( ) = 0), while the retention slope is equal to 1. Further, with an inc rainfall, the slope of the runoff increases, while, correspondingly, the slope of the r decreases so that when the slope of the runoff tends to 1, the slope of the retention b 0. In other words, the sum of the runoff slope and retention slope is always eq (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 coeff determination (R 2 ) for predicting runoff was between 0.78 and 0.96 in the Navr Darjazin basins, respectively, and are within an acceptable range.The root mean Pa(mm) The model evaluation results are shown in Table 6.As can be seen, the coefficient of determination (R 2 ) for predicting runoff was between 0.78 and 0.96 in the 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.

1.
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 (F max ) 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, F max , and S max ); 2.
Unlike the models that have fixed coefficients, in this model, all the parameters, including P a , I, F max , and S max , are defined.For this reason, the model has high sensitivity; 3.
In the present model, the method of estimating the parameters (I, F max , and S max ) is completely different from the SCS-CN method and can be easily calculated; 4.
The results of fitting the model to the rainfall-runoff data based on R 2 , 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; 5.
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

Figure 2 .
Figure 2. Schematic of the mathematical Shamohammadi model for retention.

Figure 2 .
Figure 2. Schematic of the mathematical Shamohammadi model for retention.
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 I ER values occurred in the Kasilian (8.4 mm) and Darjazin (0.5 mm) basins, respectively.The ET 0 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).

Figure 7 .
Figure 7.The land use of the Kardeh Basin.

Figure 7 .
Figure 7.The land use of the Kardeh Basin.

Figure 10 .
Figure 10.Changes in total retention and runoff compared to changes in precipitation in th basin.

Figure 10 .
Figure 10.Changes in total retention and runoff compared to changes in precipitation in the studied basin.

Table 2 .
Calibration and evaluation of the Shamohammadi model for the Emameh Basin.

ER Calibration Evaluation P a obs (mm) (P+I ER ) Q obs (mm) S
t (mm) (P a −Q) Q (mm) (P a −S

max P a F max +P a ) P a obs (mm) Q obs (mm)
Q (mm) (P a −32.72

P a 30.43+P a )
* rains that have been corrected.

Table 4 .
The land use in the basins (%).

Table 6 .
Evaluation results of the modified 3RM model.