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.

On Bayesian data assimilation for PDEs with ill-posed forward problems

A unified framework for the analysis of accuracy and stability of a class of approximate Gaussian filters for the Navier-Stokes Equations

A Dynamic Likelihood Approach to Filtering for Advection-Diffusion Dynamics

On posterior consistency of data assimilation with Gaussian process priors: the 2D Navier-Stokes equations

Bayesian inverse problems for functions and applications to fluid mechanics

Data Assimilation for Unsaturated Flow Models with Restart Adaptive Probabilistic Collocation Based Kalman Filter

Data assimilation for the stochastic Camassa-Holm equation using particle filtering: a numerical investigation

Enhancing nonlinear solvers for the Navier–Stokes equations with continuous (noisy) data assimilation

Data assimilation for the heat equation using stabilized finite element methods

Data assimilation in 2D incompressible Navier-Stokes equations, using a stabilized explicit $O(Δt)^2$ leapfrog finite difference scheme run backward in time

Bayesian inference for geophysical fluid dynamics using generative models

Adaptivity in Bayesian Inverse Finite Element Problems: Learning and Simultaneous Control of Discretisation and Sampling Errors

Enhancing nonlinear solvers for the Navier-Stokes equations with continuous (noisy) data assimilation

Preserving Nonlinear Constraints in Variational Flow Filtering Data Assimilation

Inverse Problems in a Bayesian Setting

Postprocessing Galerkin method applied to a data assimilation algorithm: a uniform in time error estimate

Data Assimilation for Combined Parameter and State Estimation in Stochastic Continuous-Discrete Nonlinear Systems

A Computational Study of a Data Assimilation Algorithm for the Two-dimensional Navier--Stokes Equations

On the stability and the uniform propagation of chaos properties of Ensemble Kalman-Bucy filters

Dynamical approximation and sensor placement for filtering problems

Generalized Bayes Approach to Inverse Problems with Model Misspecification