Experimental Investigation of Flood Energy Dissipation by Single and Hybrid Defense System

In this study, a series of laboratory experiments were conducted to investigate the energy loss through the hybrid defense system (HDS) in the order of dike, moat, and emergent vegetation in steady subcritical flow conditions. The results of HDS were compared with a single defense system (SDS) comprising only vegetation (OV). The dimensions of dike were kept constant while two different shapes (trapezoidal and rectangular) of moat were considered. The impacts of vegetation of variable thickness and density were investigated. Two combinations of HDS were investigated including the combination of dike and vegetation (DV) and the combination of dike, moat, and vegetation (DMV). The effect of backwater rise due to the vegetation, hydraulic jump formation and the impact of the arrival time of floodwater on energy dissipation were investigated. It was observed that on the upstream side of obstructions, the backwater depth increased by increasing the Froude number in both the SDS and HDS. The hydraulic jump observed in HDS was classified and the energy dissipation due to it was calculated. Under various conditions investigated in this paper, the maximum average energy dissipation was 32% in SDS and 46% in HDS. The trapezoidal moat performed better than rectangular moat as energy dissipater. The delay time was also greater with trapezoidal moat as compared to that in rectangular one. The maximum delay time was 140 s in the case of HDS. Hence, the hybrid defense system offered maximum resistance to the flow of water, thus causing a significant energy loss. For each case of SDS and HDS, empirical equations were developed by regression analysis to estimate the energy dissipation amounts.


Introduction
Heavy floods carry serious socio-economic threats. Annually, floods claim almost 20,000 lives and adversely affect at least 20 million people throughout the world [1]. The recent studies indicate that losses from flood hazards are expected to increase in the near future [2]. In various European countries, whenever there is flood in main rivers, the surrounding small rivers are also flooded at the same time creating heavy damage to properties and human lives [3]. In 2011, the Great East Tsunami in Japan arrived with a return period of 1000 years [4] which destroyed embankments, coastal forests, tsunami gates, and thousands of human lives [5][6][7]. Similarly, the other south Asian countries were also affected by heavy floods in the recent past. Pakistan faces extreme flood events almost every year. The statistical analysis shows that a financial loss of US$38.165 billion occurred, more than 12,000 causalities were reported and 616,598 km 2 land area was affected as a result of 23 major flood events from 1947 to 2015 in Pakistan [8]. The extreme flood events in India resulted in the death of several thousand people in Uttarakhand 2013, in Chennai 2015, and in Kerala 2018 [9,10]. Almost 38 districts were affected by a Water 2019, 11,1971 2 of 24 continuous spell of heavy rainfall in Bangladesh in 2017 [11]. Similarly, Indonesia is also vulnerable to flood damages [12,13].
In the last decade, many researchers investigated the energy mitigation and flow phenomena over embankments using physical models, numerical models and analytical approaches [14][15][16][17]. These methods may be costly for developing countries because of significant capital investment requirements [18]. Currently, ecological, natural methods are widely considered in many developed and underdeveloped countries [19][20][21][22][23]. Ecosystem degradation is a major cause of increasing water resources management challenges [24][25][26][27][28]. In the floodplain of rivers, vegetation may increase the hydraulic resistance by dragging and turbulence to safeguard the ecosystem [26]. Consequently, increasing tree cover is promoted as an effective means to enhance the storage capacity of catchment areas to increase the lag time, reduce discharge and flood risk [27].
In addition to the above-mentioned solutions, the hybrid or compound defense system were also studied by different researchers. Various studies consist of the energy dissipation through a compound defense system of vegetation [28], combined effects of sea embankment and coastal forest [4,29], a moat with/without a forest [30] and a double embankment system [31]. Furthermore, Zaha (2019) conducted a series of experiments to investigate a hybrid defense system consisting of a forest, a moat (depression) and an embankment under unsteady flow conditions, considering the dike as unbreakable [32]. It was noticed that in the recent catastrophic floods in various regions, soon after the breaching of embankments, the floodwater expanded into the entire floodplains within a short duration, but the barriers in the path of water caused delay in arrival of water to important places, providing a chance for adopting the preventing measures [32,33]. The field survey of different reaches of Indus basin revealed that the cross boundaries (embankments) are unprotected due to lack of proper maintenance and low density of coastal forests (Figure 1a-c). In order to improve their performance, there is a dire need to repair existing damaged embankments and construct a new sustainable system of flood protection [33].
However, the energy dissipation phenomena of a hybrid defense system and assessment of the effect of various shapes of its components still need further investigation [4,32]. Due to the hybrid defense system, the flow structure becomes more complex because of the formation of a hydraulic jump. The jump characteristics and its location should be identified to improve the efficiency of the hybrid defense system against the destructive forces [4]. In the present paper, the hydraulic jump developed in SDS and HDS was classified on the basis of its location and the Froude number. The energy dissipation due to the formation of a hydraulic jump in various cases was also calculated. In HDS, the effect of the density of vegetation and the shape of moat on total energy dissipation was analyzed and equations were developed by regression analysis to calculate energy dissipation, which has hardly been reported previously.
Therefore, the main objective of this study was to investigate the energy dissipation phenomena in the case of a single defense system and a hybrid defense system. To achieve the objective of the study, a series of experiments were conducted using laboratory models. The energy dissipation effects of each system were compared by paying special attention to the backwater rise on the upstream side of models, the formation of a hydraulic jump and the arrival time of floodwater. This study aimed to provide findings that will be useful in designing a hybrid defense system for the flood energy dissipation. The selected range of Froude numbers was based on previous data of different locations on the Indus basin, however, the results can be used worldwide for similar situations. The proposed defense system will be effective for the protection of various floodplains. The floodwater can also be stored in moat which may act as a floodwater harvesting structure. The imagery of proposed HDS at River Indus near Taunsa Barrage, Pakistan is shown in Figure 1d.

Flume Characteristics
A glass-sided open channel flume of 10-m length, 0.31-m width and 0.5-m height with a constant bed slope of 1/500 was used for the laboratory experiments with different conditions, as shown in Table 1. The field data of past 59 years regarding Froude number for Jinnah barrage (on the upstream side) and Taunsa barrage (on the downstream side) located on the Indus river were collected from the Punjab Irrigation Department, Punjab, Pakistan. These data were used to select the experimental range of initial Froude number (Fr o = V/(gh) 0.5 , where V = depth-averaged velocity (m/s), g = gravitational acceleration (m/s 2 ), and h = water depth (m)). The initial Froude number (Fr o ) was estimated with the velocity and water depth without any obstruction placed in the channel. The calculated range of Froude number at Jinnah barrage was 0.17-0. 59

Flume Characteristics
A glass-sided open channel flume of 10-m length, 0.31-m width and 0.5-m height with a constant bed slope of 1/500 was used for the laboratory experiments with different conditions, as shown in Table 1. The field data of past 59 years regarding Froude number for Jinnah barrage (on the upstream side) and Taunsa barrage (on the downstream side) located on the Indus river were collected from the Punjab Irrigation Department, Punjab, Pakistan. These data were used to select the experimental range of initial Froude number (Fro = V/(gh) 0.5 , where V = depth-averaged velocity (m/s), g = gravitational acceleration (m/s 2 ), and h = water depth (m)). The initial Froude number (Fro) was estimated with the velocity and water depth without any obstruction placed in the channel. The calculated range of Froude number at Jinnah barrage was 0.17-0.59 and at Taunsa barrage was 0.10-0. 59 Table 1. The field data of past 59 years regarding Froude number for Jinnah barrage (on the upstream side) and Taunsa barrage (on the downstream side) located on the Indus river were collected from the Punjab Irrigation Department, Punjab, Pakistan. These data were used to select the experimental range of initial Froude number (Fro = V/(gh) 0.5 , where V = depth-averaged velocity (m/s), g = gravitational acceleration (m/s 2 ), and h = water depth (m)). The initial Froude number (Fro) was estimated with the velocity and water depth without any obstruction placed in the channel. The calculated range of Froude number at Jinnah barrage was 0.17-0.59 and at Taunsa barrage was 0.10-0. 59 Table 1. The field data of past 59 years regarding Froude number for Jinnah barrage (on the upstream side) and Taunsa barrage (on the downstream side) located on the Indus river were collected from the Punjab Irrigation Department, Punjab, Pakistan. These data were used to select the experimental range of initial Froude number (Fro = V/(gh) 0.5 , where V = depth-averaged velocity (m/s), g = gravitational acceleration (m/s 2 ), and h = water depth (m)). The initial Froude number (Fro) was estimated with the velocity and water depth without any obstruction placed in the channel. The calculated range of Froude number at Jinnah barrage was 0.17-0.59 and at Taunsa barrage was 0.10-0. 59 Table 1. The field data of past 59 years regarding Froude number for Jinnah barrage (on the upstream side) and Taunsa barrage (on the downstream side) located on the Indus river were collected from the Punjab Irrigation Department, Punjab, Pakistan. These data were used to select the experimental range of initial Froude number (Fro = V/(gh) 0.5 , where V = depth-averaged velocity (m/s), g = gravitational acceleration (m/s 2 ), and h = water depth (m)). The initial Froude number (Fro) was estimated with the velocity and water depth without any obstruction placed in the channel. The

Experimental Conditions
The tree species selected for the vegetation model was the Eucalyptus species. Its tree crown was higher relative to the maximum flood level in the past, and thus, the trees can be considered as circular cylinders [29]. The average tree height in the ecological zones of Punjab, Pakistan is in the range of 7.6-14.6 m and the average trunk diameter is 0.11-0.33 m [34]. To ensure the geometric similarity, metal rods, and steel cylinders of different diameters have often been used by various researchers [35][36][37]. In this experimental work, keeping in view the above-mentioned limits of average tree diameter, metal rods of 0.003 m diameter at a scale of 1:100 were used because the wooden cylinders of such a small diameter were expected to break or bend during the experimentation. In SDS, only vegetation (OV) models with variable densities were separately placed in the flume to check the individual capability of energy dissipation for these models. W v d × 10 2 . The dn represents the cumulative diameter of trees at breast height and is defined as the product of tree diameter and the number of trees in a rectangle with a frontage of unit length along the shoreline and depth equal to the width of the forest (Wv) [37].
In the current study, dn = 180 No.cm was selected for both the intermediate and sparse vegetation cases. The water level was measured by using a point gauge at the interval of 2-3 cm. Without placing any obstruction in the flume, the discharge was measured using a flow meter against each selected

Experimental Conditions
The tree species selected for the vegetation model was the Eucalyptus species. Its tree crown was higher relative to the maximum flood level in the past, and thus, the trees can be considered as circular cylinders [29]. The average tree height in the ecological zones of Punjab, Pakistan is in the range of 7.6-14.6 m and the average trunk diameter is 0.11-0.33 m [34]. To ensure the geometric similarity, metal rods, and steel cylinders of different diameters have often been used by various researchers [35][36][37]. In this experimental work, keeping in view the above-mentioned limits of average tree diameter, metal rods of 0.003 m diameter at a scale of 1:100 were used because the wooden cylinders of such a small diameter were expected to break or bend during the experimentation. In SDS, only vegetation (OV) models with variable densities were separately placed in the flume to check the individual capability of energy dissipation for these models. W v d × 10 2 . The d n represents the cumulative diameter of trees at breast height and is defined as the product of tree diameter and the number of trees in a rectangle with a frontage of unit length along the shoreline and depth equal to the width of the forest (W v ) [37].
In the current study, d n = 180 No.cm was selected for both the intermediate and sparse vegetation cases. The water level was measured by using a point gauge at the interval of 2-3 cm. Without placing any obstruction in the flume, the discharge was measured using a flow meter against each selected initial water depth (h o ). The depth-averaged velocity was computed using the relationship between measured water depth and discharge. The dike and the moat were modeled at a scale of 1/100. The recommended height of embankment (dike) on Indus river in Pakistan is in the range of about 4-7 m high with addition to 1.2-1.8 m freeboard [38]. In the current study, the selected height of embankment (including the freeboard) was 5.5 m. Therefore, by using the scale of 1/100, the height of the embankment model was kept 5.5 cm. The dimensions of dike and moat are already mentioned in Table 1. In case of flooding, water flows from the embankment towards the floodplain in an oblique direction with respect to flow in the main channel. However, the direction of flow of water towards the weir-like obstacles/embankments in the floodplains of the river can be encountered in various orientations i.e., oblique or perpendicular depending upon the angle of placement of the obstruction [39][40][41][42]. In this study, it is assumed that the obstruction in the field will be constructed in such a way that the direction of flow of water towards the obstruction will be perpendicular. It was also considered that the slope of the floodplain was negligible provided that the banks of the river were already limited by defense work. However, it is worth mentioning here that if the obstruction is placed as per usual cases then the flow will approach the obstruction in an oblique manner and the results of this study will have to be reduced by a reduction factor depending upon the angle of obliqueness of flow. This angle will be different for different situations, and the reduction factor developed by [42][43][44] may be used to modify our results for real-life designs.
As described above, two cases of the flood energy dissipation defense system were tested. These two cases consist of a single defense system (SDS) and the hybrid defense system (HDS) as elaborated in the schematic initial water depth (ho). The depth-averaged velocity was computed using the relationship between measured water depth and discharge. The dike and the moat were modeled at a scale of 1/100. The recommended height of embankment (dike) on Indus river in Pakistan is in the range of about 4-7 m high with addition to 1.2-1.8 m freeboard [38]. In the current study, the selected height of embankment (including the freeboard) was 5.5 m. Therefore, by using the scale of 1/100, the height of the embankment model was kept 5.5 cm. The dimensions of dike and moat are already mentioned in Table 1. In case of flooding, water flows from the embankment towards the floodplain in an oblique direction with respect to flow in the main channel. However, the direction of flow of water towards the weir-like obstacles/embankments in the floodplains of the river can be encountered in various orientations i.e., oblique or perpendicular depending upon the angle of placement of the obstruction [39][40][41][42]. In this study, it is assumed that the obstruction in the field will be constructed in such a way that the direction of flow of water towards the obstruction will be perpendicular. It was also considered that the slope of the floodplain was negligible provided that the banks of the river were already limited by defense work. However, it is worth mentioning here that if the obstruction is placed as per usual cases then the flow will approach the obstruction in an oblique manner and the results of this study will have to be reduced by a reduction factor depending upon the angle of obliqueness of flow. This angle will be different for different situations, and the reduction factor developed by [42][43][44] may be used to modify our results for real-life designs.
As described above, two cases of the flood energy dissipation defense system were tested. These two cases consist of a single defense system (SDS) and the hybrid defense system (HDS) as elaborated in the schematic

Non-Dimensional Pi Groups
The scaling of experimental models to the prototype conditions is a challenging task [39]. Figure  4 shows a flow structure scheme and important parameters of the proposed HDS. Flow conditions, type of structure, vegetation arrangement and their combination are the main parameters that decide the amount of energy lost. The following dimensionless groups were developed by using well known Buckingham's pi theorem.
where Δh1 = backwater rise at the upstream of dike, ho = initial water depth without any model,  Figure 3. Schematic diagram of experimental model.

Non-Dimensional Pi Groups
The scaling of experimental models to the prototype conditions is a challenging task [39]. Figure 4 shows a flow structure scheme and important parameters of the proposed HDS. Flow conditions, type of structure, vegetation arrangement and their combination are the main parameters that decide the amount of energy lost. The following dimensionless groups were developed by using well known Buckingham's pi theorem. To minimize the sidewall effects on the flow structure the ratio between the flume width and the cylinder diameter should be greater than five [4,45]. In this study, this ratio was kept greater than 5 to minimize the sidewall effect. The backwater rise, i.e., ∆h 1 /h o , ∆h 2 /h o and total energy loss (∆E 1 /E 1 ) and (∆E 2 /E 2 ) are the function of the initial Froude number Fr o , vegetation density (B/d) and d n . Froude number is commonly used in free-surface flows, so the Reynolds number was not considered here.
The scaling of experimental models to the prototype conditions is a challenging task [39]. Figure  4 shows a flow structure scheme and important parameters of the proposed HDS. Flow conditions, type of structure, vegetation arrangement and their combination are the main parameters that decide the amount of energy lost. The following dimensionless groups were developed by using well known Buckingham's pi theorem.
where Δh1 = backwater rise at the upstream of dike, ho = initial water depth without any model,  Figure 4. Experimental model of the hybrid defense system (HDS).

Description of Energy Dissipation
The effectiveness of the hybrid defense structure was assessed by calculating the energy reduction (∆E) in two ways. One was the reduction due to the formation of a hydraulic jump within the hybrid defense system, and the other was the total energy dissipation by HDS. Energy reduction due to the hydraulic jump in Sections 2 and 3 are represented as (∆E j2 /y 1 ) and (∆E j3 /y 1 ), where, , y 1 is the minimum depth at the toe of the hydraulic jump, y 2 is the average water depth after the hydraulic jump [43]. The specific energy is defined as the energy per unit volume of water at any section of a channel with respect to the channel bed [46,47]. Therefore, where E = specific energy, y = water depth, v = velocity, and α = coefficient to account for variations in velocity, in this study a rectangular channel with a uniform cross-section and a constant discharge was taken, so the value of α is assumed to be unity as considered in the previous research [4]. The depth-averaged velocity against the known discharge (Q) was calculated using the known relationship for steady flow, i.e., V = Q/A. The relative percentage energy loss was calculated by using Equation (2) for DV and DMV models and Equation (3) for OV models.
where, ∆E 1 is the relative energy loss, E 1 = energy on the upstream side of dike in DV and DMV models, E 2 = energy on the upstream side of vegetation and E 3 = energy on the downstream side of vegetation.

Delay Time Analysis
In this study, the delay time was considered as the time required for water to reach from the inlet of 10 m flume to the crest of the crump weir (Figure 7d) installed at 800 cm from the inlet having length = 70 cm, width = 31 cm (throughout the channel width), height = 15cm, 1:5 slope on the upstream and 1:2 on the downstream side of crump weir. The time (T WM ) required to reach the water from inlet to the crest level of the crump weir was noted without placing any model in the channel and compared with the time required in SDS and HDS models.

Backwater Rise
The vegetation always provides a barrier to the flowing water which creates backwater on the upstream side [30,32,39]. In the case of OVS and OVI, the initial water depth (h o ) was raised upstream of the vegetation model, and the water surface slope was increased inside the vegetation due to the resistance offered by it. It was found that by increasing the density of vegetation from sparse to intermediate, the backwater rise increases to 26% and the surface slope inside the vegetation was also increased. The backwater rise and relative backwater rise varied linearly with increasing Froude number as shown in Figure 5a,b. However, the relative backwater rise was only slightly increased for different Froude number values as the range of the values on y-axis was lower as shown in Figure 5b. A similar trend was observed in the literature [48]. This shows that the initial flow depth had very little effect on the relative backwater rise for a given Froude number. In these figures, ∆h and h o are the water depths on the upstream side of vegetation and initial water depth without any model placed in the channel respectively.

Hydraulic Jump and Water Surface Profile Classification
A hydraulic jump is formed when supercritical flow transforms to subcritical flow [37]. Previously, various researchers studied and classified the formation of a hydraulic jump due to only vegetation [48] and due to the combined effect of dike and vegetation [4]. In the hybrid defense system, the nature of the hydraulic jump varied with changing flow and vegetation conditions [4]. It is important to control the position of a hydraulic jump to reduce the flow energy and to stabilize the hydraulic structure [49,50]. Thus, the knowledge of jump type and associated water surface profiles and their classification for flow over the hydraulic structures are essential for assessment of the mitigation effect of a defense system. The types of jumps formed in this experimental investigation Similarly, the dike and vegetation (DV) together will offer more resistance to flow and a significant backwater rise. In this study, it was observed that after placing dike on the upstream side of vegetation, the maximum bac water rise increased by 15%. The backwater rise in DVI was 8.50% greater than that in DVS. The DMV model offers more resistance to the flow of water which resulted in higher backwater rise as compared to the previous models because of the moat. By placing the moat in between dike and vegetation, backwater rise was increased by 21%. By increasing the initial Froude number, the backwater rise and relative backwater rise increases linearly in DTMVI, DRMVI and DTMVS models (Figure 5a,b). In the case of a higher density of vegetation, the increase in backwater rise is 23% more than that in the case of DTMVI. In DTMVI, this value was 12% higher than that in DRMVI model. For the lower values of Froude number, i.e., Fr o = 0.40 and 0.44, the backwater rise was slightly higher in SDS than that in HDS and for the remaining values of Froude number the backwater rise in the case of HDS was comparatively higher. It happened, because in HDS, for lower values of Fr o the depth of water overtopping the dike was low so the depth achieved in front of vegetation was also low. However, for higher values of Fr o the overtopping depth was high, due to which higher depth was achieved on the upstream side of vegetation.

Hydraulic Jump and Water Surface Profile Classification
A hydraulic jump is formed when supercritical flow transforms to subcritical flow [37]. Previously, various researchers studied and classified the formation of a hydraulic jump due to only vegetation [48] and due to the combined effect of dike and vegetation [4]. In the hybrid defense system, the nature of the hydraulic jump varied with changing flow and vegetation conditions [4]. It is important to control the position of a hydraulic jump to reduce the flow energy and to stabilize the hydraulic structure [49,50]. Thus, the knowledge of jump type and associated water surface profiles and their classification for flow over the hydraulic structures are essential for assessment of the mitigation effect of a defense system. The types of jumps formed in this experimental investigation are classified in six different categories depending upon the location of the jump, as shown in Figures 6 and 7. They represent experimental photographs captured during the experimentation. The hydraulic jump was observed between dike and vegetation model (Section 2, Figure 6) and/or on the downstream side of vegetation (Section 3, Figure 6).

Hydraulic Jump Classification in SDS
In SDS the hydraulic jump formed is classified as Type I (shown in Figure 6a Table 2. A higher backwater rise, larger undulations, and a larger undulated jump were observed in OVI as compared to OVS. The relative length of the hydraulic jump L j /y 1 was maximum against Fr o = 0.65 in OVI case. In Figure 6a h 1 and h 2 are the average depths of water on the upstream and downstream side of vegetation after the jump, y 1 = minimum depth of water in hydraulic jump, L j = length of hydraulic jump, h c = critical depth, V 1 , V 2 are the average velocities on the upstream and downstream side of vegetation patch respectively and V c = critical velocity against critical depth.

Hydraulic Jump Classification in HDS
In the case of hybrid defense system (HDS), in DVS and DVI cases, the hydraulic jump formed is classified as Type II and III, as shown in Figure 6b,c. In these figures, h 1 , h 2 and h 3 are the average water depths in front of dike, in front of vegetation, and downstream of vegetation patch, respectively, h d = the average depth of water in the middle of L dv , L o = the length of start of hydraulic jump on the downstream slope of dike from the toe of dike towards the crest, L dv = distance between dike and vegetation, Fr 1 = Froude number on the upstream of dike, Fr o , h o , V o , V 1 , V 2 , V c, and L j are already define previously. In Type II, a hydraulic jump developed in Section 2 with no jump in Section 3. However, in Type III no jump formed in Section 2 with the only jump developed in Section 3. The classification of the hydraulic jump of HDS is described below and summarized in Table 2.
In  Figure 8b.
Similarly, in DRMVI model rectangular moat was used in combination with dike and intermediate vegetation to examine the effect of the shape of moat on total energy dissipation. The profiles were the same as those in the case of DTMVI. However, larger oscillations in between dike and vegetation than DTMVI was observed due to rectangular shape, because in rectangular shape a sudden drop of water was observed due to the straight walls of moat while in trapezoidal moat due to sloping surface fewer oscillations were formed. Against the value of  Figure 8c. The variation of length of the hydraulic jump can be represented by the relative length of the hydraulic jump (L j /y 1 ) [48]. In these cases, the relative length of the hydraulic jump (L j /y 1 ) increased by increasing Fr o and the maximum relative jump was noted in the case of DTMVI as compared to DTMVS and DRMVS, as shown in Figure 9.

Energy Dissipation in Single Defense System (SDS)
The energy loss due to a hydraulic jump only and the total energy loss were calculated against the Froude number range Fro = 0.40-0.65. In this study, the total energy loss for all the cases of SDS and HDS was determined, however, the energy loss due to a hydraulic jump is explained only for HDS. For SDS, the maximum energy reduction in OVS and OVI was 28% and 36.34%, respectively, against Fro = 0.40 and it decreased with higher values of Fro. While increasing the values of Fro from 0.40 to 0.65 the energy dissipation rate decreased by about 14.25% and 33.86% in OVS and OVI, respectively. The average energy loss was 25.32% and 31.49% for OVS and OVI, respectively. The results of energy dissipation show that initial Froude number and the density of vegetation (B/d) were the factors affecting the total energy loss. Based on these two important factors regression equations were developed for all the cases of SDS and HDS. In this paper, the Equations (4)-(6) were developed using EXCEL spreadsheet and graph to get the exponent type functions. The best fit curve and 95% confidence bands were computed and plotted (Figure 10a-c) for each of Equations (4)-(6) by software built in EXCEL. Equation (4), as given below, was developed to calculate the relative percentage energy loss in the case of OVS and OVI.
By using Equation (4), the energy loss (%) in the case of OVI and OVS was calculated for a certain range of Froude number as a ready reference for future use, as shown in Figure 11a.

Energy Dissipation in Hybrid Defense System (HDS)
In HDS, the maximum energy loss (%) for DVS and DVI was 49.71% and 57.28% against Fro = 0.40. While increasing the values of Fro from 0.40 to 0.65 the energy reduction rate decreased by about 44% and 37.13% and the average energy loss was 36.63% and 41%, respectively. Equation (5) given below was developed by regression analysis as described above in Section 3.2.1, to compute the relative percentage energy reduction within the selected range of Fro. The computed values of energy loss using Equation (5) for a certain range of Fro are shown in Figure 11b.
To check the effectiveness of all the cases of HDS, i.e., DTMVS, DTMVI, and DRMVI models, the maximum energy reduction due to a hydraulic jump in Section 2 (ΔEj2), in Section 3 (ΔEj3), maximum overall energy dissipation and average energy dissipation were calculated (summary is shown in Table 3) In all these cases, DTMVI performed better than DTMVS and DRMVI regarding energy

Energy Dissipation in Single Defense System (SDS)
The energy loss due to a hydraulic jump only and the total energy loss were calculated against the Froude number range Fr o = 0.40-0.65. In this study, the total energy loss for all the cases of SDS and HDS was determined, however, the energy loss due to a hydraulic jump is explained only for HDS. For SDS, the maximum energy reduction in OVS and OVI was 28% and 36.34%, respectively, against Fr o = 0.40 and it decreased with higher values of Fr o . While increasing the values of Fr o from 0.40 to 0.65 the energy dissipation rate decreased by about 14.25% and 33.86% in OVS and OVI, respectively. The average energy loss was 25.32% and 31.49% for OVS and OVI, respectively. The results of energy dissipation show that initial Froude number and the density of vegetation (B/d) were the factors affecting the total energy loss. Based on these two important factors regression equations were developed for all the cases of SDS and HDS. In this paper, the Equations (4)-(6) were developed using EXCEL spreadsheet and graph to get the exponent type functions. The best fit curve and 95% confidence bands were computed and plotted (Figure 10a-c) for each of Equations (4)-(6) by software built in EXCEL. Equation (4), as given below, was developed to calculate the relative percentage energy loss in the case of OVS and OVI.
By using Equation (4), the energy loss (%) in the case of OVI and OVS was calculated for a certain range of Froude number as a ready reference for future use, as shown in Figure 11a.

Energy Dissipation in Hybrid Defense System (HDS)
In HDS, the maximum energy loss (%) for DVS and DVI was 49.71% and 57.28% against Fr o = 0.40. While increasing the values of Fr o from 0.40 to 0.65 the energy reduction rate decreased by about 44% and 37.13% and the average energy loss was 36.63% and 41%, respectively. Equation (5) given below was developed by regression analysis as described above in Section 3.2.1, to compute the relative percentage energy reduction within the selected range of Fr o . The computed values of energy loss using Equation (5) for a certain range of Fr o are shown in Figure 11b.
To check the effectiveness of all the cases of HDS, i.e., DTMVS, DTMVI, and DRMVI models, the maximum energy reduction due to a hydraulic jump in Section 2 (∆E j2 ), in Section 3 (∆E j3 ), maximum overall energy dissipation and average energy dissipation were calculated (summary is shown in Table 3) In all these cases, DTMVI performed better than DTMVS and DRMVI regarding energy dissipation. In the case of DTMVI, maximum ∆E j2 = 27% and ∆E j3 = 4%. The maximum total energy reduced in this case was 60% and the average energy reduced was 46%. Figure 10c shows the combine trends of DTMVS, DTMVI, and DRMVI. Based on the data of DTMVS, DTMVI, and DRMVI cases, Equation (6) was developed to compute the relative percentage energy loss. The results from Equation (6) are shown in Figure 11c.
Water 2019, 11, x FOR PEER REVIEW 16 of 25 dissipation. In the case of DTMVI, maximum ΔEj2 = 27% and ΔEj3 = 4%. The maximum total energy reduced in this case was 60% and the average energy reduced was 46%. Figure 10c shows the combine trends of DTMVS, DTMVI, and DRMVI. Based on the data of DTMVS, DTMVI, and DRMVI cases, Equation (6) was developed to compute the relative percentage energy loss. The results from Equation (6) are shown in Figure 11c.  dissipation. In the case of DTMVI, maximum ΔEj2 = 27% and ΔEj3 = 4%. The maximum total energy reduced in this case was 60% and the average energy reduced was 46%. Figure 10c shows the combine trends of DTMVS, DTMVI, and DRMVI. Based on the data of DTMVS, DTMVI, and DRMVI cases, Equation (6) was developed to compute the relative percentage energy loss. The results from Equation (6) are shown in Figure 11c.   In all the cases, the delay time decreases by increasing Fro. The relative time with respect to TWM was also maximum in the case of DTVMI, which will result in a maximum delay in flood arrival. To examine the relationship between the relative delay time for different cases of SDS and HDS, a pairwise comparison method namely least significant difference (LSD) method explained in the literature [52] was used. The calculated critical value of LSD was 0.21. So, the minimum difference between a pair of means (µ1-µ2) necessary for statistical significance was 0.21. There were 21 pairwise groups (G1 to G21). It was observed that for the difference of means of pairs G5-G6, G10-G11, and G14-G15 were greater than 0.21 means that there is a significant difference between the means of these pairs exist. OVI and DTMVI. It is concluded from G5, G6, G10, and G11 that the relative delay time of HDS (DTMVI and DRMVI) is higher than SDS (OVS and OVI). The results of G14 and G15 indicated that the relative time of delay increased by placing the moat in between the dike and vegetation of higher density. The complete pairwise comparison list is shown in Table A1 (Appendix A).
The average relative delay time and single standard error bars in either direction were represented against the relative time of delay for each case of SDS and HDS in Figure 12. The results indicated that the average relative delay time was observed to be maximum in the DTMVI model and minimum in the OVS model.

Delay in Floodwater Arrival Time and Water Level Rise
In the first step, the delay time T WM was computed without placing any obstruction model in the flume. Similarly, the delay times for a single defense system and a hybrid defense system were computed. The values of T WM were in the range of 28-80 s, against the different values of initial Froude number, whereas it was 29-115 s for OVS, 30-118 s for OVI, 32-129 s for DVS and 34-132 s for DVI. Similarly, the delay time for DTMVS, DTMVI, and DRMVI was 33-137 s, 37-140 s and 36-139 s respectively. The maximum delay time was 2 min 20 s for DTMVI against Fr o = 0.40. In all the cases, the delay time decreases by increasing Fr o . The relative time with respect to T WM was also maximum in the case of DTVMI, which will result in a maximum delay in flood arrival. To examine the relationship between the relative delay time for different cases of SDS and HDS, a pairwise comparison method namely least significant difference (LSD) method explained in the literature [52] was used. The calculated critical value of LSD was 0.21. So, the minimum difference between a pair of means (µ1-µ2) necessary for statistical significance was 0.21. There were 21 pairwise groups (G1 to G21). It was observed that for the difference of means of pairs G5-G6, G10-G11, and G14-G15 were greater than 0.21 means that there is a significant difference between the means of these pairs exist. OVI and DTMVI. It is concluded from G5, G6, G10, and G11 that the relative delay time of HDS (DTMVI and DRMVI) is higher than SDS (OVS and OVI). The results of G14 and G15 indicated that the relative time of delay increased by placing the moat in between the dike and vegetation of higher density. The complete pairwise comparison list is shown in Table A1 (Appendix A).
The average relative delay time and single standard error bars in either direction were represented against the relative time of delay for each case of SDS and HDS in Figure 12. The results indicated that the average relative delay time was observed to be maximum in the DTMVI model and minimum in the OVS model.

Hydraulic Jump Formation and Energy Dissipation in SDS
In this paper, single and hybrid defense systems were tested in relation to the backwater rise, relative energy dissipation and relative delay time of flood. In SDS, vegetation of variable densities (B/d) were considered without any other model as many researchers did in the past [47,48].  25) is rare and difficult to construct in the field. Various researchers performed experimental and numerical analysis of single defense system (dike or vegetation). The formation of undulated hydraulic jump behind the weir-like obstacles was studied experimentally and numerically to correlate the results obtained [39]. In the case of overtopping flow over weir embankment, different types of flow regimes were classified into four types [53]. Al-hashimi (2016) studied the energy dissipation and the flow characteristics over broad crested and single-step broad crested weirs and the maximum energy dissipation was 27.35% [41].
Previously, the reduction of energy was examined by single line vegetation model placed laterally across open channel flume under subcritical flow conditions and an experimental formula for the energy reduction efficiency of flood protection vegetation was proposed based on two variables i.e., Froude number and density of vegetation. In that study, it was found that a denser cylinder arrangement is more efficient in dissipating the flood energy [20]. Recently, Valyrakis (2015) carried out experimental studies to quantify the effect of increasing the density of bank vegetation on the flow velocity in the main channel and inside the vegetation [54]. Various researchers performed an experimental analysis of the formation of a hydraulic jump in an open channel flow [53][54][55]. The formation of an undulated hydraulic jump by using a sluice gate mounted in open channel flume was also reported in the literature [56][57][58]. The hydraulic jump was also classified as breaking and non-breaking undulated jump based on Froude number with aspect ratio [56] and ho/L [48], where b = width of the channel, y1 = depth of water at the toe of a hydraulic jump, ho = initial water depth and L = wavelength. In the current study, for OVS and OVI undulated hydraulic jump was observed only for the higher values of Fro (0.60, 0.63, and 0.65) and no jump was observed for lower values of Fro. The similar results were found in the literature [48]. The relative energy dissipation for a vegetation thickness of dn-180 No. cm calculated by Pasha and Tanaka

Hydraulic Jump Formation and Energy Dissipation in SDS
In this paper, single and hybrid defense systems were tested in relation to the backwater rise, relative energy dissipation and relative delay time of flood. In SDS, vegetation of variable densities (B/d) were considered without any other model as many researchers did in the past [47,48] 25) is rare and difficult to construct in the field. Various researchers performed experimental and numerical analysis of single defense system (dike or vegetation). The formation of undulated hydraulic jump behind the weir-like obstacles was studied experimentally and numerically to correlate the results obtained [39]. In the case of overtopping flow over weir embankment, different types of flow regimes were classified into four types [53]. Al-hashimi (2016) studied the energy dissipation and the flow characteristics over broad crested and single-step broad crested weirs and the maximum energy dissipation was 27.35% [41].
Previously, the reduction of energy was examined by single line vegetation model placed laterally across open channel flume under subcritical flow conditions and an experimental formula for the energy reduction efficiency of flood protection vegetation was proposed based on two variables i.e., Froude number and density of vegetation. In that study, it was found that a denser cylinder arrangement is more efficient in dissipating the flood energy [20]. Recently, Valyrakis (2015) carried out experimental studies to quantify the effect of increasing the density of bank vegetation on the flow velocity in the main channel and inside the vegetation [54]. Various researchers performed an experimental analysis of the formation of a hydraulic jump in an open channel flow [53][54][55]. The formation of an undulated hydraulic jump by using a sluice gate mounted in open channel flume was also reported in the literature [56][57][58]. The hydraulic jump was also classified as breaking and non-breaking undulated jump based on Froude number with aspect ratio [56] and h o /L [48], where b = width of the channel, y 1 = depth of water at the toe of a hydraulic jump, h o = initial water depth and L = wavelength. In the current study, for OVS and OVI undulated hydraulic jump was observed only for the higher values of Fr o (0.60, 0.63, and 0.65) and no jump was observed for lower values of Fr o . The similar results were found in the literature [48]. The relative energy dissipation for a vegetation thickness of dn-180 No. cm calculated by Pasha and Tanaka (2017) was 26-28% and 29-32% in sparse and intermediate vegetation, respectively [48]. However, in the current research, the relative energy dissipated was in the range of 24-28% and 24-36% in sparse and intermediate vegetation respectively. The range of Fr o used by Pasha and Tanaka (2017) was 0.57 to 0.73 [48] and in this study, it was 0.40 to 0.65. The back-water rise increases by increasing the density of vegetation [47]. The back-water rise was observed higher in OVI than OVS as shown in Figure 5a.
As discussed in Section 2.1.2, the overtopping flow due to high floods or embankment breaching is always in an oblique direction with respect to the flow in the main channel. However, the obstacles like roads, summer dikes, groins, or vegetation in the floodplains of a river can be encountered in various orientations. i.e., oblique or perpendicular and laid out by the floodplain management authorities. The difference between an oblique weir and plain weir is that in the case of the oblique weir, its crest makes an angle less than 90 • with the direction of flow approaching the obstruction and in the perpendicular (plain) weirs the flow is at 90 • to the weir or embankment crest [42]. Previously, some researchers used oblique weirs in their studies of open channel hydraulics [42,59,60] and other used plain or perpendicular weir/embankment [20,40,61,62]. Ali (2010) studied the vegetated weir-like structures oriented perpendicular to the flow direction to quantify the energy head loss [63]. However, it was reported in the previous studies that the flow always turns its direction when it reaches and passes the weir, towards an almost perpendicular orientation with respect to crest and the results are more deviating in the case of high obliqueness and when flow reaches an oblique weir (weir-like obstacles), it tends to change its direction normal to the weir crest [42]. For an oblique weir, the flow velocity is decomposed into two components, parallel and perpendicular to the flume axis as well as parallel and perpendicular to the weir axis. It is assumed by various authors that when the flow approaches the weir the crest parallel component does not change its magnitude, the accelerating force only acts on the velocity component perpendicular to the weir crest and the same holds for the deceleration process [43]. Therefore, to achieve the maximum possible energy dissipation similar to the results obtained in the current study, the arrangement in the real situation (in the field) the weir-like obstacles will have to be oblique at such an angle that its crest becomes perpendicular to the flow. When the flow is oblique, and weir will be oblique at the same angle the both will become perpendicular to each other, and energy dissipation will be similar to that obtained in the current research. However, for most commonly constructed obstructions, the overtopping flow is oblique to the obstruction. Therefore, the energy loss results reported in this paper will have to be reduced by a certain factor for real-life uses. These factors can be found by rigorous analysis taking into consideration the momentum and continuity equations or from previous publications [42,43,59,60]. However, the flow for a hybrid system investigated in this paper is highly complex. Therefore, some further experimental studies are needed to find a relationship between the energy loss by an obstruction across an oblique and perpendicular flow.

Hydraulic Jump Formation and Energy Dissipation in HDS
Recently, the hybrid defense system was introduced to minimize the damages of human lives and infrastructure due to extreme natural events [4], although this system was introduced in the response to the damages caused by tsunamis. However, Pasha and Tanaka (2017) investigated the energy loss through vegetation and a downstream undulated hydraulic jump for steady subcritical flow conditions with Froude number less than unity [48] and the Froude numbers of past tsunamis were estimated at 0.7-2.0 from field measurements and theory [64]. Moreover, at many locations, the tsunami's flow was subcritical and had a Froude number between 0.7 and 1 [65]. The supercritical flow generated by tsunami converted to subcritical towards the inland side as Tanaka (2013) estimated the Froude number was 0.9-0.6 at 500 to 550 m distance from the shoreline [66] and a simulation study resulted that the Froude number decreased from 1.6 near the shoreline to 0.6 at 550 m from the shoreline [29]. Thus, it is equally important to determine the energy loss through single and hybrid defense system either the flow is subcritical or supercritical.
When flow passed through the hybrid defense system, its velocity decreased due to the contraction in the flow and a decelerating region is formed due to sudden expansion on the downstream side of the vegetation model resulted from the energy head loss in the expansion region [39]. Previously, the researchers reported the damage of houses and trees due to water overtopping the embankment [55]. Thus, a secondary structure on the downstream side of the embankment is required to reduce the flow velocity. The hydraulic jump is formed after placing any structure on the downstream of embankment [50]. In this regard, a secondary structure of different layers of single-and double-layer vegetation on the downstream of the embankment was tested with supercritical flow conditions [4]. The jump classified was jump A (formed on the bed of channel) and jump B (formed on the downstream slope of the embankment). According to that study for the initial values of Fr o , Type A jump was formed on the bed of the channel in between embankment and vegetation while Type B jump was observed for the higher values of Fr o . The maximum energy dissipated due to jump was 39%, and the maximum total energy reduction was in between 27% and 54% [4]. The hydraulic jump based on the values of Froude number was also classified by 2-D numerical analysis in the previous study [67].
In supercritical flow conditions, a higher Froude number increases the length of the hydraulic jump and the behavior of an undular hydraulic jump downstream of vegetation is investigated better against a low initial Froude number as reported by different researchers [66,68]. In the current study, the performance of single and hybrid defense systems was studied against the floodwater with subcritical flow conditions. For DVS and DVI, the hydraulic jump was classified against the initial values of Fr o (0.40-0.57). However, for higher values (i.e., Fr o = 0.60, 0.63 and 0.65) the jump formed was only on the downstream of vegetation and no jump was observed between dike and vegetation. The total relative energy reduction for DVS and DVI was between 28% and 57% and the average relative energy dissipation was 37% and 41% in DVS and DVI, respectively.
The hybrid defense systems are used to enhance coastal resilience in the protection from flooding [32]. Recently, hybrid defense system comprising dike, moat, and vegetation was tested by changing the layout of the components. The fluid force, tsunami arrival time, and overflow volume for different cases were measured. The hybrid defense system in the order of vegetation, a trapezoidal moat, and an embankment performed best in relation to the overflow volume, arrival time, and fluid force [32].
In the current paper, a hybrid system in the order of dike, moat (trapezoidal and rectangular) and vegetation of variable densities were tested with subcritical flow conditions. Backwater rise, energy dissipation due to hydraulic jump, and total relative energy dissipation were measured for each case to clarify the optimal hybrid defense system. To examine the effect of the shape of the moat, rectangular and trapezoidal shapes were used. For DTMVS, DTMVI, and DRMVI weak and undular jumps of Type IV, V, and VI were observed. The maximum backwater rise was observed for DTMVI case because the denser vegetation offered more resistance to flow resulting in higher backwater rise, as shown in Figure 5a. The energy loss due to a hydraulic jump was usually low [49] and the fraction of energy dissipated due to an undular jump is less than 5% [69]. In the current study, the amount of energy dissipation due to an undular jump on the downstream side of vegetation was well below 5% (∆E j3 in Table 3). The maximum total relative energy dissipation and maximum average energy dissipation by hybrid defense system was 60% and 46% in DTMVI model. The arrival time was also highest in the DTMVI case. Contrarily, the hybrid defense system in the order of vegetation, moat and dike (where vegetation was in front) performed best in the previous study [32]. However, in this combination, in the case of a devastating flood event, the water may overtop the embankment and it will be difficult to connect this water to the river. Secondly, in that study, there was no space considered between the moat and embankment. In case of the breaching or erosion of embankment, the moat may be filled with erodible material and will not serve as a moat after its filling.
In the current research, a combination of dike, moat, and vegetation (where dike being on the riverside) was tested and in case of overtopping flow the water will come back to the moat which will further join the river by seepage action. Some water will be stored in the moat which will act as a flood harvesting structure. Moreover, the space between the downstream toe of dike and the start of moat (L dv ) was kept 15 cm to ensure the safety of the moat structure in case of breaching of dike. The rectangular moat was less effective than the trapezoidal moat based on the total energy loss and the arrival time. The sides of a rectangular hydraulic structure are also not stable. The tree breaking phenomena also occurred in some regions of Japan due to the large wave thrust of tsunami [29]. However, the floodwater has less capacity to break trees and other components of hybrid defense system [70]. Therefore, the defense system proposed in the current study will be more effective in the river environment. In the light of the results obtained in this study, a hybrid defense system in the order of dike, trapezoidal moat, and intermediate vegetation is highly recommended. The limitation of the proposed system is the space available to construct the defense system. In case of non-availability of the space, a single defense system will be preferred.

Conclusions
The flow structure with energy loss in a subcritical flow against the selected range of initial Froude numbers (Fr o = 0.40-0.65) around single defense system (SDS) and the hybrid defense system (HDS) was investigated experimentally in a flume. The main conclusions of the study are as follows: • The backwater rise is maximum for OVI and DTMVI in SDS and HDS, respectively. The backwater rise is directly proportional to the density of vegetation and the value of initial Froude number. The water surface slope also increases by increasing vegetation density. The denser and wider the vegetation, the larger is the total energy dissipation in both SDS and HDS cases.

•
In SDS only undulated hydraulic jump was observed in both OVS and OVI, resulting in a significant energy loss. In HDS both weak and undulated hydraulic jumps were formed and in the case of DTMVI, the maximum energy loss due to hydraulic jump formed in between the dike and vegetation was 27% and 4% energy was dissipated due to the formation of jump on the downstream side of vegetation. The maximum total energy reduced in this case was 60% and the average energy reduced was 46%. Similarly, in the case of DRMVI, the maximum value of energy loss due to hydraulic jump between dike and vegetation was 22% and 3% energy was dissipated due to the formation of hydraulic jump on the downstream side of vegetation. The maximum total energy reduction and average energy reduction were 60% and 43.75%, respectively. In all these cases, the rate of energy reduction due to the hydraulic jump decreases by increasing Fr o .

•
The performance of the DTMVI model to delay the arrival time of floodwater is the highest among all the models investigated in this paper • The "moat" can serve as floodwater harvesting and increasing the response time of flash floods generated from hill torrents. The shape of the moat affects the reduction of energy in general. However, its trapezoidal shape performs better than rectangular shape.
Under the above conclusions, a hybrid defense system is highly recommended where enough space is available across the river flow. The proposed system will be highly effective in various flood plains similar to that of Indus River. Recommendations: The response of the hybrid defense system with varying dimensions of dike and moat should be investigated in future and the true behavior of the hybrid defense system should be analyzed by considering the oblique flow towards the obstruction.