Metapopulation regulation acts at multiple spatial scales: Insights from a century of seabird colony census data

Abstract Density‐dependent feedback is recognized as important regulatory mechanisms of population size. Considering the spatial scales over which such feedback operates has advanced our theoretical understanding of metapopulation dynamics. Yet, metapopulation models are rarely fit to time‐series data and tend to omit details of the natural history and behavior of long‐lived, highly mobile species such as colonial mammals and birds. Seabird metapopulations consist of breeding colonies that are connected across large spatial scales, within a heterogeneous marine environment that is increasingly affected by anthropogenic disturbance. Currently, we know little about the strength and spatial scale of density‐dependent regulation and connectivity between colonies. Thus, many important seabird conservation and management decisions rely on outdated assumptions of closed populations that lack density‐dependent regulation. We investigated metapopulation dynamics and connectivity in an exemplar seabird species, the Northern gannet ( Morus bassanus ), using more than a century of census data of breeding colonies distributed across the Northeast Atlantic. We developed and fitted these data to a novel hierarchical Bayesian state‐space model, to compare increasingly complex scenarios of metapopulation regulation through lagged, local, regional, and global density dependence, as well as different mechanisms for immigration. Models with conspecific attraction fit the data better than the equipartitioning of immigrants. Considering local and regional density dependence jointly improved model fit slightly, but importantly, future colony size projections based on different mechanistic regulatory scenarios varied widely: a model with local and regional dynamics estimated a lower metapopulation capacity (645,655 Apparently Occupied Site [AOS]) and consequently higher present saturation (63%) than a model with local density dependence (1,367,352 AOS, 34%). Our findings suggest that metapopulation regulation in the gannet is more complex than traditionally assumed, and highlight the importance of using models that consider colony connectivity and regional dynamics for conservation management applications guided by precautionary principles. Our study advances our understanding of metapopulation dynamics in long‐lived colonial species and our approach provides a template for the development of metapopulation models for colonially living birds and mammals.