the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
A Fluvial Flood Risk Model for Quantifying the Benefit of Mitigation Measures under Uncertainty
Abstract. We present a dynamic probabilistic flood risk model that addresses key challenges in the implementation of integrated flood risk management. These include the need for holistic, large-scale risk assessments that adopt a system-based perspective, and a decision-making framework based on benefit-cost analysis. The proposed model allows for the explicit simulation and dynamic coupling of the flood process components, including downstream flood wave propagation and possible dike failures, in a computationally efficient and data-sparse manner. It enables the consideration of aleatory and epistemic uncertainties in a 2-level Monte Carlo framework. By separating these uncertainties, the model supports robust risk assessments and facilitates the uncertainty-aware evaluation of the benefit of mitigation measures. The model is applied to the Bavarian Danube, demonstrating its ability to estimate the flood risk reduction potential from mitigation measures.
- Preprint
(17973 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (extended)
-
RC1: 'Comment on egusphere-2025-4875', Anonymous Referee #1, 02 Dec 2025
reply
-
AC1: 'Reply on RC1', Mara Ruf, 05 Dec 2025
reply
Dear reviewer,
Thank you for the detailed and extremely valuable feedback!
We will incorporate your comments and address the questions you raised (also those directly in the manuscript) in our final version. Regarding your general comments:
1) We will make sure to point that out in the abstract and introduction of the paper.
2) Good point — We will state the limitations more clearly in the paper. The large uncertainties introduced by this simplification are, however, accounted for in the uncertainty analysis. The rationale for this model choice is that we wanted to run a very large number of simulations (which requires runtimes far shorter (~1sec per simulation) than any hydraulic model can provide, as far as we can judge) in order to holistically quantify the influence of epistemic and aleatory uncertainties (and therein also dike breaches) and to use this model in subsequent studies, f.e.x in optimization problems. Independent of that, we found your judgement very interesting: We had considered all model components to be in a roughly comparable degree of process representation…
3) Yes, we agree and struggled with this while writing the paper. On the one hand, the uncertainty and sensitivity analysis deserves a much larger share of the manuscript, since this is one of the main contributions of our work. On the other hand, the description of the flood risk model has already become quite long, and we wanted to avoid further inflating the paper. As a consequence, we chose to illustratively focus on only three epistemic uncertainties and therewith intended to show how the model — and the hierarchical structure — can incorporate and support UA and SA, without going into all modeling choices and background assumptions. I will try to shorten the flood risk model chapter (if possible) and expand the UA section.
4) In our case study, we limited the assessment to scenarios with recurrence intervals > 100 years because the mitigation measure we assess is only activated beyond that threshold. I realize this was not clearly described in the paper: we will add a description of the detention basin to the study area section and emphasize it again in the hydrological load and results sections. If other mitigation measures are evaluated, the catalogue of flood scenarios has to, of course, be extended
5) We agree and will add a short discussion in the end. Despite the more advanced approaches pursued in more recent studies (e.g. fully coupled model chains, like the benchmark study by Sairam et al. 2021), we believe that there are use cases where the development or application of these models is not feasible or, in the case of AI approaches, more useful. The model was developed in a fairly short time frame and with limited data available. Applying more advanced (and thus computationally more demanding) approaches to problem settings like feasibility studies, cost benefit analyses or optimization problems could be, out of our perspective, challenging. For these use cases, a simpler, computationally cheaper model based on surrogates can be advantageous. We will emphasize the model’s intended use but also its limitations in the paper.
Thanks again for your feedback and comments!
Citation: https://doi.org/10.5194/egusphere-2025-4875-AC1
-
AC1: 'Reply on RC1', Mara Ruf, 05 Dec 2025
reply
-
RC2: 'Comment on egusphere-2025-4875', Anonymous Referee #2, 19 Jan 2026
reply
The paper presents a flood risk model to enable quick appraisal of flood mitigation measures. The model calculates expected benefits of a measure as a reduction in flood damage with respect to the ‘no measure’ baseline. The key feature of the model is its computational efficiency, which is obtained through an implementation that allows to pre-calculate a range of systems relationships, so that run time is faster. Such computational efficiency in turn enables a comprehensive quantification of output uncertainty and sensitivity through Monte Carlo simulations. I think the paper is potentially an interesting and valuable contribution, but it needs revision in order to clarify its novelty, scalability and transferability.
General remarks
1. The Introduction does not do enough to clarify the key focus and contribution of the study.
The description of the “three key necessary developments” on page 1-2 is very broad and potentially misleading. It raises a wide set of issues but after reading the full paper I think only the first one is addressed in this work. For example, the conclusions of point 2 and 3 that “integration into official practice is still lacking” and “practical implementation often remains fragmented across various professional sectors that focus on local flood protection planning” are certainly true but there is no follow up on them in the remainder of the paper. Point 3 on page 2 talks about “cross-border approaches” and “large spatial scales” but the proposed method is demonstrated here on a 380 km section of the Danube river in Germany, which is neither “cross-border” or necessarily “large scale”. The literature review on L. 57 onwards includes studies at very different scales, from national flood risk assessment (in the UK and Germany) to regional studies within Germany or even reach-level studies, for example the 50 km reach of the Po River of Domeneghetti et al. (2013). Hence, throughout the Introduction it is unclear what the paper focus is, in terms of spatial scale.
The statement of the study contribution (L.85-90) could be better connected to the previous literature review, so to clarify how exactly the proposed model differs or advance on the methods reviewed in previous lines 57-84. A much clearer statement is given in Sec 3.6: “The developed flood risk model differs from other methodologies, such as those presented by Vorogushyn (2008) and Domeneghetti et al. (2013), primarily in that ….” and in the Conclusions: “A key element enabling this independence is the developed flood routing method, along with segment-specific lookup tables that store the relationships between inundation elevation, floodwater volume, and damages.”. It might be useful to have these or similar statements made earlier on in the Introduction and not left at such later stages into the paper.2. Section 3 is very long and it is difficult to grasp which parts are a description of approaches that have been previously used and published by others, and which are the truly novel steps. Considering that the key contribution is the way the model is implemented (=with pre-computed look-up tables instead of dynamic coupling at runtime) it seems strange that “Sec. 3.6 - Model Implementation” takes only 8 lines out of a 12-pages “Sec 3 – Flood risk model” section! Indeed, novel ideas on how to perform pre-calculations are also spread out across Secs 3.1-3.5, but they are less obvious to find out. In summary, I suggest re-organising materials in this Section so that the key differences from previous studies and novel ideas are easier to identify.
3. The transferability and scalability of the proposed approach should be discussed. Could the approach sketched out in Figure 10 be applied to the entire Danube river or Germany-wide? What would be the challenges in terms of data availability, software implementation, run-time, etc.? Could the approach be applied in other regions of the world with different climatic and landscape characteristics? What are the critical assumptions underpinning the approach that may hinder its application elsewhere?
4. I am also confused by the organisation of materials through sections 3 and 5. “Sec. 5 - Results” only includes Figures of the aggregate loss estimates and their uncertainty and sensitivity to three selected parameters. This draws the reader’s attention to the UA&SA part of the paper, which is the least “interesting” as there is no novel contribution here (all the definitions and methods used for the UA&SA have been published before, as also acknowledged in Sec. 4). If the paper’s key contribution is the computationally-efficient implementation of the flood risk model, I would have expected the Results to mainly report on the behaviour of the model, and how its predictions compare to those of a dynamically coupled but more computationally expensive model - maybe at different river sections and under different scenarios with or without a dyke breaches. I suspect Fig. 7 and. 8 partially do this, but oddly they are embedded in Section 3 (“Flood risk model”), which somehow becomes both a description of the method and (implicitly) results section. In summary, I would suggest a better organisation of the materials and expanding results that demonstrate the “success” of the proposed approach for efficient implementation of flow routing and flood risk calculations, and less on the UA&SA results – where the authors apply standard methods instead of suggesting new ones.
Other specific points
Lines 150: why was the RCP8.5 scenario chosen, and does this means that everything that follows about hazard estimation only applies under such scenario? Please clarify.
Lines 153: “The detention basin evaluated in this study (5 is designed to be activated only when the discharges exceed the threshold of a 100-year event” Something wrong with this sentence, what does “(5” mean?
Line 204: “A breach resistance value — defined as the maximum water level the segment can withstand before failure — is randomly sampled from the fragility function”
What does this mean? How is the breach resistance value sampled from the fragility curve? UnclearLine 240 onwards: Does the approach to derive D use roughness, or not? First the authors say “This study employed a DEM with a 50x50 m grid cell resolution and average, land-use-specific roughness values based on literature data” but they do not mention how in the subsequent description, and in the end they conclude “However, this approach neglects the influence of various parameters such as terrain roughness, …” Please clarify
Line 396: “The reduction of risk associated with these factors is communicated as an additional benefit of the measures implemented, but it is not included in the formal benefit-cost analysis”
What does “communicated” mean? Are these factors quantified or assessed in any other qualitative way? If not, is this point even worth raising?Line 399: “The developed flood risk model differs from other methodologies, such as those presented by Vorogushyn (2008) and Domeneghetti et al. (2013), primarily in that the flood simulation model is executed independently of a coupling with hydraulic and hydrodynamic simulation models”
Confusing. I understand from Fig.2 that the “flood simulation model” is exactly the coupling of the hydraulic (=inundation) and hydrodynamic module, so the wording “is executed independently of a coupling with…” is confusing. Maybe the authors mean something along the lines of: “The developed flood risk model differs … primarily in that the execution of the flood simulation model does not require a dynamic coupling of the hydraulic and hydrodynamic simulation modules at runtime, but rely on reading from pre-computed lookup tables of …”?Sec. 4.1 “Xa is a set of multiple hundreds of aleatory random variables, which represent the flood scenarios, the resistance of all dike segments and the breach widths of failed dike segments”. Unclear. What is the flood scenario? How many options are there? How are the resistance and breach widths sampled (according to which distribution?) and what is the space explored (I assume proportional to the number of segments but this should be clarified also to get an idea of whether using 3,600 samples is enough or too little given the dimension of the space explored!)
Table 1: How are the epistemic uncertainties sampled? For example, I suppose that “breach width PDF” or “V-E relation dike segment” are functions, so how are they sampled? How can you represent the uncertainty in a function through a “beta distribution”? Unclear
Line 444: “This is shown in Fig. 12 for a0, no measure, and the implementation of the measure under study, aM” Unclear. I guess this result is obtained for a given realisation of the epistemic uncertainty? please clarify
Line 456: “To approximate E[Ye|Xe,i] based on the samples, a one-dimensional linear regression Sis applied” – linear assumption seems ok in the three examples provided here but may not in others, and is unnecessary – it is easy to use other non-linear approximators – see for example: https://doi.org/10.1016/j.envsoft.2015.09.011
Figure 2: add legend or info in the caption about the meaning of variables Q, W, V, E, D
Figure 10: fonts are very small and hardly readable
Figure 13: what does “No numbers are shown on the x-axis due to the preliminary nature of the results.” mean!? What is the point of showing results if they are so preliminary that their physical interpretation is still impossible/unclear?
Citation: https://doi.org/10.5194/egusphere-2025-4875-RC2 -
AC2: 'Reply on RC2', Mara Ruf, 05 Feb 2026
reply
Dear Reviewer,
Thank you for your detailed and constructive feedback! Here are our answers to your comments as well as modified sections based on your improvement suggestions.
1. The Introduction does not do enough to clarify the key focus and contribution of the study.
The description of the “three key necessary developments” on page 1-2 is very broad and potentially misleading. It raises a wide set of issues but after reading the full paper I think only the first one is addressed in this work. For example, the conclusions of point 2 and 3 that “integration into official practice is still lacking” and “practical implementation often remains fragmented across various professional sectors that focus on local flood protection planning” are certainly true but there is no follow up on them in the remainder of the paper. Point 3 on page 2 talks about “cross-border approaches” and “large spatial scales” but the proposed method is demonstrated here on a 380 km section of the Danube river in Germany, which is neither “cross-border” or necessarily “large scale”. The literature review on L. 57 onwards includes studies at very different scales, from national flood risk assessment (in the UK and Germany) to regional studies within Germany or even reach-level studies, for example the 50 km reach of the Po River of Domeneghetti et al. (2013). Hence, throughout the Introduction it is unclear what the paper focus is, in terms of spatial scale. The statement of the study contribution (L.85-90) could be better connected to the previous literature review, so to clarify how exactly the proposed model differs or advance on the methods reviewed in previous lines 57-84. A much clearer statement is given in Sec 3.6: “The developed flood risk model differs from other methodologies, such as those presented by Vorogushyn (2008) and Domeneghetti et al. (2013), primarily in that ….” and in the Conclusions: “A key element enabling this independence is the developed flood routing method, along with segment-specific lookup tables that store the relationships between inundation elevation, floodwater volume, and damages.”. It might be useful to have these or similar statements made earlier on in the Introduction and not left at such later stages into the paper.
We acknowledge that the introduction was not sufficiently regarding the actual novelty/contribution of our work. To address this, we propose to rewrite the main paragraph in the introduction, reframe some other points in the introduction and clarify the spatial scale in the conclusions. You find the proposed new introduction below. We do think that our flood risk model can be used to support all three necessary developments that we outlined in the beginning of the introduction. The results we show demonstrate that the model can be used to estimate the flood risk (plus the uncertainties involved), which is necessary for point 1), but that it can also quantify the expected benefit resulting from mitigation measures, which in turn can be used in benefit-cost analyses and decision making contexts, hence point 2). It is very helpful that you pointed out that the relation between point 3) and our model is not really clear – we agree and will amend the manuscript to clarify this point. The strong advantage of our model (in our view) is, besides the ability to perform extensive uncertainty and sensitivity analyses, that it can be applied to large rivers. It is true that the case study is not very large, but we will clarify in the manuscript that the methodology can be applied to rivers which are much longer. Together with that, we will also discuss the question of how transferable the model is to other case studies (easily), and how it could, if applied to cross-border river sections, help to increase system thinking. Furthermore, we will discuss how the computational efficiency of the model can help to identify global optimal mitigation strategies (rather than, what is currently done in practice, implement local strategies on municipal/regional level).
Updated introduction, starting from line 57:
“Indispensable for these three developments are flood-risk models that can represent interactions in flood processes, provide uncertainty-aware risk estimates, and operate at river or catchment scale. In the scientific literature, different models have been presented that meet one or more of these criteria. [summary of literature without changes, lines 58-82]
Overall, this comparative analysis of existing large-scale flood models highlights the necessary trade-off between modeling complexity, spatial extent, and comprehensiveness of the uncertainty analysis. In this work we introduce a probabilistic flood-risk model that seeks an innovative balance between these three objectives. The developed model achieves computational efficiency by replacing runtime coupling of computationally expensive hydraulic and hydrodynamic simulations with pre-processed surrogates and lookup tables, while still dynamically coupling flood-process components (Fig. 2). This allows the model to represent the key process interactions, without the prohibitive runtime associated with coupled simulation models. Crucially, the model explicitly captures downstream propagation of impacts following dike failures or activation of mitigation measures via an efficient flood-routing scheme. Computational efficiency supports embedding the risk model within a two-level Monte-Carlo framework: aleatory and epistemic uncertainties are introduced and propagated separately, thereby enabling a more comprehensive uncertainty and sensitivity analysis than is typically feasible at river scales.
The reduction in modeling accuracy compared with flood risk models coupled with high-resolution hydraulic simulations is, as we argue, acceptable for many decision making contexts. Moreover, the resulting uncertainty is explicitly quantified within the uncertainty framework. The model is applied to the Bavarian Danube to evaluate the benefit of flood mitigation measures. The successful application in the case study proves that the model can (i) deliver a large-scale probabilistic risk assessment, (ii) support benefit-cost based decision making in practice by providing uncertainty-aware estimates of risk reduction, and (iii) support the identification of jointly optimized mitigation portfolios at river scale by explicitly modeling interactions of mitigation measures and downstream effects.
The remainder of the paper is organized as follows. Section 2 introduces the study area; Section 3 outlines the flood risk model components; Section 4 describes the two-level Monte-Carlo framework; and Section 5 presents results from the application of the model to the study area. The paper concludes with a final discussion and remarks in Section 6.”
2. Section 3 is very long and it is difficult to grasp which parts are a description of approaches that have been previously used and published by others, and which are the truly novel steps. Considering that the key contribution is the way the model is implemented (=with pre-computed look-up tables instead of dynamic coupling at runtime) it seems strange that “Sec. 3.6 - Model Implementation” takes only 8 lines out of a 12-pages “Sec 3 – Flood risk model” section! Indeed, novel ideas on how to perform pre-calculations are also spread out across Secs 3.1-3.5, but they are less obvious to find out. In summary, I suggest re-organising materials in this Section so that the key differences from previous studies and novel ideas are easier to identify.
We carefully reconsidered the paper’s organization and concluded that aggregating all novelties in a single section while placing the remaining model components elsewhere does not serve the presentation of the developed model well. We do not believe that a paper must always follow that structure, particularly when the aim is to present a new flood-risk model in which many components are at least partly novel or innovative. As you noted, it is not only the routing mechanism that is new, but rather a combination of different modelling approaches. Our intention with this manuscript is to provide a comprehensive overview of the model for researchers and practitioners who need a computationally inexpensive yet effective flood-risk framework that does not depend on coupling with external simulation software. For this target audience, it is important to understand the model configuration and the interactions among modules, rather than to itemize which elements are novel. Since the individual model components are described in detail in Sections 3.1–3.5, an extensive implementation description is not necessary. Moreover, aside from the eight lines you highlight, the principal content of this subsection is Figure 10, which we consider highly informative.
3. The transferability and scalability of the proposed approach should be discussed. Could the approach sketched out in Figure 10 be applied to the entire Danube river or Germany-wide? What would be the challenges in terms of data availability, software implementation, run-time, etc.? Could the approach be applied in other regions of the world with different climatic and landscape characteristics? What are the critical assumptions underpinning the approach that may hinder its application elsewhere?
We fully accept this point. We added a new paragraph to the discussion that explicitly addresses these issues and more clearly emphasizes the model’s limitations. The new paragraph is inserted at line 519:
“The flood risk model is applied to the 380 km of Bavarian Danube in Germany to assess the potential risk reduction of mitigation measures. Computational time scales approximately linearly with the number of dike segments; an extension of the model to the entire Danube, which spans roughly 2,850 km across ten countries, is thus expected to increase runtime by no more than a factor of ten. On this basis, we argue that the model can be applied in large-scale studies and support increasing system thinking. More critical for large cross-border applications are the pre-processing steps, which entail acquiring the necessary data, models and large-scale simulation outputs from multiple stakeholders. We acknowledge that this can be challenging in practice, but this challenge will be faced by any modeling approach. The extent to which uncertainties propagate and accumulate downstream over such long river stretches would require careful examination. Applying the model to a new study area follows a similar process, provided that region-specific climatic, hydrological and hydraulic simulation models are available. The pre-processing steps are repeated to capture relevant regional characteristics.
To further support system thinking, the model can be extended to include the entire catchment, rather than representing only the main river. Although tributary contributions are accounted for in the flood scenarios, potential protection failures, local mitigation measures, and other dynamic effects along the tributaries are not explicitly modeled; this limitation should be addressed in future work. Explicitly modeling tributaries in the same manner as the main river is straightforward and, if implemented with efficient parallelization, would not substantially increase computational cost. Such an extension would enable catchment-wide optimization of mitigation measures and improve the accuracy of risk estimates.
A further limitation of the model is due to the hydrological simulation of catchment runoff, which is performed during pre-processing in the hydrological load module. While the implemented design choice permits short run times of the main flood simulation module, it does not allow the assessment of interventions influencing the hydrological flood-forming processes. Consequently, many nature-based solutions cannot be evaluated using the present model implementation.
Notwithstanding, the developed model provides a solid foundation for decision-making, such as in benefit-cost analyses, which is demonstrated in a practical example at the case study. Therein, the benefit of a flood mitigation measure, along with its associated uncertainties, is quantified. In the future, the model can support assessments of various mitigation strategies, help identify high-risk areas, or optimize the combination of mitigation measures. In this way, we demonstrate how dynamic, probabilistic, and integrated flood risk management can be implemented in practice. We hope that our work encourages and supports decision-makers in moving toward a more holistic approach to flood risk management.”
4. I am also confused by the organisation of materials through sections 3 and 5. “Sec. 5 - Results” only includes Figures of the aggregate loss estimates and their uncertainty and sensitivity to three selected parameters. This draws the reader’s attention to the UA&SA part of the paper, which is the least “interesting” as there is no novel contribution here (all the definitions and methods used for the UA&SA have been published before, as also acknowledged in Sec. 4). If the paper’s key contribution is the computationally-efficient implementation of the flood risk model, I would have expected the Results to mainly report on the behaviour of the model, and how its predictions compare to those of a dynamically coupled but more computationally expensive model - maybe at different river sections and under different scenarios with or without a dyke breaches. I suspect Fig. 7 and. 8 partially do this, but oddly they are embedded in Section 3 (“Flood risk model”), which somehow becomes both a description of the method and (implicitly) results section. In summary, I would suggest a better organisation of the materials and expanding results that demonstrate the “success” of the proposed approach for efficient implementation of flow routing and flood risk calculations, and less on the UA&SA results – where the authors apply standard methods instead of suggesting new ones.
We would welcome a comparison of our results with those from more sophisticated flood-risk models; however, developing such a model for the sole purpose of comparison is beyond the scope of the present study. While it is not possible to compare the entire model as a whole, we have evaluated all model components against hydraulic or hydrodynamic simulations, available post-event data, or other established models (some of which are presented in Section 3). With respect to uncertainty and sensitivity analysis, we believe that conducting an in-depth, two-level UA/SA represents a substantive extension beyond much of the existing literature, which typically performs UA/SA but not within a two-level framework. We agree that the Results section should be revised to more clearly demonstrate the model’s strengths and capabilities; we will make these revisions in the final manuscript. Your detailed comments will also be incorporated into the revised version.
Citation: https://doi.org/10.5194/egusphere-2025-4875-AC2
-
AC2: 'Reply on RC2', Mara Ruf, 05 Feb 2026
reply
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 304 | 143 | 29 | 476 | 19 | 16 |
- HTML: 304
- PDF: 143
- XML: 29
- Total: 476
- BibTeX: 19
- EndNote: 16
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
Review of manuscript „A Fluvial Flood Risk Model for Quantifying the Benefit of Mitigation Measures under Uncertainty” by Mara Ruf et al., submitted to NHESS
The manuscript presents an interesting extension to the existing flood risk estimation and uncertainty quantification methods, with the aim of quantifying flood risk in economic terms for large river reaches. It applies a concept of simple surrogate models to the required flood model chain (rainfall-runoff simulation, river routing, dike breach/overtopping, inundation and damage estimation). It builds on the probabilistic flood risk assessment concepts developed about 20 years ago, including an uncertainty estimation using the concept of aleatory and epistemic uncertainty. The main advancement of the work is the inclusion of dike failures in a probabilistic manner and to include this in the risk and uncertainty estimation by a2-staged Monte Carlo analysis. Moreover, the work also presents the framework for a cost-benefit analysis of large-scale flood protection measures, like flood retention basins. This is possible because the model concept enables the downstream effect of flood protection measures in a river reach.
The proposed concept is valid and the work takes large scale flood estimation a step further, mainly by including probabilistic dike failures in the risk assessment. The manuscript is properly structured and the work is well presented. Overall, I recommend the publication in NHESS, but some issues should be addressed beforehand:
Next to these more general points I made some comments in the attached annotated manuscript.