An Efficient Approach to Spatiotemporal Analysis and Modeling of Air Pollution Data

A statistically efficient approach is adopted for modeling spatial time-series of large data sets. The estimation of the main diagnostic tool such as the likelihood function in Gaussian spatiotemporal models is a cumbersome task when using extended spatial time-series such as air pollution. Here, using the Innovation Algorithm, we manage to compute it for many spatiotemporal specifications. These specifications refer to the spatial periodic-trend, the spatial autoregressive moving average, the spatial autoregressive integrated and fractionally integrated moving average Gaussian models. Our method is applied to daily pollutants over a large metropolitan area like Athens. In the applied part of our paper, we first diagnose temporal and spatial structures of data using non-likelihood based criteria, such as the empirical autocorrelation and covariance functions. Second, we use likelihood and non-likelihood based criteria to select a spatiotemporal model among various specifications. Finally, using kriging we regionalize the resulting parameter estimates of the best-fitted model in space at any unmonitored location in the Athens region. The results show that a specific autoregressive integrated moving average spatiotemporal model can optimally perform in within and out of spatial sample estimation. Supplemental materials for this article are available from the journal website.

[1]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[2]  K. Mardia,et al.  Maximum likelihood estimation of models for residual covariance in spatial regression , 1984 .

[3]  Noel A Cressie,et al.  Statistics for Spatial Data. , 1992 .

[4]  Phaedon C. Kyriakidis,et al.  Stochastic modeling of atmospheric pollution: a spatial time-series framework. Part I: methodology , 2001 .

[5]  Peter J. Brockwell,et al.  Gaussian maximum likelihood estimation for ARMA models II: Spatial processes , 2006 .

[6]  Michael L. Stein,et al.  Ozone Exposure and Population Density in Harris County, Texas: Comment , 1997 .

[7]  Rong Chen,et al.  Ozone Exposure and Population Density in Harris County, Texas , 1997 .

[8]  P. Whittle ON STATIONARY PROCESSES IN THE PLANE , 1954 .

[9]  P. Guttorp,et al.  A review of statistical methods for the meteorological adjustment of tropospheric ozone , 2001 .

[10]  J. Seinfeld RETHINKING THE OZONE PROBLEM IN URBAN AND REGIONAL AIR POLLUTION , 1991 .

[11]  K. Ord Estimation Methods for Models of Spatial Interaction , 1975 .

[12]  G. Shaddick,et al.  Modelling daily multivariate pollutant data at multiple sites , 2002 .

[13]  Comment on "Ozone exposure and population density in Harris County, Texas" , 1997 .

[14]  J. Besag Statistical Analysis of Non-Lattice Data , 1975 .

[15]  Christos Zerefos,et al.  Forecasting peak pollutant levels from meteorological variables , 1995 .

[16]  J. R. Wallis,et al.  An Approach to Statistical Spatial-Temporal Modeling of Meteorological Fields , 1994 .

[17]  Gwilym M. Jenkins,et al.  Time series analysis, forecasting and control , 1972 .

[18]  Richard A. Davis,et al.  Time Series: Theory and Methods , 2013 .

[19]  Lippincott Williams Wilkins,et al.  Part I , 1997, Neurology.

[20]  J. Besag Spatial Interaction and the Statistical Analysis of Lattice Systems , 1974 .

[21]  Peter C. Young,et al.  A recursive estimation approach to the spatio-temporal analysis and modelling of air quality data , 2006, Environ. Model. Softw..