A Bayesian spatio-temporal geostatistical model with an auxiliary lattice for large datasets

Ganggang Xu, Faming Liang, Marc G. Genton

Research output: Contribution to journalArticlepeer-review

13 Scopus citations


When spatio-temporal datasets are large, the computational burden can lead to failures in the implementation of traditional geostatistical tools. In this paper, we propose a computationally efficient Bayesian hierarchical spatio-temporal model in which the spatial dependence is approximated by a Gaussian Markov random field (GMRF) while the temporal correlation is described using a vector autoregressive model. By introducing an auxiliary lattice on the spatial region of interest, the proposed method is not only able to handle irregularly spaced observations in the spatial domain, but it is also able to bypass the missing data problem in a spatio-temporal process. Because the computational complexity of the proposed Markov chain Monte Carlo algorithm is of the order O(n) with n the total number of observations in space and time, our method can be used to handle very large spatio-temporal datasets with reasonable CPU times. The performance of the proposed model is illustrated using simulation studies and a dataset of precipitation data from the coterminous United States.
Original languageEnglish (US)
Pages (from-to)61-79
Number of pages19
Issue number1
StatePublished - Jan 9 2014

ASJC Scopus subject areas

  • Statistics and Probability
  • Statistics, Probability and Uncertainty


Dive into the research topics of 'A Bayesian spatio-temporal geostatistical model with an auxiliary lattice for large datasets'. Together they form a unique fingerprint.

Cite this