Abstract:In this paper we establish a mathematical framework for a range of inverse problems for functions, given a finite set of noisy observations. The problems are hence underdetermined and are often ill-posed. We study these problems from the viewpoint of Bayesian statistics, with the resulting posterior probability measure being defined on a space of functions. We develop an abstract framework for such problems which facilitates application of an infinite-dimensional version of Bayes theorem, leads to a well-posedness result for the posterior measure (continuity in a suitable probability metric with respect to changes in data), and also leads to a theory for the existence of maximizing the posterior probability (MAP) estimators for such Bayesian inverse problems on function space. A central idea underlying these results is that continuity properties and bounds on the forward model guide the choice of the prior measure for the inverse problem, leading to the desired results on well-posedness and MAP estimators; the PDE analysis and probability theory required are thus clearly dileneated, allowing a straightforward derivation of results. We show that the abstract theory applies to some concrete applications of interest by studying problems arising from data assimilation in fluid mechanics. The objective is to make inference about the underlying velocity field, on the basis of either Eulerian or Lagrangian observations. We study problems without model error, in which case the inference is on the initial condition, and problems with model error in which case the inference is on the initial condition and on the driving noise process or, equivalently, on the entire time-dependent velocity field. In order to undertake a relatively uncluttered mathematical analysis we consider the two-dimensional Navier–Stokes equation on a torus. The case of Eulerian observations—direct observations of the velocity field itself—is then a model for weather forecasting. The case of Lagrangian observations—observations of passive tracers advected by the flow—is then a model for data arising in oceanography. The methodology which we describe herein may be applied to many other inverse problems in which it is of interest to find, given observations, an infinite-dimensional object, such as the initial condition for a PDE. A similar approach might be adopted, for example, to determine an appropriate mathematical setting for the inverse problem of determining an unknown tensor arising in a constitutive law for a PDE, given observations of the solution. The paper is structured so that the abstract theory can be read independently of the particular problems in fluid mechanics which are subsequently studied by application of the theory.

Multilevel Markov Chain Monte Carlo for Bayesian inverse problems for Navier Stokes equation with Lagrangian Observations

Motion Estimation under Location Uncertainty for Turbulent Fluid Flows

Application of Bayesian Inference to Estimate the Parameters in 2D Convection-diffusion Equation with Source

Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems

Bayesian inverse problems for functions and applications to fluid mechanics

Multilevel Monte Carlo methods for the Dean-Kawasaki equation from Fluctuating Hydrodynamics

Bayesian inversion of log-normal eikonal equations

Multiscale model reduction method for Bayesian inverse problems of subsurface flow

Forward and Inverse Uncertainty Quantification using Multilevel Monte Carlo Algorithms for an Elliptic Nonlocal Equation

Multigrid Monte Carlo Revisited: Theory and Bayesian Inference

Bayesian Parameter Inference for Partially Observed Diffusions using Multilevel Stochastic Runge-Kutta Methods

A unified Bayesian inversion approach for a class of tumor growth models with different pressure laws

MCMC for a hyperbolic Bayesian inverse problem in traffic flow modelling

Multilevel higher order Quasi-Monte Carlo Bayesian Estimation

Bayesian inverse Navier-Stokes problems: joint flow field reconstruction and parameter learning

Sequential Markov chain Monte Carlo for Lagrangian data assimilation with applications to unknown data locations

2D Backward Stochastic Navier–Stokes Equations with Nonlinear Forcing

Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation

Efficient, multimodal, and derivative-free Bayesian inference with Fisher-Rao gradient flows

Markov Chain Monte Carlo from Lagrangian Dynamics