Quasi-3D Statistical Inversion of Oceanographic Tracer Data

Document Type

Presentation Abstract

Presentation Date

3-2-2006

Abstract

A forward problem is the task of finding the solution u of a differential equation, given a set of inputs ɸ (coefficients, boundary conditions). The problem of determination of ɸ from u may be regarded as “inverse” to the one described above. Inverse problems arise in numerous fields such as general acoustics, earth sciences, algorithm development, medical imaging, etc. Inverse problems are statistical problems. The purpose is to estimate ɸ having (in general) noisy, sparse and sometimes only partial measurements of the solution u. A robust solution to an inverse problem can be obtained by introducing prior information on ɸ and modeling the measurement error.

The application we are currently working on involves estimating water velocities and mixing coefficients in a 2 km deep, rectangular region in the South Atlantic Ocean. Partial and sparse measurements of tracer concentrations (salinity, oxygen, etc.) are available. The data are filtered to eliminate outliers, then interpolated to the nearest points of a regular lattice and restricted to thin neutral density layers. The connection between velocities, diffusion coefficients, boundary conditions and tracer concentrations is made via a 3D advection-diffusion equation and a geostrophic flow model. The (un-normalized) posterior density of the parameters conditionally on the data is summarized using Markov chain Monte Carlo techniques. We reconstruct the tracer fields as well, thus, for regions where no data was available, concentrations are now estimated in a manner that is consistent with physical principles.

Additional Details

Thursday, 2 March 2006
4:10 p.m. in Math 109

This document is currently not available here.

Share

COinS