Comparing distribution of harbour porpoise using generalized additive models and hierarchical Bayesian models with integrated nested laplace approximation

• Published in: Ecological Modelling
• Main Authors: Williamson, LD, Scott, BE, Laxton, M, Illian, JB, Todd, VLG, Miller, PI, Brookes, KL
• Format: Article in Journal/Newspaper
• Language: English
• Published: Elsevier 2022
• Subjects: Gam; Laplace; Harbour porpoise; Phocoena phocoena
• Online Access: http://plymsea.ac.uk/id/eprint/9727/; http://plymsea.ac.uk/id/eprint/9727/1/Williamson_et_al_Ecol_Model_2022_porpoise_model.pdf; https://doi.org/10.1016/j.ecolmodel.2022.110011

Species Distribution Models (SDMs) are used regularly to develop management strategies, but many modelling methods ignore the spatial nature of data. To address this, we compared fine-scale spatial distribution predictions of harbour porpoise (Phocoena phocoena) using empirical aerial-video-survey data collected along the east coast of Scotland in August and September 2010 and 2014. Incorporating environmental covariates that cover habitat preferences and prey proxies, we used a traditional (and commonly implemented) Generalized Additive Model (GAM), and two Hierarchical Bayesian Modelling (HBM) approaches using Integrated Nested Laplace Approxi�mation (INLA) model-fitting methodology. One HBM-INLA modelled gridded space (similar to the GAM), and the other dealt more explicitly in continuous space using a Log-Gaussian Cox Process (LGCP). Overall, predicted distributions in the three models were similar; however, HBMs had twice the level of certainty, showed much finer-scale patterns in porpoise distribution, and identified some areas of high relative density that were not apparent in the GAM. Spatial differences were due to how the two methods accounted for autocorrelation, spatial clustering of animals, and differences between modelling in discrete vs. continuous space; consequently, methods for spatial analyses likely depend on scale at which results, and certainty, are needed. For large-scale analysis (>5–10 km resolution, e.g. initial impact assessment), there was little difference be�tween results; however, insights into fine-scale (<1 km) distribution of porpoise from the HBM model using LGCP, while more computationally costly, offered potential benefits for refining conservation management or mitigation measures within offshore developments or protected areas.