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.