arxiv: 2602.01099 · v2 · submitted 2026-02-01 · 📊 stat.AP · cs.NA· math.NA

Recognition: 2 theorem links

· Lean Theorem

Simultaneous Estimation of Seabed and Its Roughness With Longitudinal Waves

Babak Maboudi Afkham , Ana Carpio

Authors on Pith no claims yet

Pith reviewed 2026-05-16 08:47 UTC · model grok-4.3

classification 📊 stat.AP cs.NAmath.NA

keywords acoustic tomographyBayesian inversionseabed estimationroughness quantificationwave scatteringinverse problemsstatistical isotropyfractional differentiability

0 comments

The pith

An infinite-dimensional Bayesian method uses statistical isotropy and fractional differentiability to simultaneously estimate seabed topography and roughness from acoustic wave scattering.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes a Bayesian approach to acoustic seabed tomography that estimates both the seabed profile and its roughness level at the same time. It treats the inverse problem as ill-posed but makes it solvable by assuming the seabed has statistical isotropy, which lets fractional differentiability capture the roughness. This yields a numerical algorithm that not only gives estimates but also quantifies the uncertainty in those estimates. The method is tested through extensive simulations showing it works for different seabed types. If successful in practice, it would support more reliable large-scale mapping of ocean floors for exploration and safety.

Core claim

The authors claim that by leveraging the statistical isotropy of the seabed and modeling its roughness via fractional differentiability within an infinite-dimensional Bayesian framework, it is possible to simultaneously reconstruct the seabed geometry and estimate its roughness parameter from longitudinal wave scattering data, while also providing uncertainty quantification.

What carries the argument

The infinite-dimensional Bayesian posterior that incorporates a fractional regularity prior under statistical isotropy to regularize the wave scattering inverse problem.

If this is right

The algorithm provides both seabed estimates and uncertainty bounds from wave measurements.
Extensive numerical tests confirm recovery accuracy across varied seabed configurations.
The approach renders the otherwise ill-posed tomography problem computationally tractable.
It supports large-scale seabed exploration by quantifying roughness effects.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If real seabeds satisfy statistical isotropy, this method could integrate into existing sonar systems for better mapping.
Similar Bayesian regularization might apply to other inverse problems involving rough surfaces in wave propagation.
Field tests with known anisotropic seabeds would reveal the assumption's practical limits.
Extending the framework to multi-frequency waves could improve resolution in real ocean settings.

Load-bearing premise

The seabed is statistically isotropic, permitting roughness to be identified through its fractional differentiability properties.

What would settle it

If measurements from a seabed with known non-isotropic roughness produce multiple equally likely Bayesian solutions without a clear roughness estimate, the method would be falsified.

Figures

Figures reproduced from arXiv: 2602.01099 by Ana Carpio, Babak Maboudi Afkham.

**Figure 8.** Figure 8: We notice that the trace of s, in Figure 8a, shows some correlation between the successive samples, after the warm-up phase. We report that the effective sample size [53], estimated using the ArviZ Python package [46], is approximately 10 for this chain, indicating that roughly 10 effectively independent samples can be extracted [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗

read the original abstract

This paper introduces an infinite-dimensional Bayesian framework for acoustic seabed tomography, leveraging wave scattering to simultaneously estimate the seabed and its roughness. Tomography is considered an ill-posed problem where multiple seabed configurations can result in similar measurement patterns. We propose a novel approach focusing on the statistical isotropy of the seabed. Utilizing fractional differentiability to identify seabed roughness, the paper presents a robust numerical algorithm to estimate the seabed and quantify uncertainties. Extensive numerical experiments validate the effectiveness of this method, offering a promising avenue for large-scale seabed exploration.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Bayesian joint recovery of seabed shape and roughness from waves via isotropy and fractional priors, with synthetic validation but limited robustness checks.

read the letter

The paper sets up an infinite-dimensional Bayesian framework to recover both seabed topography and its roughness parameter at once from acoustic wave scattering data. It treats roughness as a random field whose statistical isotropy and fractional differentiability supply the regularization needed to make the ill-posed tomography problem tractable, then runs a numerical algorithm that produces posterior estimates and uncertainty measures. The synthetic experiments show the method can reconstruct both quantities when the data are generated from isotropic fields. That joint formulation and the explicit use of isotropy to control roughness are the genuinely new pieces. The setup is internally consistent and avoids obvious circularity by grounding the prior in measurable wave data. The main soft spot is that the isotropy assumption carries most of the load. If real seabeds deviate even modestly from isotropy, the posterior may spread or concentrate incorrectly, and the reported experiments do not appear to test that case. Without those checks the uncertainty quantification remains plausible only under ideal conditions. This is useful reading for people already working on Bayesian inverse problems for wave equations or geophysical tomography. A reader looking for a concrete way to handle joint shape-and-roughness recovery would find the formulation worth examining. It is worth sending to peer review because the core idea is fresh and the numerical evidence, while narrow, is reproducible enough to merit referee scrutiny.

Referee Report

2 major / 1 minor

Summary. This paper introduces an infinite-dimensional Bayesian framework for acoustic seabed tomography that leverages wave scattering measurements to simultaneously estimate the seabed profile and its roughness. The approach addresses the ill-posed nature of the inverse problem by assuming statistical isotropy of the seabed and employing fractional differentiability to characterize roughness, with a numerical algorithm proposed and validated through extensive numerical experiments.

Significance. If the central claims hold, the work would contribute a Bayesian method for uncertainty quantification in infinite-dimensional seabed tomography, potentially aiding large-scale marine exploration. The combination of wave scattering data with fractional regularity under isotropy offers a way to regularize non-uniqueness, though its practical impact depends on validation beyond synthetic isotropic cases.

major comments (2)

[Framework and assumptions] The statistical isotropy assumption is load-bearing for tractability via fractional differentiability, yet the manuscript provides no explicit tests showing posterior concentration when isotropy holds or how uncertainties degrade under mild violations (common in real seabeds). This directly affects the robustness claim for the algorithm.
[Numerical experiments] Numerical experiments are described as validating effectiveness, but if limited to synthetic isotropic fields, they do not address whether the method remains reliable when the isotropy assumption is violated, weakening support for the central claim of a robust algorithm with reliable uncertainty quantification.

minor comments (1)

[Abstract] The abstract would benefit from inclusion of specific performance metrics (e.g., error norms or posterior variance values) from the experiments to substantiate the 'robust' descriptor.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address the major points regarding the isotropy assumption and numerical experiments below, with planned revisions to clarify scope and limitations.

read point-by-point responses

Referee: [Framework and assumptions] The statistical isotropy assumption is load-bearing for tractability via fractional differentiability, yet the manuscript provides no explicit tests showing posterior concentration when isotropy holds or how uncertainties degrade under mild violations (common in real seabeds). This directly affects the robustness claim for the algorithm.

Authors: We agree that the isotropy assumption is central to enabling fractional differentiability in the infinite-dimensional Bayesian setting. The framework is derived and analyzed under this assumption, with numerical results demonstrating posterior concentration and uncertainty quantification when isotropy holds. The manuscript does not claim robustness to violations, which are indeed common in real seabeds. We will revise the text to explicitly qualify all claims as holding under the isotropy assumption and add a dedicated paragraph discussing potential degradation under mild anisotropy as a limitation and avenue for future work. revision: partial
Referee: [Numerical experiments] Numerical experiments are described as validating effectiveness, but if limited to synthetic isotropic fields, they do not address whether the method remains reliable when the isotropy assumption is violated, weakening support for the central claim of a robust algorithm with reliable uncertainty quantification.

Authors: The experiments are intentionally restricted to synthetic data generated from isotropic fields to validate the method within its modeling assumptions. This provides evidence for effectiveness and uncertainty quantification under the stated conditions. We acknowledge that tests under violated isotropy would be needed to support broader robustness claims. We will revise the abstract, results, and conclusions sections to remove unqualified references to 'robustness' and instead specify validation under isotropy, while noting the limitation for non-isotropic cases. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation relies on external wave scattering data and standard Bayesian updating

full rationale

The paper's central framework applies an infinite-dimensional Bayesian inversion to acoustic wave scattering measurements, using the statistical isotropy assumption to enable fractional differentiability as a roughness regularizer. No self-definitional reductions, fitted inputs renamed as predictions, or load-bearing self-citations appear in the derivation chain. The approach treats isotropy as an explicit modeling assumption rather than deriving it from the target estimates, and numerical experiments are described as validation on synthetic data without tautological closure. The claim remains independent of its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the domain assumption of statistical isotropy to enable roughness characterization via fractional differentiability; no free parameters or invented entities are explicitly introduced in the abstract.

axioms (1)

domain assumption The seabed exhibits statistical isotropy.
Invoked to identify seabed roughness using fractional differentiability and to regularize the ill-posed inverse problem.

pith-pipeline@v0.9.0 · 5380 in / 1041 out tokens · 46017 ms · 2026-05-16T08:47:50.985233+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We identify the roughness of an isotropic seabed by means of its level of fractional differentiability... η ∼ N(m, C_s) ... λ_j = ((1/ℓ)^2 + j^2)^{-2s}
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

statistical isotropy of the seabed... fractional differentiability to identify seabed roughness

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages · 2 internal anchors

[1]

Abugattas, A

[1]C. Abugattas, A. Carpio, E. Cebri ´an, and G. Oleaga,Quantifying uncertainty in inverse scattering problems set in layered environments, Applied Mathematics and Computation, 500 (2025), p. 129453, https://doi.org/https://doi.org/10.1016/j.amc.2025.129453. [2]R. A. Adams,Sobolev Spaces, (1975). [3]R. A. Adams and J. J. Fournier,Sobolev spaces, vol. 140,...

work page doi:10.1016/j.amc.2025.129453 2025
[2]

ESTIMATION OF SEABED AND ROUGHNESS WITH LONGITUDINAL WAVES29 [4]B. M. Afkham, K. Knudsen, A. K. Rasmussen, and T. Tarvainen,A bayesian approach for consistent reconstruction of inclusions, Inverse Problems, 40 (2024), p. 045004. [5]B. M. Afkham, N. A. B. Riis, Y. Dong, and P. C. Hansen,Inferring object boundaries and their roughness with uncertainty quant...

work page doi:10.1007/s10851-024-01207-9 2024
[3]

Bonnel, A

[9]J. Bonnel, A. Vardi, J. Leonard, and S. Dosso,From geoacoustic inversion to seabed tomog- raphy using a distributed network of sources and receivers, The Journal of the Acoustical Society of America, 157 (2025), pp. A196–A196, https://doi.org/10.1121/10.0037871. [10]L. Borcea, J. Garnier, A. V. Mamonov, and J. Zimmerling,Reduced order model approach fo...

work page doi:10.1121/10.0037871 2025
[4]

Carpio, E

[15]A. Carpio, E. Cebri ´an, and A. Guti ´errez,Object based Bayesian full-waveform inversion for shear elastography, Inverse Problems, 39 (2023), p. 075007, https://doi.org/10.1088/ 1361-6420/acd5f8, https://dx.doi.org/10.1088/1361-6420/acd5f8. [16]A. Carpio and M.-L. Rap ´un,Multifrequency topological derivative approach to inverse scat- tering problems...

work page doi:10.1088/1361-6420/acd5f8 2023
[5]

doi: 10.1007/978-3-319-12385-1

[20]M. Dashti and A. M. Stuart,The Bayesian Approach to Inverse Problems, Springer Inter- national Publishing, 2017, pp. 311–428, https://doi.org/10.1007/978-3-319-12385-1

work page doi:10.1007/978-3-319-12385-1 2017
[6]

[21]S. E. Dosso and C. W. Holland,Bayesian inversion of seabed reflection data, in Acoustic Sensing Techniques for the Shallow Water Environment, A. Caiti, N. R. Chapman, J.- P. Hermand, and S. M. Jesus, eds., Springer, Dordrecht, 2006, https://doi.org/10.1007/ 978-1-4020-4386-4

work page 2006
[7]

Duane, A

[22]S. Duane, A. D. Kennedy, B. J. Pendleton, and D. Roweth,Hybrid monte carlo, Physics letters B, 195 (1987), pp. 216–222. [23]R. M. Dudley,Real Analysis and Probability, CRC Press,

work page 1987
[8]

[24]M. M. Dunlop, M. A. Iglesias, and A. M. Stuart,Hierarchical Bayesian level set inversion, Statistics and Computing, 27 (2017), pp. 1555–1584. [25]L. Edsberg,Introduction To Computation and Modeling For Differential Equations, John Wiley & Sons, 2 ed.,

work page 2017
[9]

[26]J. Fan, I. Gijbels, T.-C. Hu, and L.-S. Huang,A study of variable bandwidth selection for local polynomial regression, Statistica Sinica, (1996), pp. 113–127. [27]J. Feldman,Equivalence and perpendicularity of Gaussian processes, Pacific J. Math, 8 (1958), pp. 699–708. [28]N. P. Fofonoff and R. Millard Jr,Algorithms for the computation of fundamental ...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/670067 1996
[10]

30B. M. AFKHAM, A. CARPIO [32]Z. Gelbaum,White noise representation of gaussian random fields, arXiv preprint arXiv:1201.5635, (2012). [33]J. Goodman and J. Weare,Ensemble samplers with affine invariance, Communications in applied mathematics and computational science, 5 (2010), pp. 65–80. [34]P. Grisvard,Elliptic problems in nonsmooth domains, (1985). [3...

work page internal anchor Pith review Pith/arXiv arXiv 2012
[11]

[38]M. A. Iglesias, Y. Lu, and A. M. Stuart,A Bayesian level set method for geometric inverse problems, Interfaces and free boundaries, 18 (2016), pp. 181–217. [39]D. Jackson and M. Richardson,High-frequency seafloor acoustics, Springer Science & Busi- ness Media,

work page 2016
[12]

[40]F. B. Jensen, W. A. Kuperman, M. B. Porter, H. Schmidt, and A. Tolstoy,Computa- tional ocean acoustics, vol. 2011, Springer,

work page 2011
[13]

Jenserud and S

[41]T. Jenserud and S. Ivansson,Measurements and modeling of effects of out-of-plane rever- beration on the power delay profile for underwater acoustic channels, IEEE Journal of Oceanic Engineering, 40 (2015), pp. 807–821. [42]H. Jia and X. Li,Underwater reverberation suppression based on non-negative matrix factori- sation, Journal of Sound and Vibration...

work page 2015
[14]

Karamehmedovi ´c, P.-E

[44]M. Karamehmedovi ´c, P.-E. Hansen, and T. Wriedt,An efficient rough-interface scattering model for embedded nano-structures, Thin Solid Films, 541 (2013), pp. 51–56. [45]D. Koller and N. Friedman,Probabilistic Graphical Models: Principles and Techniques, MIT press,

work page 2013
[15]

Kumar, C

[46]R. Kumar, C. Carroll, A. Hartikainen, and O. Martin,Arviz a unified library for ex- ploratory analysis of bayesian models in python, Journal of Open Source Software, 4 (2019), p. 1143, https://doi.org/10.21105/joss.01143, https://doi.org/10.21105/joss.01143. [47]L. D. Landau, L. Pitaevskii, A. M. Kosevich, and E. M. Lifshitz,Theory of elasticity: volu...

work page doi:10.21105/joss.01143 2019
[16]

[49]J. L. Lions and E. Magenes,Problemes aux limites non-homogenes et applications, 1 (1968). [50]X. Lurton, G. Lamarche, C. Brown, V. Lucieer, G. RIce, A. Schimel, and T. Weber, Backscatter measurements by seafloor-mapping sonars, Guidelines and recommendations, 200 (2015). [51]A. Quarteroni,Numerical Models for Differential Problems, Springer Milan, 201...

work page doi:10.1007/978-88-470-5522-3 1968
[17]

[54]G. O. Roberts and J. S. Rosenthal,Optimal scaling of discrete approximations to Langevin diffusions, Journal of the Royal Statistical Society: Series B (Statistical Methodology), 60 (1998), pp. 255–268. [55]L. Roininen, J. M. Huttunen, and S. Lasanen,Whittle-Mat´ ern priors for Bayesian statis- tical inversion with applications in electrical impedance...

work page 1998
[18]

Schr ¨oder, A

[57]S. Schr ¨oder, A. Duparr´e, L. Coriand, A. T ¨unnermann, D. H. Penalver, and J. E. Har- vey,Modeling of light scattering in different regimes of surface roughness, Optics express, 19 (2011), pp. 9820–9835. [58]C. Sherlock and G. Roberts,Optimal scaling of the random walk metropolis on elliptically symmetric unimodal targets,

work page 2011
[19]

[59]A. M. Stuart,Inverse problems: a Bayesian perspective, Acta numerica, 19 (2010), pp. 451–

work page 2010
[20]

Suuronen, N

[60]J. Suuronen, N. K. Chada, and L. Roininen,Cauchy Markov random field priors for Bayesian inversion, Statistics and computing, 32 (2022), p

work page 2022
[21]

Tromp,Seismic wavefield imaging of earth’s interior across scales, Nature Reviews Earth & Environment, 1 (2020), pp

[62]J. Tromp,Seismic wavefield imaging of earth’s interior across scales, Nature Reviews Earth & Environment, 1 (2020), pp. 40–53. [63]C. Tsogka and G. C. Papanicolaou,Time reversal through a solid–liquid interface and super-resolution, Inverse problems, 18 (2002), p

work page 2020
[22]

[64]P. C. Valentine, B. J. Todd, and V. E. Kostylev,Classification of marine sublittoral habitats, with application to the northeastern north application to the northeastern north America region, in American Fisheries Society Symposium, vol. 41, 2005, pp. 183–200. [65]H. Wang, W. Chen, X. Li, Q. Liang, X. Qin, and J. Li,CUG-STCN: A seabed to- pography cla...

work page doi:10.1016/j.jag.2025.104383 2005