the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 11 Mar 2022
Multi-dimensional hydrological-hydraulic model with variational data assimilation for river networks and floodplains
Léo Pujol
Pierre-André Garambois
Jérôme Monnier
Abstract. This contribution presents a novel multi-dimensional (multi-D) hydraulic-hydrological numerical model with variational data assimilation capabilities. It allows multi-scale modeling over large domains, combining in situ observations with high-resolution hydro-meteorology and satellite data. The multi-D hydraulic model relies on the 2D shallow water equations solved with a 1D2D adapted single finite volume solver. 1Dlike reaches are built through meshing methods that cause the 2D solver to degenerate into 1D. They are connected to 2D portions that act as local zooms, for modeling complex flow zones such as floodplains and confluences, via 1Dlike-2D interfaces. An existing parsimonious hydrological model, GR4H, is implemented and coupled to the hydraulic model. The forward-inverse multi-D computational model is successfully validated on academic and real cases of increasing complexity, including using the second order scheme version. Assimilating multiple observations of flow signatures leads to accurate inferences of multi-variate and spatially distributed parameters among bathymetry-friction, upstream/lateral hydrographs, and hydrological model parameters. This notably demonstrates the possibility for information feedback towards upstream hydrological catchments, that is backward hydrology. A 1Dlike model of part of the Garonne river is built and accurately reproduces flow lines and propagations of a 2D reference model. A multi-D model of the complex Adour basin network, inflowed by the semi-distributed hydrological model, is built. High resolution flow simulations are obtained on a large domain, including fine zooms on floodplains, with a relatively low computational cost since the network contains mostly 1Dlike reaches. The current work constitutes an upgrade of the Dassflow computational platform. The adjoint of the whole tool chain is obtained by automatic code differentiation.
-
Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
-
Preprint
(3724 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(3724 KB) - BibTeX
- EndNote
- Final revised paper
Journal article(s) based on this preprint
Léo Pujol et al.
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2022-10', Anonymous Referee #1, 18 Apr 2022
General Comments
The manuscript presents the development of a multi-dimensional hydrological-hydraulic model. The model is calibrated/optimized using the variational data assimilation approach. Enhancing the computational efficiency of multi-dimensional river routing is important in the field of hydraulic modeling. Furthermore, it is important to investigation on state-parameter estimation utilizing the assimilation approach to identify the optimal parameters for obtaining better estimates of the physical variables of hydrodynamics. The authors describe a scientifically sound approach for performing multi-dimensional hydrologic-hydraulic modeling and calibrating model parameters using variational data assimilation. However, the manuscript lacks the reasoning and purpose for multi-dimensional hydrologic-hydraulic modeling over 2D modeling of the entire river length (a lot of commercial software are available for such 2D modelling). On the other hand, it is unclear how the combination of VDA with multi-dimensional modeling will increase the capacities of estimating physical variables. The use of several acronyms and mathematical formulae with no physical relevance impeded the intelligibility of the work considerably. Some specific and technical comments are provided to enhance the text that will be published in Geoscientific Model Development.
Specific Comments
- In the introduction, the authors did not clearly mention the motivation and the objective of the study. The authors mentioned “high resolution accuracy and fast computation times” but they introduced many studies on the matter in the next paragraph (L70-79). Authors should present more focused science questions.
- It is a bit confusing why the authors combined the multi-dimensional modelling with multi-source data assimilation methods over full 2-D modelling.
- The authors should explain the reason for using variational data assimilation method over ensemble data assimilation methods.
- It is not clear from the text in which temporal scale the parameter optimization is performed?
- What are the difference between parameter optimization (i.e., VDA) between 1D and 2D cases.
- Can the methods shown in the manuscript be applied for all spatial scales? A discussion of the spatial resolutions 1D or 2D river segments are needed for fully utilization of the methods developed in this manuscript.
- In the synthetic experiment, the authors discuss only scenarios with no lateral flow (e.g., surface and subsurface runoff) but it is better to have some discussion with lateral flow case.
- Many of the mathematical equations found in the main text are repeated in the appendices. So, I would like to suggest the authors to use the equations in the appendices to help them explain the main text more clearly. Authors can reduce the number of repeating equations by doing so.
- Is it possible for the model to modify the dimensionality on a temporal scale? When a flood occurs in one river reach, the flood is simulated using a 2D mesh, while in other cases, a 1D method is employed. If this is true, how will the model determine the flooding times?
- In section 3.3.1, the authors present a observing system simulation experiment (twin experiment) where a virtual observation is assumed. When they used virtual observations to calibrate the model parameters, they assumed the observations are available in all the river pixels and all the time. The availability of observations for all the river reaches in all the time may not be reasonable. To assess the validity of methods the authors should test more realistic scenario by assuming either spatial or temporal discontinuity.
Technical Corrections
- L60: What does “precipiton” mean
- L73 It is not easy to guess local 2D ‘zooms’. Please elaborate on it.
- L83 SWE is not defined before
- Tabel 1: what is “sources available” better to explain it in the caption
- L113: What do SW means, “shallow water”?
- L129: ]0,T] reads [0,T]
- Eq (1): doesn’t U, F, G, Sg, and Sf be introduced?
- Section 2.2.4 does not contain in any title
- L201-204: do xi and t refers to location and time, respectively?
- L211: What are Qr and Qd
- Eq(9): latter part of Eq(9) is missing.
- Figure 6: The authors can include the x, y axis in the panel (b). Is it possible to show an example of WSE variation at the 1D-2D mesh boundary? Display the legend in the center panel (blue/red lines on the hydrograph). This statement applies to all similar figures.
- L332: What is the reference discharge for calculation of NSE.
- Figure 12: Explain the area zoomed in caption, what are the two-river section shown in the focused area, explain them in the caption.
- L415: What is SMS meshing tool mean?
- L499: What is PET standard for?
- L500: What is SMASH?
Citation: https://doi.org/10.5194/egusphere-2022-10-RC1 -
RC2: 'Comment on egusphere-2022-10', Anonymous Referee #2, 19 Apr 2022
This study proposes an integrated hydrological and multi-dimensional hydraulic modeling approach that is capable of handling multi-variate optimization problems of high dimension using multi-source data. The new multi-D hydraulic computational model was coupled to the formerly developed hydrological model (GR4H) in a semi-distributed setup, called DassFlow2D-V3. The topic is of high importance to the hydraulic and hydrology community, particularly under the massively growing high-resolution data or products obtained by remote sensing. However, the novelty of the work seems to be exaggerated. The basis of the proposed hydrologic-hydraulic model was already developed by Monnier et al. (2016) and Santos et al. (2018) and as the authors acknowledge, this work presents an upgrade to the above setting with respect to a new multi-dimensional hydraulic computational model. Furthermore, although the authors have described the capabilities of similar models (L50-70), the need for and the striking advantage of the proposed framework over the competing models have not been demonstrated. The manuscript also lack the underlying science questions that need to be outlined clearly. Further technical and editorial comments are listed below to consider before manuscript can be published in Geoscientific Model Development.
Major comments
- The material presented in section 2 is hard to follow in many parts: too many acronyms and multiple cross-references to other sections not presented, yet. Also, the figures are not referred to as they appear in the order presented. The authors are invited to carefully review the manuscript for a clearer and smoother presentation.
- L215-217: It seems only some parameters were calibrated in the integrated model, but the authors should describe the reason behind this choice.
- L240: Why did not you consider the square root of the current objective functions to make their unit tangible and comparable to the unit of the estimating variables, e.g., Q?
- Section 2.4: It is not clear what the implication of “variational” is in the VDA framework.
- L334-338: You have repeated this experiment setting for at least three times in the manuscript. The same issue is seen in other parts. The authors are highly recommended to avoid repeating the same material, but with different toning, as well as to confine the results section to what are really the results. Currently, the results sections includes material related to the details of different experiments that should be stated in section 2.
- I am maybe missing something, but from the results it looks like the proposed new modeling framework does not that remarkable advantage in comparison to the formerly developed models of the same purpose. The authors should clearly highlight the distinct advantages of the proposed framework based on the reported results.
- To evaluate the accuracy and efficiency of the proposed model, the authors should expand their test cased to real-world river basins of small of medium size (< 1,000 km2) and compare the reproduced hydrographs with observations at multiple points across the basin. Currently, it is really difficult to judge about the applicability of the model as well as the relative advantages relative to the other competing models.
- The authors should also discuss the computational cost of the proposed model. Obviously, this issue would be interesting in case of real-world medium size basins, not the very simple virtual experiments considered in this study. In summary, the author should discuss if their model has the potential of application in the flood forecasting systems at a country or continent scale.
Minor comments
- Figure 1 needs improvement, and a reader cannot easily grasp the illustrated conceptual meshing approach.
- L7: virtual experiments instead of “academic”?
- L42: “a” key to…
- L84: SWE not defined before in text.
- Table 1: The citations are numbered in this table, while they are not labeled in the reference list.
- L109: VDA not defined before.
- Section 2.2.4 lacks a title.
- L230: the aim here is to … please fix. The manuscript needs several other English proofing instances that the authors should take care of them.
- Figure 6: What is the blue time series in the second panel?
Citation: https://doi.org/10.5194/egusphere-2022-10-RC2 - AC1: 'Response to reviewers', Léo Pujol, 10 Jun 2022
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2022-10', Anonymous Referee #1, 18 Apr 2022
General Comments
The manuscript presents the development of a multi-dimensional hydrological-hydraulic model. The model is calibrated/optimized using the variational data assimilation approach. Enhancing the computational efficiency of multi-dimensional river routing is important in the field of hydraulic modeling. Furthermore, it is important to investigation on state-parameter estimation utilizing the assimilation approach to identify the optimal parameters for obtaining better estimates of the physical variables of hydrodynamics. The authors describe a scientifically sound approach for performing multi-dimensional hydrologic-hydraulic modeling and calibrating model parameters using variational data assimilation. However, the manuscript lacks the reasoning and purpose for multi-dimensional hydrologic-hydraulic modeling over 2D modeling of the entire river length (a lot of commercial software are available for such 2D modelling). On the other hand, it is unclear how the combination of VDA with multi-dimensional modeling will increase the capacities of estimating physical variables. The use of several acronyms and mathematical formulae with no physical relevance impeded the intelligibility of the work considerably. Some specific and technical comments are provided to enhance the text that will be published in Geoscientific Model Development.
Specific Comments
- In the introduction, the authors did not clearly mention the motivation and the objective of the study. The authors mentioned “high resolution accuracy and fast computation times” but they introduced many studies on the matter in the next paragraph (L70-79). Authors should present more focused science questions.
- It is a bit confusing why the authors combined the multi-dimensional modelling with multi-source data assimilation methods over full 2-D modelling.
- The authors should explain the reason for using variational data assimilation method over ensemble data assimilation methods.
- It is not clear from the text in which temporal scale the parameter optimization is performed?
- What are the difference between parameter optimization (i.e., VDA) between 1D and 2D cases.
- Can the methods shown in the manuscript be applied for all spatial scales? A discussion of the spatial resolutions 1D or 2D river segments are needed for fully utilization of the methods developed in this manuscript.
- In the synthetic experiment, the authors discuss only scenarios with no lateral flow (e.g., surface and subsurface runoff) but it is better to have some discussion with lateral flow case.
- Many of the mathematical equations found in the main text are repeated in the appendices. So, I would like to suggest the authors to use the equations in the appendices to help them explain the main text more clearly. Authors can reduce the number of repeating equations by doing so.
- Is it possible for the model to modify the dimensionality on a temporal scale? When a flood occurs in one river reach, the flood is simulated using a 2D mesh, while in other cases, a 1D method is employed. If this is true, how will the model determine the flooding times?
- In section 3.3.1, the authors present a observing system simulation experiment (twin experiment) where a virtual observation is assumed. When they used virtual observations to calibrate the model parameters, they assumed the observations are available in all the river pixels and all the time. The availability of observations for all the river reaches in all the time may not be reasonable. To assess the validity of methods the authors should test more realistic scenario by assuming either spatial or temporal discontinuity.
Technical Corrections
- L60: What does “precipiton” mean
- L73 It is not easy to guess local 2D ‘zooms’. Please elaborate on it.
- L83 SWE is not defined before
- Tabel 1: what is “sources available” better to explain it in the caption
- L113: What do SW means, “shallow water”?
- L129: ]0,T] reads [0,T]
- Eq (1): doesn’t U, F, G, Sg, and Sf be introduced?
- Section 2.2.4 does not contain in any title
- L201-204: do xi and t refers to location and time, respectively?
- L211: What are Qr and Qd
- Eq(9): latter part of Eq(9) is missing.
- Figure 6: The authors can include the x, y axis in the panel (b). Is it possible to show an example of WSE variation at the 1D-2D mesh boundary? Display the legend in the center panel (blue/red lines on the hydrograph). This statement applies to all similar figures.
- L332: What is the reference discharge for calculation of NSE.
- Figure 12: Explain the area zoomed in caption, what are the two-river section shown in the focused area, explain them in the caption.
- L415: What is SMS meshing tool mean?
- L499: What is PET standard for?
- L500: What is SMASH?
Citation: https://doi.org/10.5194/egusphere-2022-10-RC1 -
RC2: 'Comment on egusphere-2022-10', Anonymous Referee #2, 19 Apr 2022
This study proposes an integrated hydrological and multi-dimensional hydraulic modeling approach that is capable of handling multi-variate optimization problems of high dimension using multi-source data. The new multi-D hydraulic computational model was coupled to the formerly developed hydrological model (GR4H) in a semi-distributed setup, called DassFlow2D-V3. The topic is of high importance to the hydraulic and hydrology community, particularly under the massively growing high-resolution data or products obtained by remote sensing. However, the novelty of the work seems to be exaggerated. The basis of the proposed hydrologic-hydraulic model was already developed by Monnier et al. (2016) and Santos et al. (2018) and as the authors acknowledge, this work presents an upgrade to the above setting with respect to a new multi-dimensional hydraulic computational model. Furthermore, although the authors have described the capabilities of similar models (L50-70), the need for and the striking advantage of the proposed framework over the competing models have not been demonstrated. The manuscript also lack the underlying science questions that need to be outlined clearly. Further technical and editorial comments are listed below to consider before manuscript can be published in Geoscientific Model Development.
Major comments
- The material presented in section 2 is hard to follow in many parts: too many acronyms and multiple cross-references to other sections not presented, yet. Also, the figures are not referred to as they appear in the order presented. The authors are invited to carefully review the manuscript for a clearer and smoother presentation.
- L215-217: It seems only some parameters were calibrated in the integrated model, but the authors should describe the reason behind this choice.
- L240: Why did not you consider the square root of the current objective functions to make their unit tangible and comparable to the unit of the estimating variables, e.g., Q?
- Section 2.4: It is not clear what the implication of “variational” is in the VDA framework.
- L334-338: You have repeated this experiment setting for at least three times in the manuscript. The same issue is seen in other parts. The authors are highly recommended to avoid repeating the same material, but with different toning, as well as to confine the results section to what are really the results. Currently, the results sections includes material related to the details of different experiments that should be stated in section 2.
- I am maybe missing something, but from the results it looks like the proposed new modeling framework does not that remarkable advantage in comparison to the formerly developed models of the same purpose. The authors should clearly highlight the distinct advantages of the proposed framework based on the reported results.
- To evaluate the accuracy and efficiency of the proposed model, the authors should expand their test cased to real-world river basins of small of medium size (< 1,000 km2) and compare the reproduced hydrographs with observations at multiple points across the basin. Currently, it is really difficult to judge about the applicability of the model as well as the relative advantages relative to the other competing models.
- The authors should also discuss the computational cost of the proposed model. Obviously, this issue would be interesting in case of real-world medium size basins, not the very simple virtual experiments considered in this study. In summary, the author should discuss if their model has the potential of application in the flood forecasting systems at a country or continent scale.
Minor comments
- Figure 1 needs improvement, and a reader cannot easily grasp the illustrated conceptual meshing approach.
- L7: virtual experiments instead of “academic”?
- L42: “a” key to…
- L84: SWE not defined before in text.
- Table 1: The citations are numbered in this table, while they are not labeled in the reference list.
- L109: VDA not defined before.
- Section 2.2.4 lacks a title.
- L230: the aim here is to … please fix. The manuscript needs several other English proofing instances that the authors should take care of them.
- Figure 6: What is the blue time series in the second panel?
Citation: https://doi.org/10.5194/egusphere-2022-10-RC2 - AC1: 'Response to reviewers', Léo Pujol, 10 Jun 2022
Peer review completion
Journal article(s) based on this preprint
Léo Pujol et al.
Léo Pujol et al.
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 439 | 135 | 17 | 591 | 7 | 6 |
- HTML: 439
- PDF: 135
- XML: 17
- Total: 591
- BibTeX: 7
- EndNote: 6
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(3724 KB) - Metadata XML