Skip to content

Research at St Andrews

Non-stationary Gaussian models with physical barriers

Research output: Contribution to journalArticle

Open Access Status

  • Embargoed (until 18/01/20)


Haakon Bakka, Jarno Vanhatalo, Janine B. Illian, Daniel Simpson, Haavard Rue

School/Research organisations


The classical tools in spatial statistics are stationary models, like the Matérn field. However, in some applications there are boundaries, holes, or physical barriers in the study area, e.g. a coastline, and stationary models will inappropriately smooth over these features, requiring the use of a non-stationary model.

We propose a new model, the Barrier model, which is different from the established methods as it is not based on the shortest distance around the physical barrier, nor on boundary conditions. The Barrier model is based on viewing the Matérn correlation, not as a correlation function on the shortest distance between two points, but as a collection of paths through a Simultaneous Autoregressive (SAR) model. We then manipulate these local dependencies to cut off paths that are crossing the physical barriers. To make the new SAR well behaved, we formulate it as a stochastic partial differential equation (SPDE) that can be discretised to represent the Gaussian field, with a sparse precision matrix that is automatically positive definite.

The main advantage with the Barrier model is that the computational cost is the same as for the stationary model. The model is easy to use, and can deal with both sparse data and very complex barriers, as shown in an application in the Finnish Archipelago Sea. Additionally, the Barrier model is better at reconstructing the modified Horseshoe test function than the standard models used in R-INLA.



Original languageEnglish
JournalSpatial Statistics
VolumeIn press
Early online date18 Jan 2019
Publication statusE-pub ahead of print - 18 Jan 2019

    Research areas

  • Archipelago, Barriers, Coastline problem, INLA, Spatial statistics, Stochastic partial differential equations

Discover related content
Find related publications, people, projects and more using interactive charts.

View graph of relations

Related by journal

  1. Level set Cox processes

    Hildeman, A., Bolin, D., Wallin, J. & Illian, J. B., Dec 2018, In : Spatial Statistics. 28, p. 169-193

    Research output: Contribution to journalArticle

  2. A changepoint analysis of spatio-temporal point processes

    Altieri, L., Scott, E. M., Cocchi, D. & Illian, J. B., 2015, In : Spatial Statistics.

    Research output: Contribution to journalArticle

ID: 257463684