Geoscientific Model Development (Aug 2022)

Multi-dimensional hydrological–hydraulic model with variational data assimilation for river networks and floodplains

  • L. Pujol,
  • P.-A. Garambois,
  • J. Monnier,
  • J. Monnier

DOI
https://doi.org/10.5194/gmd-15-6085-2022
Journal volume & issue
Vol. 15
pp. 6085 – 6113

Abstract

Read online

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 hydrometeorology and satellite data. The multi-D hydraulic model relies on the 2D shallow-water equations solved with a 1D–2D adapted single finite-volume solver. One-dimensional-like 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 1D-like–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 virtual 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 and lateral hydrographs and hydrological model parameters. This notably demonstrates the possibility for information feedback towards upstream hydrological catchments, that is, backward hydrology. A 1D-like 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, with inflow from 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 1D-like 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.