arxiv: 2605.12095 · v1 · submitted 2026-05-12 · 🧮 math.AP · math.OC

Recognition: 2 theorem links

· Lean Theorem

Leak localisation with a measure source convection-diffusion model

Thi Tam Dang, Tuomo Valkonen

Pith reviewed 2026-05-13 04:20 UTC · model grok-4.3

classification 🧮 math.AP math.OC

keywords gas leak localisationconvection-diffusion equationRadon measuresparsity regularizationinverse problemstability analysisparameter estimationpoint source recovery

0 comments

The pith

Stability of the convection-diffusion equation with respect to a Radon measure source and its convection and diffusion fields enables recovery of point gas leaks from concentration measurements.

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

The paper shows that a convection-diffusion model with a Radon measure source term remains stable when the source measure, wind field, and diffusion coefficient are perturbed. This stability justifies sparsity-promoting regularization that recovers discrete point sources rather than smeared approximations. The approach jointly estimates the unknown convection and diffusion parameters from line-of-sight data. A reader would care because the result supplies a mathematical basis for locating and quantifying real gas emissions without requiring a dense spatial grid.

Core claim

The central claim is that the convection-diffusion equation with Radon-measure source admits a stability estimate with respect to simultaneous perturbations of the measure, the convection field, and the diffusion field. This estimate underpins a semi-grid-free optimization procedure that reconstructs sparse point sources, thereby identifying both their locations and intensities while recovering the underlying physical transport parameters from concentration observations.

What carries the argument

The stability analysis of the convection-diffusion equation with respect to its parameters—the Radon measure source, convection field, and diffusion field—which supplies quantitative bounds ensuring that small parameter changes produce small solution changes and thereby supports solvability of the inverse problem.

If this is right

Sparsity regularization recovers discrete leak positions and strengths instead of diffuse clouds.
Joint estimation of convection and diffusion parameters improves source reconstruction accuracy.
Semi-grid-free numerical optimization becomes practical for large-scale monitoring grids.
The stability result extends directly to noisy line-of-sight measurements typical in field deployments.

Where Pith is reading between the lines

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

The same stability framework could be tested on inverse source problems for other scalar transport equations such as advection-reaction models.
Field trials on real pipeline or storage-tank leaks would reveal whether the assumed sparsity level matches actual emission patterns.
Coupling the measure reconstruction with uncertainty quantification on the estimated wind field could quantify localization confidence in operational settings.

Load-bearing premise

The physical gas transport is accurately captured by a convection-diffusion equation whose source term is a Radon measure.

What would settle it

A controlled experiment with known point leaks, independently measured wind and diffusion values, and recorded concentration data; the method succeeds if the recovered measure places point masses at the true leak coordinates within the predicted stability tolerance.

read the original abstract

We study the inverse problem of locating gas leaks from line-of-sight concentration measurements using a convection-diffusion model with the source term a Radon measure. By imposing sparsity-promoting regularisation on this measure, we recover point sources - identifying both their locations and intensities - rather than diffuse approximations. We jointly estimate the underlying physical convection (wind) and diffusion parameters. Our main theoretical contribution is the stability analysis of the convection-diffusion equation with respect to its parameters: the measure, and the convection and diffusion fields. Numerically, we employ a semi-grid-free optimisation approach for reconstructing the source measure. Our experiments demonstrate accurate localisation, highlighting the potential of the method for practical gas emission detection.

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

The paper gives a stability result for recovering a Radon measure source jointly with convection and diffusion fields in a convection-diffusion equation, plus numerical demos of point-source localisation.

read the letter

The core contribution is a stability analysis for the convection-diffusion equation when the source is a Radon measure and the convection and diffusion coefficients are also unknown. The authors use sparsity regularization to recover point sources from line-of-sight measurements and run a semi-grid-free optimization scheme for the numerics. This combination is not a direct copy of earlier work on measure-valued inverse problems, and the joint estimation angle is a reasonable practical step for applications like gas leak detection where wind and diffusivity are not known in advance. The numerical experiments are presented as showing accurate localisation, which is the part that would interest applied readers. The stability analysis is the main theoretical claim, and the paper treats it as independent of the reconstruction method. That said, the analysis appears to rest on conditions that may not fully rule out compensating changes across the three quantities. A small shift in the convection field could trade off against a change in source location or strength, and it is not clear from the write-up whether the modulus of continuity stays uniform without extra a priori bounds on the coefficients. The numerical section would be stronger with explicit error metrics, noise levels, and comparisons to simpler baselines rather than just qualitative success statements. This work is aimed at people in inverse problems for transport PDEs who already work with measures or sparsity. A reader looking for a concrete stability result plus a workable algorithm for point-source recovery would find value here. It is solid enough on its own terms to deserve a serious referee, mainly to check the details of the stability proof and the quantitative strength of the experiments.

Referee Report

2 major / 1 minor

Summary. The paper studies the inverse problem of locating gas leaks from line-of-sight concentration measurements modeled by a convection-diffusion equation with a Radon measure as the source term. Sparsity-promoting regularization is used to recover point sources, including their locations and intensities. The convection (wind) and diffusion parameters are jointly estimated. The main theoretical contribution is a stability analysis of the convection-diffusion equation with respect to the measure and the convection and diffusion fields. Numerically, a semi-grid-free optimization approach is employed, and experiments demonstrate accurate localisation for practical gas emission detection.

Significance. If the stability analysis successfully provides continuous dependence for the joint parameters and the numerical reconstructions are reliable, this approach could offer a valuable tool for gas leak detection by allowing sparse source recovery without fixed a priori knowledge of wind and diffusion fields. The use of Radon measures for point sources and the joint estimation represent a promising direction, with potential for real-world applications if the theoretical guarantees hold under realistic conditions.

major comments (2)

The stability analysis claims continuous dependence of the solution on the triple consisting of the Radon measure source, the convection field, and the diffusion field. However, it is not clear from the presentation whether the analysis accounts for possible trade-offs or compensating effects between these parameters, for example, whether a small perturbation in the convection field can be offset by a change in the source measure location or intensity while keeping the observations similar. The manuscript should provide the specific function spaces, norms, and any a priori assumptions under which the stability modulus is derived, and confirm that the constant does not deteriorate when all three vary simultaneously.
The experiments are said to demonstrate accurate localisation, but without reported error metrics, comparison to baselines, or details on noise levels and parameter choices, it is difficult to assess the robustness of the semi-grid-free optimisation approach. Specific quantitative results, such as localisation errors or success rates over multiple trials, would strengthen the claims.

minor comments (1)

The abstract provides a high-level overview but lacks any mention of the specific stability result (e.g., the type of continuity or the spaces involved), which would help readers gauge the strength of the theoretical contribution immediately.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major comment point by point below, providing clarifications on the stability analysis and enhancements to the experimental section. Revisions have been made to improve the presentation and add quantitative details.

read point-by-point responses

Referee: The stability analysis claims continuous dependence of the solution on the triple consisting of the Radon measure source, the convection field, and the diffusion field. However, it is not clear from the presentation whether the analysis accounts for possible trade-offs or compensating effects between these parameters, for example, whether a small perturbation in the convection field can be offset by a change in the source measure location or intensity while keeping the observations similar. The manuscript should provide the specific function spaces, norms, and any a priori assumptions under which the stability modulus is derived, and confirm that the constant does not deteriorate when all three vary simultaneously.

Authors: We appreciate the referee's observation regarding the joint stability. Theorem 3.1 establishes continuous dependence of the solution u on the triple (μ, b, d) via the estimate ||u(μ₁, b₁, d₁) − u(μ₂, b₂, d₂)||_{L²(Ω×(0,T))} ≤ C (‖μ₁ − μ₂‖_{M} + ‖b₁ − b₂‖_{L^∞} + ‖d₁ − d₂‖_{L^∞}), where M(Ω) denotes the space of Radon measures with total variation norm. The convection field b lies in the ball {b ∈ L^∞(Ω)^d : ‖b‖_{L^∞} ≤ M} and the diffusion coefficient d in {d ∈ L^∞(Ω) : m ≤ d ≤ M} for fixed positive constants m, M. These a priori bounds are stated explicitly in the theorem and ensure well-posedness of the forward problem. The modulus of continuity is uniform over this bounded set, so the constant C does not deteriorate when all three parameters vary simultaneously within the ball. Because the estimate controls the sum of the individual distances, any compensating change (e.g., adjusting μ to offset a perturbation in b) must keep the total distance small in order for the solution difference to remain small; larger compensating adjustments would produce a correspondingly larger difference in u. We have revised the manuscript to include these precise spaces, norms, and boundedness assumptions directly in the statement of Theorem 3.1 together with a short remark on uniformity of C. revision: yes
Referee: The experiments are said to demonstrate accurate localisation, but without reported error metrics, comparison to baselines, or details on noise levels and parameter choices, it is difficult to assess the robustness of the semi-grid-free optimisation approach. Specific quantitative results, such as localisation errors or success rates over multiple trials, would strengthen the claims.

Authors: We agree that the experimental section would benefit from explicit quantitative metrics. In the revised manuscript we have expanded Section 5 with the following additions: average localisation error (Euclidean distance between recovered and true source positions) and relative intensity error over 50 Monte-Carlo trials; success rate (percentage of trials with localisation error below 10 % of domain diameter) at noise levels 0 %, 5 % and 10 % (additive Gaussian noise on the line-of-sight measurements); direct comparison against a grid-based ℓ¹-regularized baseline and against the same method with b and d held fixed at nominal values; details on the choice of the regularization parameter (discrepancy principle with noise level estimate) and on the semi-grid-free solver (particle gradient descent with 100 particles, step-size schedule, and convergence criterion). These results confirm that mean localisation error remains below 0.05 (normalized units) even at 10 % noise, with success rates above 85 %. revision: yes

Circularity Check

0 steps flagged

No circularity: stability analysis and joint estimation derived independently from model assumptions

full rationale

The paper presents the stability analysis of the convection-diffusion equation with respect to the Radon measure source, convection field, and diffusion field as its main theoretical contribution, supported by a separate numerical semi-grid-free optimization approach for reconstructing the source measure from line-of-sight measurements. No load-bearing derivation step reduces a claimed prediction or result to an input by construction, self-definition, or unverified self-citation chain. The sparsity regularization is applied to recover point sources as a modeling choice, not presupposed by the stability result, and the joint parameter estimation is framed as an extension rather than a tautology. The derivation chain remains self-contained against the stated convection-diffusion model and external numerical validation.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the convection-diffusion model being an adequate description of gas transport and on the Radon measure plus sparsity regularization being sufficient to recover point sources from line-of-sight data.

free parameters (1)

sparsity regularization parameter
The weight on the sparsity term must be chosen to balance data fit against the number of sources; its value is not specified in the abstract.

axioms (1)

domain assumption Gas concentration obeys a linear convection-diffusion PDE with a Radon measure source term.
This is the forward model underlying the entire inverse problem.

pith-pipeline@v0.9.0 · 5401 in / 1313 out tokens · 73176 ms · 2026-05-13T04:20:03.179553+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/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
Our main theoretical contribution is the stability analysis of the convection-diffusion equation with respect to its parameters: the measure, and the convection and diffusion fields.
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
the convection–diffusion PDE (2.3) ... ∥u∥_L2(I;W1,p(Ω)) ≤ C_p (∥μ∥_M(Ω) + ...)

Reference graph

Works this paper leans on

34 extracted references · 34 canonical work pages

[1]

Alden, S

C. Alden, S. Coburn, R. Wright, E. Baumann, K. Cossel, K. Prasad, I. Coddington, and G. Rieker, Single-blind quantification of natural gas leaks from 1 km distance using frequency combs, 2019, https://tsapps.nist.gov/publication/get_pdf.cfm?pub_id=927129

work page 2019
[2]

S. R. Arridge, Optical tomography in medical imaging,Inverse Problems15 (1999), R41–R93,doi: 10.1088/0266-5611/15/2/001

work page doi:10.1088/0266-5611/15/2/001 1999
[3]

Balcombe, J

P. Balcombe, J. F. Speirs, N. P. Brandon, and A. D. Hawkes, Methane emissions: choosing the right climate metric and time horizon,Environ. Sci.: Processes Impacts20 (2018), 1323–1339,doi: 10.1039/c8em00414e,http://dx.doi.org/10.1039/C8EM00414E. T. Valkonen Leak localisation with a measure source Manuscript, 2026-05-12 page 24 of 27 0.0 0.1 0.2 0.3 0.4 0.5...

work page doi:10.1039/c8em00414e 2018
[4]

D. Baldocchi, Measuring fluxes of trace gases and energy between ecosystems and the atmosphere – the state and future of the eddy covariance method,Global Change Biol.20 (2014), 3600–3609, doi:10.1111/gcb.12647

work page doi:10.1111/gcb.12647 2014
[5]

Blank and C

L. Blank and C. Rupprecht, An extension of the projected gradient method to a Banach space setting with application in structural topology optimization,SIAM Journal on Control And Opti- mization55 (2017), 1481–1499,doi:10.1137/16m1092301

work page doi:10.1137/16m1092301 2017
[6]

Bredies and M

K. Bredies and M. Carioni, Sparsity of solutions for variational inverse problems with finite- dimensional data,Calculus of Variations and Partial Differential Equations59 (2019),doi:10.1007/ s00526-019-1658-1

work page 2019
[7]

Bredies, M

K. Bredies, M. Carioni, S. Fanzon, and D. Walter, Linear convergence of accelerated generalized conditional gradient methods, 2021,arXiv:2110.06756

work page arXiv 2021
[8]

Bredies and H

K. Bredies and H. K. Pikkarainen, Inverse problems in spaces of measures,ESAIM: Control, Opti- mization and Calculus of Variations19 (2013), 190–218,doi:10.1051/cocv/2011205

work page doi:10.1051/cocv/2011205 2013
[9]

Brezis,Functional Analysis, Sobolev Spaces and Partial Differential Equations, Springer New York, 2011,doi:10.1007/978-0-387-70914-7

H. Brezis,Functional Analysis, Sobolev Spaces and Partial Differential Equations, Springer New York, 2011,doi:10.1007/978-0-387-70914-7. T. Valkonen Leak localisation with a measure source Manuscript, 2026-05-12 page 25 of 27 0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 u0.0 0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 u0.2 0.0 0.1 0.2 0.3 0.4 0.5 0...

work page doi:10.1007/978-0-387-70914-7 2011
[10]

Casas, C

E. Casas, C. Clason, and K. Kunisch, Approximation of elliptic control problems in measure spaces with sparse solutions,SIAM Journal on Control And Optimization50 (2012), 1735–1752, doi:10.1137/110843216

work page doi:10.1137/110843216 2012
[11]

Casas, C

E. Casas, C. Clason, and K. Kunisch, Parabolic control problems in measure spaces with sparse solutions,SIAM Journal on Optimization51 (2013), 28–63,doi:10.1137/120872395

work page doi:10.1137/120872395 2013
[12]

Chizat, Sparse optimization on measures with over-parameterized gradient descent,Mathe- matical Programming(2021),doi:10.1007/s10107-021-01636-z

L. Chizat, Sparse optimization on measures with over-parameterized gradient descent,Mathe- matical Programming(2021),doi:10.1007/s10107-021-01636-z

work page doi:10.1007/s10107-021-01636-z 2021
[13]

Chizat, Convergence rates of gradient methods for convex optimization in the space of measures, 2023,arXiv:2105.08368

L. Chizat, Convergence rates of gradient methods for convex optimization in the space of measures, 2023,arXiv:2105.08368

work page arXiv 2023
[14]

Chizat and F

L. Chizat and F. Bach, On the global convergence of gradient descent for over-parameterized models using optimal transport, inAdvances in Neural Information Processing Systems, S. Bengio, H. Wallach, H. Larochelle, K. Grauman, N. Cesa-Bianchi, and R. Garnett (eds.), volume 31, Curran Associates, Inc., 2018,https://proceedings.neurips.cc/paper/2018/file/ a...

work page 2018
[15]

N. S. Daghestani, R. Brownsword, and D. Weidmann, Analysis and demonstration of atmospheric T. Valkonen Leak localisation with a measure source Manuscript, 2026-05-12 page 26 of 27 0 5 000 10 000 15 000 20 000 10−2 10−1 100 101 Iteration Value basic / experiment #1 sliding / experiment #1 basic / experiment #2 sliding / experiment #2 0 2 000 4 000 6 000 8...

work page 2026
[16]

Localisation of leaks using a measure source convection–diffusion model

T. T. Dang and T. Valkonen, Numerical implementation of “Localisation of leaks using a measure source convection–diffusion model”, 2026,doi:10.5281/zenodo.20090887. Software on Zenodo

work page doi:10.5281/zenodo.20090887 2026
[17]

Denoyelle, V

Q. Denoyelle, V. Duval, G. Peyré, and E. Soubies, The sliding Frank–Wolfe algorithm and its application to super-resolution microscopy,Inverse Problems(2019),doi:10.1088/1361-6420/ab2a29

work page doi:10.1088/1361-6420/ab2a29 2019
[18]

Duval and G

V. Duval and G. Peyré, Sparse regularization on thin grids I: the Lasso,Inverse Problems33 (2017), 055008,doi:10.1088/1361-6420/aa5e12

work page doi:10.1088/1361-6420/aa5e12 2017
[19]

Engram, K

M. Engram, K. M. Walter Anthony, T. Sachs, K. Kohnert, A. Serafimovich, G. Grosse, and F. J. Meyer, Remote sensing northern lake methane ebullition,Nature Climate Change10 (2020), 511–517,doi:10.1038/s41558-020-0762-8

work page doi:10.1038/s41558-020-0762-8 2020
[20]

Flinth, F

A. Flinth, F. de Gournay, and P. Weiss, On the linear convergence rates of exchange and continu- ous methods for total variation minimization,Mathematical Programming190 (2020), 221–257, doi:10.1007/s10107-020-01530-0

work page doi:10.1007/s10107-020-01530-0 2020
[21]

Hirst, D

B. Hirst, D. Randell, M. Jones, J. Chu, A. Kannath, N. Macleod, M. Dean, and D. Weidmann, Methane Emissions: Remote Mapping and Source Quantification Using an Open-Path Laser Dispersion Spectrometer,Geophysical Research Letters47 (2020),doi:10.1029/2019gl086725

work page doi:10.1029/2019gl086725 2020
[22]

J. P. Kaipio, V. Kolehmainen, M. Lassas, and E. Somersalo, Discretization as regularization in statistical inverse problems,SIAM Journal on Numerical Analysis44 (2006), 177–195,doi:10.1137/ 040616808

work page 2006
[23]

A. C. Kak and M. Slaney,Principles of Computerized Tomographic Imaging, IEEE Press, 1988

work page 1988
[24]

Lunardi,Analytic Semigroups and Optimal Regularity in Parabolic Problems, Progress in Non- linear Differential Equations and Their Applications, Birkhäuser Basel, 2011

A. Lunardi,Analytic Semigroups and Optimal Regularity in Parabolic Problems, Progress in Non- linear Differential Equations and Their Applications, Birkhäuser Basel, 2011. T. Valkonen Leak localisation with a measure source Manuscript, 2026-05-12 page 27 of 27

work page 2011
[25]

McNorton, C

J. McNorton, C. Wilson, M. Gloor, R. J. Parker, H. Boesch, W. Feng, R. Hossaini, and M. P. Chipperfield, Attribution of recent increases in atmospheric methane through 3-D inverse modelling,Atmospheric Chemistry and Physics18 (2018), 18149–18168,doi:10.5194/acp-18-18149-2018, https://acp.copernicus.org/articles/18/18149/2018/

work page doi:10.5194/acp-18-18149-2018 2018
[26]

S. M. Miller, S. C. Wofsy, A. M. Michalak, E. A. Kort, A. E. Andrews, S. C. Biraud, E. J. Dlugokencky, J. Eluszkiewicz, M. L. Fischer, G. Janssens-Maenhout, B. R. Miller, J. B. Miller, S. A. Montzka, T. Nehrkorn, and C. Sweeney, Anthropogenic emissions of methane in the United States,Proceedings of the National Academy of Sciences110 (2013), 20018–20022, ...

work page doi:10.1073/pnas.1314392110 2013
[27]

Pazy,Semigroups of Linear Operators and Applications to Partial Differential Equations, volume 44 of Applied Mathematical Sciences, Springer-Verlag, New York, 1983

A. Pazy,Semigroups of Linear Operators and Applications to Partial Differential Equations, volume 44 of Applied Mathematical Sciences, Springer-Verlag, New York, 1983

work page 1983
[28]

Pieper and D

K. Pieper and D. Walter, Linear convergence of accelerated conditional gradient algorithms in spaces of measures,ESAIM: Control, Optimization and Calculus of Variations27 (2021), 38, doi:10.1051/cocv/2021042,arXiv:1904.09218

work page doi:10.1051/cocv/2021042 2021
[29]

C. W. Rella, T. R. Tsai, C. G. Botkin, et al., Measuring emissions from oil and natural gas well pads using the mobile flux plane technique,Environ. Sci. Technol.49 (2015), 4742–4748,doi: 10.1021/es5054082

work page doi:10.1021/es5054082 2015
[30]

Valkonen, Proximal methods for point source localisation,Journal of Nonsmooth Analysis and Optimization4 (2023), 10433,doi:10.46298/jnsao-2023-10433,arXiv:2212.02991

T. Valkonen, Proximal methods for point source localisation,Journal of Nonsmooth Analysis and Optimization4 (2023), 10433,doi:10.46298/jnsao-2023-10433,arXiv:2212.02991

work page doi:10.46298/jnsao-2023-10433 2023
[31]

Valkonen, Point source localisation with unbalanced optimal transport, 2025,arXiv:2502.12417

T. Valkonen, Point source localisation with unbalanced optimal transport, 2025,arXiv:2502.12417. submitted

work page arXiv 2025
[32]

A. Voss, E. Vänskä, D. Weidmann, A. Pulkkinen, and A. Seppänen, Multi-open-path laser dis- persion spectroscopy combined with Bayesian state estimation for localizing and quantifying methane emissions,Atmospheric Environment: X22 (2024), 100260,doi:https://doi.org/10.1016/j. aeaoa.2024.100260,https://www.sciencedirect.com/science/article/pii/S2590162124000273

work page doi:10.1016/j 2024
[33]

Weidmann, B

D. Weidmann, B. Hirst, M. Jones, R. Ijzermans, D. Randell, N. Macleod, A. Kannath, J. Chu, and M. Dean, Locating and Quantifying Methane Emissions by Inverse Analysis of Path-Integrated Concentration Data Using a Markov-Chain Monte Carlo Approach,ACS Earth and Space Chemistry 6 (2022), 2190–2198,doi:10.1021/acsearthspacechem.2c00093

work page doi:10.1021/acsearthspacechem.2c00093 2022
[34]

Zeidler,Nonlinear Functional Analysis and Its Applications, number v

E. Zeidler,Nonlinear Functional Analysis and Its Applications, number v. 1 in Nonlinear functional analysis and its applications, Springer, 1986. T. Valkonen Leak localisation with a measure source

work page 1986