In ﬂ uences of the Runo ﬀ Partition Method on the Flexible Hybrid Runo ﬀ Generation Model for Flood Prediction

: The partition of surface runo ﬀ and in ﬁ ltration is crucial in hydrologic modeling. To improve the ﬂ ood prediction, we designed four strategies to explore the in ﬂ uences of the runo ﬀ partition method on the ﬂ exible hybrid runo ﬀ generation model. The runo ﬀ partition strategies consist of a hydrological model without the runo ﬀ partition module, a two-source runo ﬀ partition method, an improved two-source runo ﬀ partition method considering the heterogeneity of the subsurface topography and land cover, and a three-source runo ﬀ partition method. The Xin’anjiang hydrological model was used as the modeling framework to simulate a six-hourly stream ﬂ ow for the Xun River watershed in Shaanxi Province, China. And the saturation-excess runo ﬀ generation and in ﬁ ltration-excess runo ﬀ generation mechanisms were combined to construct the ﬂ exible hybrid runo ﬀ generation model. The performances of the four strategies were compared and analyzed based on the continuous ﬂ ow discharge as well as the ﬂ ood events. The runo ﬀ components analysis method was used to test the model’s conformity with the reality of the watershed. The results showed that the three-source runo ﬀ partition method was not applicable to the ﬂ exible hybrid runo ﬀ generation model because it overestimated the surface runo ﬀ and almost ignored the subsurface storm ﬂ ow runo ﬀ . The improved two-source runo ﬀ partition method outperformed the others as it considered the heterogeneity of the watershed.


Introduction
A hybrid runoff generation process pattern comprised of multiple mechanisms can often happen in semi-arid, semi-humid, and mountainous watersheds due to the heterogeneity of meteorological factors and underlying surface conditions (i.e., rainfall, land covers, soil types, etc.) [1][2][3][4]. Runoff generated by the integration of the three components, including the subsurface stormflow runoff, saturation-excess runoff, and infiltration-excess runoff generation, is known as hybrid runoff, which leads to rapid flood occurrences and high flood peak discharges, and thus makes hydrological forecast even more challenging [5,6].
Hydrological modeling using conceptualized hybrid runoff generation mechanisms attracted lots of attention to solve the limitations of single runoff generation mechanisms [7]. Lots of conceptual mixed runoff generation models have been developed, for instance, the vertically mixed runoff generation model [8], the XAJ-Green-Ampt model [7], as well as the variable infiltration capacity (VIC) runoff generation model [9], and so on. These hybrid runoff generation models were constructed in accordance with the vertical combination of saturation-excess and infiltration-excess modules or in accordance with the spatial combination of saturation-excess and infiltration-excess modules [2]. The current mixed runoff generation models can acquire a good performance in the semihumid and semi-arid regions but also have some shortcomings. For instance, they assumed that a fixed runoff generation mechanism dominated each natural watershed. This assumption could not reflect the heterogeneity of the dominant runoff generation mechanisms within a basin [5,10,11]. To this end, the study of hydrological modeling has recently progressed toward developing and applying flexible models, which provided a new method for solving flood forecasting issues in semi-arid, semi-humid, and mountainous areas [5].
The flexible hydrological modeling framework allowed for the selection of alternative representations of runoff processes and various combinations of both linear and nonlinear components [12]. This kind of model can be regulated timely to construct a suitable model structure for a particular basin [5]. A lot of research has been carried out recently on developing flexible framework models. For example, Liu et al. [5] developed a novel flexible hybrid runoff generation modeling framework appropriate for hydrological modeling in semi-arid and semi-humid areas, which is known as the spatial combination computing models for runoff generation (SCCMs). Huang et al. [12] discussed the performances of four traditional hydrological models and compared them with those of four flexible models for semi-arid environments. Yi et al. [2] proposed an improved flexible hybrid runoff generation strategy and compared its accuracy with that of four traditional flexible strategies. The fundamentals of the flexible hybrid runoff generation models discussed above were the combinations of three single runoff generation mechanisms, sometimes one or two.
In the land surface hydrological processes, different water balance components are closely related and interact with each other. An inappropriate separation of runoff components directly or indirectly influences the simulations of other water balance components [13]. Runoff partition methods can also be significant for flexible runoff generation modeling [14]. While runoff generation dominates the simulated discharge volumes, runoff separation dominates the simulated hydrograph shapes. It is generally accepted that the hydrograph is comprised of different components with various response times, such as a fast or a slow runoff component [15][16][17][18]. The fast runoff may further be separated into surface flow and subsurface stormflow, and the slow runoff represents subsurface flow [19,20]. In addition, the same runoff component may originate from different runoff generation mechanisms [2,21]. For instance, surface runoff may originate from saturation excess or infiltration excess [22]. This issue can also be studied by using the runoff partition methods.
Several runoff partition methods have been commonly adopted in previous studies. The XAJ model proposed a runoff partition method according to the steady infiltration rate, known as the two-source XAJ model [23]. However, the infiltration can be spatially unevenly distributed due to the heterogeneity of the subsurface topography and land covers, etc. [14]. In addition, the separated surface runoff does not consider the subsurface stormflow. Kirkby [24] supposed that the third runoff component, subsurface stormflow, ought to be regarded as infiltration rate, which varies via the soil layer. The upper soil layer often has better permeability than the lower soil layer leading to a relatively impermeable layer in between, which is known as the interface. Infiltrating water from the upper soil layer will flow laterally, driven by topography-i.e., subsurface stormflow, which is now considered in the three-source XAJ model [25][26][27]. The surface runoff, subsurface stormflow, and subsurface runoff are separated using the free water storage reservoir in the three-source XAJ model [28]. Kling and Nachtnebel [15] presented a simple method for the regional estimation of runoff separation parameters of a spatially distributed monthly water balance model. Qi et al. [14] compared four different surface runoff and infiltration partition methods based on a Richards-equation-based SWAT model (RSWAT) to understand their impacts on watershed modeling. Pelletier and Andréassian [29] proposed a novel hydrograph separation method that is based only on quantitative streamflow data and climate descriptors and does not require a priori physical parametrization. However, there are hardly any reports on the influence of runoff separation methods on the flexible hybrid runoff generation model for flood prediction.
The objective of this study was, therefore, to probe into the influence of the runoff partition method on the flexible hybrid hydrological model. The major contributions of this study are presented as follows. (1) The saturation excess runoff generation and infiltration excess runoff generation mechanisms were combined to construct the flexible hybrid runoff generation models, and the XAJ model was adopted as the hydrological modeling framework. (2) An improved two-source runoff partition method considering the heterogeneity of the land cover was proposed. (3) Performances of various runoff partition methods, which consist of the two-source method, improved two-source method, and three-source method, were compared and discussed. Finally, the Xun River basin was chosen as a case study, and the influence of the runoff partition method on the flexible hybrid runoff generation model was investigated.
The paper is arranged as follows: In Section 2, the flexible hybrid runoff generation models and three runoff partition methods are described, and the parameter calibration method and evaluation criteria are introduced. In Section 3, the data and study area are described. In Section 4, the performances of the continuous flow discharge and flood events and the analysis of the runoff components are demonstrated. In Section 5, the model performance and applicability, as well as the nonlinear components of the four strategies, are discussed. Finally, the overall conclusions are drawn.

Construction of Flexible Hybrid Runoff Generation Method
At present, the well-developed and widely utilized conceptual hydrological model for saturation-excess runoff is in a mature state. However, there is still a need for a more precise comprehension of the hybrid runoff generation, particularly the challenge of distinguishing between saturation-excess and infiltration-excess runoff that coexist within a watershed in semi-humid regions. The existing prevalent hybrid runoff models are constructed in accordance with the vertical integration of saturation-excess and infiltration-excess modules, such as the Sacramento model [30].
In this study, to address this issue, we constructed a flexible hybrid runoff generation module by incorporating the vertical integration of saturation-excess and infiltrationexcess modules, referring to Yi et al. [2], Huang et al. [12], and Bao and Zhao [31]. A depiction of the constructed flexible hybrid runoff generation module is shown in Figure  1, and the model fluxes are shown in Table 1.  The hybrid runoff generation model (as shown in Figure 1) simulates the surface runoff and subsurface runoff processes. The surface runoff will be generated when the precipitation rate exceeds the infiltration rate. The surface runoff generation is simulated by combining an infiltration equation with the parabolic infiltration capacity distribution curve. The infiltration capacity distribution curve is similar to the storage capacity distribution curve within the XAJ model, which can be displayed by Equations (1) and (2).
where pi A is the partial area where the infiltration capacity is not more than or equal to capacity; and f is the areal mean infiltration rate, which can be computed by Green-Ampt formula, Horton formula, and Philip infiltration formula, and so on [32][33][34]. In the present study, f was computed by using the Green-Ampt formula [35], which is presented as where f is the infiltration rate (mm·h −1 ); F is the depth of the cumulative infiltration (mm); K is the saturated hydraulic conductivity (mm·h −1 );  is the wetting front capillary pressure head (mm); and   is the change of soil moisture content across the wetting front. Subsurface runoff is simulated by adopting the parabolic storage capacity distribution curve [2,12], which is given by Equation (3). When the soil moisture reaches or exceeds the field capacity, subsurface runoff occurs.
where ps A is the partial pervious area in which the tension water storage capacity is less than or equal to the value WM, which is the tension water capacity at a point, ranging from 0 to a maximum WMM; p A is the pervious area of the watershed; and 2 B is the exponential of distribution of the tension water capacity.

Construction of Runoff Partition Method
The flexible hybrid runoff generation module established in Section 2.1 consists of the surface and subsurface runoffs. However, it ignores the subsurface stormflow, which is ubiquitous, especially in steep and humid watersheds. According to Huang et al. [12] and Yi et al. [2], the soil profile was separated into the upper and lower layers, and the subsurface stormflow in their model was produced between the upper and lower soil layers. Most studies have shown that subsurface stormflow is a saturated (or nearsaturated) water flow phenomenon. To this end, the runoff partition method separated the subsurface stormflow from the saturation excess. In the present study, the runoff partition methods of the two-source and three-source XAJ models were introduced to the flexible runoff generation module. In addition, an improved two-source partition method was proposed to consider the heterogeneity of the subsurface topography and land covers.

Two-Source Runoff Partition Method
According to the original formulation in the two-source XAJ model, runoff was separated into two components using Horton's concept of a final, constant infiltration rate ( c f ). To this end, the steady infiltration rate was used to separate the subsurface stormflow runoff and the subsurface runoff in this study. The flow hydrograph at the outlet of the watershed was composed of the infiltration-excess surface runoff, the saturation-excess subsurface stormflow runoff, and subsurface runoff. The subsurface stormflow runoff and subsurface runoff can be expressed by where PE is the rainfall that exceeds evaporation; sg R is the saturation-excess subsurface runoff; and si R is the saturation-excess subsurface stormflow runoff.

Improved Two-Source Runoff Partition Method
The traditional two-source partition method did not consider the heterogeneity of the subsurface topography and land covers. To this end, we proposed to introduce a parabolic distribution curve to characterize the uneven distribution of the steady infiltration rate between the upper and lower soil layers. The parabolic distribution curve was analogous to the tension water storage capacity distribution curve of the XAJ model, which can be given as where FM is the steady infiltration rate at a point in a watershed, ranging from 0 to a maximum FMM, and 3 B is the exponential of distribution of the steady infiltration rate.
Subsequently, the subsurface stormflow runoff and subsurface runoff can be calculated according to Equation (7) and are expressed by where FR, equaling to g R F A , is the proportion of the runoff-producing region over the whole watershed, and AU is the vertical coordinate corresponding to FM.

Three-Source Runoff Partition Method
According to the principle of hillslope hydrology, the XAJ model was adjusted to a three-source partition method to compute the runoff generation in 1992 [23]. A conceptual structure of a free water reservoir is considered based on the vertical distribution of soil moisture. More details can be found in Chen et al. [36] and Yi et al. [2]. To this end, the three-source runoff partition method was introduced to separate the total runoff of the saturation excess ( g R ). Using the free water reservoir, the depth of the total runoff is divided into three different components, which consist of the surface runoff ( ss R ), subsurface stormflow runoff ( si R ), and subsurface runoff ( sg R ). It is noteworthy that the surface runoff consists of both the saturation excess and infiltration excess in the flexible hybrid runoff generation models.

Hydrological Modeling Framework
In the present study, the XAJ model was utilized as the hydrological modeling framework. The XAJ model includes four modules, which are the evapotranspiration, runoff generation, runoff partition, and runoff routing modules. The saturation-excess runoff generation module of the XAJ model was substituted with the hybrid runoff generation module, and the runoff partition module was eliminated, which is named Strategy P1. Then, the runoff partition module was replaced with the two-source runoff partition module (Strategy P2), the improved two-source runoff partition module (Strategy P3), and the three-source runoff partition module (Strategy P4), respectively. The evapotranspiration and runoff routing modules remain unaltered with the current XAJ model. For detailed information on the parameters of the XAJ model, refer to Yi et al. [2]. The parameters of the runoff partition methods in the four strategies are presented in Table 2. Specifically, the runoff components of the four runoff partition strategies are illustrated in Table 3. Strategy P1 is designed to have no runoff partition module, which only consists of surface and subsurface runoff. Strategies P2, P3, and P4 consist of all three runoff components (i.e., the surface runoff, subsurface stormflow runoff, and subsurface runoff), and the differences lie in that Strategies P2 and P3 do not consider the attenuation through storage caused by the watershed, while Strategy P4 does. In addition, Strategy P4 considers the surface runoff due to the saturation-excess runoff generation mechanism.

Model Calibration and Evaluation
To optimize the parameters of the flexible hybrid runoff generation model, we used the Shuffled Complex Evolution Algorithm (SCE-UA) [28,37]. For parameter calibration, a composite objective function composed of three measurement indexes was used [2,38,39]. The composite objective function and three metrics are given as Several criteria were adopted for the evaluation of the model performance; these criteria consist of the ENS, the Root Mean Square Error (RMSE), the relative flood peak error (Qp), the relative flood peak error (Wp), as well as the flood peak time error (Tp), which can be given as

Study Area and Data
As an important tributary of the Han River, the Xun River is situated in Shaanxi Province, China. Its length is 218 km. The basin area is 6448 km 2 . The average slope of the whole watershed amounts to 2.9‰. The Xun River basin is situated in mountainous terrain with significant variations in topographic relief, in which the runoff generation heterogeneity characteristics of the watershed should be taken into account in the processes of hydrological modeling. The Xun River basin is located in the transition zone between the warm temperate and the northern subtropical regions. It experiences low rainfall in winter, with an average temperature of 7.8 °C. In contrast, it receives abundant rainfall during summer with an average temperature of 24.8 °C. The average annual rainfall is 798 mm in the entire basin. The flood season occurs from July to October, which accounts for 60% of the total annual runoff.
As depicted in Figure 2, the basin contains one flow station (the Xiangjiaping Station) and six meteorological stations. As the outlet of the Xun River basin, the Xiangjiaping hydrological station has an annual average discharge of 63.4 m 3 /s. The highest recorded flood peak discharge at the basin outlet stands at 6090 m 3 /s. Floods in the Xun River basin are typically characterized by short durations and high flow peaks, posing significant threats to downstream areas. Consequently, relying solely on a single runoff generation mechanism for watershed flood prediction in this basin is insufficient, and hybrid runoff generation mechanisms are to be constructed urgently. The meteorological stations provided the six-hourly rainfall and evapotranspiration data spanning from 2010 to 2021. Additionally, we collected the six-hourly runoff data at the Xiangjiaping Station during the same period. To calibrate and validate the constructed model, we specifically chose the annual flow series during the flood seasons. The period from 2010 to 2017 is used for calibration, and the period from 2018 to 2021 is used for validation. However, we excluded the flood season of 2016 due to a significant error in the measured flow discharge. In addition, we identified 23 flood events with flood peak discharge exceeding 500 m 3 /s. Both the calibration and validation periods aligned with the continuous flow discharge data. Moreover, we calculated the antecedent precipitation by using the daily recession coefficient of the water storage, taking into account the initial condition.

Model Calibration and Performance Evaluation for the Continuous Flow Discharge
This study used the XAJ model as the hydrological forecasting framework, replacing the saturation-excess runoff generation module with a hybrid runoff generation module. We designed four runoff partition strategies to explore the effects of the runoff partition method over the flexible hybrid runoff generation model for flood prediction. The other modules remained unaltered within the original XAJ model. The continuous discharge processes in the flood season from 2010 to 2017 were adopted to calibrate the established hydrological model, while the rest were used for validation. Several evaluation indices, including ENS, EKG, and RMSE, were adopted to evaluate the performances of different runoff partition strategies, as listed in Table 4. Figure 3 shows the discharge hydrographs of the four runoff partition methods, in which the calibration years are 2010, 2011, and 2017, and the validation years are 2019, 2020, and 2021.   Figure 3 that Strategy P1 is the poorest in performance. The performances of Strategies P2, P3, and P4 are almost consistent, and the performances of Strategies P3 and P4 are slightly better than that of Strategy P2. However, model complexity control is necessary as the model with more parameters can have higher precision in most cases [40]. The number of parameters of partition methods in Strategies P2, P3, and P4 are 1, 2, and 4, respectively. To this end, it is difficult to determine which strategy can reflect the actual conditions of the natural watershed. Figure 4 presents the density distribution of the measured and simulated discharge hydrographs in order to demonstrate the performances of different strategies more visually. Figure 4 plots the density distribution curves for the measured discharge and simulated discharge hydrographs of the four strategies. The results show that the density curves of the four strategies are similar to that of the measured discharge besides 2020. Similar to the results in Table 4, the performances of the simulated results in 2011, 2019, and 2021 are the best. As shown in Figure 4a,b, the results demonstrate that the simulated results are lower than the measured discharge when the actual discharge is high. The errors in 2020 are large, possibly because the discharge in 2020 is lower than those of the other years. In addition, it can also be observed from Figure 4 that the density distribution curves of Strategy P2 are more consistent with that of the measured discharge.

Performances of Simulated Results of Flood Events
Depending on the application, hydrological modeling can be event-based or continuous [41]. Individual flood events are simulated by event-based models. Mostly, we focus on flood events because of their extensive damage. To this end, 23 flood events were selected from the continuous flow discharge process from 2010 to 2021 to further analyze the performances of flood events.
It is noteworthy that the hydrological models' parameters for both the flood events and the continuous flow discharge were calibrated together. The discharge peak, discharge volume, discharge process, and occurrence time of discharge peak are the four essential elements to describe the discharge hydrograph, and the evaluation index distributions are plotted in Figure 5.  Strategies P3 and P4 outperformed the other strategies, followed by Strategy P2. Specifically, the ENS of Strategies P3 and P4 have a more centralized distribution, demonstrating that the models are more robust than others. The Tp distributions of the four strategies are almost the same. The Qp and Wp distributions of Strategies P3 and P4 are slightly better than those of Strategies P1 and P2. Overall, the hydrological models' performances with runoff partition modules are better than Strategy P1. Although Strategy P4 introduced four parameters to describe the runoff partition module, it did not perform as satisfactorily as expected. Therefore, quantitative studies are needed to determine the differences in the four runoff partition strategies and further explore their influences on the performances of the flexible hybrid runoff generation model.

Comparisons of Runoff Components for the Four Runoff Partition Strategies
In order to explore the quantitative influences of the runoff partition method on the flexible hybrid runoff generation model for flood prediction, we analyze the three runoff components of the four runoff partition strategies. Figure 6 plots the runoff components computed by the four runoff partition strategies for the flow discharge processes in 2010, 2011, 2017, 2019, 2020, and 2021.
It can be observed from Figure 6 that the proportions of each runoff component for the four strategies are significantly different, while the proportion of subsurface runoff for all strategies is high. Since the runoff partition module was not considered, Strategy P1 consists of only surface and subsurface runoffs, and its results are different from other strategies in that the subsurface runoff occupies about 80% of the total runoff. Strategies P2, P3, and P4 contain all three runoff components, while the proportion of the subsurface stormflow runoff of Strategy P4 is tiny, with surface runoff and subsurface runoff each occupying about 50%. Strategy P2 has the highest percentage of subsurface stormflow runoff of all strategies. Based on our previous research [2], the Xun River basin has abundant subsurface stormflow runoff, accounting for about 30% of the total runoff. We noted that the results of component ratios of Strategy P3 are the most consistent with our previous research.

Discussion on Model Performance and Applicability
It is noteworthy that although the ratios of runoff components are significantly different, Strategies P3 and P4 show almost consistent performances, which indeed raised an issue on the parameter equifinality [42].
We tried to analyze the differences between each strategy theoretically. For Strategies P2 and P3, the two-source runoff partition methods are used to divide the infiltration into subsurface stormflow and subsurface runoffs. The differences lie in that Strategy P3 considers the heterogeneity of the steady infiltration of the watershed, and the results show that Strategy P3 performed a lot better than Strategy P2. However, it is noteworthy that Strategy P3 performs consistently with Strategy P4, even though Strategy P4 has a larger number of parameters. Generally speaking, complex models are more likely to achieve better simulation results than simple models, but this does not mean that complex models can accurately reflect the actual hydrological conditions of a watershed. To this end, the physical mechanism of Strategy P4 deserves further clarification and exploration.
The three-source runoff partition method adopted in Strategy P4 has been widely used in the XAJ hydrological model [36,43,44]. Unlike the flexible hybrid runoff generation model used in this study, the XAJ model operates based on the saturationexcess runoff generation mechanism, resulting in surface runoff when there is no available soil moisture storage [45]. This is significantly different from the constructed flexible hybrid runoff generation model used in this study, in which the surface runoff occurs when the precipitation rate exceeds the saturated hydraulic conductivity of the surface soil [2]. When Strategy P4 is adopted to separate the excess rainfall, the surface runoff consists of both the saturation-excess and infiltration-excess runoffs.
In addition, it is still not specific whether it is reasonable to consider both the saturated-and infiltrated-surface runoffs in a hydrological model. Bao et al. proposed an improved Green-Ampt infiltration equation, which is related to the storage of soil moisture [35,46]. Since the infiltration equation considers the influence of the soil moisture storage, the infiltration-excess surface runoff computed using the infiltration equation can be explained by the saturation-excess runoff generation and the infiltration-excess runoff generation mechanisms [47,48]. For this reason, the saturation-excess surface runoff should not be recalculated doubly.
In summary, we designed four strategies to investigate the influences of the runoff partition method on the flexible hybrid runoff generation model for flood prediction. The simulation results demonstrated that the runoff partition method was critical to the flexible runoff generation model. Contrary to common perceptions, Strategy P4 did not perform exceptionally well, and it could not reflect the real conditions of the reality of the watershed. Strategy P4 overestimated the surface runoff and almost ignored the subsurface stormflow runoff. Strategies P2 and P3 showed good performances and, considering the heterogeneity of steady infiltration, could improve the simulation of the discharge hydrographs at the outlet of the watershed.

Discussion on the Nonlinear Components of the Four Strategies
As we discussed in the Introduction Section, the flexible hybrid runoff hydrological models allowed for the selection of alternative representations of runoff processes and various combinations of both linear and nonlinear components. The nonlinear components of each sub-modules can significantly improve the performance of the hydrological model [49,50]. There are 2, 2, 3, and 3 nonlinear components of the four strategies, respectively. In Strategies P1 and P2, the key nonlinear components are the distribution curves of infiltration capacity as well as soil moisture storage capacity in a parabolic manner, and their shapes are determined by two shape parameters, B1 and B2, respectively. In Strategies P3, we proposed a third nonlinear component, namely the steady infiltration capacity distribution curve over the interface, to divide the infiltration from the upper soil into the subsurface stormflow runoff and the subsurface runoff. Simultaneously, the shape of the curve was characterized by a shape parameter, B3. In Strategy P4, the free water storage reservoir was introduced to separate the three runoff components from infiltration.
It is unreasonable to compare models with different degrees of flexibility, and it is obvious that the hydrological model with a higher degree of flexibility is more likely to perform better than those with a low degree of flexibility [12]. This conclusion is consistent with the results in Sections 4.1 and 4.2. The high flexibility of those nonlinear components in Strategies P3 and P4 may answer why both the models have better performances than Strategies P1 and P2 in the study case. For Strategies P3 and P4, the performances of these two strategies show similar results, while the runoff components are significantly different. This phenomenon has been discussed in Section 5.1.
The above discussions demonstrate that enough key nonlinear components are necessary to accommodate complex and diverse hydrological environments. Simultaneously, Huang et al. [12] argued that the nonlinear components ought to possess enough flexibility. In addition, we found that the complex model does not necessarily reflect the actual condition of the real watershed, and the improvement of the model's performances can be due to the complexity of the modeling structure or the increase in the number of parameters.

Conclusions
We constructed the flexible hybrid runoff generation models by combining the saturation-excess runoff generation and infiltration-excess runoff generation mechanisms and used the XAJ model as the hydrological modeling framework in this study. We designed four strategies to investigate the influences of the runoff partition method on the flexible hybrid runoff generation model for flood prediction. Performances of various runoff partition methods, which consist of the two-source, improved two-source, and three-source methods, were compared and discussed based on the continuous flow discharge and flood events. In addition, the model performances and applicability, as well as the nonlinear components of the four strategies, were discussed. The main conclusions are summarized as follows: (1) Strategy P3 and P4 outperform other strategies, followed by Strategy S2. And Strategy P1 with no runoff partition module cannot reflect the actual conditions of the watershed; (2) Although the performance of Strategy P4 is good, it is not applicable to the flexible hybrid runoff generation models because it overestimates the surface runoff and almost ignores the subsurface stormflow runoff; (3) Strategy P3 is a compromise between strategies P2 and P4. It retains the advantages of the free reservoir in Strategy P4 and considers the heterogeneity of the watershed; (4) The runoff partition method is of great influence on the performances of the flexible hybrid runoff generation model. Given that most watersheds are dominated by a mixed runoff generation mechanism rather than a single runoff generation mechanism, our study has great practical value for hydrological modeling and flood prediction. Data Availability Statement: Data will be made available on request.