Comparing estimates of abundance trends and distribution shifts using single- and multispecies models of fishes and biogenic habitat

Abstract

Several approaches have been developed over the last decade to simultaneously estimate distribution or density for multiple species (e.g. “joint species distribution” or “multispecies occupancy” models). However, there has been little research comparing estimates of abundance trends or distribution shifts from these multispecies models with similar single-species estimates. We seek to determine whether a model including correlations among species (and particularly species that may affect habitat quality, termed “biogenic habitat”) improves predictive performance or decreases standard errors for estimates of total biomass and distribution shift relative to similar single-species models. To accomplish this objective, we apply a vector-autoregressive spatio-temporal (VAST) model that simultaneously estimates spatio-temporal variation in density for multiple species, and present an application of this model using data for eight US Pacific Coast rockfishes (Sebastes spp.), thornyheads (Sebastolobus spp.), and structure-forming invertebrates (SFIs). We identified three fish groups having similar spatial distribution (northern Sebastes, coastwide Sebastes, and Sebastolobus species), and estimated differences among groups in their association with SFI. The multispecies model was more parsimonious and had better predictive performance than fitting a single-species model to each taxon individually, and estimated fine-scale variation in density even for species with relatively few encounters (which the single-species model was unable to do). However, the single-species models showed similar abundance trends and distribution shifts to those of the multispecies model, with slightly smaller standard errors. Therefore, we conclude that spatial variation in density (and annual variation in these patterns) is correlated among fishes and SFI, with congeneric fishes more correlated than species from different genera. However, explicitly modelling correlations among fishes and biogenic habitat does not seem to improve precision for estimates of abundance trends or distribution shifts for these fishes.