Departmental Papers (MEAM)

Document Type

Conference Paper

Date of this Version

May 2006


Postprint version. Published in Lecture Notes in Computer Science, Volume 3993, Computational Science – ICCS 2006 6th International Conference, Proceedings, Part III, 2006, pages 481-488.
Publisher URL:


We are interested in a DDDAS problem of localization of airborne contaminant releases in regional atmospheric transport models from sparse observations. Given measurements of the contaminant over an observation window at a small number of points in space, and a velocity field as predicted for example by a mesoscopic weather model, we seek an estimate of the state of the contaminant at the begining of the observation interval that minimizes the least squares misfit between measured and predicted contaminant field, subject to the convection-diffusion equation for the contaminant. Once the "initial" conditions are estimated by solution of the inverse problem, we issue predictions of the evolution of the contaminant, the observation window is advanced in time, and the process repeated to issue a new prediction, in the style of 4D-Var. We design an appropriate numerical strategy that exploits the spectral structure of the inverse operator, and leads to efficient and accurate resolution of the inverse problem. Numerical experiments verify that high resolution inversion can be carried out rapidly for a well-resolved terrain model of the greater Los Angeles area.



Date Posted: 30 November 2006

This document has been peer reviewed.