Identification of a Point of Release by Use of Optimally Weighted Least Squares

This paper addresses the parametric inverse problem of locating the point of release of atmospheric pollution. A finite set of observed mixing ratios is compared, by use of least squares, with the analogous mixing ratios computed by an adjoint dispersion model for all possible locations of the release. Classically, the least squares are weighted using the covariance matrix of the measurement errors. However, in practice, this matrix cannot be determined for the prevailing part of these errors arising from the limited representativity of the dispersion model. An alternative weighting proposed here is related to a unified approach of the parametric and assimilative inverse problems corresponding, respectively, to identification of the point of emission or estimation of the distributed emissions. The proposed weighting is shown to optimize the resolution and numerical stability of the inversion. The importance of the most common monitoring networks, with point detectors at various locations, is stressed as a misleading singular case. During the procedure it is also shown that a monitoring network, under given meteorological conditions, itself contains natural statistics about the emissions, irrespective of prior assumptions.

[1]  Dusanka Zupanski,et al.  Model Error Estimation Employing an Ensemble Data Assimilation Approach , 2006 .

[2]  John H. Seinfeld,et al.  Inverse air pollution modelling of urban-scale carbon monoxide emissions , 1995 .

[3]  Marc Bocquet,et al.  Diagnosis and impacts of non-Gaussianity of innovations in data assimilation , 2010 .

[4]  Hubert B. Keller,et al.  Source localization by spatially distributed electronic noses for advection and diffusion , 2005, IEEE Transactions on Signal Processing.

[5]  J. Issartel,et al.  Emergence of a tracer source from air concentration measurements, a new strategy for linear assimilation , 2005 .

[6]  R. Daley The Effect of Serially Correlated Observation and Model Error on Atmospheric Data Assimilation , 1992 .

[7]  Marc Bocquet,et al.  Data assimilation for short-range dispersion of radionuclides: An application to wind tunnel data , 2006 .

[8]  Christopher M. Danforth,et al.  Accounting for Model Errors in Ensemble Data Assimilation , 2009 .

[9]  Janusz A. Pudykiewicz,et al.  APPLICATION OF ADJOINT TRACER TRANSPORT EQUATIONS FOR EVALUATING SOURCE PARAMETERS , 1998 .

[10]  James A. Hansen Accounting for Model Error in Ensemble-Based State Estimation and Forecasting , 2002 .

[11]  Maithili Sharan,et al.  An inversion technique for the retrieval of single-point emissions from atmospheric concentration measurements , 2009, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[12]  G. Milliken,et al.  On Necessary and Sufficient Conditions for Ordinary Least Squares Estimators to Be Best Linear Unbiased Estimators , 1984 .

[13]  K. Mardia Characterizations of Directional Distributions , 1975 .

[14]  Alberto Carrassi,et al.  Accounting for Model Error in Variational Data Assimilation: A Deterministic Formulation , 2010 .

[15]  K. S. Rao Source estimation methods for atmospheric dispersion , 2007 .

[16]  R. Fisher Dispersion on a sphere , 1953, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[17]  M. Singh,et al.  A mathematical model for the dispersion of air pollutants in low wind conditions , 1996 .

[18]  Maithili Sharan,et al.  An inversion technique to retrieve the source of a tracer with an application to synthetic satellite measurements , 2007, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.