Spatio-temporal modeling of particulate matter concentration through the SPDE approach

Michela Cameletti*, Finn Lindgren, Daniel Simpson, Håvard Rue

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

236 Scopus citations


In this work, we consider a hierarchical spatio-temporal model for particulate matter (PM) concentration in the North-Italian region Piemonte. The model involves a Gaussian Field (GF), affected by a measurement error, and a state process characterized by a first order autoregressive dynamic model and spatially correlated innovations. This kind of model is well discussed and widely used in the air quality literature thanks to its flexibility in modelling the effect of relevant covariates (i. e. meteorological and geographical variables) as well as time and space dependence. However, Bayesian inference-through Markov chain Monte Carlo (MCMC) techniques-can be a challenge due to convergence problems and heavy computational loads. In particular, the computational issue refers to the infeasibility of linear algebra operations involving the big dense covariance matrices which occur when large spatio-temporal datasets are present. The main goal of this work is to present an effective estimating and spatial prediction strategy for the considered spatio-temporal model. This proposal consists in representing a GF with Matérn covariance function as a Gaussian Markov Random Field (GMRF) through the Stochastic Partial Differential Equations (SPDE) approach. The main advantage of moving from a GF to a GMRF stems from the good computational properties that the latter enjoys. In fact, GMRFs are defined by sparse matrices that allow for computationally effective numerical methods. Moreover, when dealing with Bayesian inference for GMRFs, it is possible to adopt the Integrated Nested Laplace Approximation (INLA) algorithm as an alternative to MCMC methods giving rise to additional computational advantages. The implementation of the SPDE approach through the R-library INLA ( is illustrated with reference to the Piemonte PM data. In particular, providing the step-by-step R-code, we show how it is easy to get prediction and probability of exceedance maps in a reasonable computing time.

Original languageEnglish (US)
Pages (from-to)109-131
Number of pages23
JournalAStA Advances in Statistical Analysis
Issue number2
StatePublished - Apr 2013


  • Covariance functions
  • Gaussian Markov random fields
  • Gaussian fields
  • Hierarchical models
  • Integrated Nested Laplace Approximation

ASJC Scopus subject areas

  • Analysis
  • Statistics and Probability
  • Modeling and Simulation
  • Social Sciences (miscellaneous)
  • Economics and Econometrics
  • Applied Mathematics


Dive into the research topics of 'Spatio-temporal modeling of particulate matter concentration through the SPDE approach'. Together they form a unique fingerprint.

Cite this