A Bayesian space-time model for discrete spread processes on a lattice

In this article we present a Bayesian Markov model for investigating environmental spread processes. We formulate a model where the spread of a disease over a heterogeneous landscape through time is represented as a probabilistic function of two processes: local diffusion and random-jump dispersal. This formulation represents two mechanisms of spread which result in highly peaked and long-tailed distributions of dispersal distances (i.e., local and long-distance spread), commonly observed in the spread of infectious diseases and biological invasions. We demonstrate the properties of this model using a simulation experiment and an empirical case study - the spread of mountain pine beetle in western Canada. Posterior predictive checking was used to validate the number of newly inhabited regions in each time period. The model performed well in the simulation study in which a goodness-of-fit statistic measuring the number of newly inhabited regions in each time interval fell within the 95% posterior predictive credible interval in over 97% of simulations. The case study of a mountain pine beetle infestation in western Canada (1999-2009) extended the base model in two ways. First, spatial covariates thought to impact the local diffusion parameters, elevation and forest cover, were included in the model. Second, a refined definition for translocation or jump-dispersal based on mountain pine beetle ecology was incorporated improving the fit of the model. Posterior predictive checks on the mountain pine beetle model found that the observed goodness-of-fit test statistic fell within the 95% posterior predictive credible interval for 8 out of 10 years. The simulation study and case study provide evidence that the model presented here is both robust and flexible; and is therefore appropriate for a wide range of spread processes in epidemiology and ecology.