Estimation of Reservoir Sediment Flux through Bottom Outlet with Combination of Numerical and Empirical Methods

: Sediment deposition issues for reservoirs are important in Taiwan because the severe deposition could excessively decrease the reservoir lifecycle. Extreme storm events usually can carry a massive amount of sediment into reservoirs, and deposition will happen unless the incoming material can pass through sluice gates. When it comes with high concentration, the density current ﬂow is prone to be generated, and the bottom outlets are the most e ﬀ ective sluice gate to release the sediment. In order to improve the sediment release e ﬃ ciency, an accurate estimation of arriving concentration and time of the density current can be useful for the reservoir management. This study develops a two-stage approach which combines a numerical model (SRH2D) and the modiﬁed Rouse equation to predict the sediment ﬂux of the reservoir. The numerical model was veriﬁed and applied to establish the relation between inﬂow and dam face concentration. The modiﬁed Rouse equation then adopted this relation to estimate the proper exponential parameter. As a result, the sediment ﬂux amount at each bottom outlet can be accurately predicted by this equation. With this means, an early warning system can be established for reservoir operation, which can improve release e ﬃ ciency during typhoons.


Introduction
The water crisis is a highly-discussed issue in the world because of increasing population, agriculture irrigation, and flooding issues. In addition, the extreme flooding and sedimentation events caused by the climate change have frequently occurred in the recent years [1][2][3]. Water resource management has become a concerning topic for government agencies. Although the reservoir is the proper solution to store water resources, its capacity may rapidly decrease due to severe sedimentation. According to the statistics, the number of reservoirs losing storage capacity is more than that of those increasing, and the dead storage of reservoirs significantly increases year by year [4,5]. The reservoir deposition, therefore, should be relieved to extend the lifecycle of the reservoir. This urgent water scarcity issue happens not only globally but locally [6,7]; for instance, Taiwan suffers from reservoir deposition and faces the difficulty of water resource management [8]. In order to maintain the storage capacity for the existing reservoirs, several desilting methods were conducted by the Water Resource Agency in Taiwan, including mechanical removal, dredging, and sediment bypass tunnels [9,10]. However, the operational duration is limited for mechanical removal and dredging, and these methods also affect reservoir efficiency. Besides, the transport of sediments will not only be expensive but also cause environmental problems such as dust and traffic issues. The bypass tunnel is the better can be established between the inflow information and the density current, and applied to the empirical-based method. Next, the Rouse equation [33] was modified to be embedded in the second stage by using the numerical outcome and calculating the height of the muddy lake. The corresponding parameter of the Rouse equation was decided by the numerical results and the modified equation was conducted to predict the total amount of release sediment from bottom outlets. The field measurement of the sediment concentration is compared to the result of this twostage approach. As a result, the method can reflect the physical mechanism of the density current flow. The early warning system, therefore, can be established based on this modified equation as the sediment concentration at the reservoir entrance is obtained. Moreover, this two-stage approach can be applied to predict the future trend of the density current and the sediment amount released at the dam face.

Study Site
The original total storage capacity of reservoirs is 2.7 billion m 3 in Taiwan. Sediment deposition occupied more than 30% of the total storage capacity, and the storage was reduced by 0.9 billion to 1.9 billion m 3 after several decades (Huang et al. [34]). The total deposition among the top five reservoirs (Zengwen Reservoir, Feitsui Reservoir, Shimen Reservoir, Deji Reservoir, and Wushe

Study Site
The original total storage capacity of reservoirs is 2.7 billion m 3 in Taiwan. Sediment deposition occupied more than 30% of the total storage capacity, and the storage was reduced by 0.9 billion to 1.9 billion m 3 after several decades (Huang et al. [34]). The total deposition among the top five reservoirs (Zengwen Reservoir, Feitsui Reservoir, Shimen Reservoir, Deji Reservoir, and Wushe Reservoir) has already increased to 0.57 billion m 3 (Yu [35]). It is greatly damaging to the function of the reservoir. Density current flow in a reservoir still cannot fully be captured by numerical modeling even though the sediment flux and efficiency are dependent on it. The study site is the Shimen Reservoir, which is the ranked third largest storage capacity in Taiwan, and is located in North Taiwan. The first storage occurred in 1963, and it employed an earth and rock-fill dam type with a height of 133.1 m, intercepting the Dahan River. The serious deposition described above will not only affect the water storage ability but also dam safety management. Therefore, the focus of this study is how to release as much sediment as possible in the flood period. The estimated release efficiency is presented in the following discussion section. The relative location of the Shimen Reservoir in Taiwan is shown in Figure 2. Reservoir) has already increased to 0.57 billion m 3 (Yu [35]). It is greatly damaging to the function of the reservoir. Density current flow in a reservoir still cannot fully be captured by numerical modeling even though the sediment flux and efficiency are dependent on it. The study site is the Shimen Reservoir, which is the ranked third largest storage capacity in Taiwan, and is located in North Taiwan. The first storage occurred in 1963, and it employed an earth and rock-fill dam type with a height of 133.1 m, intercepting the Dahan River. The serious deposition described above will not only affect the water storage ability but also dam safety management. Therefore, the focus of this study is how to release as much sediment as possible in the flood period. The estimated release efficiency is presented in the following discussion section. The relative location of the Shimen Reservoir in Taiwan is shown in Figure 2. Numerous typhoons have passed Shimen Reservoir in the past five decades, and the bed elevation was raised by inflow sediment over time (see Figure 3). The original bed elevation (1963) at the dam front was 139.20 m above sea level (a.s.l.) and was raised by more than 30 m after several operations. In the current situation, the PRO (Permanent River Outlet) and PPI (Power Plant Intake) are located at the bottom outlet, and this study emphasizes these two bottom outlets in order to investigate the sediment release amount, and further improve the prediction of the release efficiency. Numerous typhoons have passed Shimen Reservoir in the past five decades, and the bed elevation was raised by inflow sediment over time (see Figure 3). The original bed elevation (1963) at the dam front was 139.20 m above sea level (a.s.l.) and was raised by more than 30 m after several operations. In the current situation, the PRO (Permanent River Outlet) and PPI (Power Plant Intake) are located at the bottom outlet, and this study emphasizes these two bottom outlets in order to investigate the sediment release amount, and further improve the prediction of the release efficiency.  Data from three typhoons from 2015 to 2016 were collected as shown in Table 1. As a result, the total release efficiency was between 19.92 to 36.27%. Typhoon Soudelor was adopted to be the calibrated case because the peak discharge was indicated to be the highest value, 5634 m 3 /s, among these three typhoons. In addition, the bottom outlets, PRO and PPI, showed higher release efficiency in comparison with the other outlets. How to well-estimate the sediment movement, then wellpredict the sediment release efficiency of the PRO and PPI is the target of this study and will be discussed below.  Figure 4 shows the simulation area. The inflow boundary is located at Lofu (section No. 32), and the outflow boundaries are situated in front of the Shimen Dam, including the spillway, diversion tunnel, canal, PPI, and PRO. Next, a calculation mesh consisting of 10,230 nodes with a resolution from 5 × 5 to 50 × 50 m is shown in Figure 4. The field bed elevation was measured in February 2015, and the initial simulation elevation was generated, as shown in the same Figure 4. Data from three typhoons from 2015 to 2016 were collected as shown in Table 1. As a result, the total release efficiency was between 19.92 to 36.27%. Typhoon Soudelor was adopted to be the calibrated case because the peak discharge was indicated to be the highest value, 5634 m 3 /s, among these three typhoons. In addition, the bottom outlets, PRO and PPI, showed higher release efficiency in comparison with the other outlets. How to well-estimate the sediment movement, then well-predict the sediment release efficiency of the PRO and PPI is the target of this study and will be discussed below.  Figure 4 shows the simulation area. The inflow boundary is located at Lofu (section No. 32), and the outflow boundaries are situated in front of the Shimen Dam, including the spillway, diversion tunnel, canal, PPI, and PRO. Next, a calculation mesh consisting of 10,230 nodes with a resolution from 5 × 5 to 50 × 50 m is shown in Figure 4. The field bed elevation was measured in February 2015, and the initial simulation elevation was generated, as shown in the same Figure 4.

Governing Equations for Density Current Flow
A 2-D layer-averaged turbidity current model, SRH2D, based on the finite-volume method is presented to solve the governing equations [31]; the interaction with suspended sediment and water on the reservoir bed is considered. With the adoption of SRH2D, several processes should be taken into consideration. First, the model is coupled with a sediment formula to simulate the turbidity current flow pattern. Next, as the turbidity current reaches the dam front, the releasing relationship of each sluice gate is used. The governing equation is described below: where h = current thickness; t = time; x and y = x -direction and y -direction in Cartesian coordinates; u and v = layer-averaged velocity in x -direction and y -direction; D U = average velocity defined as 2 2 u v + ; w e = dimensionless entrainment coefficient, as defined in Equation (2).
This equation is in line with the research of Parker et al. [24].
where i R = bulk Richardson number, the relationship between i R and r F (Froude number) is ; g = acceleration of gravity; t C = total suspended sediment concentration defined as k k C  ; k C = layer-averaged volumetric concentration of the th k sediment size class.
The momentum equations are presented in Equations (4) and (5).

Governing Equations for Density Current Flow
A 2-D layer-averaged turbidity current model, SRH2D, based on the finite-volume method is presented to solve the governing equations [31]; the interaction with suspended sediment and water on the reservoir bed is considered. With the adoption of SRH2D, several processes should be taken into consideration. First, the model is coupled with a sediment formula to simulate the turbidity current flow pattern. Next, as the turbidity current reaches the dam front, the releasing relationship of each sluice gate is used. The governing equation is described below: where h = current thickness; t = time; x and y = x-direction and y-direction in Cartesian coordinates; u and v = layer-averaged velocity in x-direction and y-direction; U D = average velocity defined as √ u 2 + v 2 ; e w = dimensionless entrainment coefficient, as defined in Equation (2). This equation is in line with the research of Parker et al. [24].
where R i = bulk Richardson number, the relationship between R i and F r (Froude number) is R i = 1/F 2 r ; g = acceleration of gravity; C t = total suspended sediment concentration defined as k C k ; C k = layer-averaged volumetric concentration of the kth sediment size class.
The momentum equations are presented in Equations (4) and (5). In the above equations, T xx , T xy , T yy , = dispersion terms defined as (6); R = ρ D /ρ w − 1 = submerged specific gravity of sediment in the turbidity current; ρ D = density of sediment; ρ w = density of ambient water; Z = current top elevation; r w = friction between upper layer water and the turbidity current; ρ = mixture density; τ x and τ y = bed shear stresses in x-direction and y-direction.
Equation (6) is calculated with the Boussinesq formulation (Lai and Greimann [36]) where υ = kinematic viscosity of water; υ t = turbulent eddy viscosity. In addition, a turbulence model, also known as a depth-averaged parabolic model, is presented to calculate the turbulent eddy viscosity. The equation of the parabolic model is shown in Equation (7).
In the above equations, C P = constant and ranges from 0.05-1.00; V * = bed frictional velocity.
In Equations (4) and (5), r w is defined as the friction velocity components, for instance, u * and v * = shear velocities in the x-direction and y-direction. These terms can be written as Equations (8) and (9).
In the equations above, C f = drag coefficient. It is considered to be the total drag friction including both the bed and interfacial drag friction.
Equations (10) and (11) are the sediment concentration equations, which rely on the law of conservation of mass, and could be represented as what follows: In the above formulas, the right hand side is determined as the erosion and deposition term, in which υ k = fall velocity of the kth sediment size class; p k = volume fraction of the kth sediment size class; E k = erosion rate potential; C bk = near-bed concentration of the kth size class; γ = porosity of bed sediment; The relationship between C bk and the shape factor of sediment particle (r 0k ) is C bk = r 0k × C k (Garcia [37]). The shape factor can be computed by Equation (11): where d k = diameter of sediment size k; d gm = geometric mean diameter.

Variation of Muddy Lake
When the density current moves to the dam face, the sluice gate can be operated to release the sediment. However, the bottom outlet sometimes transports only part of inflow material due to the huge flood discharge. As the inflow is larger than the outflow discharge, the density current gives form to a muddy lake. The interface of the clear and turbid water may rise suddenly due to the variation of the inflow hyper-concentrated fluid by time series. The turbid water will diffuse from the lower to the higher water level and finally deposit on the bottom of the reservoir bed if it is not to released. In order to understand the variation of the muddy lake, the relationship of the discharge and storage capacity can be adapted to calculate the raised height of the muddy lake by time series, which is shown in Equations (13) and (14).
H ml t = height of muddy lake at corresponding time t and further conducted as the input parameter for the Equation (17); v ml t = velocity of muddy lake movement; H c t = height of clear lake at corresponding time t; v c t = velocity of clear lake movement. Here, V ml t can be calculated by Equation (15) as: where Q i = inflow discharge and represents as turbidity water; Q o = outflow discharge through different sluice gates. A = storage area in different elevation. Q o (H ml t ) can be distinguished by using Equation (16).
where Q PRO = outflow discharge of the bottom sluice gate (PRO) with a foundation sill located at 169.5 m a.s.l.; Q PPI = outflow discharge of the power plant intake (PPI) with a foundation sill located at 173 m a.s.l.; Q Canal = outflow discharge of the Canal with a foundation sill located at 192.5 m a.s.l.; Q DT = outflow discharge of the diversion tunnel with a foundation sill located at 220.0 m a.s.l.; Q SP = outflow discharge of the spillway with a foundation sill located at 235.0 m a.s.l.

Concentration Calculation
The modified Rouse equation [33] was adapted to predict the sediment release flux at bottom outlets of the Shimen Reservoir. The description of the predicted equation is shown in Figure 5 and Equation (17): where C a = concentration at distances a; a = distances between reservoir bed and the centerline of each bottom outlets; C = concentration at distances H ml t ; D = water depth; e = exponent of concentration distribution, and is the critical parameter to predict the sediment flux from the bottom outlets accurately. In this study, e needs to be calibrated, and the proper range will be presented.
In Equation (17), a is a fixed value and does not need to be adjusted. Next, D and H ml t can be obtained and calculated from the observation data and Equations (13) and (14). In addition, C is the average concentration of the density current. Namely, e is the singular major parameter which needs to be calibrated. The vertical concentration profile of the density current can be determined by the exponent value. Finally, the variation of C a by time series can be calculated rapidly and used to further estimate the released amount of sediment combining the height of aspiration, which will be explained in the next section.
where a C = concentration at distances a ; a = distances between reservoir bed and the centerline of each bottom outlets; C = concentration at distances ml t H ; D = water depth; e = exponent of concentration distribution, and is the critical parameter to predict the sediment flux from the bottom outlets accurately. In this study, e needs to be calibrated, and the proper range will be presented. In Equation (17), a is a fixed value and does not need to be adjusted. Next, D and ml t H can be obtained and calculated from the observation data and Equations (13) and (14). In addition, C is the average concentration of the density current. Namely, e is the singular major parameter which needs to be calibrated. The vertical concentration profile of the density current can be determined by the exponent value. Finally, the variation of a C by time series can be calculated rapidly and used to further estimate the released amount of sediment combining the height of aspiration, which will be explained in the next section.

Height of Aspiration
When the density current reaches the dam face, the sluice gate should be operated to release the sediment. Figure 6 shows two different releasing conditions, the differences of (a) and (b) are the height of the density current. In case (a), the sluice gate owns the higher elevation and the clear water can be released with the turbid water at the same time. On the other hand, the released fluid in case (b) is mainly density current.
Fan [38] exploited a laboratory experiment to determine the flow pattern around the orifice.
p Fr the densimetric Froude number, was used to explain the aspiration height. In case (a), the equations used are shown as Equations (18) and (19); in case (b), the equation is shown as Equation (20).

Height of Aspiration
When the density current reaches the dam face, the sluice gate should be operated to release the sediment. Figure 6 shows two different releasing conditions, the differences of (a) and (b) are the height of the density current. In case (a), the sluice gate owns the higher elevation and the clear water can be released with the turbid water at the same time. On the other hand, the released fluid in case (b) is mainly density current.

Research Results
Given that the arrival time and the sediment concentration of the density current can be simulated with the model, it is possible to estimate the sediment release efficiency and the sediment amount of PRO. The transport duration of the density current from the entrance to the outlet was simulated and compared to the measurement during the typhoon (see Figure 7). Comparison of the Fan [38] exploited a laboratory experiment to determine the flow pattern around the orifice. Fr p the densimetric Froude number, was used to explain the aspiration height. In case (a), the equations used are shown as Equations (18) and (19); in case (b), the equation is shown as Equation (20).
Q = total discharge through orifice; h L = height of aspiration as in definition sketch. The aspiration height can be estimated by rearranging Equation (20).

Research Results
Given that the arrival time and the sediment concentration of the density current can be simulated with the model, it is possible to estimate the sediment release efficiency and the sediment amount of PRO. The transport duration of the density current from the entrance to the outlet was simulated and compared to the measurement during the typhoon (see Figure 7). Comparison of the measured and simulated results showed that the model could capture the density movement. The travel velocity of the simulated and measured data is close along the Shimen Reservoir. The difference of the arrival time at the dam face between observed and simulated is 0.22 h, which is good enough for the gate operation.  In addition, S24 is located at 9100 m, S20 is located at 7700 m, S15 is located at 5800 m, S12 is located at 4200 m, and S07 is located at 1200 m).  In addition, S24 is located at 9100 m, S20 is located at 7700 m, S15 is located at 5800 m, S12 is located at 4200 m, and S07 is located at 1200 m). Figure 8 shows the comparison of measured and simulated peak concentration in Typhoon Soudelor. The measured and simulated results own the same peak concentration value, 18,287 ppm, at the location 16,200 m. Each comparison point along Shimen Reservoir is 9100 m, 5800 m, 4200 m, and 1200 m. As a result, the measured data is missed at the location of 9100 and 5800 m. It shows a decreased trend from 16,200 m to 1200 m. In the other hands, the simulation shows the continuity decline trend from 16,177 to 5037 ppm. Above all, the simulated result fitted the measured data at the dam front. The values were 5623, and 6387 ppm. In other words, the simulated outcome shows the close decay trend by comparing it with the measured data. Figure 7. Arrival time at sections during Typhoon Soudelor (Lofu, at the location of 16,200 m is defined to the inflow boundary. In addition, S24 is located at 9100 m, S20 is located at 7700 m, S15 is located at 5800 m, S12 is located at 4200 m, and S07 is located at 1200 m). Figure 8 shows the comparison of measured and simulated peak concentration in Typhoon Soudelor. The measured and simulated results own the same peak concentration value, 18,287 ppm, at the location 16,200 m. Each comparison point along Shimen Reservoir is 9100 m, 5800 m, 4200 m, and 1200 m. As a result, the measured data is missed at the location of 9100 and 5800 m. It shows a decreased trend from 16,200 m to 1200 m. In the other hands, the simulation shows the continuity decline trend from 16,177 to 5037 ppm. Above all, the simulated result fitted the measured data at the dam front. The values were 5623, and 6387 ppm. In other words, the simulated outcome shows the close decay trend by comparing it with the measured data.   Figure 9 shows the released sediment through the PRO and PPI during Typhoon Soudelor. A more significant error for the discharge sediment appeared after the first 6 and 14 h of the PRO and PPI. The reason for this could be related to the outflow elevation. In the Shimen Reservoir, the PRO and PPI are the bottom outlets located at 169.5 and 173 m. The model shows a reasonable trend, but the sediment discharge of bottom outlets can be further improved. The total outflow amounts of sediment were 103,323 and 81,715, and 50,882 and 63,527 m 3 of the PRO and PPI in the observation and simulation respectively. SRH2D shows that the limitation of the 2-D layer-averaged model and the sediment flux at the bottom outlet is difficult to simulate. Therefore, a couple of numerical and empirical method should be developed to raise the prediction accuracy to establish a robust early-warning system of the reservoir operation. sediment were 103,323 and 81,715, and 50,882 and 63,527 m 3 of the PRO and PPI in the observation and simulation respectively. SRH2D shows that the limitation of the 2-D layer-averaged model and the sediment flux at the bottom outlet is difficult to simulate. Therefore, a couple of numerical and empirical method should be developed to raise the prediction accuracy to establish a robust early-warning system of the reservoir operation.

Discussion
First, the decay rate of sediment concentration should be determined. Three typhoons were applied and shown in Figure 10 Figure 10 represents all of the simulations in different events and investigates the relationship between dimensionless distance and decay rate. The results show that the peak concentration decreases along the distance, and the

Discussion
First, the decay rate of sediment concentration should be determined. Three typhoons were applied and shown in Figure 10. The decay rates of observation and simulation are close. Values for Typhoons Soudelor, Dujuan, and Megi are 0.35 and 0.31; 0.31 and 0.33; and 0.32 and 0.31. The simulation reports a consistent and stable trend along the reservoir. Figure 10 represents all of the simulations in different events and investigates the relationship between dimensionless distance and decay rate. The results show that the peak concentration decreases along the distance, and the relationship between the two is highly positive. In addition, a regression formula, Equation (21), can be applied to estimate the peak concentration of different locations as: where D R = decay rate of the sediment concentration and M D = dimensionless distance. relationship between the two is highly positive. In addition, a regression formula, Equation (21), can be applied to estimate the peak concentration of different locations as: Where R D = decay rate of the sediment concentration and D M = dimensionless distance. The relationship between the inflow peak discharge and travel velocity is shown in Figure 11. The simulated data include Typhoons Soudelor, Dujuan, and Megi. The highest peak discharge was observed during Typhoon Soudelor, the middle value was seen during Typhoon Megi, and the event undergoing the shortest peak discharge was Typhoon Dujuan. For Soudelor and Megi, the results show that the travel velocity is proportional to the increase in peak discharge. In addition, the peak discharge for Typhoon Dujuan was the shortest, and its travel velocity was ranked third. As a result, the relationship between peak discharge and travel velocity presents a highly positive correlation. Equation (22) can be used to estimate the travel velocity. The relationship between the inflow peak discharge and travel velocity is shown in Figure 11. The simulated data include Typhoons Soudelor, Dujuan, and Megi. The highest peak discharge was observed during Typhoon Soudelor, the middle value was seen during Typhoon Megi, and the event undergoing the shortest peak discharge was Typhoon Dujuan. For Soudelor and Megi, the results show that the travel velocity is proportional to the increase in peak discharge. In addition, the peak discharge for Typhoon Dujuan was the shortest, and its travel velocity was ranked third. As a result, the relationship between peak discharge and travel velocity presents a highly positive correlation. Equation (22) can be used to estimate the travel velocity.
V T = dimensionless travel velocity of peak concentration of density current, and Q p = peak dimensionless discharge of density current. According to Figures 10 and 11, the inflow concentration graph can shift to dam face, and to be the outflow concentration C shown in Equation (17), and the sediment flux of bottom outlets is calculated by the same Equation (17). This study focuses on the release sediment amount of the PRO and PPI, and the release condition was mainly related to Figure 6b, aspiration of turbid water from the muddy lake through a bottom outlet. The predicted sediment flux of the PRO (Figure 12a) presented a close trend to the measured data before the 22nd hour. After that time, the predicted sediment flux showed a slowly increased pattern than that of the measured data between 23rd and 33rd hour. Then, the measurement indicates a lower increased trend than that of the prediction. Finally, the prediction presents a higher accuracy than the simulation by comparing with the measurement. Figure 12b shows the predicted result of the PPI, and it presents similarly to the PRO. The prediction indicates a close value to the measured data after 42nd hour. Moreover, it performs a more reasonable result than the simulation. This study focuses on the release sediment amount of the PRO and PPI, and the release condition was mainly related to Figure 6b, aspiration of turbid water from the muddy lake through a bottom outlet. The predicted sediment flux of the PRO (Figure 12a) presented a close trend to the measured data before the 22nd hour. After that time, the predicted sediment flux showed a slowly increased pattern than that of the measured data between 23rd and 33rd hour. Then, the measurement indicates a lower increased trend than that of the prediction. Finally, the prediction presents a higher accuracy than the simulation by comparing with the measurement. Figure 12b shows the predicted result of the PPI, and it presents similarly to the PRO. The prediction indicates a close value to the measured data after 42nd hour. Moreover, it performs a more reasonable result than the simulation. Table 2 is the comparison of the measured, simulated, and predicted release amount of PRO and PPI during these three typhoons. In Soudelor case, the performance of the prediction is better than the simulation mentioned before. Base on the same calculation procedure, the verified typhoons, Dujuan, and Megi, also show higher accuracy than the simulation. Above all, the prediction of the three typhoons shows a lower error than the simulated result.   Figure 11. Relationship between peak discharge and travel velocity. The accurately predicted ability relies on the critical parameter, e , which presented in the equation (17). A reasonable and suitable value will improve the accuracy of the predicted release amount in different typhoons. Figure 13 shows the relationship between e and the ratio of the fluid ( i Q ) and sediment ( s Q ). The scale of different typhoons including the inflow discharge and sediment will affect the e value, and Figure 13 shows the recommendation range (0.2-0.5) in the future  Table 2 is the comparison of the measured, simulated, and predicted release amount of PRO and PPI during these three typhoons. In Soudelor case, the performance of the prediction is better than the simulation mentioned before. Base on the same calculation procedure, the verified typhoons, Dujuan, and Megi, also show higher accuracy than the simulation. Above all, the prediction of the three typhoons shows a lower error than the simulated result. The error percentage of the PRO of prediction and simulation is 1.89% and 20.91%, 5.58% and 52.75%, 1.45%, and 51.96% respectively. The error percentage of the PPI of prediction and simulation is 14.03 and 24.85, 5.98 and 17.89, 8.09 and 23.31, respectively. Furthermore, the total release amount of Soudelor, Dujuan, and Megi of the measurement, simulation, and prediction is 192,577, 167,459, and 194,083 from PRO; and 165,145, 174,634, 171,352 m 3 from PPI. A significant result presents that the prediction is a robust approach which can grasp the release pattern from bottom outlets because of the mismatch between measurement and prediction is only 0.78% and 5.75%. The accurately predicted ability relies on the critical parameter, e, which presented in the equation (17). A reasonable and suitable value will improve the accuracy of the predicted release amount in different typhoons. Figure 13 shows the relationship between e and the ratio of the fluid (Q i ) and sediment (Q s ). The scale of different typhoons including the inflow discharge and sediment will affect the e value, and Figure 13 shows the recommendation range (0.2-0.5) in the future scenario.

Conclusion
This study develops a two-stage approach which combines numerical and empirical methods to estimate the sediment flux released through bottom outlets simply with the sediment concentration at the entrance in the reservoir. A 2-D model, SRH2D, was adopted to simulate the flow pattern of the density current in the reservoir and verified with three typhoons. The research result shows the 2-D layer-averaged model was limited to well describe the movement of the current in the reservoir. The sediment flux released through the bottom outlet could not be captured with the numerical model but this was improved with a modified equation based on the Rouse equation. The equation was modified by using the estimation of the muddy lake thickness and the aspiration height to

Conclusions
This study develops a two-stage approach which combines numerical and empirical methods to estimate the sediment flux released through bottom outlets simply with the sediment concentration at the entrance in the reservoir. A 2-D model, SRH2D, was adopted to simulate the flow pattern of the density current in the reservoir and verified with three typhoons. The research result shows the 2-D layer-averaged model was limited to well describe the movement of the current in the reservoir. The sediment flux released through the bottom outlet could not be captured with the numerical model but this was improved with a modified equation based on the Rouse equation. The equation was modified by using the estimation of the muddy lake thickness and the aspiration height to calibrate the critical parameter, e. The prediction value of this two-stage approach is shown to be highly accurate and efficient when compared to the numerical simulation.
The advantage of this two-stage approach is that it is user-friendly and quickly accessible. The predicted release amount from the bottom outlet can be grasped when the inflow sediment concentration and the parameter e, are determined. In other words, the reasonable future variation of the density current can be presented before 8 to 10 h. The modified empirical approach can be combined with real-time monitoring to further improve the early-warning ability.

Conflicts of Interest:
There is no conflict of interest. Notation a distances between reservoir bed and the centerline of bottom outlet A storage area in different elevation C concentration at distances y t C a concentration at distances a C bk near-bed concentration of the kth size class C f drag coefficient C P constant and ranges from 0.05-1.00 concentration of the peak discharge Q SP outflow discharge of Spillway r 0k shape factor of sediment particle r w friction between upper ambient water and the turbidity current R current specific gravity R i bulk Richardson number T xx , T xy , T yy depth-averaged stresses due to turbulence and dispersion t time u layer-averaged velocity in x-direction u * shear velocity in the x-direction U D average velocity v layer-averaged velocity in y-direction v c t velocity of clear lake movement v ml t velocity of muddy lake movement V T dimensionless travel velocity of peak concentration of density current v * shear velocities in the y-direction V * bed frictional velocity x x-direction in Cartesian coordinate y y-direction in Cartesian coordinate Z current top elevation Z b bed elevation ρ mixture density ρ D density of sediment ρ w density of ambient water τ x bed shear stresses in x-direction τ y bed shear stresses in y-direction υ kinematic viscosity of water υ k fall velocity of the kth sediment size class υ t turbulent eddy viscosity γ porosity of bed sediment