The Addicks and Barker Reservoirs, built in the forties, are located in Houston and collect precipitation and run-off from upstream areas to reduce flood risks along Buffalo Bayou to protect downtown Houston. During Hurricane Harvey (August 25 - August 30, 2017), the precipitation reached a new record of 910 mm [36.2 inches] in a 4 day period in Houston. The gates of Addicks and Barker Reservoirs were opened during the night of 27-28 August which led to major damages due to downstream flooding. Besides, non-government owned land upstream was flooded due to high water levels in the reservoirs. In this report, new design water levels for Addicks and Barker Reservoir are calculated based on inflowing discharge into the reservoirs and precipitation directly onto the reservoirs, including data of Hurricane Harvey. These calculated design water levels are compared with the critical water levels calculated based on the failure mechanisms of the dams. This study shows that the original design water level of the dams, based on the Probable Maximum Flood, are 2.83 m and 1.01 m higher than the critical water level for which failure of the dams can occur due to piping for Addicks and Barker Reservoir. However, the maximum allowed water level which is currently maintained by the United State Army Corps of Engineers, is 2.19 m and 2.46 m below the calculated critical water level. During Hurricane Harvey, these maximum allowed water levels were exceeded with 3.46 m and 1.93 m. The damage of residential properties upstream and downstream of the reservoirs are minimized based on the distribution of excess volume from the inflow of creeks and precipitation onto the reservoirs. The ratio of the amount of volume which should remain upstream of the dams and the volume discharged into the Buffalo Bayou is calculated for every considered event with its duration and return period. The ratio of Addicks Reservoir is the dominant ratio, which should be used for both reservoirs. Run-off alone already produces damage, especially for the 12h and 24h precipitation, so the Addicks and Barker Reservoirs should not release discharge into the Buffalo Bayou for small durations. For events with a longer duration, it would cause less damage to open the outlets of the reservoirs than to keep them closed. However, if the water level in the reservoir exceeds the critical water level for piping, it is advised to discharge more to the downstream area to prevent breaching of the dams. Since the critical water level is reached for approximately 25% of the events at Addicks Reservoir, mitigations against piping should be taken to improve the minimization of damage. For Barker Reservoir, the critical water level is not reached in the optimization. During big events, people living upstream will be more affected by the flooding than people living downstream since this optimization is based on the damage minimization of residential properties.

The number of fatalities due to a potential flood event is traditionally determined utilizing 'mortality functions’. Data of recent large-scale flooding in the Netherlands are not available since the Netherlands was successful in flood prevention. Therefore, only data from the last coastal flood event in 1953 with 1795 direct fatalities are available. The mortality functions are empirical relationships to provide mortality as a function of three explicit flood characteristics, namely water depth, flow velocity, and water level rise rate. Many more factors are included implicitly since the functions were derived from 1953 data. These underlying factors are thus based on the circumstances of the coastal flooding in 1953 and might not be representative anymore for future flood events elsewhere in the Netherlands. The three flood characteristics in current flood risk assessments are determined by means of coarse flood simulations. Since modern software is becoming more advanced, more detailed flood simulations are becoming possible. Therefore, the applicability of the mortality functions needs to be studied if finer model resolutions are used. This report presents the case study of river area the Bommelerwaard in which the validity of the 1953-based functions, possibilities for alternative functions, and finer model resolutions in hydrodynamic models are tested and analyzed with regards to their impact on flood fatality risk. A hydrodynamic model is developed using the new software program D-Flow Flexible Mesh which is able to apply finer resolutions at locations that require more detail. The different model resolutions that are tested are 100m and 25m, and 5m for the area close to the breach. The flood simulations with these model resolutions resulted in similar outcomes for the number of estimated fatalities in this case study. Overall, the 100m model is preferred because it is sufficiently able to indicate the dangerous locations, provide the order of magnitude of the flood characteristics, while it demands short computation times and matches the level of detail of the data of 1953. However, it is recommended to model the area around the breach (‘breach zone’) with finer model resolutions because the resulting higher local peak velocities are relevant for potential building collapse. For the areas around obstacles and underpasses, it is also recommended to use finer resolutions or to make use of 1D objects or fixed weirs. This study concluded that finer model resolutions at dangerous locations have an impact on the individual risk value of the neighbourhood, and this can have consequences for the maximum individual risk value and thus the overall safety standard of a large dike ring. Furthermore, the case study illustrated that compartment dikes have a significant impact on the local mortality because of the high water level rise rates just upstream. It is recommended to look into possibilities to reduce this high local mortality rate and hence, individual risk, for example by optimizing the location and number of compartment dikes or exploring the effects of openings in the dikes. Moreover, this study identified the discussion points in the current Dutch loss of life approach by a literature study, knowledge of recent flood events abroad, and loss of life approaches internationally. Alternative mortality functions are proposed based on the literature and analyzed through sensitivity analyses in the case study. It is recommended to substantiate and take into account the factors water arrival time, improved building characteristics, and age in the loss of life approach. Preventive evacuation is already taken into account in this approach, but in addition, water arrival time can be included by means of fleeing. This study shows that water arrival time has a great effect on the number of fatalities because some areas have relatively large arrival times and this enables inhabitants to flee the area. Emergency response is thereby of crucial importance. Also in 1953 this factor proved to be relevant. The improved building characteristics compared to 1953 are shown to have a limited impact on the absolute number of fatalities in this case study but it reduced the maximum value of the individual risk and is thus of relevance, especially for dike ring areas with large water depths (>2.1m) and high rise rates (>0.5m/h). Moreover, this study underlines the vulnerability of the elderly during flood events. Since the age distribution has shifted since 1953 and significantly more elderly are present in society nowadays, it is relevant to take this explicitly into account. This case study shows that correcting for age can have a significant impact on the number of fatalities. The impact on the individual risk is limited, but this depends on the spatial distribution of the elderly and should be further analyzed. Finally, the individual risk is sensitive to the configuration of the neighbourhoods. It is therefore also recommended to look into more robust approaches to determine the individual risk.