Sparse optimization for inverse problems in atmospheric modelling

We consider inverse problems in atmospheric modelling represented by a linear system which is based on a source-receptor sensitivity matrix and measurements. Instead of using the ordinary least squares, we add a weighting matrix based on the topology of measurement points and show the connection with Bayesian modelling. Since the source-receptor sensitivity matrix is usually ill-conditioned, the problem is often regularized, either by perturbing the objective function or by modifying the sensitivity matrix. However, both these approaches may be heavily dependent on specified parameters. To ease this burden, we propose to use techniques looking for a sparse solution with a small number of positive elements. Finally, we compare all these methods on the European Tracer Experiment (ETEX) data where there is no apriori information apart from the release position and some measurements. Inverse problems in atmospheric modelling are investigated.Sparse optimization techniques are introduced and discussed.Natural restrictions on the solution (nonnegativity of the release) are considered.The methods are compared on the European Tracer Experiment (ETEX) data.Free Matlab codes are provided.

[1]  Gerhard Wotawa,et al.  Xenon-133 and caesium-137 releases into the atmosphere from the Fukushima Dai-ichi nuclear power plant: determination of the source term, atmospheric dispersion, and deposition , 2011 .

[2]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[3]  S. Trini Castelli,et al.  Atmospheric tracer experiment uncertainties related to model evaluation , 2014, Environ. Model. Softw..

[4]  Heather Savoy,et al.  Software framework for inverse modeling and uncertainty characterization , 2015, Environ. Model. Softw..

[5]  Ian G. Enting,et al.  Inverse problems in atmospheric constituent transport , 2002 .

[6]  Scot M. Miller,et al.  Atmospheric inverse modeling with known physical bounds: an example from trace gas emissions , 2013 .

[7]  R. Tibshirani.  Regression Shrinkage and Selection via the Lasso , 1996 .

[8]  V. Šmídl,et al.  Efficient Sequential Monte Carlo Sampling for Continuous Monitoring of a Radiation Situation , 2014, Technometrics.

[9]  Xiaojin Zheng,et al.  Recent Advances in Mathematical Programming with Semi-continuous Variables and Cardinality Constraint , 2013 .

[10]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[11]  Kaisa Miettinen,et al.  Nonlinear multiobjective optimization , 1998, International series in operations research and management science.

[12]  A. Stohl,et al.  Validation of the lagrangian particle dispersion model FLEXPART against large-scale tracer experiment data , 1998 .

[13]  Joseph H. A. Guillaume,et al.  Characterising performance of environmental models , 2013, Environ. Model. Softw..

[14]  Michael Elad,et al.  From Sparse Solutions of Systems of Equations to Sparse Modeling of Signals and Images , 2009, SIAM Rev..

[15]  Martin Vetterli,et al.  A robust method for inverse transport modeling of atmospheric emissions using blind outlier detection , 2014 .

[16]  Yaakov Tsaig,et al.  Fast Solution of $\ell _{1}$ -Norm Minimization Problems When the Solution May Be Sparse , 2008, IEEE Transactions on Information Theory.

[17]  慧 廣瀬.  A Mathematical Introduction to Compressive Sensing , 2015 .

[18]  Stephen P. Boyd,et al.  Enhancing Sparsity by Reweighted ℓ1 Minimization , 2007, 0711.1612.

[19]  John C. Gille,et al.  Inverse modeling of carbon monoxide surface emissions using Climate Monitoring and Diagnostics Laboratory network observations , 2002 .

[20]  Ivan Dokmanic,et al.  The Fukushima inverse problem , 2013, 2013 IEEE International Conference on Acoustics, Speech and Signal Processing.

[21]  Enrico Pisoni,et al.  A new approach to design source-receptor relationships for air quality modelling , 2015, Environ. Model. Softw..

[22]  Tobias Achterberg,et al.  SCIP: solving constraint integer programs , 2009, Math. Program. Comput..

[23]  Michael J. Friedel,et al.  Hybrid modeling of spatial continuity for application to numerical inverse problems , 2013, Environ. Model. Softw..

[24]  Jonathan Currie,et al.  Opti: Lowering the Barrier Between Open Source Optimizers and the Industrial MATLAB User , 2012 .

[25]  S. Cohn,et al.  Ooce Note Series on Global Modeling and Data Assimilation Construction of Correlation Functions in Two and Three Dimensions and Convolution Covariance Functions , 2022 .

[26]  Nancy Nichols,et al.  Correlated observation errors in data assimilation , 2008 .

[27]  R. Tibshirani,et al.  Regression shrinkage and selection via the lasso: a retrospective , 2011 .

[28]  Marc Bocquet,et al.  Inverse modelling of atmospheric tracers: non-Gaussian methods and second-order sensitivity analysis , 2008 .

[29]  Julio Lumbreras,et al.  Advancements in the design and validation of an air pollution integrated assessment model for Spain , 2014, Environ. Model. Softw..

[30]  Boris Kompare,et al.  Environmental Modelling & Software , 2014 .

[31]  Dimitris Bertsimas,et al.  Algorithm for cardinality-constrained quadratic optimization , 2009, Comput. Optim. Appl..

[32]  S. Foucart,et al.  Sparsest solutions of underdetermined linear systems via ℓq-minimization for 0 , 2009 .

[33]  Guillermo Sapiro,et al.  Online Learning for Matrix Factorization and Sparse Coding , 2009, J. Mach. Learn. Res..

[34]  Holger Rauhut,et al.  A Mathematical Introduction to Compressive Sensing , 2013, Applied and Numerical Harmonic Analysis.

[35]  J. Whitaker,et al.  Distance-dependent filtering of background error covariance estimates in an ensemble Kalman filter , 2001 .

[36]  Petra Seibert,et al.  Source-receptor matrix calculation with a Lagrangian particle dispersion model in backward mode , 2004 .

[37]  Arthur Albert,et al.  Regression and the Moore-Penrose Pseudoinverse , 2012 .

[38]  Yuantao Gu,et al.  The Convergence Guarantees of a Non-Convex Approach for Sparse Recovery , 2012, IEEE Transactions on Signal Processing.

[39]  Martin Vetterli,et al.  Outlier removal for improved source estimation in atmospheric inverse problems , 2014, 2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[40]  Kerstin Stebel,et al.  Determination of time- and height-resolved volcanic ash emissions and their use for quantitative ash dispersion modeling: the 2010 Eyjafjallajökull eruption , 2011 .

[41]  Stephen P. Boyd,et al.  Graph Implementations for Nonsmooth Convex Programs , 2008, Recent Advances in Learning and Control.