Using 2D HEC-RAS Modeling and Embankment Dam Break Scenario for Assessing the Flood Control Capacity of a Multi-Reservoir System (NE Romania)

: Using hydraulic modeling techniques (e.g., one-dimensional/two-dimensional (1D/2D) hydraulic modeling, dam break scenarios) for extracting the ﬂood settings is an important aspect of any action plan for dam failure (APDF) and ﬂood mitigation strategy. For example, the ﬂood hydraulic models and dam break scenario generated based on light detection and ranging (LiDAR)-derived digital elevation models (DEMs) and processed in the dedicated geographic information systems (GIS) and hydraulic modeling software (e.g., HEC-RAS—Hydrologic Engineering Center River Analysis System, developed by USACE HEC, Davis, CA, USA) can improve the ﬂood hazard maps in case of potentially embankment dam failure. In this study, we develop a small-scale conceptual approach using 2D HEC-RAS software according to the three embankment dam break scenarios, LiDAR data (0.5 m spatial resolution), and 2D hydraulic modeling for the Ba¸seu multi-reservoir system which belongs to the Ba¸seu River (NE Romania) including R1—Cal Alb reservoir, R2—Movileni reservoirs, R3—T ă t ă r ă ¸seni reservoirs, R4—Negreni reservoirs, and R5—H ă ne¸sti reservoirs. In order to test the ﬂood control capacity of the Bas , eu multi-reservoir system, the Cal Alb (R1) dam break scenario (piping failure) was taken into account. Three 2D stream ﬂow modeling conﬁgurations based on R1 inﬂow rate with a 1% (100 year), 0.5% (500 year), and 0.1% (1000 year) recurrence interval and the water volume which can be accumulated with that speciﬁc inﬂow rate (1% = 10.19 × 10 6 m 3 ; 0.5% = 12.39 × 10 6 m 3 ; 0.1% = 17.35 × 10 6 m 3 ) were computed. The potential ﬂood wave impact was achieved on the basis of different ﬂood severity maps (e.g., ﬂood extent, ﬂood depth, ﬂood velocity, ﬂood hazard) generated for each recurrence interval scenario and highlighted within the built-up area of 27 settlements (S1–S27) located downstream of R1. The results showed that the multi-reservoir system of Bas , eu River has an important role in ﬂood mitigation and contributes to the APDF in the context of climate change and the intensiﬁcation of hydrological hazard manifestation in northeastern Romania.


Introduction
The location of human settlements near water sources (e.g., lotic ecosystems-rivers, streams; lentic ecosystems-lakes, ponds, marshes) has proven to be a very important aspect of habitation practice since ancient times [1][2][3]. In this context, when the water supply was limited during drought periods, the water reserve necessity led to the construction In Romania, according to 1422/192 Ordinance from 16 May 2012, which is additional to 2007/60/CE Directive, the owners or the administrators of a reservoir with a dam height over 10 m and/or with water volume more than 10 million m 3 have the obligation to make an APDF [29]. Furthermore, in the same ordinance, it is specified that an APDF can be made for reservoirs with a dam height lower than 10 m and/or a water volume > 10 million m 3 only if a hypothetical dam failure can cause significant material damage or human casualties. However, with more than 1000 registered dams in Romania (e.g., 250 dams classified as big dams, 85 classified as embankment big dams), only few of them have an APDF. In this case is also the Cal Alb reservoir (R1) located in the upper sector of the Başeu multi-reservoir system (Figure 1a). Similar to other embankment dams in the country which were built up during the communist period (before 1989), the Cal Alb reservoir has not yet received an APDF, even if it accomplishes all the 1422/192 Ordinance requirements (e.g., dam height more than 10 m, water volume greater than 10 million m 3 , frequent dam faced infiltrations). Thus, given the existence of a precedent on the integrity of the dam body (e.g., piping process), the fact that the Cal Alb reservoir provides the water resource for an important community in the Botos , ani County, the fact that downstream of Cal Alb reservoir are located 25 rural settlements (S1-S10; S12-S23; S25-S27) and two urban settlements (S11-Săveni and S24-S , tefănes , ti), we consider that the Cal Alb reservoir provides the optimal conditions to test the integrity of the dam body and the capacity of the multi-reservoir system (R2, R3, R4 and R5) located downstream to take over and reduce the flood wave produced by a potential dam break scenario [30] (Figure 1b).  [31]. The abbreviations within the (b) map indicate the location of R1-R5 reservoirs and dams built up for protection of S1-S27 settlements.
In this study, in order to test the flood control capacity of the Bașeu multi-reservoir system, we generated the first dam break scenario for the Cal Alb (R1) reservoir using a LiDAR-derived DEM with 0.5 m spatial resolution and 2D HEC-RAS hydraulic modeling [30][31][32][33][34][35]. The different scenarios were based on three R1 inflow rates with different recur-  [31]. The abbreviations within the (b) map indicate the location of R1-R5 reservoirs and dams built up for protection of S1-S27 settlements.
In this study, in order to test the flood control capacity of the Bas , eu multi-reservoir system, we generated the first dam break scenario for the Cal Alb (R1) reservoir using a LiDAR-derived DEM with 0.5 m spatial resolution and 2D HEC-RAS hydraulic modeling [30][31][32][33][34][35]. The different scenarios were based on three R1 inflow rates with different recurrence intervals (1% or 100 years; 0.5% or 500 years; 0.1% or 1000 years) and the water volume which can be accumulated with that specific inflow rate (1% = 10.19 × 10 6 m 3 ; 0.5% = 12.39 × 10 6 m 3 ; 0.1% = 17.35 × 10 6 m 3 ) [7]. We used 2D HEC-RAS hydraulic modeling instead of 1D HEC-RAS hydraulic modeling because the 1D method cannot calculate the propagation of the multidirectional flood wave [31,[36][37][38]. On the other hand, the 2D hydraulic models are the most appropriate solution in order to simulate the flood wave lateral diffusion [39][40][41][42][43]. Therefore, in order to accomplish the main objective of the study, a 2D hydraulic break scenario was created using HEC-RAS v.5.0.7 software [23][24][25]. Considering the recent infiltration history of Cal Alb (R1) dam, a piping failure mode was used [7,44]. Additionally, on the basis of the adapted Federal Emergency Management Agency (FEMA) methodology on hazard classification [45], after the Australian Institute for Disaster Resilience (AIDR) hazard classification [46], the flood severity assessment was done for all 27 settlements potentially affected by R1 dam failure. The results provide for the first time the flood control capacity of the Bas , eu multi-reservoir system in case of Cal Alb (R1) failure and contribute to the APDF in the northeastern region of Romania.

Study Area
Başeu river basin (96,770 ha) is located in northeastern Romania, in the administrative territory of Botos , ani County. With a total length of 118 km, the Başeu River is one of the most important tributaries of Prut River (natural border between Romania and Republic of Moldova) (Figure 1a) [47][48][49]. The lowest altitude is found on the common floodplain of Bas , eu and Prut River (54 m above sea level (a.s.l.)) and the highest altitude is on the top of Pădurea Cristesti Hill (318 m a.s.l.) [7,30]. The geological conditions are characterized by the presence of recent alluvial deposits (upper Holocene) accumulated over 200-300 m sandy clay deposits (Miocene). Climate conditions control 60% of the water flow rate, often overtaken in transition seasons and in periods with maximum rainfall. The average air temperature values are between 7.7 • C and 9.6 • C. The annual average precipitation is between 450 and 624 mm, with peak values occurring in the high basin areas. The groundwater contributes between 30% and 40% of the average annual flow. The multi-annual average flow rate at S , tefănes , ti station (south of the basin) is 1.6 m 3 /s, with fluctuations between 0.1 m 3 /s (2016) and 6.9 m 3 /s (1969). The historical flow was recorded in July 1969 and was 330 m 3 /s. Negative hydrological events associated with the maximum discharge occurred in 1973 (100 m 3 /s) and 2005 (124 m 3 /s) [30].
The Başeu multi-reservoir system consists of five permanent and nonpermanent multi-reservoirs areas (total water surface 953.3 ha) from upstream to downstream as follows: permanent multi-reservoirs-Cal Alb (R1) and Negreni (R4); nonpermanent multi-reservoirs-Movileni (R2), Tătărăşeni (R3), and Hănes , ti (R5) (Table 1; Figure 1b). According to Romanian legislation regarding the management and operating rules of dam reservoirs, permanent reservoirs are filled up with water throughout the operation period, and the nonpermanent reservoirs are maintained empty for flood mitigation. In the last decade, within Başeu river basin, the nonpermanent reservoirs were used more for fish farming instead of flood control. Moreover, more than 50 polders or lateral structures (total surface of 241.8 ha), which ensure the protection of 1071 ha (690 ha of agricultural land, 320 ha of fish farming, 61 ha of urban area), were built up in the Başeu floodplain over time. Downstream of Cal Alb (R1) reservoir are located 27 settlements (S1-S27) with a total number of 33,869 inhabitants (S1-S10; S12-S23; S25-S27-rural settlements with 19,094 inhabitants; S11 and S24-urban settlements with 14,775 inhabitants) ( Table 2; Figure 1b). The total built-up area is 3014 ha (rural built-up area-2613 ha; urban built-up area-401 ha) with an average population density of 1123 inhabitants/km 2 (rural popula-tion density-731 inhabitants/km 2 ; urban population density-3685 inhabitants/km 2 ) [30] ( Figure 1b).
Within the Bas , eu hydrographic basin, the Cal Alb reservoir (R1) (water surface-174.5 ha at normal water level (NWL)) is placed in the upper sector of Başeu River (northeastern part of the basin) (Figure 2a), at 900 m downstream of the confluence with the Ciolac Brook. The river basin of the R1 reservoir has a total surface of 18,300 ha, with an altitude between 318 m (Pădurea Cristesti Hill) and 123 m (R1 water mirror level). The average altitudinal class is between 125 and 200 m (>50%). The embankment dam was built between 1971 and 1973 through decree No. 118 of 28 May 1971 [50]. According to the Romanian National Classification (RNC) of reservoirs and lakes, the Cal Alb reservoir is classified as category B (special importance) and as the third class of importance [51]. The embankment dam of R1 reservoir (Figure 2b) was built up transversely on the Başeu riverbed axis. The main feature characteristics of the R1 dam are as follows: the longitudinal section is trapezoidal where the base has a length of 295.07 m and the dam crest width does not exceed 5 m; the maximum height of the dam is 14.5 m (10.22 m without dam crest); the upstream slope report coefficient is 1:3.5 due to the concrete slab wall for protection, and the downstream slope is 1:3 being equipped with a grass layer for protection; the dam crest elevation is 134.13 m a.s.l. and the dam base elevation is 123.91 m a.s.l. (Figure 2c). The dam body was filled with local materials (>95% clay) and isolated with a 10 cm concrete layer on the water side and with a 20 cm grass layer on the opposite side. Based on the lithological profiles, two types of geological formations are highlighted: the geological foundation, made up of gray clays deposits, and the superficial layer (Başeu alluvial deposits) with a thickness between 6.5 m and 8.5 m Table 1. Characteristics of the main dams within Bas , eu river basin according to Romanian National Classification (RNC) [51].  foundation, made up of gray clays deposits, and the superficial layer (Başeu alluvial deposits) with a thickness between 6.5 m and 8.5 m   This status was acquired because the R1 reservoir ensures the water supply for irrigation of 1060 ha of agricultural land and also ensures 88 ha of water surface for fish farming. With a dam height of 11 m and a total storage volume of 25 million m 3 , Negreni reservoir (R4) is the second permanent reservoir downstream of Cal Alb reservoir. The nonpermanent reservoirs are R2, R3, and R5. These reservoirs have a dam height between 4 and 5 m and a storage water volume between 2.5 and 6 million m 3 (Table 1).

Reservoir
Currently, the R1 reservoir and dam are under the management of the Prut-Bîrlad Water Basin Administration (PBWBA) [50,51].
In 1973, less than a year after the construction of the R1 reservoir, there were first infiltrations occurring through the dam foundation and embankment body with immediate effects downstream where the entire area became unstable due to the swamping process. According to the technical reports made at that time by the authorities, due to the presence of the aquifer layer at just 1 m below the dam foundation combined with the high permeability of the lithological deposits, which allowed a good circulation of groundwater, all of these affected the dam body integrity. However, the dam operated for over 35 years with these issues until the period between 2009 and 2013 when a dam stabilization project was carried out.      [51] in 2015, the northeastern territory of Romania was airborne scanned using LiDAR technology ( Figure 3a). Within this project, the detailed flight level on the Başeu multi-reservoir system (study area) was classified as the first (A) level of importance-flights in the area adjacent to watercourses and urban areas. Accordingly, using the Leica ALS60 Airborne Laser Scanner (ALS) [52], a detailed elevation point cloud (.las files) of the topographic surface was achieved. After this stage, the obtained .las files were classified and vegetation areas were extracted [53], and the individual raster files (.tiff files) with an area of 0.25 km 2 and a spatial resolution of 0.5 m/pixel were generated. The LiDAR-derived DEMs obtained via the ArcGIS processing of more than 2000 tiles (.tiff files) were unified to a spatial resolution of 0.5 m/pixel. In this way, we obtained a high-accuracy DTM which included detailed topographic information related to hydro-technical structures (e.g., dams, reservoirs, polders, lateral structures) which considerably improved the 2D hydraulic model. Furthermore, to capture all the topographic details given by the presence of new constructions (post 2015) within the built-up area such as houses or other types of buildings (e.g., attachment buildings, administrative buildings, industrial buildings), they were integrated with the DTM for higher accuracy of the 2D hydraulic model [30]. The building database was obtained by on-screen digitizing techniques using orthophotos collected in 2018. To each building, an average height of 6 m was assigned, in order to rasterize and join to the final LiDAR-derived DTM [54].

Development of Bathymetric Model
In order to generate the 2D hydraulic models, we decided to integrate the R1 bathymetric model within the LiDAR-derived DTM. The bathymetric data for the R1 reservoir were obtained through the SMIS-CSNR No. 17945 project results [51]. The bathymetric measurements were made in the summer of 2015 using a motorboat equipped with a singlebeam echo-sounder. More than 40,000 depth point networks were acquired in two different coordinate systems: WGS 84-Global Datum, EPSG: 4326 and Stereo 1970-National Datum, EPSG: 31700. After processing the depth points, a bathymetric model with a spatial resolution of 0.5 m/pixel was obtained [55,56]. Unfortunately, from all reservoirs within the study area, only the R1 reservoir has bathymetric measurements, whereas, for R2-R5 reservoirs, we used the normal water level extracted from LiDAR-derived DTM. However, for the final DTM used in the 2D hydraulic modeling of the R1 dam failure, the bathymetric data were integrated. In this context, a raster dataset was generated in order to join to the LiDAR-derived DTM and the R1 bathymetric model.

Hydrological Data
The Başeu river sector between the source and the confluence with the Cal Alb (R1) reservoir has a length of 17 km and ensures an average multiannual inflow rate of the 0.365 m 3 /s. In 2014, the official operating rules (OORs) for R1 reservoir were drawn up by PBWBA. According to the OORs of R1 reservoir, the inflow rate and the water volume for the three scenarios intervals (0.1% or 1000 year; 0.5% or 500 year; 1% or 100 year) were extracted ( Table 3). The recurrence intervals were calculated by PBWBA on the basis of hydrological data obtained from the two hydrometric stations located downstream of the R1 reservoir. The hydrological data series were between 31 years (1969-2000, Havârna gauging station) and 53 years (1966-2019 S , tefănes , ti gauging station). The flow sizing rate Q2% (50-year recurrence interval) was 145 m 3 /s. The flow sizing rate Q0.5% + 20% (500-year recurrence interval + 20%) was 270 m 3 /s. Thus, to ensure the life of the aquatic ecosystem for R1 reservoir, an inflow rate of 0.015 m 3 /s is needed. Table 3. The Q max (m 3 /s) inflow rates and volume of water contained (million m 3 ) estimated by PBWBA for 1% (100-year), 0.5% (500-year), and 0.1% (1000-year) recurrence intervals for Cal Alb (R1) reservoir.

2D HEC-RAS Modeling
The 2D hydraulic modeling requires an impressive and accurate volume of information regarding the spatial data (e.g., terrain model, bathymetric model, roughness coefficient) and hydrological data. According to these requirements, the HEC-RAS software [23][24][25] has the ability to simulate an inflow flood using three different methods: one-dimensional (1D) unsteady flow model [57], two-dimensional (2D) unsteady flow model [58], and level pool flow model [59]. Therefore, due to the high density of human settlements downstream of the R1 reservoir in conjunction with the availability of all necessary data (e.g., LiDAR-derived DTM, bathymetry data, hydrological data, and roughness information), we considered the 2D hydraulic modeling as the most appropriate method for testing the flood control capacity of Başeu multi-reservoir system. Furthermore, due to the high density of settlements (rural and urban areas) downstream of R1 reservoir, we chose to represent the built-up area by integrating each building in the LiDAR-derived DTM for a high accuracy of the 2D HEC-RAS hydraulic model. This process was detailed in Section 2.1.1 [51,54]. The advantage of using 2D modeling is that this method can capture the changes in water elevation when a flood wave arrives, as well as when the water elevation starts decreasing due to the dam break [60,61].
Overall, the HEC-RAS software can perform a 2D unsteady flow routing with two different equations: Equation (1), 2D full Saint Venant (also called full momentum equation) [23,25,31,62]; Equations (2) and (3), 2D Diffusion wave equation (which is set as the default in HEC-RAS software) [23,25]. In the process of choosing the proper equation for the breach scenario, two different hydraulic plans with Equation (1) and with Equations (2) and (3) were created, and the results were compared. If there are no significant differences between the two models, the user can proceed with the diffusion wave equation [23]. In our case, the two models were similar, and we decided to complete the analysis using the default equation (2D diffusion wave equation) which assures a greater stability property and a faster computational time.

∂ζ ∂t
∂p ∂t ∂q ∂t where h is the water depth (m), p and q are the specific flow in the xand y-directions (m 2 ·s −1 ), ζ is the surface elevation (m), g is the acceleration due to gravity (m·s −2 ), n is the Manning resistance, ρ is the water density (kg·m −3 ), τ xx , τ yy , and τ xy are the components of the effective shear stress, and f is the Coriolis (s −1 ). When the diffusive wave is selected, the inertial terms of the momentum equations are neglected (Equations (2) and (3)).

Dam Break Scenario
On the basis of historical data obtained from the dam breaks around the world, scientists identified two main modes of dam failure: piping and overtopping [63][64][65]. Although information on dam breaks is numerous, no pattern has yet been identified regarding the shape of the breach, its location within the dam, or the time formation. To estimate breach characteristics, four methods can be used: (i) comparative analysis of the analyzed dam with another similar historical dam failure [23,66]; (ii) regression equations developed from historical dam failures [23,44,67]; (iii) using the velocity vs. erosion rates [23]; (iv) using computer models and principles of hydraulics [23,59,[68][69][70][71]. In this case, we used regression equations to estimate the dam breach characteristics [44].
For this study, we assumed that a historical flood can occur in the river basin of the Cal Alb reservoir. The inflow for each scenario would create a dam breach due to the piping process. The flood wave created by the dam break would lead to the filling of the downstream reservoirs. Due to the location of the reservoirs, polders, and lateral structures downstream of the R1, the flood wave severity would be diminished. The flood wave would overlap with normal level water for each reservoir located downstream of the R1. In order to run the unsteady flow analysis with the breach scenario [15,23,25,31,72], six steps were followed: (a) the storage area extent was imported, (b) the dam characteristics were introduced in order to create the dam, (c) the breach parameters were calculated, (d) the subgrid model (mesh) was created, (e) the roughness data were introduced, and (f) the inflow data were set and the computational time was calculated [23,25].
The first step was to import the storage area extent in the HEC-RAS software. The following characteristics of the dam were introduced: dam length-295.07 m; dam crest width-5 m; dam height-14.5 m. Once the dam was created, the following parameters were introduced to breach the dam: center station, final bottom width, final bottom elevation, left-and right-side slope, breach weir coefficient, and piping coefficient. According to the equations developed by Froehlich [44], who studied 74 dam breaks, we estimated two parameters: average breach width using Equation (4) and breach formation time according to Equation (5). On the basis of the average breach width, average slide slope, and height of the final breach (Figure 3c), we calculated the breach bottom width using Equation (6).
where B ave is the average breach width (m), K 0 is constant (1.0 for piping failure), V w is the reservoir volume at time of failure (m 3 ), h b is the height of the final breach (m), g is the gravitational acceleration (9.80665 m·s −2 ), t f is the breach formation time (s), W b is the final bottom width, H is the left slope side (0.7 for piping failure), and V is the right slope side (1.0 for piping failure). According to Equation (4), the average breach width is 55.6 m; according to Equation (5), the breach formation time is 1.5 h; according to Equation (6), the breach bottom width is 42.6 m. The connection between the storage area and the dam was realized in the moment when the break line for the dam was drawn [23,73]. The next step was to create the subgrid model (mesh). The subgrid models represent a polygonal network where each cell has information about the underlying terrain. In the 2D model, each computational cell is similar to each cross-section from the 1D model, and a property table is created for each cell [25,31]. The table contains information about the roughness, the elevation-volume relationship, area, maximum elevation, minimum elevation, average elevation, etc. On the basis of these data obtained from the underlying terrain, the flow between the cells is computed. To obtain an accurate representation of the flood wave, a proper cell size must be chosen [74,75]. On the one hand, a large cell may lead to loss of details about the flow moving; however, on the other hand if the cell size is too small, this could cause computational instability [23] (Figure 3d). For the current study, the 2D flow area was created on the basis of the floodplain extent with a surface of 160 km 2 and more than 2,500,000 computational cells with an average size of 64 m 2 . The 2D flow area was directly connected with the dam. Regarding the Manning's roughness values, i.e., the study area being a typical plain with a low urbanization process, the roughness of the land-use categories was as follows: degraded land-0.03; lakes and reservoirs-0.04; shrubbery-0.04; streams-0.04; wetlands-0.07; exploitation roads-0.08; railroad-0.08; forest vegetation-0.10; orchard-0.01; unproductive land-0.01; attachment buildings-0.015; county roads-0.015; houses-0.015; local roads-0.015; main streets-0.015; national roads-0.015; secondary streets-0.015; arable land-0.035; grassland-0.035; vineyard-0.035; yards-0.035 [23,[76][77][78].
For the unsteady flow data, boundary conditions and initial conditions were set. As initial conditions, the elevation versus volume curve was used for the storage area as input. For each scenario, we assumed that the reservoir has a water volume of 17.35 million m 3 for the 0.1% recurrence interval (which corresponds to 315 m 3 /s flow rate), 12.39 million m 3 for the 0.5% recurrence interval (which corresponds to 225 m 3 /s flow rate), and 10.16 million m 3 for the 1% recurrence interval (which corresponds to 185 m 3 /s flow rate) (Figure 3e,f). For each reservoir located downstream of the R1, the initial conditions were set as a partial fill with water because these reservoirs are used for fish farming. A percentage of the total storage capacity of 28.8% was used for R2, 63.3% for R3, 41.2% for R4, and 37% for R5. As boundary conditions, a normal depth (the energy slope of the riverbed) was set for the 2D flow area [31].
The last step within the dam break modeling process was to calculate the computational time to ensure the stability of the hydraulic model. The main effect of a wrong computational time is the instability of the model [23]. Explained in more detail, if the time step is too big, in addition to the instability of the model, attenuation of the peak can occur, and, if the time step is too small, the computational period will be very long. However, in order to estimate the computational time, the Courant condition was used (Equation (7)) [23,31].
where C is the Courant number, ∆T is the time step (s), ∆x is the distance step in m (average two-dimensional cell size), V w is the flood wave speed (m/s), dQ is the change in discharge over a short time interval (Q 2 − Q 1 ), dA is the change in cross-section area over a short time interval (A 2 − A 1 ). According to Equation (7), a time step of 10 s was used to run the model.

Flood Hazard Assessment
Worldwide, flood hazard classification studies have been performed [31,45,46,62,79], where the most common hazard classification is based on the flood depth and flood extent [79]. In this study, we generated the flood hazard maps on the basis of 2D HEC-RAS modeling and Cal Alb embankment dam break multi-scenarios via the depth × velocity (D × V) raster exported from the RAS Mapper module. In order to generate the hazard severity classes, the flood D × V was classified according to AIDR [46] in six categories: H1 (D × V ≤ 0.3 m 2 /s), H2 (D × V range between >0.3 m 2 /s and ≤0.6 m 2 /s), H3 (D × V range between >0.6 m 2 /s and ≤1.2 m 2 /s), H4 (D × V range between >1.2 m 2 /s and ≤2 m 2 /s), H5 (D × V range between >2 m 2 /s and ≤4 m 2 /s), and H6 (D × V >4 m 2 /s) ( Table 4) [46].

Flood Pattern
On the basis of the 24 h breach simulation period, using the RAS Mapper interface [25], we managed to export, as raster and vector, two types of data for each break scenario (Figure 3g). To analyze the impact of the flood created by the hypothetical dam breach, depth and velocity (raster files) were needed. To analyze the extent of the flood, we exported the inundation boundary (vector file) from RAS Mapper [25,80] Figure S3, Supplementary Materials). Considering the flood extent in each dam break scenario, we can say that the multi-reservoir system proved to be effective in diminishing the flood extent in the case of two scenarios (1% and 0.5%). In the case of the third scenario (0.1%), the multi-reservoir system did not cope with the flow rate discharged upstream, and the flood extent was over 58.2% (with 1659 ha more) higher than the 1% scenario and 56.8% (with 1621 ha more) higher than the 0.5% scenario. Regarding the built-up area affected by floods, the multi-reservoir system led to the protection of 22 settlements (of which two are urban settlements) in the case of the 0.5% and 1% scenarios and just four rural settlements in the case of the 0.1% scenario.   For the first scenario (1%), the built-up area of five settlements is potentially affected. A total surface of 20.1 ha is potentially affected. In terms of land-use categories, the most potentially affected categories (within the built-up area) are the arable land (12.8 ha), wetlands (0.75 ha), grassland (1.99 ha), forest vegetation (1.78 ha, yards (1.58 ha), and buildings (0.45 ha-25 houses and 28 attachment buildings). Other categories have a potentially affected surface below 0.2 ha (e.g., streams, local roads, orchards, national roads, degraded land, shrubbery, vineyard). Referring to the affected area of the settlements, Havârna (S1) is the most affected settlement with 10.44 ha (Table 5), which is also the first settlement situated downstream of Cal Alb reservoir. A number of 25-50 inhabitants can be potentially affected by the flood wave. The second most affected settlement is Balint , i (S4) with 5.79 ha and 25-50 inhabitants potentially affected. The least affected settlement is Tătărăs , eni (S3) with 0.44 ha and fewer than five inhabitants potentially affected. Other affected built-up areas are Gârbeni (S2-2.67 ha) and Negreni (S6-0.79 ha). The potential number of affected inhabitants for these two settlements is between 5 and 50. According to the second scenario (0.5%), a built-up area of 21.74 ha (five settlements) is potentially affected and the result is similar to the first scenario (0.1%). The most affected landuse categories (within the built-up area) are arable land (13.6 ha), grassland (2.29 ha), forest vegetation (1.9 ha), yards (1.8 ha), wetlands (0.78 ha), and buildings (0.5 ha-28 houses and 32 attachment buildings). Other categories have a potentially affected surface below 0.2 ha (streams, lakes and reservoirs, local roads, county roads, orchards, degraded land, shrubbery, vineyard). Havârna (S1) remains the most affected settlement with 11.2 ha and 50-75 inhabitants potentially affected (Table 6)  In the third scenario (0.1%), a built-up area of 87.2 ha is potentially affected. A total of 23 settlements are at risk to be affected by the potential flood wave ( Table 6). The most affected land-use categories (within the built-up area) over 10 ha are arable land (44.6 ha), grassland (12.3 ha), and yards (10.5 ha). Land-use categories with a potential affected surface between 1 and 10 ha are forest vegetation (7.8 ha), buildings (3.4 ha-179 houses, 194 attachment buildings, two industrial buildings), local roads (2.3 ha), and wetlands (1.82 ha). Land-use categories with above 1 ha of potentially affected surface are streams, lakes and reservoirs, exploitation roads, county roads, national roads, orchards, secondary streets, unproductive land, degraded land, shrubbery, and vineyard ( Table 6). The most potentially affected settlements remain those located in the vicinity of the dam. Havârna (S1-25.78 ha) and Balint , i (S3-15.4 ha) settlements are the most affected with a potential number of affected inhabitants of 150-200. Other settlements that have a significant surface with high potential to be affected are Tătărăs , eni (S3-5.6 ha), Balint , i (S4-9.45 ha), S , tiubieni (S7-4.72 ha), and S , tefănes , ti (S24-4.1), with a number of potentially affected inhabitants between 5 and 100. The least affected settlements (less than 1 ha) are Niculcea (S5-0.08 ha), Sârbi (S14-0.7 ha), Moara Jorii (S18-0.01 ha), Păun (S21-0.76 ha), Năstase (S22-0.01 ha), and Stânca ((S23-0.71 ha), with a number of potentially affected inhabitants between 1 and 25. Due to the topographical location on the first terrasse of Bas , eu River and the high flood control capacity of Hănes , ti reservoir, four settlements (Slobozia-S16, Mihălăs , eni-S19, Negres , ti-S20, and Românes , ti-S27) are fully protected from the flood wave. According to all three flood scenarios developed in the R1 dam break scenario, the first four settlements (S1-Havârna; S2-Gârbeni; S3-Tătărăs , eni; S4-Balint , i) have the highest potential of being affected by the flood wave. The high flood control capacity of the multi-reservoir system is proven by the small number of settlements (four from 27 settlements) that are subject to a potential flood event.

Flood Depth
The second parameter that was calculated and generated by HEC-RAS software was flood depth, which was calculated for each cell on the basis of the underlying terrain, taking into consideration the maximum flood depth of the computational cell no matter the time when the maximum values were registered [25]. The maximum flood depth for the first scenario is 7.4 m and corresponds to the main channel of the Bas , eu River ( Figure S4, Supplementary Materials). In terms of affected buildings, 75.5% (40 buildings) of them can be potentially affected by a flood depth that does not exceed 1 m, and 24.5% (13 buildings) are potentially affected by a flood depth between 1 and 2 m (Table 7). In the case of the second scenario, the maximum depth is 7.6 m ( Figure S5, Supplementary Materials). The result is very similar to the first scenario, even in flood depth classes and affected buildings. A percentage of 71.7% of the buildings are potentially affected by a flood depth that does not exceed 1 m, 25% by a flood depth between 1 and 2 m, and, unlike the first scenario, two buildings (3.3%) are affected by a flood depth between 2 and 3 m ( Table 7). In the case of the third scenario, a maximum depth of 9.95 m ( Figure 5, Figure S6, Supplementary Materials) was registered, and 46.9% (175 buildings) of buildings are potentially affected by a flood depth that does not exceed 1 m, 22.8% (85 buildings) by a flood depth between 1 and 2 m, 15% (56 buildings) by a flood depth between 2 and 3 m, 9.4% (35 buildings) by a flood depth between 3 and 4 m, 2.7% (10 buildings) by a flood depth between 4 and 5 m, and only 3.2% (12 buildings) are potentially affected by a flood that exceed 5 m (S1-S5 settlements) ( Table 7). computed for each dam break scenario with 1% (100-year), 0.5% (500-year), and 0.1% (1000-year) recurrence intervals.

Flood Velocity
The velocity parameter was calculated in the same way as the flood depth [25]. For each computational cell, the maximum value of velocity from the 24 h simulation time was registered. For all the computed scenarios, the flood velocity exceeds 7 m/s: 7.3 m/s in the case of the first and second scenario and 9.1 m/s in the case of the third scenario ( Figure 6). The high values of the flood velocity can be attributed to the geomorphologic and vegetation conditions which lead to low roughness values, accentuated by the narrow areas of the drainage channel. In the case of the first scenario, 88.7% (47 buildings) are potentially affected by a flood velocity that does not exceed 1 m/s, and 11.3 (six buildings) are potentially affected by a flood velocity between 1 and 2 m/s (Table 8, Figure S7, Supplementary Materials). In the case of the second scenario, 86.7% (52 buildings) are potentially affected by a flood velocity that does exceed 1 m/s, and 13.3% (eight buildings) are potentially affected by a flood velocity between 1 and 2 m/s (Table 8, Figure S8, Supplementary Materials). In the case of the third scenario, 69.9% of the total buildings (257 buildings) are potentially affected by a flood velocity that does not exceed 1 m/s, 28.2% (105 buildings) are potentially affected by a flood velocity between 1 and 2 m/s, and only 2.9% (11 buildings) are potentially affected by a flood velocity between 2 and 3 m/s (Table 8, Figure S9, Supplementary Materials). Table 8. Number of buildings with a high potential to be affected and the flood wave velocity computed for each dam break scenario with 1% (100-year), 0.5% (500-year), and 0.1% (1000-year) recurrence intervals.

Discussion
Permanent and nonpermanent reservoirs have played an important role in protecting human settlements which are located in floodplains, especially rural settlements where subsistence farming is still practiced [47][48][49]. The intensification of extreme hydrological events (e.g., floods, flash floods, overtopping) and uncontrolled urbanization within flood hazard areas led to the point where the reservoirs became a key factor in defense against flood hazard [31]. The newly identified purpose of the reservoirs, in addition to those attributed to agricultural, social, and economic activities, makes reservoirs important strategic structures. In addition to the many benefits that dams bring to human communities located in the proximity of the reservoirs, their management also poses the risk that these water storages will yield to certain special conditions. The failure of an embankment dam, depending on its characteristics (e.g., surface area of the reservoir, water storage volume, dam characteristics) can generate flood waves that are significantly higher in terms of water volume release and water velocity than floods waves triggered by natural conditions [4] and, consequently, they can lead to more significant economic damages and human loses [4,5]. However, the flood and flood wave generated by the dam failure are categorized as accidental floods [2,5].
Therefore, complex and in-depth studies on accidental floods, especially those caused by embankment dam failure, are now mandatory due to the dynamics of hydroclimatic conditions and the destructive impact that such events can have on areas downstream of the reservoir. In this context, the awareness and perception of the risk that these hydrological phenomena induce to human society have led to the adoption of strategies to reduce and prevent the negative effects that may occur through their manifestation [48]. In Romania, these strategies should be implemented at national, regional, and even local level for the most important river basins and reservoirs, but they are not implemented in all the cases (e.g., R1-Cal Alb embankment dam and reservoir).
The present study complemented the national and local flood defense plans, given that the two permanent reservoirs on the Bas , eu River (R1-Cal Alb reservoir and R4-Negreni reservoir) did not have an APDF. At the same time, the importance of conducting such a study is given by the fact that, at national and regional level, no unanimously accepted methodology is developed to assess the risk of accidental floods caused by the failure of a dam. Also, it is the first study of this type applied to the Başeu river basin, i.e., the R1-R5 multi-reservoir system.
For high accuracy of the hydraulic model, a digital surface model (DSM) would have been needed. The advantage of using a DSM is that the houses, attachment buildings, road network, or other artificial structures are integrated in the hydraulic modeling. The presence of these solid objects in the DSM can cause a different flood behavior, as well as change the flood parameters (flood extent, flood depth, flood velocity) and the number of the potential affected buildings. Due to the unavailability of this dataset, the buildings and the road network were used to increase the accuracy level of the DTM. The availability of a LiDAR model with a resolution of 0.5 m made it possible to create a complex hydraulic model that simulated the maximum flood wave and its overlap with the local topographic (e.g., old drain channels, abandoned meanders). On the basis of LiDAR-derived DTM (0.5 m), a bathymetrical model (0.5 m), hydrological data (correlation between R1 inflow with water volume and water level for three different recurrence intervals), and built-up area and land use data from Bas , eu floodplain, we developed and applied a methodology which can be easily adapted also by regional authorities (e.g., PBWBA) and used within flood mitigation strategy.
In order to create the hydraulic model, a polygonal network with a computation point spacing of 7.5 m was created. For each computational cell, the LiDAR terrain characteristics were captured in the hydraulic model and a detailed elevation-volume relationship was developed by the HEC-RAS software. According to the underlying terrain process used in the HEC-RAS modeling, a large computation cell can be used without decreasing the accuracy of the model [25]. An important role in hydraulic modeling is played by the roughness coefficient. The flood parameters, especially flood velocity, depend on the roughness coefficient. Between roughness coefficient and flood velocity exists an inversely proportional relationship. For example, an increase in the roughness coefficient can cause a decrease in the flood velocity. Thus, for this study, we updated and optimized the roughness coefficient according to the spatial distribution of land-use categories, which is imperative and recommended for high-accuracy hydraulic modeling [30,31]. Moreover, by creating a complex hydraulic model and simulating a piping breach in the body of the R1-Cal Alb reservoir, the flood control capacity of the R1-R5 multi-reservoir system was assessed. Water depth levels, flood wave velocities, and flood extent were extracted and calculated on the basis of the results obtained from the RAS Mapper module. In the future, we want to redo the 0.1% (1000 year) recurrence interval scenario, while replacing the LiDAR DTM with LiDAR DSM in order to simulate the presence of artificial structures (e.g., roads, houses, attachment buildings).
The efficiency of the R1-R5 multi-reservoir system in the present study was also proven by the study [30] where the flood risk on the Bas , eu River was analyzed by creating a 1D steady flow hydraulic model. A part of the previous study was a scenario generated by the introduction of hydrotechnical works in the 1D hydraulic model, but without the dam failure. Comparing the results of the two studies, the neighboring settlements (S1) Havârna, (S2) Gârbeni, (S3) Tătărăs , eni, (S4) Balint , i, and (S6) Negreni of the R1-Cal Alb reservoir, R2-Movileni reservoir, R3-Tătărăşeni reservoir, and R4-Negreni reservoir showed the most important damages regarding the potentially flooded surface. Starting with R4-Negreni reservoir, the flood wave generated by the failure of the Cal Alb dam is diminished, and the drainage is made through the riverbed and floodplain. No builtup area represents a potentially flooded area. Downstream of R5-Hăneşti reservoir, the drainage is made exclusively through the riverbed, only in the area of (S23) Stânca, (S24) S , tefănes , ti, (S25) Bădiut , i, (S26) Bobules , ti, and (S27) Românes , ti settlements (where the altitude drops to 84 m), and the flood wave covers the floodplain to a width of about 1-2 km. None of the settlements at the confluence with the Prut River (S23-Stânca, S24-S , tefănes , ti, S25-Bădiut , i, S26-Bobules , ti, and S27-Românes , ti) present potentially flooded areas.

•
Combining 2D hydraulic modeling of the R1 dam break scenario (piping failure) with high-density LiDAR data (0.5 m spatial resolution) and local hydrological parameters (correlation between R1 inflow (m 3 /s) with water volume (m 3 ) and water level (m)) proved to be an efficient method to improve the action plan for dam failure (APDF) and flood mitigation strategy within the Başeu multi-reservoir system. • The multi-scenario approach using the inflow rate with 1% (100 year), 0.5% (500 year), and 0.1% (1000 year) recurrence intervals allowed the testing of the flood control capacity of Başeu multi-reservoir system according to R1 water volume which can cause a flood event in the case of dam failure. Accordingly, the first two scenarios (100 year and 500 year) indicate that only first four settlements (S1-Havârna, S2-Gârbeni, S3-Tătărăs , eni, S4-Balint , i) located downstream of R1 are potentially affected by floods due to the location between R1 and R2 reservoirs, and only in case of 1000-year scenario are all 27 settlements potentially affected.

•
The 2D hydraulic models were exported into a set of flood hazard parameters (e.g., flood extent, flood depth, flood velocity) and can be used for improve the flood hazard maps and answer real questions regarding the flood hazard threat at the local level in case of a dam failure.
Supplementary Materials: The following are available online at https://www.mdpi.com/2073-444 1/13/1/57/s1: Figure S1. Flood extent within built-up area of 27 settlements located downstream of the Cal Alb reservoir computed for dam break scenarios with the 1% (100 year) recurrence interval; Figure S2. Flood extent within built-up area of 27 settlements located downstream of the Cal Alb reservoir computed for dam break scenarios with the 0.5% (500 year) recurrence interval; Figure S3. Flood extent within built-up area of 27 settlements located downstream of the Cal Alb reservoir computed for dam break scenarios with the 0.1% (1000 year) recurrence interval; Figure S4. Data Availability Statement: Data sharing not applicable.

Conflicts of Interest:
The authors declare no conflict of interest. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.