References

ACP

Atmospheric Chemistry and Physics

ACP

1680-7324

Copernicus GmbH

Göttingen, Germany

10.5194/acp-8-7697-2008

Aerosol model selection and uncertainty modelling by adaptive MCMC technique

Laine

¹ Tamminen

Finnish Meteorological Institute, Helsinki, Finland

19 12 2008

8 24 7697 7707

This is an open-access article ditributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

This article is available from http://www.atmos-chem-phys.net/8/7697/2008/acp-8-7697-2008.html

The full text article is available as a PDF file from http://www.atmos-chem-phys.net/8/7697/2008/acp-8-7697-2008.pdf

We present a new technique for model selection problem in atmospheric remote sensing. The technique is based on Monte Carlo sampling and it allows model selection, calculation of model posterior probabilities and model averaging in Bayesian way. <br><br> The algorithm developed here is called Adaptive Automatic Reversible Jump Markov chain Monte Carlo method (AARJ). It uses Markov chain Monte Carlo (MCMC) technique and its extension called Reversible Jump MCMC. Both of these techniques have been used extensively in statistical parameter estimation problems in wide area of applications since late 1990's. The novel feature in our algorithm is the fact that it is fully automatic and easy to use. <br><br> We show how the AARJ algorithm can be implemented and used for model selection and averaging, and to directly incorporate the model uncertainty. We demonstrate the technique by applying it to the statistical inversion problem of gas profile retrieval of GOMOS instrument on board the ENVISAT satellite. Four simple models are used simultaneously to describe the dependence of the aerosol cross-sections on wavelength. During the AARJ estimation all the models are used and we obtain a probability distribution characterizing how probable each model is. By using model averaging, the uncertainty related to selecting the aerosol model can be taken into account in assessing the uncertainty of the estimates.

References 1

Bertaux, J L., Kyrölä, E., and Wehr, T.: Stellar Occultation Technique for Atmospheric Ozone Monitoring: GOMOS on Envisat, Earth Observation Quarterly, 67, 17–20, 2000.

ESA 2007: GOMOS Product Handbook Issue 3.0, European Space Agency, http://envisat.esa.int/dataproducts/, 2007.

Gamerman, D.: Markov Chain Monte Carlo – Stochastic simulation for Bayesian inference, Chapman & Hall, 245 pp., 1997.

Gelman, A., Roberts, G O., and Gilks, W R.: Efficient Metropolis Jumping Rules, Bayesian Statistics, 5, 599–607, 1996.

Green, P J.: Reversible jump Markov chain Monte Carlo computation and Bayesian model determination, Biometrika, 82, 711–732, 1995.

Green, P J.: Trans-dimensional Markov chain Monte Carlo, in: Highly Structured Stochastic Systems, edited by: Green, P J., Hjort, N L., and Richardson, S., no 27, Oxford Statistical Science Series, pp. 179–198, Oxford University Press, 2003.

Haario, H., Saksman, E., and Tamminen, J.: An adaptive Metropolis algorithm, Bernoulli, 7, 223–242, 2001.

Haario, H., Laine, M., Lehtinen, M., Saksman, E., and Tamminen, J.: MCMC methods for high dimensional inversion in remote sensing, J. Roy. Stat. Soc. B., 66, 591–607, \doi10.1111/j.1467-9868.2004.02053.x, 2004.

Haario, H., Laine, M., Mira, A., and Saksman, E.: DRAM: Efficient adaptive MCMC, Stat. Comput., 16, 339–354, \doi10.1007/s11222-006-9438-0, 2006.

Hastie, D.: Towards Automatic Reversible Jump Markov Chain Monte Carlo, Ph.D. thesis, University of Bristol Department of Mathematics, 2005.

Kaipio, J P. and Somersalo, E.: Computational and Statistical Methods for Inverse Problems, Springer, 353 pp., 2004.

Laine, M.: MCMC toolbox for Matlab website, http://www.helsinki.fi/~mjlaine/mcmc/, 2008.

Mira, A.: On Metropolis-Hastings algorithms with delayed rejection, Metron, LIX, 231–241, 2001.

Spiegelhalter, D J., Best, N G., Carlin, B P., and van~der Linde, A.: Bayesian measures of model complexity and fit, J. Roy. Stat. Soc. B., 64, 583–639, 2002.

Tamminen, J.: Adaptive Markov chain Monte Carlo algorithms with geophysical applications, Finnish Meteorological Institute Contributions~47, Finnish Meteorological Institute, Helsinki, http://urn.fi/URN:ISBN:952-10-2016-4, 2004.

Tamminen, J. and Kyrölä, E.: Bayesian solution for nonlinear and non-Gaussian inverse problem by Markov chain Monte Carlo method, J. Geophys. Res., 106, 14 377–14 390, 2001.

Vanhellemont, F., Fussen, D., Dodion, J., Bingen, C., and Mateshvili, N.: Choosing a suitable analytical model for aerosol extinction spectra in the retrieval of UV/visible satellite occultation measurements, J. Geophys. Res., 111, 1–8, doi:10.1029/2005JD0069412006.

Veihelmann, B., Levelt, P. F., Stammes, P., and Veefkind, J. P.: Simulation study of the aerosol information content in OMI spectral reflectance measurements, Atmos. Chem. Phys., 7, 3115–3127, 2007.