Methods
We collated records of P. muralis presence through field surveys and a citizen science campaign. We used presence-only models to predict climate suitability at a national scale (5km resolution), and fine-scale habitat suitability at the local scale (2m resolution). We then integrated local models into an individual-based modelling platform to simulate population dynamics and forecast range expansion for 10 populations in heterogeneous landscapes.