A Direct Approach for Detection of Bottom Topography in Shallow Water
Lamsahel Noureddine, Carole Rosier
TL;DR
The paper tackles the problem of reconstructing underwater bottom topography from surface measurements in a 1D open-channel flow. It develops a direct inverse approach derived from the shallow-water equations, requiring the free surface and its first two time derivatives at a single time, plus upstream discharge and bed information, and solves for the bed via a second-order finite-difference discretization with a Lipschitz stability guarantee. The method yields analytic solutions in steady flows and demonstrates high accuracy and noise robustness across subcritical, transcritical, and supercritical regimes in extensive numerical tests. This work provides a fast, data-efficient tool for bathymetry estimation with rigorous stability analysis, potentially extendable to deeper-water models and no-stagnation settings.
Abstract
We propose a fast, stable, and direct analytic method to detect underwater channel topography from surface wave measurements, based on one-dimensional shallow water equations. The technique requires knowledge of the free surface and its first two time derivatives at a single instant $t^{\star}$ above the fixed, bounded open segment of the domain. We first restructure the forward shallow water equations to obtain an inverse model in which the bottom profile is the only unknown, and then discretize this model using a second-order finite-difference scheme to infer the floor topography. We demonstrate that the approach satisfies a Lipschitz stability and is independent of the initial conditions of the forward problem. The well-posedness of this inverse model requires that, at the chosen measurement time $t^{\star}$, the discharge be strictly positive across the fixed portion of the open channel, which is automatically satisfied for steady and supercritical flows. For unsteady subcritical and transcritical flows, we derive two empirically validated sufficient conditions ensuring strict positivity after a sufficiently large time. The proposed methodology is tested on a range of scenarios, including classical benchmarks and different types of inlet discharges and bathymetries. We find that this analytic approach yields high approximation accuracy and that the bed profile reconstruction is stable under noise. In addition, the sufficient conditions are met across all tests.
