Development of a New Event-Based Rainfall-Runoff Equation Based on Average Rainfall Intensity During an Event

Event-based rainfall-runoff models are practical tools commonly used to predict catchments’ response to a rainfall event. However, one of the main concerns is that the characteristics of rain events are neglected in the model development. This paper develops a novel event-based rainfall-runoff equation to incorporate rainfall characteristics into account. The performance of the new equation is evaluated based on the root mean square error, Nash–Sutcliffe efficiency coefficient, and per cent bias for 13,339 rainfall-runoff events between 2005 and 2020 over 23 catchments across New Zealand and Australia with oceanic, mediterranean, tropical, subtropical, and semiarid climates. Compared to the previous event-based models, the new equation shows an improvement in runoff estimation in almost all case studies. Furthermore, considering the new equation is simple, efficient, and takes the rain event duration into account, the new equation has the potential to become a robust alternative method to the conventional curve number method in hydrological engineering projects.


Introduction
The hydrological response of a catchment to a rainfall event is commonly known as a rainfall-runoff process [1]. A rainfall-runoff process often consists of complex physical procedures, i.e. interception, depression storage, evapotranspiration, infiltration, subsurface flow, groundwater flow, overland flow, and channel flow [2]. Previous studies indicate that rainfall characteristics such as rainfall rate, duration, and distribution, besides the catchment characteristics such as size, soil type, land use, and slope, can impact the occurrence and volume of runoff [3][4][5][6]. In addition, considering the physical condition of catchments is not homogenous, each catchment may respond differently to a rainfall event [1].
Several physically based mathematical expressions and computer simulation models have been developed that fully or partially represent the hydrological processes. Sometimes, infiltration equations are implemented to predict runoff from a rainfall event. The most common infiltration equations cited in the literature for rainfall-runoff modelling are Green-Ampt [7], Philip [8], Horton [8], Holtan [9], and Kostiakov [10] methods. In addition, sometimes, distributed surface water models are combined with distributed groundwater models to provide a better estimation [11][12][13][14][15][16][17][18][19]. However, these models are highly complex, data-hungry, and often require several spatial and temporal parameters, which in most cases are not readily available [20].
A simple alternative to physically based and continuous models is to estimate runoff from an individual rain event, often known as event-based models [21]. For example, the runoff coefficient has been used to estimate runoff from a rainfall event for many decades, which assumes that the depth of direct runoff is a percentage of the rainfall depth: where C is the runoff coefficient, P is the total runoff depth, and Q is the total runoff depth in each rain event [22][23][24][25], or sometimes Q is the total direct runoff depth for each event [26][27][28][29], where "depth" refers to the total volume of runoff or direct runoff for each event over the catchment area. The terminology of "runoff coefficient" can be confusing as it

3
is not consistent throughout the literature [30]. For example, the parameter called response factor [26], hydrologic response [27], runoff ratio [23], water yield [31], conversion efficiency [24], direct runoff response ratio [5], and runoff coefficient [6,22,28]. The other definition of runoff coefficients also can come from the rational method, which claims the peak flow is proportional to rainfall intensity for a given catchment [32]. For clarification, in the current study, the runoff coefficient means the total direct runoff over total precipitation for each rain event.
Even though the concept of estimating runoff for each event from runoff is attractive, the runoff coefficient is highly variable as the runoff coefficient is proportional to catchment and rainfall characteristics [3,5].
The SCS-CN method originally comes from a lumpedbased approach that calculates the total direct runoff from a storm event [58]. For deriving the SCS-CN expression, the proportionality between retention and runoff is assumed to be: where F is the actual retention, S is the maximum potential retention, and Q is the total direct runoff during a rainfall event. Assuming: and subtracting a certain amount of rainfall, referred to as initial abstraction from the total runoff depth P and solving Eq. 2 for Q yields: where I a is the initial abstraction. Commonly, a linear correlation is assumed between S and I a , where λ is the initial abstraction ratio. λ is considered 0.2 in the original formula [40,59,60]. However, the background of this assumption is known [58].
In the SI units, where CN value is defined by hydrological soil group, land use, hydrological surface condition, and soil moisture condition [41]. Theoretically, the CN may vary between 0 and 100. From a hydrologic engineering perspective, the SCS-CN method is simple, transparent, appealing, and only requires one empirical parameter, the CN. The SCS is popular because SCS-CN databases can easily link to distributed soil and vegetation layers stored within a GIS [34]. Hawkins [41] reported that there is no alternative with the benefits of the curve number method available, and within its group, the SCS-CN method is monotypic. However, besides many advantages, too many ambiguities raise questions about this method's accuracy [33,34,41,42]. For example, Eq. 2, as the central assumption of the SCS-CN method, has no physical justification [34].
Another concern is that the rainfall characteristics, such as rain duration, are neglected in the equation [61]. For example, the SCS-CN method does not count for the differences between 40 cm rainfall in 1 day and 10 days, while infiltration and runoff would be considerably different [32,58]. However, several studies statistically indicated that the runoff is proportional to the rainfall duration [3][4][5].
This study aims to introduce a novel, robust, and simple equation for estimating runoff from a rain event by taking the rain event duration into consideration. The new equation only has one empirical parameter. However, compared to the other eventbased models, it more precisely represents the physics of the rainfall-runoff process by incorporating rainfall characteristics.

Hydro-Meteorological Data
For this study, 23 catchments ranging from 0.3 to 564 km 2 are selected from catchments in New Zealand and Australia. The sites in New Zealand have an oceanic climate, with cool but not cold winters and warm summers. In contrast, the case studies in Australia have a range of climates, including Oceanic, Mediterranean, tropical, subtropical, and semi-arid climates. Commonly, no snow is recorded in any sites during the study periods.
Rain events are defined based on the minimum-intereventtime method [62], in which the beginning and end of each event are defined by a minimum rainless period of at least 24 h. The available data set is divided into the calibration and validation periods. In total, 13,339 events are observed across 23 study areas. In all cases, rainfalls are recorded either at the outlet or in the vicinity of the outlet. A summary of catchment size, climate type, calibration, and validation periods and the number of events are listed in Table 1.

Baseflow Separation
The quick flow (direct runoff), which is the direct response to a catchment to rainfall, is separated from the base flow, water seeping into the stream from groundwater. Therefore for the development of a rainfall-runoff model often, the quick flow is separated from the runoff. The BFI + , which is a powerful and user-friendly software with eleven methods for baseflow and quick flow separation methods, is selected [63,64]. The one-parameter recession digital filter method in the BFI + software is used for baseflow quick-flow separation for the current study. As an example of baseflow separation, the base flow separation for the Opanuku river is shown in Fig. 1. Then, the direct runoff depth for each event is calculated by dividing the total direct runoff for each event over the catchment area.

Development of the New Equation
The event average rainfall rate is calculated for each event: where i av is the average rainfall rate (mm/h) in each event, and t is the event duration (h). The existing correlation between runoff coefficient (Q/P) and average rainfall rate [1] has been suggested in previous studies. The correlation between the runoff coefficient and average rainfall rate in Opanuku and Pukekohe Ngakoroa rivers is illustrated in Fig. 2.
Considering the duration varies in each event, the proportionality between runoff and average rainfall rates (I av ) is not substantial.
Intensity-duration-frequency (IDF) curves were implemented to enhance the correlation by unifying rain event durations. The IDF curves concept was established in early 1930, and many mathematical expressions for IDF curves were suggested subsequently [32].
IDF curve can be expressed in a general form: where t is the rain storm duration, T is the storm event return period, f (T) is a function of the return period, and g(t) is a function of the storm duration t. f (T) and g(t) are unique functions for each site and can be estimated directly from previous rainfall records. For this study, following a method suggested by [32,65], the annual maximum rainfall for each site was rank-ordered, and the return period for each annual maximum rainfall value was estimated. Then, the rank-ordered rainfall value was divided by the event duration for each return period, resulting in the average intensity plotted versus the duration. As an example, the IDF curves of Pukekohe for return periods ranging from 0.01 to 100 years are given in Fig. 3.
For storms with a return period of T, Eq. 8 gives the following relation between the equivalent rainfall rate in 24 h and the rainfall rate in each rainstorm with a duration of t [32]: where i 24 is the equivalent rainfall rate in 24 (h).
Rearranging for i 24 yields: Assuming a power function form for g(t) after Koutsoyiannis et al. [66]: where a and n are local constant variables obtained from the past rainfall records for each site.
Therefore, Substituting g(t) and g (24) in Eq. 10 yields: Then, Eq. 13 is used to calculate i 24 for all events. With the new adjustment on rainfall rate, the runoff coefficient shows a stronger correlation to i 24 than i av in all case studies.
For all sites, the n value is calculated using the IDF curves. For example, from the rainfall record of the Opanuku and Pukekohe Ngakoroa sites, the IDF curves are generated, and then the n value is calculated as 0.37 and 0.56, respectively. In the next step, i 24 is calculated, and the relationships between the runoff coefficient and i 24 are plotted in Fig. 4.
(11) g(t) = at n ,  Fig. 5 compare the correlation (R 2 ) of Q/P vs i av and Q/P vs i 24 for all case studies. Considering all points in the scatterplot located above the 1:1 line indicates that using the IDF curves concept to unify the event durations effectively improved the runoff coefficient and rainfall rate significantly. A comparison between medians in the boxplot shows that the median of R 2 is increased by around 100% from R 2 = 0.22 in Q 0 /P − i av to R 2 = 0.44 in Q 0 /P − i 24 .
Then, it is assumed that the runoff coefficient (Q/P) is directly proportional to i 24 : Then, by multiplying the right side of the proportionality by a constant, α, the proportionality is replaced by an equation: where α (h/mm) is a local constant and varies for each site.
Substituting i 24 from Eq. 13 in Eq. 15 yields: From this point, Eq. 17 is known as the "new equation" in this paper. Considering quick flow data is used to develop the new equation, the concept of initial abstraction is already included in the P/Q correlation.

Results and Discussion
For all case studies, α and CN values were calculated using the least squares fitting technique to minimise the sum of the squared differences between the observed data where Q 0 is the observed runoff depth, and Q C is the calculated runoff depth from the new equation and SCS-CN method. The optimised CN and α values from the calibration process are listed in Table 2. Also, as some examples, the observed and estimated runoff by the new equation and SCS-CN method versus the total rainfall for each event for Opanuku and Pukekohe Ngakoroasites are illustrated in Fig. 6.
The SCS-CN method predicts only a single value for each rainfall event regardless of the rainfall characterises. However, the new equation estimates various runoffs for each rainfall event depending on the rainfall durations. For example, looking at three extreme events around P = 250 mm in Pukekohe Ngakoroa (Fig. 5), the SCS-CN method predicts 71 mm for the three events. In contrast, the new equation predicts 42, 51, and 82 mm, respectively. Similarly, three different runoffs, 29, 36 and 83 mm, are also observed at the outlet.
As a way of comparison, the estimated runoff by the new equation and SCS-CN method is plotted against the observed runoff for the Opanuku river in Fig. 7.
The R 2 in the SCS-CN plot is around 23% lower than the R 2 in the new equation (from R 2 = 0.79 to 0.64). Also, the slope of the linear trendline in the new equation is closer to the 1:1 line than the SCS-CN. So, overall, the new equation shows a better performance in the runoff depth estimation than the SCS-CN method in the Opanuku river for the study period.

The Performance Evaluation Strategies
Three statistical criteria were used to evaluate the new equation and SCS-CN performances.

Root Mean Square Error (RMSE)
where m is the total number of events. RMSE = 0 indicates the perfect model fit, and the closer the RMSE values are to zero, the better the model performance.

Nash-Sutcliffe Model Efficiency (NSE)
The Nash-Sutcliffe model efficiency coefficient [67] that is often used for evaluation of the performance of hydrologic models [68,69] is calculated for both methods: where Q O is the mean of observed flow depth, and NSE is the Nash-Sutcliffe model efficiency coefficient. The NSE values range from − ∞ to 1, where 1 shows a perfect prediction, whereas NSE ≤ 0 suggests that a simple average of observed values is better than the model estimations. Therefore, the model performance is unacceptable.

Percent Bias (PBIAS)
PBIAS measures the average tendency of the predicted flow depth to be larger or smaller than the observed data [70].
A low-magnitude value indicates an accurate model simulation where 0 is a perfect fit. Positive values suggest model underestimates, and negative values indicate model overestimates.

Calibration
RMSE, NSE, and PBIAS were calculated for the calibration periods for all case studies. A summary of the calibrated values (CN and α) and performance evaluation of the new equation and SCS-CN models in the calibration process are listed in Table 2.
Overall, the new equation shows a better performance compared to the SCS-CN model. For example, the average NSE of the new equation for 23 case studies is 0.74, which is 25% higher than the NSE of the SCS-CN. Likewise, the average RMSE improved to 14%. Scatterplots in Fig. 8 Figure 9 shows a scatterplot that compares the α and CN for 23 case studies.
The scatterplot shows a relatively weak to moderate correlation between α value and CN with an R 2 of 0.42. It means the higher the α value in a catchment, the less  Pukekohe Ngakoroa Opanuku rainfall infiltrates the groundwater system and generates more runoff.

Validation
For validation, the calculated values for CN and α from the calibration process were directly used in the SCS-CN and the new equation on the validation data sets, defined in Table 1. A summary of NSE, RMSE, and PBIAS for the validation process is listed in Table 3.
A comparison between the calibration and validation performances shows that in both models, all performance indicators are slightly dropped in the calibration process. There are some cases where both models perform exceptionally well, for example, in Rutherford Creek and Elizabeth Valley case studies. But there are some cases in which both models failed to predict runoff depth, for instance, Wakefield and West Hoe Heights.    Figure 10 shows no correlation between the new equation performance (NSE) and the α value. Considering that the α value represents the capacity of catchments to generate runoff, the lack of a correlation means the performance of the new equation is independent of the α value and can be used in a wide range of catchments.
In general, it can be concluded that the new equation is a robust alternative to the conventional models in oceanic climate conditions. This research also studied a few catchments in Mediterranean, tropical, subtropical, and semiarid classmates. However, no advantage or disadvantage is observed in the performance of the new equation in different climates. However, the numbers of samples in each climate are relatively small. Therefore, further investigation with a broader range of climates is required to conclude whether the new equation performs reasonably in all climate conditions.

Discussion on the Core Assumption of the New Equation and the SCS-CN Method
The proportionality of Q/S vs F/S (Eq. 2) and Q/S vs i 24 (Eq. 14) are the central assumptions for SCS-CN and the new equation, respectively. In order to investigate the validity of these assumptions, the R 2 of Q 0 /P vs F/S, and Q 0 /P vs i 24 for all case studies are calculated. Equations 3, 6, 8, and 13 are used to calculate F, S, i av , and i 24 , respectively. A scatterplot and a boxplot in Fig. 11 compare the R 2 of Q 0 /P vs F/S and i 24 for all case studies.
From the boxplot, the median R 2 of Q 0 /P vs i 24 is 0.44, which is 101% higher than the R 2 of Q 0 /P vs F/S. Also, a one-to-one compression in the scatterplot shows that all points are above the 1:1 line in the advantage of the new equation. Therefore, it can be concluded that the central assumption of the new equation (Eq. 14) is more reliable than the central assumption of the SCS-CN method (Eq. 2).

Conclusion
A new empirical event-based rainfall-runoff equation was developed. The IDF curve concept was implemented to introduce a 24-h event average rainfall (i 24 ) for each rain event. An acceptable correlation between the runoff coefficient (Q/P) and i 24 was observed. Based on this correlation, a new equation was developed. The new equation and SCS-CN equation were used to estimate runoff for 23 catchments across New Zealand and Australia with oceanic and semi-arid climates ranging from 0.3 to 584 km 2 . Statistical indicators such as NSE, RMSE, and PBIAS were used to compare the predictive skill of the new equation and SCS-CN. It was shown that the new equations perform better in both calibration and validation processes. Therefore, it can be concluded that the new equation is a robust alternative to the conventional event-based model in oceanic and semi-arid climate conditions. Also, there is a possibility that the new equation will be valid in other climate conditions. However, it needs further investigation.  Author Contribution Not applicable.
Funding Open Access funding enabled and organized by CAUL and its Member Institutions.
Data Availability All data that support the findings of this study are available from public websites.

Declarations
Ethical Approval Not applicable.

Conflict of Interest
The author declares no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.