Multispecies integrated population model reveals bottom‐up dynamics in a seabird predator–prey system

Multifaceted density dependence: Social structure and seasonality effects on Serengeti lion demography

Interactions between density and environmental conditions have important effects on vital rates and consequently on population dynamics and can take complex pathways in species whose demography is strongly influenced by social context, such as the African lion, Panthera leo. In populations of such species, the response of vital rates to density can vary depending on the social structure (e.g. effects of group size or composition). However, studies assessing density dependence in populations of lions and other social species have seldom considered the effects of multiple socially explicit measures of density, and-more particularly for lions-of nomadic males. Additionally, vital-rate responses to interactions between the environment and various measures of density remain largely uninvestigated. To fill these knowledge gaps, we aimed to understand how a socially and spatially explicit consideration of density (i.e. at the local scale) and its interaction with environmental seasonality affect vital rates of lions in the Serengeti National Park, Tanzania. We used a Bayesian multistate capture-recapture model and Bayesian generalized linear mixed models to estimate lion stage-specific survival and between-stage transition rates, as well as reproduction probability and recruitment, while testing for season-specific effects of density measures at the group and home-range levels. We found evidence for several such effects. For example, resident-male survival increased more strongly with coalition size in the dry season compared with the wet season, and adult-female abundance affected subadult survival negatively in the wet season, but positively in the dry season. Additionally, while our models showed no effect of nomadic males on adult-female survival, they revealed strong effects of nomads on key processes such as reproduction and takeover dynamics. Therefore, our results highlight the importance of accounting for seasonality and social context when assessing the effects of density on vital rates of Serengeti lions and of social species more generally.

Bottom-up effects drive the dynamic of an Antarctic seabird predator-prey system

Understanding how populations respond to variability in environmental conditions and interspecific interactions is one of the biggest challenges of population ecology, particularly in the context of global change. Although many studies have investigated population responses to climate change, very few have explicitly integrated interspecific relationships when studying these responses. In this study, we aimed to understand the combined effects of interspecific interactions and environmental conditions on the demographic parameters of a prey-predator system of three sympatric seabird populations breeding in Antarctica: the south polar skua (Catharacta maccormicki) and its two main preys during the breeding season, the Adélie penguin (Pygoscelis adeliae) and the emperor penguin (Aptenodytes forsteri). We built a two-species integrated population model (IPM) with 31 years of capture-recapture and count data and provided a framework that made it possible to estimate the demographic parameters and abundance of a predator-prey system in a context where capture-recapture data were not available for one species. Our results showed that predator-prey interactions and local environmental conditions differentially affected south polar skuas depending on their breeding state of the previous year. Concerning prey-predator relationships, the number of Adélie penguin breeding pairs showed a positive effect on south polar skua survival and breeding probability, and the number of emperor penguin dead chicks showed a positive effect on the breeding success of south polar skuas. In contrast, there was no evidence for an effect of the number of south polar skuas on the demography of Adélie penguins. We also found an important impact of sea ice conditions on both the dynamics of south polar skuas and Adélie penguins. Our results suggest that this prey-predator system is mostly driven by bottom-up processes and local environmental conditions.

A multispecies hierarchical model to integrate count and distance-sampling data

Integrated community models-an emerging framework in which multiple data sources for multiple species are analyzed simultaneously-offer opportunities to expand inferences beyond the single-species and single-data-source approaches common in ecology. We developed a novel integrated community model that combines distance sampling and single-visit count data; within the model, information is shared among data sources (via a joint likelihood) and species (via a random-effects structure) to estimate abundance patterns across a community. Parameters relating to abundance are shared between data sources, and the model can specify either shared or separate observation processes for each data source. Simulations demonstrated that the model provided unbiased estimates of abundance and detection parameters even when detection probabilities varied between the data types. The integrated community model also provided more accurate and more precise parameter estimates than alternative single-species and single-data-source models in many instances. We applied the model to a community of 11 herbivore species in the Masai Mara National Reserve, Kenya, and found considerable interspecific variation in response to local wildlife management practices: Five species showed higher abundances in a region with passive conservation enforcement (median across species: 4.5× higher), three species showed higher abundances in a region with active conservation enforcement (median: 3.9× higher), and the remaining three species showed no abundance differences between the two regions. Furthermore, the community average of abundance was slightly higher in the region with active conservation enforcement but not definitively so (posterior mean: higher by 0.20 animals; 95% credible interval: 1.43 fewer animals, 1.86 more animals). Our integrated community modeling framework has the potential to expand the scope of inference over space, time, and levels of biological organization, but practitioners should carefully evaluate whether model assumptions are met in their systems and whether data integration is valuable for their applications.

Integrated community models: A framework combining multispecies data sources to estimate the status, trends and dynamics of biodiversity

Data deficiencies among rare or cryptic species preclude assessment of community-level processes using many existing approaches, limiting our understanding of the trends and stressors for large numbers of species. Yet evaluating the dynamics of whole communities, not just common or charismatic species, is critical to understanding and the responses of biodiversity to ongoing environmental pressures. A recent surge in both public science and government-funded data collection efforts has led to a wealth of biodiversity data. However, these data collection programmes use a wide range of sampling protocols (from unstructured, opportunistic observations of wildlife to well-structured, design-based programmes) and record information at a variety of spatiotemporal scales. As a result, available biodiversity data vary substantially in quantity and information content, which must be carefully reconciled for meaningful ecological analysis. Hierarchical modelling, including single-species integrated models and hierarchical community models, has improved our ability to assess and predict biodiversity trends and processes. Here, we highlight the emerging 'integrated community modelling' framework that combines both data integration and community modelling to improve inferences on species- and community-level dynamics. We illustrate the framework with a series of worked examples. Our three case studies demonstrate how integrated community models can be used to extend the geographic scope when evaluating species distributions and community-level richness patterns; discern population and community trends over time; and estimate demographic rates and population growth for communities of sympatric species. We implemented these worked examples using multiple software methods through the R platform via packages with formula-based interfaces and through development of custom code in JAGS, NIMBLE and Stan. Integrated community models provide an exciting approach to model biological and observational processes for multiple species using multiple data types and sources simultaneously, thus accounting for uncertainty and sampling error within a unified framework. By leveraging the combined benefits of both data integration and community modelling, integrated community models can produce valuable information about both common and rare species as well as community-level dynamics, allowing for holistic evaluation of the effects of global change on biodiversity.

Integrated Population Models: Achieving Their Potential

Precise and accurate estimates of abundance and demographic rates are primary quantities of interest within wildlife conservation and management. Such quantities provide insight into population trends over time and the associated underlying ecological drivers of the systems. This information is fundamental in managing ecosystems, assessing species conservation status and developing and implementing effective conservation policy. Observational monitoring data are typically collected on wildlife populations using an array of different survey protocols, dependent on the primary questions of interest. For each of these survey designs, a range of advanced statistical techniques have been developed which are typically well understood. However, often multiple types of data may exist for the same population under study. Analyzing each data set separately implicitly discards the common information contained in the other data sets. An alternative approach that aims to optimize the shared information contained within multiple data sets is to use a “model-based data integration” approach, or more commonly referred to as an “integrated model.” This integrated modeling approach simultaneously analyzes all the available data within a single, and robust, statistical framework. This paper provides a statistical overview of ecological integrated models, with a focus on integrated population models (IPMs) which include abundance and demographic rates as quantities of interest. Four main challenges within this area are discussed, namely model specification, computational aspects, model assessment and forecasting. This should encourage researchers to explore further and develop new practical tools to ensure that full utility can be made of IPMs for future studies.

Pathways to global-change effects on biodiversity: new opportunities for dynamically forecasting demography and species interactions

In structured populations, persistence under environmental change may be particularly threatened when abiotic factors simultaneously negatively affect survival and reproduction of several life cycle stages, as opposed to a single stage. Such effects can then be exacerbated when species interactions generate reciprocal feedbacks between the demographic rates of the different species. Despite the importance of such demographic feedbacks, forecasts that account for them are limited as individual-based data on interacting species are perceived to be essential for such mechanistic forecasting-but are rarely available. Here, we first review the current shortcomings in assessing demographic feedbacks in population and community dynamics. We then present an overview of advances in statistical tools that provide an opportunity to leverage population-level data on abundances of multiple species to infer stage-specific demography. Lastly, we showcase a state-of-the-art Bayesian method to infer and project stage-specific survival and reproduction for several interacting species in a Mediterranean shrub community. This case study shows that climate change threatens populations most strongly by changing the interaction effects of conspecific and heterospecific neighbours on both juvenile and adult survival. Thus, the repurposing of multi-species abundance data for mechanistic forecasting can substantially improve our understanding of emerging threats on biodiversity.

Spatial consistency in drivers of population dynamics of a declining migratory bird

Many migratory species are in decline across their geographical ranges. Single-population studies can provide important insights into drivers at a local scale, but effective conservation requires multi-population perspectives. This is challenging because relevant data are often hard to consolidate, and state-of-the-art analytical tools are typically tailored to specific datasets. We capitalized on a recent data harmonization initiative (SPI-Birds) and linked it to a generalized modelling framework to identify the demographic and environmental drivers of large-scale population decline in migratory pied flycatchers (Ficedula hypoleuca) breeding across Britain. We implemented a generalized integrated population model (IPM) to estimate age-specific vital rates, including their dependency on environmental conditions, and total and breeding population size of pied flycatchers using long-term (34-64 years) monitoring data from seven locations representative of the British breeding range. We then quantified the relative contributions of different vital rates and population structure to changes in short- and long-term population growth rate using transient life table response experiments (LTREs). Substantial covariation in population sizes across breeding locations suggested that change was the result of large-scale drivers. This was supported by LTRE analyses, which attributed past changes in short-term population growth rates and long-term population trends primarily to variation in annual survival and dispersal dynamics, which largely act during migration and/or nonbreeding season. Contributions of variation in local reproductive parameters were small in comparison, despite sensitivity to local temperature and rainfall within the breeding period. We show that both short- and long-term population changes of British breeding pied flycatchers are likely linked to factors acting during migration and in nonbreeding areas, where future research should be prioritized. We illustrate the potential of multi-population analyses for informing management at (inter)national scales and highlight the importance of data standardization, generalized and accessible analytical tools, and reproducible workflows to achieve them.