Comparing spatiotemporal species distribution models: A case study of a Scotian Shelf sea cucumber (Cucumaria frondosa)

AbstractNumerous spatiotemporal species distribution modeling frameworks are now available to the ecological practitioner. This study compared three such frameworks accessible in the R programming language: generalized additive models with spatiotemporal smooths as implemented by mgcv, spatiotemporal generalized linear mixed models based on nearest neighbor Gaussian processes as implemented by starve, and spatiotemporal generalized linear mixed models based on the stochastic partial differential equations approach as implemented by sdmTMB. The primary focus was to compare the inferences obtained from applying these frameworks to the case study of the orange‐footed sea cucumber, Cucumaria frondosa, on the Scotian Shelf off Nova Scotia, Canada. Each model was fit to catch data (2000–2019) from Fisheries and Oceans Canada's annual Research Vessel and Snow Crab surveys. Environmental covariates were sourced from high‐resolution data layers, including physical oceanographic, bathymetric, and seafloor morphometric datasets. The three models captured variability in sea cucumber distribution that would have been overlooked without a spatiotemporal approach. Although their predictions were similar, including within C. frondosa spatial reserves, the models provided different inferences regarding covariate effects. This suggests that while practitioners primarily interested in mapping species distributions need only apply the most familiar framework, those most concerned with identifying predictive environmental covariates may benefit from comparing the output from multiple approaches. Employing multiple approaches can also serve as a validation technique.