Simplified Simulation Method for Flood-Induced Bend Scour — A Case Study Near the Shuideliaw Embankment on the Cho-Shui River

The modeling of flood-induced bend scour near embankment toes can provide important information for river engineering, embankment safety warnings, and emergency action management. During the rainy seasons, short-term general scour and bend scour are the most common causes for the failure of reinforced concrete embankments in Taiwan. To gain a deeper understanding of the scouring process near levee foundations, this study proposed a straightforward and practical method for bend scour simulation. The proposed simulation method is subdivided into three stages: two-dimensional flow simulation, general scour estimation, and bend scour estimation. A new bend-scour computation equation is proposed and incorporated into a two-dimensional hydraulic finite-volume model for simulating the evolution of bend scour depth around embankment toes. The proposed method is applied to simulate the temporal evolution of bend scouring near the Shuideliaw Embankment on the Cho-Shui River in Taiwan, where serious failure occurred during the June 2012 monsoon. Field data were gathered using the numbered-brick technique at the Shuideliaw Embankment to demonstrate the accuracy of the proposed method. The results of the bend scour simulations compared reasonably well with field measurements, indicating close agreement in terms of water levels and bend scour depths near the Shuideliaw Embankment. The proposed method was found to quickly estimate the maximum short-term general scour and bend scour depths for further enhancement of the safety of the embankment toe.


Introduction
Embankments are natural or man-made hydraulic structures used to prevent inundation disasters in adjacent flood-prone areas.In general, embankment failure occurs during a flood as the result of several possible factors: (1) overtopping of the embankment crest; (2) scouring at the toe; (3) piping or internal erosion; and (4) toe or slope sliding of the embankment [1][2][3][4].These potential mechanisms can result in the breaching or failure of embankments, which can occur as a single mechanism or a combination of different mechanisms.The failure of embankments due to multiple mechanisms may cause widespread flooding and result in significant loss of life.In Taiwan, the Water Resource Agency (WRA) [5] reported that approximately 43% of embankment failures are associated with scouring at the toe.On 8 August 2009, embankment failures due to Typhoon Morakot caused catastrophic inundation and serious property damage in Southern Taiwan.The Shuideliaw Embankment on the Cho-Shui Water 2017, 9, 324 2 of 19 River was damaged during the June 2012 monsoon and Typhoon Saola (August 2012), with peak discharges of 6120 and 10,600 m 3 /s, respectively.Moreover, the adjacent irrigation channels were damaged during these flood events.
To mitigate potential embankment hazards, the prediction of flood-induced embankment failure is a basic necessity.A number of numerical models have been developed over the past 50 years to predict the embankment erosion and failure process [6][7][8][9][10][11][12][13].These models can be divided into two major categories: (1) parametric models and (2) physically based models.The parametric model uses regression analysis to establish statistically based equations that can be used to estimate embankment breaching width and discharge.The physically based model uses sediment transport formulas to estimate the rate of embankment erosion according to a reasonable approximation of the evolution of embankment breaching.According to a review of the literature, these models can be used to simulate embankment breach formation, growth, and closure.Although many published studies have investigated earthen embankment erosion and failure processes, there has been little focus on reinforced concrete embankments.In Taiwan, an extreme flood can erode the riverbed near embankment toes and result in their failure.Hence, scour has been the major safety concern for reinforced concrete embankment structures.
The more commonly used approach adopts numerical methods to simulate the scouring around cross-river structures.Several studies have adopted numerical methods to investigate the scour process around bridge piers [14][15][16][17][18] and spur dikes [19][20][21].Time-consuming three-dimensional (3D) numerical models produce sharp resolutions of 3D flow fields and local scour around hydraulic structures under laboratory conditions.Moreover, these 3D models have not been as successful in practical applications, where more field data are required.By contrast, two-dimensional (2D) numerical models can provide accurate hydraulic characteristics at the locations far from the local structures.Additionally, from a practical perspective, 2D models can provide fundamental approach flow information to estimate both general and local scour depths.This enables the combination of 2D models with reliable scour-depth equations for the simulation of total scour processes [22].
As the majority of embankment failure in Taiwan is due to toe scouring, a new practical method is required that explicitly relates to embankment toe scour.The scouring process around the embankment toe in a river bend leads to bend scouring at the toe of the outer embankment, in which riverbed variation during typhoon periods results in the undermining of embankment structures [23].Many researchers proposed formulas to calculate the bend scour depth [23].However, when applied to a field river, these equations might not be suited because the sediment transport characteristics in an alluvial river are different from those in an experimental flume.Thus, a reliable method is needed for estimating the maximum scour depth in a river bend in order to design appropriate embankment protection or establish an embankment warning system.In addition, according to the literature review, most previous studies have focused on embankment erosion and earthen embankment failure.Few investigators have reported the simulations of embankment toe scour in a river bend, especially for reinforced concrete embankment structures.
Therefore, this paper follows the idea of Hong et al. [22] to propose a new modeling methodology for simulating the temporal evolution of scour depth at embankment toes in river bend reaches.The difference between the previous study [22] and the presented study is mainly in the investigated river structures.The previous study [22] investigated the bridge scour, and the present study focuses on the bend scour at the embankment toe.Due to the different physical mechanism of bridge and bend scours, a new bend scour computation equation is proposed and integrated with a 2D shallow water flow model based on the finite-volume model [24][25][26], leading to a new and reliable method, which is the main contribution of this paper.Its practical applicability is demonstrated by a case study of the Cho-Shui River, in which the Shuideliaw Embankment is located in a gravel-bed river bend.Due to river bend effects, bend scour occurs at the outer embankment, i.e., the riverside of the Shuideliaw Embankment.Hence, this paper investigates flood-induced embankment toe scour by using the field data collected with the numbered-brick method.The proposed model is then verified through simulations of bend scour depth at the Shuideliaw Embankment under two flood events.Comparisons with three existing bend scour equations [27][28][29] are made between the simulated and field-measured maximum scour depths [30,31].The practicability of the method is further confirmed through its application to the embankment failure event that occurred during Typhoon Saola in 2012.Based on nine scenario simulations, a new practical embankment safety curve representing the scoured depth-discharge relationship is finally proposed for use in river management.

Field Data Collection and Formula Derivation
As described in Section 1, several bend scour equations [23] were proposed based on the experimental results in the bend flume.These equations might not be suited to estimating bend scour in a field river.Therefore, a reliable bend scour computation equation is needed for investigating the present study site.Field data collection and equation derivation for bend scour estimation are described in this section.
Field measurement of embankment erosion or scour is very difficult, especially during high flow conditions.Flow conditions during floods are extremely dangerous, complicating the measurement of riverbed variation such as riverbed scour.In this study, 12 sets of bend scour data were collected to develop a new bend scour computation equation for rivers in Taiwan.Specifically, 10 and two sets of scour data were collected from the basins of the Cho-Shui and Da-Chia Rivers, respectively, for the database.Two sets of data collected in the Cho-Shui River basin were near the Shuideliaw Embankment and an irrigation channel directly behind the embankment.Therefore, the scour data collected near the Shuideliaw Embankment can be used to accurately evaluate the performance of the proposed method.

Site Description
The field study site was chosen at the middle reach of the longest river in Taiwan, Cho-Shui River, from Chi-Chi Weir to Formosa Freeway No. 3 Bridge on the Cho-Shui River.A map of the Cho-Shui River basin and an aerial photograph near the Shuideliaw Embankment are shown in Figures 1  and 2, respectively.The Cho-Shui River upstream of the Chi-Chi Weir, which was constructed to provide water for irrigation and industrial consumption, has a watershed area of 2034 km 2 .The aerial photograph in Figure 2a shows two concave embankments located near the Mingchu Bridge: one on the right-hand side upstream (but facing downstream) of the Mingchu Bridge, and the other on the left-hand side downstream of the bridge, where the Shuideliaw Embankment is located.bend scour equations [27][28][29] are made between the simulated and field-measured maximum scour depths [30,31].The practicability of the method is further confirmed through its application to the embankment failure event that occurred during Typhoon Saola in 2012.Based on nine scenario simulations, a new practical embankment safety curve representing the scoured depth-discharge relationship is finally proposed for use in river management.

Field Data Collection and Formula Derivation
As described in Section 1, several bend scour equations [23] were proposed based on the experimental results in the bend flume.These equations might not be suited to estimating bend scour in a field river.Therefore, a reliable bend scour computation equation is needed for investigating the present study site.Field data collection and equation derivation for bend scour estimation are described in this section.
Field measurement of embankment erosion or scour is very difficult, especially during high flow conditions.Flow conditions during floods are extremely dangerous, complicating the measurement of riverbed variation such as riverbed scour.In this study, 12 sets of bend scour data were collected to develop a new bend scour computation equation for rivers in Taiwan.Specifically, 10 and two sets of scour data were collected from the basins of the Cho-Shui and Da-Chia Rivers, respectively, for the database.Two sets of data collected in the Cho-Shui River basin were near the Shuideliaw Embankment and an irrigation channel directly behind the embankment.Therefore, the scour data collected near the Shuideliaw Embankment can be used to accurately evaluate the performance of the proposed method.

Site Description
The field study site was chosen at the middle reach of the longest river in Taiwan, Cho-Shui River, from Chi-Chi Weir to Formosa Freeway No. 3 Bridge on the Cho-Shui River.A map of the Cho-Shui River basin and an aerial photograph near the Shuideliaw Embankment are shown in Figures 1 and 2, respectively.The Cho-Shui River upstream of the Chi-Chi Weir, which was constructed to provide water for irrigation and industrial consumption, has a watershed area of 2034 km 2 .The aerial photograph in Figure 2a shows two concave embankments located near the Mingchu Bridge: one on the right-hand side upstream (but facing downstream) of the Mingchu Bridge, and the other on the left-hand side downstream of the bridge, where the Shuideliaw Embankment is located.To prevent embankment erosion, the Fourth River Management Office constructed the Cho-Shui and Shuideliaw Embankments on the right-and left-hand sides (facing downstream) near the Mingchu Bridge, respectively.The Mingchu Bridge and Formosa Freeway No. 3 Bridge are 6.5 and 8.5 km, respectively, downstream of the Chi-Chi Weir.According to a field survey by the Water Resources Agency, the mean channel slope in the middle reach of the Cho-Shui River is 0.006, and there is no tributary between the Chi-Chi Weir and the Formosa Freeway No. 3 Bridge; therefore, the flow discharge at Mingchu Bridge and Shuideliaw Embankment was assumed to be the flow released from the Chi-Chi Weir.Because the distance between the field site and the river mouth is more than 45 km, the tidal effect does not influence the flows at the Mingchu Bridge and Shuideliaw Embankment.
Bed samples obtained in 2013 in the vicinity of the Shuideliaw Embankment show that the sediment was composed of non-cohesive alluvial gravels with various sands.Figure 3 shows the sediment particle size distribution near the Shuideliaw Embankment, clearly demonstrating that the distribution of particle size is bimodal.A survey of the bed materials revealed that the returned sediment sizes are of approximately d16 = 1.7 mm, median size d50 = 108 mm, and d84 = 329 mm.To prevent embankment erosion, the Fourth River Management Office constructed the Cho-Shui and Shuideliaw Embankments on the right-and left-hand sides (facing downstream) near the Mingchu Bridge, respectively.The Mingchu Bridge and Formosa Freeway No. 3 Bridge are 6.5 and 8.5 km, respectively, downstream of the Chi-Chi Weir.According to a field survey by the Water Resources Agency, the mean channel slope in the middle reach of the Cho-Shui River is 0.006, and there is no tributary between the Chi-Chi Weir and the Formosa Freeway No. 3 Bridge; therefore, the flow discharge at Mingchu Bridge and Shuideliaw Embankment was assumed to be the flow released from the Chi-Chi Weir.Because the distance between the field site and the river mouth is more than 45 km, the tidal effect does not influence the flows at the Mingchu Bridge and Shuideliaw Embankment.
Bed samples obtained in 2013 in the vicinity of the Shuideliaw Embankment show that the sediment was composed of non-cohesive alluvial gravels with various sands.Figure 3

Field Data
Su and Lu [32] proposed a highly efficient numbered-brick method to measure the general scour in natural channels.Lu [31] adopted this direct observation method for measuring the maximum bend scour depth near the Shuideliaw Embankment.In the present study, bend scour field data were collected as listed in Table 1.All of the study sites were associated with the scour measurements near the embankment toe and all of the scour depths were measured using the same numbered-brick method.The peak flow discharge Qp ranged from 164 to 5309 m 3 /s and the flow intensity (unit peak flow discharge qp) ranged from 1.31 to 25.37 m 2 /s.The channel slope s0 was also an important factor, with a range from 0.0063 to 0.0153, and the median particle size d50 of the bed materials was from 3.43 to 168 mm.The characteristic parameter for bend scour is the outer radius of the bend R0, which was recorded to range from 169 to 800 m.Finally, the measured maximum scour depth in a river bend was from 0.22 to 5.23 m.With these measured data, a reliable bend scour computation equation was developed, as described in Section 2.3.Furthermore, the measured data for two flood events, Typhoon Trami and Typhoon Usagi, were adopted in the present study to verify the accuracy of the proposed method, as described in Section 4.2. Figure 4 shows the flow hydrographs released from the Chi-Chi Weir for Typhoon Trami (August 2013) and Typhoon Usagi (September 2013).

Field Data
Su and Lu [32] proposed a highly efficient numbered-brick method to measure the general scour in natural channels.Lu [31] adopted this direct observation method for measuring the maximum bend scour depth near the Shuideliaw Embankment.In the present study, bend scour field data were collected as listed in Table 1.All of the study sites were associated with the scour measurements near the embankment toe and all of the scour depths were measured using the same numbered-brick method.The peak flow discharge Q p ranged from 164 to 5309 m 3 /s and the flow intensity (unit peak flow discharge q p ) ranged from 1.31 to 25.37 m 2 /s.The channel slope s 0 was also an important factor, with a range from 0.0063 to 0.0153, and the median particle size d 50 of the bed materials was from 3.43 to 168 mm.The characteristic parameter for bend scour is the outer radius of the bend R 0 , which was recorded to range from 169 to 800 m.Finally, the measured maximum scour depth in a river bend was from 0.22 to 5.23 m.With these measured data, a reliable bend scour computation equation was developed, as described in Section 2.3.Furthermore, the measured data for two flood events, Typhoon Trami and Typhoon Usagi, were adopted in the present study to verify the accuracy of the proposed method, as described in Section 4.2. Figure 4 shows the flow hydrographs released from the Chi-Chi Weir for Typhoon Trami (August 2013) and Typhoon Usagi (September 2013).

Field Data
Su and Lu [32] proposed a highly efficient numbered-brick method to measure the general scour in natural channels.Lu [31] adopted this direct observation method for measuring the maximum bend scour depth near the Shuideliaw Embankment.In the present study, bend scour field data were collected as listed in Table 1.All of the study sites were associated with the scour measurements near the embankment toe and all of the scour depths were measured using the same numbered-brick method.The peak flow discharge Qp ranged from 164 to 5309 m 3 /s and the flow intensity (unit peak flow discharge qp) ranged from 1.31 to 25.37 m 2 /s.The channel slope s0 was also an important factor, with a range from 0.0063 to 0.0153, and the median particle size d50 of the bed materials was from 3.43 to 168 mm.The characteristic parameter for bend scour is the outer radius of the bend R0, which was recorded to range from 169 to 800 m.Finally, the measured maximum scour depth in a river bend was from 0.22 to 5.23 m.With these measured data, a reliable bend scour computation equation was developed, as described in Section 2.3.Furthermore, the measured data for two flood events, Typhoon Trami and Typhoon Usagi, were adopted in the present study to verify the accuracy of the proposed method, as described in Section 4.2. Figure 4 shows the flow hydrographs released from the Chi-Chi Weir for Typhoon Trami (August 2013) and Typhoon Usagi (September 2013).

Proposed Bend Scour Computation Equation
According to Melville and Coleman [23], bend scour may occur in bend reaches and is categorized as short-term general scour.Several empirical equations for estimating bend scour have been developed in the literature [23].Three previous empirical equations are used in this study, and the corresponding schematic definition terms are shown in Figure 5.

Proposed Bend Scour Computation Equation
According to Melville and Coleman [23], bend scour may occur in bend reaches and is categorized as short-term general scour.Several empirical equations for estimating bend scour have been developed in the literature [23].Three previous empirical equations are used in this study, and the corresponding schematic definition terms are shown in Figure 5. First, Galay et al. [27] developed an empirical equation for rigid (non-eroding outer bank) bend of 60° in a gravel-bed river reach: where bs h denotes the maximum flow depth in the bend cross-section B'B, u h is the water depth in the upstream straight reach A'A, W denotes the water surface width, and c R is the centerline radius of the bend.Next, using the scour data from 70 bends in a reach of the Red River between Arkansas and Louisiana in the USA, Thorne [28] developed the following regression equation: Notably, Equation ( 2) is limited to / c R W greater than 2.0 based on the scour dataset.
However, a best-fit regression equation was found by the U.S. Army Corps of Engineering [29]: 2.57 0.36 ln These three equations were developed using the bend factors c R and W as the primary parameters.However, on the basis of the bend scour data collected in this study, the sediment size as well as the standard deviation of particle size distribution was found to significantly affect the bend scour depth.Moreover, the composite of the sediment around the embankment toe also affects the scouring processes.For example, the scouring process would be different under the same approach flow conditions with non-cohesive or cohesive sediments.In this study, the field data First, Galay et al. [27] developed an empirical equation for rigid (non-eroding outer bank) bend of 60 • in a gravel-bed river reach: where h bs denotes the maximum flow depth in the bend cross-section B'B, h u is the water depth in the upstream straight reach A'A, W denotes the water surface width, and R c is the centerline radius of the bend.Next, using the scour data from 70 bends in a reach of the Red River between Arkansas and Louisiana in the USA, Thorne [28] developed the following regression equation: Notably, Equation ( 2) is limited to R c /W greater than 2.0 based on the scour dataset.However, a best-fit regression equation was found by the U.S. Army Corps of Engineering [29]: These three equations were developed using the bend factors R c and W as the primary parameters.However, on the basis of the bend scour data collected in this study, the sediment size as well as the standard deviation of particle size distribution was found to significantly affect the bend scour depth.Moreover, the composite of the sediment around the embankment toe also affects the scouring processes.For example, the scouring process would be different under the same approach flow conditions with non-cohesive or cohesive sediments.In this study, the field data collected in Table 1 were non-cohesive sediments.To focus on the field site we investigated, only non-cohesive sediments were considered in developing bend scour formula.Therefore, this study proposes a new computation equation.The characteristic parameters of bend scour estimation can be defined as: where d bs is the bend scour depth; q is the unit discharge; ρ is the density of fluid; g is the gravitational acceleration; s 0 is the channel slope; d 50 is the median particle diameter; σ g is the geometric standard deviation of the particle size distribution; ρ s is the density of sediment particles; and R 0 (= R c + 0.5 W) is the outer radius of the bend.Based on Buckingham's π theorem, bend scour depth can be written in the following dimensionless form: In the present study, multiple regression analysis is then conducted with the measured field data listed in Table 1: which has a coefficient (R-squared value) of 0.924, indicating that the short-term bend scour depths computed by Equation ( 6) fit the measured data well, as shown in Figure 6.It should be noted that only 10 sets of data were used for calibration, while the remaining two sets of data were used for verification, as indicated in Figure 6.Bend scour depth was found to be directly proportional to unit discharge, channel slope, and water surface width, and inversely proportional to the centerline radius of the bend.
Water 2017, 9, 324 8 of 19 collected in Table 1 were non-cohesive sediments.To focus on the field site we investigated, only non-cohesive sediments were considered in developing bend scour formula.Therefore, this study proposes a new computation equation.The characteristic parameters of bend scour estimation can be defined as: where dbs is the bend scour depth; q is the unit discharge; ρ is the density of fluid; g is the gravitational acceleration; 0 s is the channel slope; d50 is the median particle diameter; σg is the geometric standard deviation of the particle size distribution; ρs is the density of sediment particles; and R0 (= Rc + 0.5 W) is the outer radius of the bend.Based on Buckingham's π theorem, bend scour depth can be written in the following dimensionless form: ( ) In the present study, multiple regression analysis is then conducted with the measured field data listed in Table 1: which has a coefficient (R-squared value) of 0.924, indicating that the short-term bend scour depths computed by Equation ( 6) fit the measured data well, as shown in Figure 6.It should be noted that only 10 sets of data were used for calibration, while the remaining two sets of data were used for verification, as indicated in Figure 6.Bend scour depth was found to be directly proportional to unit discharge, channel slope, and water surface width, and inversely proportional to the centerline radius of the bend.

Method for Embankment Toe Scour Simulation
For an embankment located on a river bend, bend scour can lead to increased scour at the embankment toe.The physically based numerical models of bend scour simulation, as discussed in Section 1, present too much uncertainty due to the many parameters used in modeling.Therefore, a simple method with few parameters is most appropriate.Following the numerical framework proposed by Hong et al. [22], a new composite method based on the combination of flow, general, and bend scour computation equations is proposed herein to estimate the evolution of scour depth at an embankment toe on a river bend.Overall, the first step is to estimate the evolution of general scour depth by combining a 2D finite-volume model with a general scour computation equation.The next step is to modify the approach flow depth and velocity upstream of the embankment toe.The last step is to use the revised approach flow conditions for determining the bend scour depth evolution.The simulation procedure for embankment toe scour is presented in detail in Section 3.3.

2D Finite-Volume Hydraulic Model
The 2D model [24][25][26] is based on a finite-volume model in the unstructured arbitrarily shaped mesh system.The governing equations for 2D shallow water flow, including the continuity and momentum equations, can be expressed as: in which where Q is the vector of conserved variables; F I and G I are the inviscid flux vectors in the xand y-directions, respectively; F V and G V are the viscous flux vectors in the xand y-directions, respectively; S is the source term; h is the water depth; u and v are the depth-averaged velocity components in the xand y-directions, respectively; ρ is the density of flow; T xx , T xy , and T yy are the depth-averaged turbulent stresses; g is the gravitational acceleration; s 0x and s 0y are the bed slopes in the xand y-directions, respectively; and s fx and s fy are the friction slopes in the xand y-directions, respectively.The finite-volume model was adopted to solve the 2D governing equations, leading to the following basic vector equation: in which A is the area of the cell; m is the index that represents the side of the cell; M is the total number of the sides for the cell; T(θ) −1 is the inverse of the rotation matrix corresponding to the m side; θ is the angle between the outward unit vector n and the x-axis; n is the outward unit vector normal to the boundary of the control volume; L m is the length of the m side for the cell; represents the inviscid numerical flux; F V (Q) denotes the viscous numerical flux; and S is the integral form of the source terms.To resolve discontinuous shock waves or hydraulic jumps, many different types of approximate Riemann solvers have been developed for estimating the inviscid numerical flux F I (Q).This paper employs an upstream flux-splitting finite-volume scheme [24] to obtain the inviscid numerical flux.In addition, the Jacobian transformation method [25] is used to estimate the viscous numerical flux.
Water 2017, 9, 324 10 of 19 To avoid the generation of numerical errors due to unphysical high velocities at wet/dry fronts, the wet/dry boundary treatment of the model is improved herein.This paper introduces a positive tolerance depth h tol to deal with wet/dry boundary tracking: 1.
The cell interface is the dry edge, when h L ≤ h tol and h R ≤ h tol , in which h L and h R are the left and right water depths at the center of the cell interface, respectively.There is no flux estimation.

2.
The cell interface is the wet edge, when h L > h tol and h R > h tol .The upstream flux-splitting finite-volume scheme coupled with the hydrostatic reconstruction method [24] is adopted to estimate the well-balanced numerical fluxes.

3.
The cell interface is the partially wet edge with flux, when h L > h tol , h R ≤ h tol , and , where z b is the bed elevation.The momentum flux is set at 0 and the mass flux (hu) LR at cell interface LR is estimated as: (hu) LR = 1.42δη |δη|.

4.
The cell interface is the partially wet edge without flux, when h L > h tol , h R ≤ h tol , and According to the bed elevation condition, there is no flux across the cell interface.Thus, no flux estimation is required.
A detailed description of the model can be found in the literature [24][25][26].In practice, the presented 2D model adopts a well-balanced upwind numerical scheme coupled with a robust wet/dry algorithm.Hence, the model is a useful hydraulic modeling tool to simulate 2D hydraulic flows involving irregular bed topography.

General Scour Computation Equation
Many of the equations for scour estimation (general, contraction, or local scour) were developed according to the analysis of data from laboratory research.However, these equations derived from laboratory studies may not be suitable for field application of scouring.To overcome this problem, Hong et al. [22] proposed a general scour computation equation.This reliable equation is based on the results of general scour depth measurements through the numbered-brick method.The dimensionless equation can be expressed as: in which d gs denotes the short-term general scour depth.The accuracy of Equation ( 15) has been verified to give satisfactory results when estimating the short-term general scour depth in the Cho-Shui River [22].

Method for Simulating Bend Scour Depth Evolution
For sediment transport modeling, the type of soil is important, as discussed in Abou-Seida et al. [33].However, the study site is located in a gravel river with non-cohesive size classes, as shown in Figure 3. Therefore, the effect of cohesive sediment transport on bend scour is not considered herein.Consequently, the proposed Equations ( 6) and ( 15) are only suitable for general and bend scour estimation of non-cohesive sediment, respectively.
As shown in Figure 7, the algorithm for estimating bend scour depth evolution is summarized as follows: 1.
Inputting flow hydrographs (simulation time step is 1 h), setting model parameters (i.e., numerical time step and Manning roughness coefficients), and then simulating the 2D flow field in the bend reach through the proposed 2D finite-volume hydraulic model (Section 3.1).At the end of each simulation time step, the water depth and velocity in each computational cell are outputted.2.
After 2D flow simulation, one obtains the approach flow conditions (i.e., entire hydrographs) at the specific computational cell, including water depth, velocities, water surface width, and the centerline radius of the bend.

3.
For each simulation time step (1 h), one substitutes the approach flow conditions into the general scour computation equation (i.e., Equation ( 15)) to achieve the estimated general scour depth.
After the entire general scour simulation, the evolution of the estimated general scour depth is thus obtained.4.
For each simulation time step (1 h), one re-evaluates the water depth by adding the water depth and the estimated general scour depth to obtain the revised water depth: h = h + d gs . 5.
For each simulation time step (1 h), one re-evaluates velocity through unit flow discharge divided by revised water depth and obtains the revised velocity as v = q/ h. 6.
For each simulation time step (1 h), one substitutes the revised flow conditions (i.e., h u = h) into Equations ( 1)-( 3) and ( 6) to obtain the bend scour depth.After entire bend scour simulation, the evolution of the estimated bend scour depth is thus obtained.
1. Inputting flow hydrographs (simulation time step is 1 h), setting model parameters (i.e., numerical time step and Manning roughness coefficients), and then simulating the 2D flow field in the bend reach through the proposed 2D finite-volume hydraulic model (Section 3.1).At the end of each simulation time step, the water depth and velocity in each computational cell are outputted.2. After 2D flow simulation, one obtains the approach flow conditions (i.e., entire hydrographs) at the specific computational cell, including water depth, velocities, water surface width, and the centerline radius of the bend.3.For each simulation time step (1 h), one substitutes the approach flow conditions into the general scour computation equation (i.e., Equation ( 15)) to achieve the estimated general scour depth.
After the entire general scour simulation, the evolution of the estimated general scour depth is thus obtained.4. For each simulation time step (1 h), one re-evaluates the water depth by adding the water depth and the estimated general scour depth to obtain the revised water depth: For each simulation time step (1 h), one re-evaluates velocity through unit flow discharge divided by revised water depth and obtains the revised velocity as / v q h =   .6.For each simulation time step (1 h), one substitutes the revised flow conditions (i.e., u h h =  ) into Equations ( 1)-( 3) and ( 6) to obtain the bend scour depth.After entire bend scour simulation, the evolution of the estimated bend scour depth is thus obtained.

General scour Bend scour
Figure 7.The algorithm of the proposed method for estimating the evolution of bend scour depth.
The output variable of Equations ( 1)-( 3) is the maximum flow depth (i.e., bs h ), whereas the output estimated value of Equation ( 6) is the bend scour depth (i.e., bs d ).In this paper, the bend scour depth is selected to be a common basis for comparison.As shown in Figure 5, the bend scour depth equals maximum flow depth minus revised water depth.Based on this relationship ( bs bs d h h = −  ), the maximum flow depths estimated by the three equations were transferred to the bend scour depths.Then, the transferred bend scour depths were compared with the bend scour depths estimated by the proposed method.Consequently, the results discussed in Section 4.2 are on the basis of the same physical variables (i.e., computed and measured bend scour depths).The output variable of Equations ( 1)-( 3) is the maximum flow depth (i.e., h bs ), whereas the output estimated value of Equation ( 6) is the bend scour depth (i.e., d bs ).In this paper, the bend scour depth is selected to be a common basis for comparison.As shown in Figure 5, the bend scour depth equals maximum flow depth minus revised water depth.Based on this relationship (d bs = h bs − h), the maximum flow depths estimated by the three equations were transferred to the bend scour depths.Then, the transferred bend scour depths were compared with the bend scour depths estimated by the proposed method.Consequently, the results discussed in Section 4.2 are on the basis of the same physical variables (i.e., computed and measured bend scour depths).

Verifications of 2D Finite-Volume Model
This section verifies the accuracy of the finite-volume model for 2D flow field simulations.Figure 8 shows the simulation domain of the study channel bend, namely a cross-section and LiDAR (0.5 m × 0.5 m) bed elevation data.The study reach was approximately 12 km long between cross section 117 (the upstream inflow boundary) and cross section 90 (the downstream outflow boundary).The field-measured water level data at the Mingchu gauging station was used for Manning roughness coefficient calibration and model verification.Toe scouring at the Shuideliaw Embankment usually results from heavy rainfall during typhoon events, and two such flood events are analyzed herein: Typhoon Trami (August 2013) and Typhoon Usagi (September 2013).
The Manning roughness coefficient has a significant effect on flow and scour modeling.For the commonly used approach, the Manning roughness coefficient was initially given in the entire simulation domain and should be calibrated based on a reasonable range.According to Chow [34], the range of Manning roughness coefficient in a natural river is from 0.025 to 0.055.Then, a reasonable value must be initially given.In this paper, the Manning roughness coefficient was initially estimated following the method in [35].As described in [35], an averaged value of the Manning roughness coefficients estimated from three formulas (Lane formula by n = 0.015d 75 1/6 , Einstein formula by n = 0.0132d 65 1/6 , Strickler formula by n = 0.015d m 1/6 ) is suggested for use in a gravel-bed river.The unit of sediment particle size in these three formulas is millimeters.Accordingly, using these three formulas one obtains 0.038, 0.032, and 0.034, respectively.Then, an average value of 0.035 was achieved and was initially given in the entire simulation domain.For Typhoon Trami, the overestimated error of water-level hydrograph was obtained using the value of 0.035.Thus, a further adjustment of value of Manning roughness coefficient is needed.The purpose of the adjustment is to make a reasonable match between the simulated and measured water-level hydrographs.After several adjustments, the best choice for calibration of Manning roughness coefficient was finally determined at a value of 0.025 and this value is used to simulate 2D flow field for Typhoon Usagi.

Verifications of 2D Finite-Volume Model
This section verifies the accuracy of the finite-volume model for 2D flow field simulations.Figure 8 shows the simulation domain of the study channel bend, namely a cross-section and LiDAR (0.5 m × 0.5 m) bed elevation data.The study reach was approximately 12 km long between cross section 117 (the upstream inflow boundary) and cross section 90 (the downstream outflow boundary).The field-measured water level data at the Mingchu gauging station was used for Manning roughness coefficient calibration and model verification.Toe scouring at the Shuideliaw Embankment usually results from heavy rainfall during typhoon events, and two such flood events are analyzed herein: Typhoon Trami (August 2013) and Typhoon Usagi (September 2013).
The Manning roughness coefficient has a significant effect on flow and scour modeling.For the commonly used approach, the Manning roughness coefficient was initially given in the entire simulation domain and should be calibrated based on a reasonable range.According to Chow [34], the range of Manning roughness coefficient in a natural river is from 0.025 to 0.055.Then, a reasonable value must be initially given.In this paper, the Manning roughness coefficient was initially estimated following the method in [35].As described in [35], an averaged value of the Manning roughness coefficients estimated from three formulas (Lane formula by ) is suggested for use in a gravel-bed river.The unit of sediment particle size in these three formulas is millimeters.Accordingly, using these three formulas one obtains 0.038, 0.032, and 0.034, respectively.Then, an average value of 0.035 was achieved and was initially given in the entire simulation domain.For Typhoon Trami, the overestimated error of water-level hydrograph was obtained using the value of 0.035.Thus, a further adjustment of value of Manning roughness coefficient is needed.The purpose of the adjustment is to make a reasonable match between the simulated and measured water-level hydrographs.After several adjustments, the best choice for calibration of Manning roughness coefficient was finally determined at a value of 0.025 and this value is used to simulate 2D flow field for Typhoon Usagi.The flow discharge hydrographs released from Chi-Chi Weir, as shown in Figure 4, were used for the subcritical inflow boundary conditions.The left and right sides along the simulation domain were considered the land boundaries and the time step was 0.5 s.The numerical accuracy of the 2D model was evaluated by comparing the simulated results with the measured water levels using the following three criteria: The flow discharge hydrographs released from Chi-Chi Weir, as shown in Figure 4, were used for the subcritical inflow boundary conditions.The left and right sides along the simulation domain were considered the land boundaries and the time step was 0.5 s.The numerical accuracy of the 2D model was evaluated by comparing the simulated results with the measured water levels using the following three criteria: in which Eη p is the peak water level error; ET p represents the error of time to peak water level; R 2 denotes the coefficient of determination; η sim p and η mea p denote the simulated and measured peak water levels, respectively; T sim p and T mea p are the simulated and measured time to peak water level, respectively; η sim i and η mea i are the simulated and measured water levels at each time; and η sim and η mea denote the average simulated and measured water levels.The results for the two examined flood events, for which the Eη p results were 0.13% and 0.46%, respectively, are provided in Table 2.Because the overall performance of Eη p is less than 0.5%, the results demonstrate that the model can provide high numerical accuracy for the peak water level in modeling 2D flow with variable bed topography.In addition, the average value of R 2 coefficient for two flood events is 0.88, close to the best value of 1.0, indicating that the proposed 2D model can achieve reasonable performance for water-level hydrograph.Figure 9 depicts the simulated and the measured water levels at Mingchu gauging station, which are compared by calibrating the Manning roughness coefficients (0.025).The simulated results reveal close agreement with the measured data at the peak flow.For the periods before and after the peak flow, the model overestimates the water level.The error between the simulated and measured data may be due to the adopted assumptions of a "fixed" bed and the constant given value of Manning roughness coefficient.Careful attention should be paid to the resolution of the entire bend scour depth evolution, and the relevant uncertainty of Manning roughness coefficient should be considered in future studies.In the present study, these assumptions do not affect the maximum bend scour depths based on the acceptable results listed in Table 3.In addition, Figure 10 displays the simulated results for velocity contours under peak flood conditions.Reasonable results are obtained using the 2D model, indicating that the velocities adjacent to the outer Shuideliaw Embankment increase significantly with peak discharge.Due to the river bend effects, high-velocity floods can attack and erode the riverbank, meaning that the outer embankment is the most vulnerable area.The results also show that the proposed 2D model is accurate for modeling flood flow with the irregular bed topography and complex geometry in a channel bend.

Verifications of Proposed Method for Bend Scour Depth
The results presented in Section 4.1 demonstrate that the proposed finite-volume hydraulic model determines accurate approach flow conditions.In this section, the evolution simulations of the bend scour depth obtained through the proposed approach are presented.
The approach flow conditions for the estimation of bend scour depth include the water depth, velocity, water surface width, and centerline radius of the bend.The proposed finite-volume hydraulic model was adopted to compute the required approach flow conditions for two typhoon events.The evolution of bend scour depth was then computed using the proposed method as well as the three existing equations (Galay et al. [27], Thorne [28], and U.S. Army Corps of Engineering [29]) identified earlier.Comparisons between the measured data and the computed results are displayed in Figure 11.Globally, all simulated scour depth variations yield the proper rising and falling limbs of the hydrograph, which are similar to those displaying in flow hydrographs, as shown in Figure 4. To further discuss the differences among the four methods, the relative error is used and defined as: where d sim bs and d mea bs represent the simulated and measured maximum bend scour depth, respectively.Table 3 compares the simulated results with the measured data, including the maximum bend scour depth and the relative error.The results indicate that the existing empirical equations by Thorne and the U.S. Army Corps of Engineering tend to overestimate scour depth.However, the proposed method herein provides a good fit compared with the measured maximum bend scour depth, and is thus the most accurate of all the tested methods.

Practical Embankment Safety Curve
A rational and practical tool to determine embankment safety in real time is of great importance for decisions made during typhoon periods.Based on the proposed method, this study attempts to provide a quantitative curve to evaluate the influence of discharge on bend scour depth at the embankment toe.The numerical simulations are performed with nine different discharges of 1.11-, 2-, 5-, 10-, 20-, 25-, 50-, 100-, and 200-year return period floods.According to the results described in the previous section, the proposed bend scour computation equation was verified to be suitable for simulating embankment toe scour in this case study.Thus, the bend scour depth at the peak flow condition obtained from the proposed equation is utilized to develop the embankment safety curve.Figure 12a shows the results of numerical experiments, indicating that the bend scour depth increases significantly with an increase in the discharge.Using statistical regression analysis results in the following regression equation: d bs = 0.0005Q.On the basis of this equation, the bend scour depth at the embankment toe downstream of the Chi-Chi Weir in Cho-Shui River can be easily estimated in real time for a given flow discharge.In addition, the relationship between the scoured bed level and the corresponding flow discharge is added to Figure 12a.Moreover, the location of Shuideliaw Embankment and adjacent bed change at cross section 105 are presented in Figure 12b.Figure 12b indicates that the toe of embankment foundation is 139.0 m and the bed level before the flood at Shuideliaw Embankment is 140.5 m.During typhoon periods, the bed level is eroded due to bend scour.Based on the proposed embankment safety curve, river managers can easily know whether the scoured bed level is smaller than the toe of the embankment foundation or not.After the embankment safety curve is established, its accuracy must be verified.From a process standpoint, the safety of an embankment is controlled by the relationship between the resisting and driving forces of the embankment.For reinforced concrete embankments, the driving forces can represent the hydraulic scouring around the embankment toe, and the resisting force can stand for the depth of the embankment foundation.During the flood induced by a monsoon in June 2012, the bend scour depth near Shuideliaw Embankment was larger than the depth of the embankment foundation.Hence, the embankment failed due to the embankment toe erosion caused by flooding with a peak discharge of 6122 m 3 /s.By using the proposed equation 0.0005 bs d Q = , when discharge is equal to 6122 m 3 /s, the total bend scour depth is estimated at 3.06 m, which is larger than the depth of the embankment foundation (3.0 m).Therefore, the results demonstrate that the proposed simplified equation could be used to correctly assess embankment toe safety.To further provide advantageous assessment information, the discharge warning value can be achieved from this equation fairly simply.Moreover, the depth of the embankment foundation can be substituted into this equation, resulting in 6000 m 3 /s waring discharge.The results herein suggest that the study embankment will be unstable and may further fail if discharge from After the embankment safety curve is established, its accuracy must be verified.From a process standpoint, the safety of an embankment is controlled by the relationship between the resisting and driving forces of the embankment.For reinforced concrete embankments, the driving forces can represent the hydraulic scouring around the embankment toe, and the resisting force can stand for the depth of the embankment foundation.During the flood induced by a monsoon in June 2012, the bend scour depth near Shuideliaw Embankment was larger than the depth of the embankment foundation.Hence, the embankment failed due to the embankment toe erosion caused by flooding with a peak discharge of 6122 m 3 /s.By using the proposed equation d bs = 0.0005Q, when discharge is equal to 6122 m 3 /s, the total bend scour depth is estimated at 3.06 m, which is larger than the depth of the embankment foundation (3.0 m).Therefore, the results demonstrate that the proposed simplified equation could be used to correctly assess embankment toe safety.To further provide advantageous assessment information, the discharge warning value can be achieved from this equation fairly simply.Moreover, the depth of the embankment foundation can be substituted into this equation, resulting in 6000 m 3 /s waring discharge.The results herein suggest that the study embankment will be unstable and may further fail if discharge from the Chi-Chi Weir is greater than 6000 m 3 /s.By contrast, the numerical experiments demonstrated that the maximum bend scour depth near the Shuideliaw Embankment is greater than the designed foundation depth, because the flow discharge is greater than 6000 m 3 /s (which is between the flood discharges for the two-and five-year return periods, similar to the bankfull discharge).This similarity between the flow and bankfull discharge denotes a significant change in the river bed, and indicates why embankment failure usually occurs during bankfull discharge.
Water 2017, 9, 324 17 of 19 the Chi-Chi Weir is greater than 6000 m 3 /s.By contrast, the numerical experiments demonstrated that the maximum bend scour depth near the Shuideliaw Embankment is greater than the designed foundation depth, because the flow discharge is greater than 6000 m 3 /s (which is between the flood discharges for the two-and five-year return periods, similar to the bankfull discharge).This similarity between the flow and bankfull discharge denotes a significant change in the river bed, and indicates why embankment failure usually occurs during bankfull discharge.In summary, the proposed embankment safety curve provides a quick and advantageous assessment for river managers to use when making management decisions associated with embankment stability or safety warning.

Conclusions
Scour at the embankment toe is a complex fluvial process, especially in bend reaches.A new composite method for simulating general and bend scour depths was therefore proposed in this paper.This composite method consists of three components: a 2D finite-volume model, a general scour computation equation, and a new bend scour computation equation.An assumption linking the flow and bend scour evolutions has been presented, and it allows the proposed method to simulate both general and bend scour evolutions in a river bend with complex bed topography under typhoon-induced flood conditions.The 2D flow fields and the temporal evolutions of bend scour depths in a bend reach of the Cho-Shui River were simulated in detail for two typhoon flood events, and the simulated results were compared to field data measured using the numbered-brick method.
The primary conclusions of this study are as follows: (1) the new practical method is capable of simulating the temporal evolution process caused by the development of bend scour at  In summary, the proposed embankment safety curve provides a quick and advantageous assessment for river managers to use when making management decisions associated with embankment stability or safety warning.

Conclusions
Scour at the embankment toe is a complex fluvial process, especially in bend reaches.A new composite method for simulating general and bend scour depths was therefore proposed in this paper.This composite method consists of three components: a 2D finite-volume model, a general scour computation equation, and a new bend scour computation equation.An assumption linking the flow and bend scour evolutions has been presented, and it allows the proposed method to simulate both general and bend scour evolutions in a river bend with complex bed topography under typhoon-induced flood conditions.The 2D flow fields and the temporal evolutions of bend scour depths in a bend reach of the Cho-Shui River were simulated in detail for two typhoon flood events, and the simulated results were compared to field data measured using the numbered-brick method.
The primary conclusions of this study are as follows: (1) the new practical method is capable of simulating the temporal evolution process caused by the development of bend scour at embankment toes; (2) the simulated results are generally in agreement with the field-measured data, including water level data at the Mingchu gauging station and maximum bend scour depths at Shuideliaw Embankment; (3) the proposed embankment safety curve enables river managers to make rational decisions when responding to embankment emergencies during the typhoon season.
The new method has been demonstrated to accurately simulate bend scour evolution based on the separating assumption of scour components.However, further testing must be conducted in comparison with more field-measured evolution of embankment toes, specifically considering the uncertainty of the sediment deposition process and the Manning roughness coefficient.In addition, our method can be integrated into a flood warning system that can be used for automated real-time forecasts of water levels and embankment toe scouring.The system provides information for effective and timely decision-making on hazard mitigation action prior to expected critical situations.

Figure 1 .
Figure 1.Topographic map of study site in the Cho-Shui River basin.Figure 1. Topographic map of study site in the Cho-Shui River basin.

Figure 1 .
Figure 1.Topographic map of study site in the Cho-Shui River basin.Figure 1. Topographic map of study site in the Cho-Shui River basin.

Figure 2 .
Figure 2. Aerial view of the study site before (a) and after (b) Typhoon Saoloa.The location of the scour measurement is also illustrated (c).

Figure 2 .
Figure 2. Aerial view of the study site before (a) and after (b) Typhoon Saoloa.The location of the scour measurement is also illustrated (c).

Figure 3 .
Figure 3. Size distribution of sediment particles at the study site where the survey was conducted in 2013.

Figure 3 .
Figure 3. Size distribution of sediment particles at the study site where the survey was conducted in 2013.

Figure 3 .
Figure 3. Size distribution of sediment particles at the study site where the survey was conducted in 2013.

Figure 7 .
Figure 7.The algorithm of the proposed method for estimating the evolution of bend scour depth.

Figure 8 .
Figure 8.The simulation domain in a gravel-bed bend river reach of the Cho-Shui River

Figure 8 .
Figure 8.The simulation domain in a gravel-bed bend river reach of the Cho-Shui River

Figure 9 .Figure 10 .
Figure 9. Comparisons of simulated and measured water level hydrographs for (a) Typhoon Trami and (b) Typhoon Usagi.

Figure 11 .
Figure 11.The simulated time variations of bend scour depths using four approaches for (a) Typhoon Trami and (b) Typhoon Usagi.

Figure 11 .
Figure 11.The simulated time variations of bend scour depths using four approaches for (a) Typhoon Trami and (b) Typhoon Usagi.

Figure 12 .
Figure 12.The relationship between the bend scoured depth, scoured bed level, and flow discharge at Chi-Chi Weir (a).The location of Shuideliaw Embankment and adjacent bed change at cross section 105 are also illustrated (b).

Figure 12 .
Figure 12.The relationship between the bend scoured depth, scoured bed level, and flow discharge at Chi-Chi Weir (a).The location of Shuideliaw Embankment and adjacent bed change at cross section 105 are also illustrated (b).

Table 1 .
Field data collection for short-term bend scour induced by floods.

Table 2 .
Simulated results for two flood events using three criteria.

Table 3 .
The simulated maximum bend scour depths and the relative errors by four approaches for Typhoon Trami and Typhoon Usagi.