Journal Home > Volume 1 , Issue 4

Numerical simulations are widely used as a predictive tool to better understand complex air flows and pollution transport on the scale of individual buildings, city blocks, and entire cities. To improve prediction for air flows and pollution transport, we propose a Variational Data Assimilation (VarDA) model which assimilates data from sensors into the open-source, finite-element, fluid dynamics model Fluidity. VarDA is based on the minimization of a function which estimates the discrepancy between numerical results and observations assuming that the two sources of information, forecast and observations, have errors that are adequately described by error covariance matrices. The conditioning of the numerical problem is dominated by the condition number of the background error covariance matrix which is ill-conditioned. In this paper, a preconditioned VarDA model is presented, it is based on a reduced background error covariance matrix. The Empirical Orthogonal Functions (EOFs) method is used to alleviate the computational cost and reduce the space dimension. Experimental results are provided assuming observed values provided by sensors from positions mainly located on roofs of buildings.


menu
Abstract
Full text
Outline
About this article

Effective Variational Data Assimilation in Air-Pollution Prediction

Show Author's information Rossella Arcucci( )Christopher PainYi-Ke Guo
Data Science Institute, Department of Computing, Imperial College London, London, SW7 2AZ, United Kingdom.
Department of Earth Science & Engineering, Imperial College London, London, SW7 2AZ, United Kingdom.

Abstract

Numerical simulations are widely used as a predictive tool to better understand complex air flows and pollution transport on the scale of individual buildings, city blocks, and entire cities. To improve prediction for air flows and pollution transport, we propose a Variational Data Assimilation (VarDA) model which assimilates data from sensors into the open-source, finite-element, fluid dynamics model Fluidity. VarDA is based on the minimization of a function which estimates the discrepancy between numerical results and observations assuming that the two sources of information, forecast and observations, have errors that are adequately described by error covariance matrices. The conditioning of the numerical problem is dominated by the condition number of the background error covariance matrix which is ill-conditioned. In this paper, a preconditioned VarDA model is presented, it is based on a reduced background error covariance matrix. The Empirical Orthogonal Functions (EOFs) method is used to alleviate the computational cost and reduce the space dimension. Experimental results are provided assuming observed values provided by sensors from positions mainly located on roofs of buildings.

Keywords: big data, preconditioning, data assimilation, reduced order space

References(22)

[1]
E. Kalnay, Atmospheric Modeling, Data Assimilation and Predictability. Cambridge University Press, 2003.
DOI
[2]
D. M. Baker, W. Huang, Y. R. Guo, J. Bourgeois, and Q. N. Xiao, Three-dimensional variational data assimilation system for MM5: Implementation and initial results, Mon. Wea. Rev., vol. 132, pp. 897-914, 2004.
DOI
[3]
E. Andersson, J. Haseler, P. Undén, P. Courtier, G. Kelly, D. Vasiljevic, C. Brancovic, C. Cardinali, C. Gaffard, A. Hollingsworth, et al., The ECMWF implementation of three dimensional variational assimilation (3DVar), Part III: Experimental results, Quarterly Journal Royal Met. Society, vol. 124, pp. 1831-1860, 1998.
[4]
R. Arcucci, L. D’Amore, L. Carracciuolo, G. Scotti, and G. Laccetti, A decomposition of the tikhonov regularization functional oriented to exploit hybrid multilevel parallelism, International Journal of Parallel Programming, vol. 45, no. 5, pp. 1214-1235, 2017.
[5]
J. N. Nichols, Mathematical concepts in data assimilation, in Data Assimilation, W. Lahoz, B. Khattatov, and R. Menard, eds. Springer, 2010.
[6]
J. S. A. Haben, Conditioning and preconditioning of the minimisation problem in variational data assimilation, PhD dissertation, University of Reading, UK, 2011.
DOI
[7]
A. M. Moore, H. G. Arango, G. Broquet, B. S. Powell, A. T. Weaver, and J. Zavala-Garay, The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part I - System overview and formulation, J. Progress in Oceanography, vol. 91, pp. 34-49, 2011.
[8]
C. C. Pain, A. P. Umpleby, C. R. E. De Oliveira, and A. J. H. Goddard, Tetrahedral mesh optimisation and adaptivity for steady-state and transient finite element calculations, Computer Methods in Applied Mechanics and Engineering, vol. 190, pp. 3771-3796, 2001.
[9]
D. R. Davies, C. R. Wilson, and S. C. Kramer, Fluidity: A fully unstructured anisotropic adaptive mesh computational modeling framework for geodynamics, Geochemistry, Geophysics, Geosystems, vol. 12, no. 6, 2011.
[10]
D. C. Liu and J. Nocedal, On the limited memory BFGS method for large scale optimization, Mathematical Programming, vol. 45, pp. 503-528, 1989.
[11]
J. Nocedal, R. H. Byrd, P. Lu, and C. Zhu, L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization, ACM Transactions on Mathematical Software, vol. 23, pp. 550-560, 1997
[12]
R. Arcucci, L. D’Amore, J. Pistoia, R. Toumi, and A. Murli, On the variational data assimilation problem solving and sensitivity analysis, Journal of Computational Physics, vol. 335, pp. 311-326, 2017.
[13]
E. N. Lorenz, Empirical orthogonal functions and statistical weather prediction, Statistical Forecasting Project, MIT, Sci. Rep. No. 1, 1956.
[14]
D. G. Cacuci, I. M. Navon, and M. Ionescu-Bujor, Computational Methods for Data Evaluation and Assimilation. CRC Press, 2013.
[15]
C. Hansen, Rank-Deficient and Discrete Ill-Posed Problems, Numerical Aspects of Linear Inversion. SIAM, 1998.
DOI
[16]
P. C. Hansen, J. G. Nagy, and D. P. O’Leary, Deblurring Images: Matrices, Spectra, and Filtering. SIAM, 2006.
DOI
[17]
J. P. Courtier, A strategy for operational implementation of 4D-VAR, using an incremental approach, Q. J. R. Meteorol Soc., vol. 120, 1367-1387, 1994.
[18]
H. K. Engl, M. Hanke, and A. Neubauer, Regularization of Inverse Problems. Kluwer, 1996.
DOI
[19]
A. C. Lorenc, Development of an operational variational assimilation scheme, Journal of the Meteorological Society of Japan, vol. 75, pp. 339-346, 1997.
[20]
E. Aristodemou, T. Bentham, C. Pain, and A. Robins, A comparison of mesh-adaptive LES with wind tunnel data for flow past buildings: Mean flows and velocity fluctuations, Atmospheric Environment Journal, vol. 43, pp. 6238-6253, 2009.
[21]
R. Ford, C. C. Pain, A. J. H. Goddard, C. R. E. De Oliveira, and A. P. Umpleby, A nonhydrostatic finite-element model for three-dimensional stratified oceanic flows. Part I: Model formulation, Monthly Weather Review, vol. 132, pp. 2816-2831, 2004.
[22]
Imperial College London AMCG, Fluidity Manual v4.1.12, https://figshare.com/articles/Fluidity_Manual/1387713, 2015.
Publication history
Copyright
Acknowledgements
Rights and permissions

Publication history

Received: 16 March 2018
Accepted: 20 March 2018
Published: 02 July 2018
Issue date: December 2018

Copyright

© The author(s) 2018

Acknowledgements

This work was supported by the EPSRC Grand Challenge grant "Managing Air for Green Inner Cities" (MAGIC) EP/N010221/1.

Rights and permissions

Return