A Semi-analytic but Biased Uncertainty Assessment Method using Sample Extensions, Analysed for Nonlinear Travel Time Tomography
Abstract
Many geophysical problems can be cast as inverse problems that estimate a set of parameter values from observed data. Within a Bayesian framework, solutions to such problems are described probabilistically by the so-called posterior probability distribution functions (pdf's). To obtain robust inference results often requires millions of model parameter value samples to be drawn, and their simulation to be performed; this is a computationally expensive procedure. We investigate the concept of sample extensions as a means to improve efficiency when solving fully nonlinear inverse problems. A sample's extension is defined as the set of models or parameter values whose forward function values are directly accessible from a sample for which the forward function has already been evaluated, obviating the need for additional forward function evaluations. In a specific case of first-arrival travel time calculations used in seismic travel time tomography, we apply sample extensions to obtain a continuous region with non-zero hypervolume within parameter space, across all of which the forward function values are known given only a single forward simulation. We devise a deterministic sampling technique that identifies the most informative extensions by solving an optimisation problem. In an illustrative tomographic example that involves a single travel time datum, we find 51 optimal samples, and use them to construct an analytic approximation to the Bayesian posterior pdf. Additionally, we propose an extensions-based algorithm for real-world tomography scenarios and apply it to a synthetic 2D example. This study highlights two fundamental problems that make the method inefficient: (1) limited hypervolumes of extensions and (2) neglecting parameter correlations to simplify analytic calculations. Finding solutions to these problems defines possible directions for future research.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.