the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Multi-scale hydraulic graph neural networks for flood modelling
Abstract. Deep learning-based surrogate models represent a powerful alternative to numerical models for speeding up flood mapping while preserving accuracy. In particular, solutions based on hydraulic-based graph neural networks (SWE-GNN) enable transferability to domains not used for training and allow including physical constraints. However, these models are limited due to four main aspects. First, they cannot model rapid differences in flow propagation speeds; secondly, they can face instabilities during training when using a large number of layers, needed for effective modelling; third, they cannot accommodate time-varying boundary conditions; and fourth, they require initial conditions from a numerical solver. To address these issues, we propose a multi-scale hydraulic-based graph neural network (mSWE-GNN) that models the flood at different resolutions and propagation speeds. We include time-varying boundary conditions via ghost cells, which enforce the solution at the domain's boundary and drop the need of a numerical solver for the initial conditions. To improve generalization over unseen meshes and reduce the data demand, we use invariance principles and make the inputs independent from coordinates' rotations. Numerical results on dike-breach floods show that the model predicts the full spatio-temporal simulation of the flood over unseen irregular meshes, topographies, and time-varying boundary conditions, with mean absolute errors in time of 0.05 m for water depths and 0.003 m2 s−1 for unit discharges. We further corroborate the mSWE-GNN in a realistic case study in The Netherlands and show generalization capabilities with only one fine-tuning sample, with mean absolute errors of 0.12 m for water depth, critical success index for a water depth threshold of 0.05 m of 87.68 %, and speed-ups of over 700 times. Overall, the approach opens up several avenues for probabilistic analyses of realistic configurations and flood scenarios.
- Preprint
(22813 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 01 Nov 2024)
-
RC1: 'Comment on egusphere-2024-2621', Anonymous Referee #1, 27 Sep 2024
reply
This paper presents a deep-learning based approach for 2D surface flood modelling that builds on a previously existing model, SWE-GNN. It proposes an alternative model that overcomes the need of a numerical solver to determine initial conditions and includes other novelties to improve the speed and generalization of the model. The paper demonstrates that the model can benefit from fine-tuning to generalize to new case studies. In my opinion, the research is interesting and is very well presented. I am therefore recommending the paper for publication, after minor revisions.
General comments
The paper proposes a model for flood modelling and is applied only to a specific case of floods; dike-breach floods. It would be nice to have a section in the discussion on how the model could apply to other kind of floods, such as fluvial and pluvial floods for example.
Specific comment
L84 – L85: For reproducibility, I would suggest clarifying how a mesh is classified as being too small and how the resolution of the fine mesh is chosen.
L93: It might be better to define what epsilon is here rather than in L104.
Eq. 3 and 4: I think that it might be useful to have the clarification of what h_di and h_si are.
L190: I would suggest replacing ‘training simulations’ by ‘training data’ here for clarity.
L193: You should explain that water level is the sum of the water depth and elevation of the cell as this is might not be straightforward.
Section 2.5: Clarify that O is the output dimension.
Table 1: To which resolution of the mesh does the edge length refer?
Section 3.1: I believe a concise explanation of how the Manning’s n coefficients are defined for the synthetic dataset would be helpful. I’m curious to know if these coefficients are spatially variable as the model might face challenges in extracting meaningful insights from them if they were spatially uniform.
L270-275: You provide a nice and clear explanation of CSI. However, you consider 2 thresholds (0.05 m and 0.3 m) while you explain in L273-275 that the flooded and non-flooded cells are distinguished from another in TP, FP and FN. However, this isn’t entirely accurate as cells with water depths below 0.3 m could still be considered as flooded. I would try to be more accurate and write for example: ‘TP are true positives, i.e. number of cells where both model and simulations predict water levels above the threshold value’. The same applies for the description of FP and FN.
L285: I would suggest adding a sentence stating that you are assuming the enhanced SWE-GNN also outperforms the other models, given that the original SWE-GNN outperforms them and despite the added modifications.
L305: Do you know why the MAE of h increases? Could it be due to error propagation throughout the simulation? If you have any insights into the cause of this increase, it might be worth adding it.
Fig. 9 and Fig. 10: Consider using different color scales as it is difficult to visualize the differences in water levels as it is. Also, could you clarify which mesh (e.g. finest mesh) is shown in the figures?
Fig. 10: It is really nice that you detail what the negative and positive values correspond to in the captions of Fig. 9 and Fig. 10. However, I would recommend revising the caption in Fig. 10 to enhance clarity, e.g. ‘positive values indicate that the model estimates later arrival times than the numerical simulation, while negative values indicate that the model predicts earlier arrival times’
Technical corrections
Throughout the manuscript (e.g., L14, L51, L223, caption of Fig.5): Uncapitalize ‘the’ in ‘The Netherlands’
L178: ‘via edges directed towards’ rather than ‘via directed edges towards’
L292: Clarify ‘is comparatively faster than the numerical model’
Throughout the manuscript: Add a spacing between the numbers and the ‘m’ for meters
L302-303 and 305: Replace the four ‘in correspondence of’ with ‘occur at’, ‘occur simultaneously to’, ‘near’ and ‘at’ respectively
Fig. 8: You might want to keep the notation of CSI0.05 and CSI0.3 like in the manuscript, i.e. add the ‘m’
Table 2 and 3: Write the units as [10-2 m] and [10-2 m2/s]
L330: Typo in ‘dataset’
L365: Typo in ‘dependent’
L393: ‘analyze’ rather than ‘analyse’ to keep consistency with the previous sentence
L466: Remove ‘the’ in ‘comes from the an increase’
Citation: https://doi.org/10.5194/egusphere-2024-2621-RC1
Data sets
Raw datasets for paper "Multi-scale hydraulic graph neural networks for flood modelling" Roberto Bentivoglio https://doi.org/10.5281/zenodo.13326595
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
297 | 61 | 66 | 424 | 2 | 1 |
- HTML: 297
- PDF: 61
- XML: 66
- Total: 424
- BibTeX: 2
- EndNote: 1
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1