Software for `Reaction-diffusion models of invasive tree pest spread: quantifying the expansion of oak processionary moth in the UK'

<p dir="ltr">This software allows the user to recreate the figures from the preprint, "Reaction-diffusion models of invasion tree pest spread: quantifying the expansion of oak processionary moth in the UK" [https://arxiv.org/abs/2509.14166]. This includes simulating the models and conducting parameter inference.</p><p dir="ltr">UK woodlands, forests, and urban treescapes are under threat from invasive species, exacerbated by climate change, trade, and transport. Invasive tree pests debilitate their host and disrupt forest ecosystems, thus it is imperative to quantitatively model and predict their spread. Addressing this, we represent the spatial distribution of the pest as a population density field which evolves according to a spatiotemporal reaction-diffusion equation. We solve this intractable system of equations numerically and, from the solution, we determine first arrival times of the pest at locations in the field. The adopted model permits us to obtain the expansion rate of pest spread directly from the model parameters, which we infer in the Bayesian paradigm, using a Markov chain Monte Carlo scheme. We apply our framework to the ongoing spread of oak processionary moth in the UK, an outbreak which continues to grow despite management efforts. We demonstrate that our approach effectively captures the spread of the pest and that this has occurred at a non-constant expansion rate. The proposed framework is a powerful tool for quantitatively modelling the spread of an invasive tree pest and could underpin future prediction and management approaches.</p>