Robustly estimating the demographic contribution of immigration: Simulation, sensitivity analysis and seals.

Identifying important demographic drivers of population dynamics is fundamental for understanding life-history evolution and implementing effective conservation measures. Integrated population models (IPMs) coupled with transient life table response experiments (tLTREs) allow ecologists to quantify the contributions of demographic parameters to observed population change. While IPMs can estimate parameters that are not estimable using any data source alone, for example, immigration, the estimated contribution of such parameters to population change is prone to bias. Currently, it is unclear when robust conclusions can be drawn from them. We sought to understand the drivers of a rebounding southern elephant seal population on Marion Island using the IPM-tLTRE framework, applied to count and mark-recapture data on 9500 female seals over nearly 40 years. Given the uncertainty around IPM-tLTRE estimates of immigration, we also aimed to investigate the utility of simulation and sensitivity analyses as general tools for evaluating the robustness of conclusions obtained in this framework. Using a Bayesian IPM and tLTRE analysis, we quantified the contributions of survival, immigration and population structure to population growth. We assessed the sensitivity of our estimates to choice of multivariate priors on immigration and other vital rates. To do so we make a novel application of Gaussian process priors, in comparison with commonly used shrinkage priors. Using simulation, we assessed our model's ability to estimate the demographic contribution of immigration under different levels of temporal variance in immigration. The tLTRE analysis suggested that adult survival and immigration were the most important drivers of recent population growth. While the contribution of immigration was sensitive to prior choices, the estimate was consistently large. Furthermore, our simulation study validated the importance of immigration by showing that our estimate of its demographic contribution is unlikely to result as a biased overestimate. Our results highlight the connectivity between distant populations of southern elephant seals, illustrating that female dispersal can be important in regulating the abundance of local populations even when natal site fidelity is high. More generally, we demonstrate how robust ecological conclusions may be obtained about immigration from the IPM-tLTRE framework, by combining sensitivity analysis and simulation.

Ninety years of change, from commercial extinction to recovery, range expansion and decline for Antarctic fur seals at South Georgia.

With environmental change, understanding how species recover from overharvesting and maintain viable populations is central to ecosystem restoration. Here, we reconstruct 90 years of recovery trajectory of the Antarctic fur seal at South Georgia (S.W. Atlantic), a key indicator species in the krill-based food webs of the Southern Ocean. After being harvested to commercial extinction by 1907, this population rebounded and now constitutes the most abundant otariid in the World. However, its status remains uncertain due to insufficient and conflicting data, and anthropogenic pressures affecting Antarctic krill, an essential staple for millions of fur seals and other predators. Using integrated population models, we estimated simultaneously the long-term abundance for Bird Island, northwest South Georgia, epicentre of recovery of the species after sealing, and population adjustments for survey counts with spatiotemporal applicability. Applied to the latest comprehensive survey data, we estimated the population at South Georgia in 2007-2009 as 3,510,283 fur seals [95% CI: 3,140,548-3,919,604] (ca. 98% of global population), after 40 years of maximum growth and range expansion owing to an abundant krill supply. At Bird Island, after 50 years of exponential growth followed by 25 years of slow stable growth, the population collapsed in 2009 and has thereafter declined by -7.2% [-5.2, -9.1] per annum, to levels of the 1970s. For the instrumental record, this trajectory correlates with a time-varying relationship between coupled climate and sea surface temperature cycles associated with low regional krill availability, although the effects of increasing krill extraction by commercial fishing and natural competitors remain uncertain. Since 2015, fur seal longevity and recruitment have dropped, sexual maturation has retarded, and population growth is expected to remain mostly negative and highly variable. Our analysis documents the rise and fall of a key Southern Ocean predator over a century of profound environmental and ecosystem change.

Partitioning variance in population growth for models with environmental and demographic stochasticity.

How demographic factors lead to variation or change in growth rates can be investigated using life table response experiments (LTRE) based on structured population models. Traditionally, LTREs focused on decomposing the asymptotic growth rate, but more recently decompositions of annual 'realized' growth rates using 'transient' LTREs have gained in popularity. Transient LTREs have been used particularly to understand how variation in vital rates translate into variation in growth for populations under long-term study. For these, complete population models may be constructed to investigate how temporal variation in environmental drivers affect vital rates. Such investigations have usually come down to estimating covariate coefficients for the effects of environmental variables on vital rates, but formal ways of assessing how they lead to variation in growth rates have been lacking. We extend transient LTREs to further partition the contributions from vital rates into contributions from temporally varying factors that affect them. The decomposition allows one to compare the resultant effect on the growth rate of different environmental factors, as well as density dependence, which may each act via multiple vital rates. We also show how realized growth rates can be decomposed into separate components from environmental and demographic stochasticity. The latter is typically omitted in LTRE analyses. We illustrate these extensions with an integrated population model (IPM) for data from a 26 years study on northern wheatears (Oenanthe oenanthe), a migratory passerine bird breeding in an agricultural landscape. For this population, consisting of around 50-120 breeding pairs per year, we partition variation in realized growth rates into environmental contributions from temperature, rainfall, population density and unexplained random variation via multiple vital rates, and from demographic stochasticity. The case study suggests that variation in first year survival via the unexplained random component, and adult survival via temperature are two main factors behind environmental variation in growth rates. More than half of the variation in growth rates is suggested to come from demographic stochasticity, demonstrating the importance of this factor for populations of moderate size.

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.

Recovery of critically endangered Nassau grouper (Epinephelus striatus) in the Cayman Islands following targeted conservation actions.

Many large-bodied marine fishes that form spawning aggregations, such as the Nassau grouper (Epinephelus striatus), have suffered regional overfishing due to exploitation during spawning. In response, marine resource managers in many locations have established marine protected areas or seasonal closures to recover these overfished stocks. The challenge in assessing management effectiveness lies largely in the development of accurate estimates to track stock size through time. For the past 15 y, the Cayman Islands government has taken a series of management actions aimed at recovering collapsed stocks of Nassau grouper. Importantly, the government also partnered with academic and nonprofit organizations to establish a research and monitoring program (Grouper Moon) aimed at documenting the impacts of conservation action. Here, we develop an integrated population model of 2 Cayman Nassau grouper stocks based on both diver-collected mark-resight observations and video censuses. Using both data types across multiple years, we fit parameters for a state-space model for population growth. We show that over the last 15 y the Nassau grouper population on Little Cayman has more than tripled in response to conservation efforts. Census data from Cayman Brac, while more sparse, show a similar pattern. These findings demonstrate that spatial and seasonal closures aimed at rebuilding aggregation-based fisheries can foster conservation success.

Age-specific survival rates, causes of death, and allowable take of golden eagles in the western United States.

In the United States, the Bald and Golden Eagle Protection Act prohibits take of golden eagles (Aquila chrysaetos) unless authorized by permit, and stipulates that all permitted take must be sustainable. Golden eagles are unintentionally killed in conjunction with many lawful activities (e.g., electrocution on power poles, collision with wind turbines). Managers who issue permits for incidental take of golden eagles must determine allowable take levels and manage permitted take accordingly. To aid managers in making these decisions in the western United States, we used an integrated population model to obtain estimates of golden eagle vital rates and population size, and then used those estimates in a prescribed take level (PTL) model to estimate the allowable take level. Estimated mean annual survival rates for golden eagles ranged from 0.70 (95% credible interval = 0.66-0.74) for first-year birds to 0.90 (0.88-0.91) for adults. Models suggested a high proportion of adult female golden eagles attempted to breed and breeding pairs fledged a mean of 0.53 (0.39-0.72) young annually. Population size in the coterminous western United States has averaged ~31,800 individuals for several decades, with λ = 1.0 (0.96-1.05). The PTL model estimated a median allowable take limit of ~2227 (708-4182) individuals annually given a management objective of maintaining a stable population. We estimate that take averaged 2572 out of 4373 (59%) deaths annually, based on a representative sample of transmitter-tagged golden eagles. For the subset of golden eagles that were recovered and a cause of death determined, anthropogenic mortality accounted for an average of 74% of deaths after their first year; leading forms of take over all age classes were shooting (~670 per year), collisions (~611), electrocutions (~506), and poisoning (~427). Although observed take overlapped the credible interval of our allowable take estimate and the population overall has been stable, our findings indicate that additional take, unless mitigated for, may not be sustainable. Our analysis demonstrates the utility of the joint application of integrated population and prescribed take level models to management of incidental take of a protected species.

Can attraction to and competition for high-quality habitats shape breeding propensity?

In many animal species, sexually mature individuals may skip breeding opportunities despite a likely negative impact on fitness. In spatio-temporally heterogeneous environments, habitat selection theory predicts that individuals select habitats where fitness prospects are maximized. Individuals are attracted to high-quality habitat patches where they compete for high-quality breeding sites. Since failures in contests to secure a site may prevent individuals from breeding, we hypothesized that attraction to and competition for high-quality habitats could shape breeding propensity. Under this hypothesis, we predicted the two following associations between breeding propensity and two key population features. (1) When mean habitat quality in the population increases in multiple patches such that availability of high-quality sites increases across the population, the resulting decrease in competition should positively affect breeding propensity. (2) When the number of individuals increases in the population, the resulting increase in competitors should negatively affect breeding propensity (negative density dependence). Using long-term data from kittiwakes Rissa tridactyla, we checked the prerequisite of prediction (1), that availability of high-quality sites is positively associated with current mean habitat quality in the population (represented by breeding success). We then applied integrated population modelling to quantify annual fluctuations in population mean breeding success, breeding propensity and number of individuals by breeding status (pre-breeders, breeders, skippers and immigrants), and tested our predictions. Our results showed that breeding propensity acts as an important driver of population growth. As expected, breeding propensity was positively associated with preceding mean habitat quality in the population, and negatively with the number of competitors. These relationships varied depending on breeding status, which likely reflects status dependence in competitive ability. These findings highlight the importance of competition for high-quality breeding sites in shaping breeding propensity. Thereby, we draw attention towards alternative and complementary explanations to more standard considerations regarding the energetic cost of reproduction, and point to possible side effects of habitat selection behaviours on individual life histories and population dynamics.

Density-dependence produces spurious relationships among demographic parameters in a harvested species.

Harvest of wild organisms is an important component of human culture, economy, and recreation, but can also put species at risk of extinction. Decisions that guide successful management actions therefore rely on the ability of researchers to link changes in demographic processes to the anthropogenic actions or environmental changes that underlie variation in demographic parameters. Ecologists often use population models or maximum sustained yield curves to estimate the impacts of harvest on wildlife and fish populations. Applications of these models usually focus exclusively on the impact of harvest and often fail to consider adequately other potential, often collinear, mechanistic drivers of the observed relationships between harvest and demographic rates. In this study, we used an integrated population model and long-term data (1973-2016) to examine the relationships among hunting and natural mortality, the number of hunters, habitat conditions, and population size of blue-winged teal Spatula discors, an abundant North American dabbling duck with a relatively fast-paced life history strategy. Over the last two and a half decades of the study, teal abundance tripled, hunting mortality probability increased slightly ( < 0.02 ), and natural mortality probability increased substantially ( > 0.1 ) at greater population densities. We demonstrate strong density-dependent effects on natural mortality and fecundity as population density increased, indicative of compensatory harvest mortality and compensatory natality. Critically, an analysis that only assessed the relationship between survival and hunting mortality would spuriously indicate depensatory mortality due to multicollinearity between abundance, natural mortality and hunting mortality. Our findings demonstrate that models that only consider the direct effect of hunting on survival or natural mortality can fail to accurately assess the mechanistic impact of hunting on population dynamics due to multicollinearity among demographic drivers. This multicollinearity limits inference and may have strong impacts on applied management actions globally.

Partial and complete dependency among data sets has minimal consequence on estimates from integrated population models.

Integrated population models (IPMs) are widely used to combine disparate data sets in joint analysis to better understand population dynamics and provide guidance for conservation activities. An often-cited assumption of IPMs is independence among component data sets within the combined likelihood. Dependency among data sets should lead to underestimation of variance and bias because individuals contribute data to more than one data set. In practice, studied individuals often occur in multiple data sets in IPMs (i.e., overlap), which is one way for the independence assumption to be violated. Such cases have the potential to dissuade practitioners and limit application of IPMs to solve emerging ecological problems. We assessed precision and bias of demographic rates estimated from IPMs using a complete gradient (0-100%) of overlap among data sets, wide ranges in demographic rates (e.g., survival 0.1-0.8) and sample sizes (100-1,200 individuals) and variable data sources. We compared results from our simulations with those from IPMs constructed using empirical data on tree swallows (Tachycineta bicolor) where data sets either had complete overlap or included different individuals. Contrary to previous investigators, we found no substantive bias or uncertainty in any demographic rate from IPMs derived from data sets with complete overlap. While variability in demographic rates was greater at low sample sizes (i.e., low capture, recapture, and survey probabilities), there were negligible differences in the posterior mean or root mean square error of demographic rates among IPMs with strong dependence vs. complete independence among data sets. Our simulations suggest IPMs can be designed using only capture-recapture data or harvest and capture-recovery data where population estimates are obtained from the same data as survival and productivity data. While we encourage researchers to carefully consider the modeling approach best suited for their data sets, our results suggest that dependence among data sets does not generally compromise IPM estimates. Thus, violation of the independence assumption should not dissuade researchers from the application of IPMs in ecological research.

Decomposing fecundity and evaluating demographic influence of multiple broods in a migratory bird.

Relevance of breeding season fecundity as a driver of population dynamics has been highlighted by many studies. Despite that, knowledge about how brood type specific (i.e. first, second or replacement) fecundity affects demography of multiple-brooded species is limited. In fact, estimation of brood type specific fecundity is often challenging due to imperfect detection of nesting attempts. We examined the demographic contribution and the feedback on population density of different components of fecundity, along with other vital rates, in a facultative multiple-brooded migratory bird. We used a novel formulation of a fecundity model that allows incorporating reproductive data for which information on the type of brood was unknown in some cases, and embedded it into an integrated population model (IPM) to obtain consensual estimates of all demographic rates, including brood type specific fecundities, reproductive success probabilities and proportion of breeding pairs that performed a second or replacement brood. We then conducted transient life table response experiments on IPM estimates to account for non-stationary environments. We applied the model to two 20-year datasets collected in a Swiss and a German local population of wrynecks Jynx torquilla. Brood type specific fecundities and temporal patterns of brood type specific probabilities of success, number of successful and unsuccessful first broods, probability of starting a second or a replacement brood and proportion of pairs that performed a second or a replacement brood differed between the two populations. However, changes in immigration rate and apparent survival were the dominant contributors to temporal variation and large sequential changes in realized population growth rates in both populations. In the Swiss population we also found that second brood fecundity declined when population size increased. Our study provides insight into the reproductive processes that affect population dynamics and mediate density-dependent fecundity in a migratory bird. In addition, the analytical approach proposed can be used in other studies of multiple-brooded species to maximize the use of available fecundity data through the estimation of unknown brood types, thus favouring a better understanding of the demographic contribution of brood type specific fecundity.

Interrelated impacts of climate and land-use change on a widespread waterbird.

Together climate and land-use change play a crucial role in determining species distribution and abundance, but measuring the simultaneous impacts of these processes on current and future population trajectories is challenging due to time lags, interactive effects and data limitations. Most approaches that relate multiple global change drivers to population changes have been based on occurrence or count data alone. We leveraged three long-term (1995-2019) datasets to develop a coupled integrated population model-Bayesian population viability analysis (IPM-BPVA) to project future survival and reproductive success for common loons Gavia immer in northern Wisconsin, USA, by explicitly linking vital rates to changes in climate and land use. The winter North Atlantic Oscillation (NAO), a broad-scale climate index, immediately preceding the breeding season and annual changes in developed land cover within breeding areas both had strongly negative influences on adult survival. Local summer rainfall was negatively related to fecundity, though this relationship was mediated by a lagged interaction with the winter NAO, suggesting a compensatory population-level response to climate variability. We compared population viability under 12 future scenarios of annual land-use change, precipitation and NAO conditions. Under all scenarios, the loon population was expected to decline, yet the steepest declines were projected under positive NAO trends, as anticipated with ongoing climate change. Thus, loons breeding in the northern United States are likely to remain affected by climatic processes occurring thousands of miles away in the North Atlantic during the non-breeding period of the annual cycle. Our results reveal that climate and land-use changes are differentially contributing to loon population declines along the southern edge of their breeding range and will continue to do so despite natural compensatory responses. We also demonstrate that concurrent analysis of multiple data types facilitates deeper understanding of the ecological implications of anthropogenic-induced change occurring at multiple spatial scales. Our modelling approach can be used to project demographic responses of populations to varying environmental conditions while accounting for multiple sources of uncertainty, an increasingly pressing need in the face of unprecedented global change.

Climate-driven carry-over effects negatively influence population growth rate in a food-caching boreal passerine.

Understanding how events throughout the annual cycle are linked is important for predicting variation in individual fitness, but whether and how carry-over effects scale up to influence population dynamics is poorly understood. Using 38 years of demographic data from Algonquin Provincial Park, Ontario, and a full annual cycle integrated population model, we examined the influence of environmental conditions and density on the population growth rate of Canada jays (Perisoreus canadensis), a resident boreal passerine that relies on perishable cached food for over-winter survival and late-winter breeding. Our results demonstrate that fall environmental variables, most notably the number of freeze-thaw events, carried over to influence late-winter fecundity, which, in turn, was the main vital rate driving population growth. These results are consistent with the hypothesis that warmer and more variable fall conditions accelerate the degradation of perishable stored food that is relied upon for successful reproduction. Future warming during the fall and winter may compromise the viability of cached food that requires consistent subzero temperatures for effective preservation, potentially exacerbating climate-driven carry-over effects that impact long-term population dynamics.

Why we should care about movements: Using spatially explicit integrated population models to assess habitat source-sink dynamics.

Assessing the source-sink status of populations and habitats is of major importance for understanding population dynamics and for the management of natural populations. Sources produce a net surplus of individuals (per capita contribution to the metapopulation > 1) and will be the main contributors for self-sustaining populations, whereas sinks produce a deficit (contribution < 1). However, making these types of assessments is generally hindered by the problem of separating mortality from permanent emigration, especially when survival probabilities as well as moved distances are habitat-specific. To address this long-standing issue, we propose a spatial multi-event integrated population model (IPM) that incorporates habitat-specific dispersal distances of individuals. Using information about local movements, this IPM adjusts survival estimates for emigration outside the study area. Analysing 24 years of data on a farmland passerine (the northern wheatear Oenanthe oenanthe), we assessed habitat-specific contributions, and hence the source-sink status and temporal variation of two key breeding habitats, while accounting for habitat- and sex-specific local dispersal distances of juveniles and adults. We then examined the sensitivity of the source-sink analysis by comparing results with and without accounting for these local movements. Estimates of first-year survival, and consequently habitat-specific contributions, were higher when local movement data were included. The consequences from including movement data were sex specific, with contribution shifting from sink to likely source in one habitat for males, and previously noted habitat differences for females disappearing. Assessing the source-sink status of habitats is extremely challenging. We show that our spatial IPM accounting for local movements can reduce biases in estimates of the contribution by different habitats, and thus reduce the overestimation of the occurrence of sink habitats. This approach allows combining all available data on demographic rates and movements, which will allow better assessment of source-sink dynamics and better informed conservation interventions.

Density-dependent population dynamics of a high Arctic capital breeder, the barnacle goose.

Density regulation of the population growth rate occurs through negative feedbacks on underlying vital rates, in response to increasing population size. Here, we examine in a capital breeder how vital rates of different life-history stages, their elasticities and population growth rates are affected by changes in population size. We developed an integrated population model for a local population of Svalbard barnacle geese, Branta leucopsis, using counts, reproductive data and individual-based mark-recapture data (1990-2017) to model age class-specific survival, reproduction and number of individuals. Based on these estimates, we quantified the changes in demographic structure and the effect of population size on age class-specific vital rates and elasticities, as well as the population growth rate. Local density regulation at the breeding grounds acted to reduce population growth through negative effects on reproduction; however, population size could not explain substantial variation in survival rates, although there was some support for density-dependent first-year survival. With the use of prospective perturbation analysis of the density-dependent projection matrix, we show that the elasticities to different vital rates changed as population size increased. As population size approached carrying capacity, the influence of reproductive rates and early-life survival on the population growth rate was reduced, whereas the influence of adult survival increased. A retrospective perturbation analysis revealed that density dependence resulted in a positive contribution of reproductive rates, and a negative contribution of the numbers of individuals in the adult age class, to the realised population growth rate. The patterns of density dependence in this population of barnacle geese were different from those recorded in income breeding birds, where density regulation mainly occurs through an effect on early-life survival. This indicates that the population dynamics of capital breeders, such as the barnacle goose, are likely to be more reproduction-driven than is the case for income breeders.

Metapopulation dynamics of roseate terns: Sources, sinks and implications for conservation management decisions.

Habitat management to restore or create breeding sites may allow metapopulations to increase in size and reduce the risk of demographic stochasticity or disasters causing metapopulation extinction. However, if newly restored or created sites are of low quality, they may act as sinks that draw individuals away from better quality sites to the detriment of metapopulation size. Following intensive conservation effort, the metapopulation of roseate tern (Sterna dougallii) in NW Europe is recovering from a large crash in numbers, but most former colonies remain unoccupied and hence are potential targets for restoration. To inform conservation efforts, we studied the dynamics of this metapopulation with a multistate integrated population model to assess each of the three main colonies for important demographic contributors to population growth rate, source/sink status and possible density dependence. All three study colonies are managed for roseate terns (and other tern species) in similar ways, but the demographic processes vary considerably between colonies. The largest colony is a source involved in almost all dispersal, and its growth is determined by survival rates and productivity. Productivity and juvenile apparent survival at the largest colony appear to be density-dependent. Although the mechanisms are unclear, this may provide an increasing impetus for emigration of recruits to other colonies in future. The smallest of the three colonies is a sink, relying on immigration for its growth. Simulation models suggest the metapopulation would be c. 10% larger in the absence of dispersal to the sink colony. This work indicates that, due to variable site quality, aims to enhance both distribution and size of metapopulations may be mutually exclusive. In this case, before future attempts to encourage recolonisation of former sites, assessments of site suitability should be undertaken, focusing on food availability and isolation from predators to maximise the likelihood of attaining levels of productivity and survival that avoid creation of a sink population to the detriment of the overall metapopulation size.

Immature survival, fertility, and density dependence drive global population dynamics in a long-lived species.

Disentangling the influence of demographic parameters and the role of density dependence on species' population dynamics is a challenge, especially when fractions of the population are unobservable. Additionally, due to the difficulty of gathering data at large spatial scales, most studies ignore the global dynamic of a species, which would integrate local heterogeneity dynamics and remove the noise of dispersal. We developed an integrated population model (IPM) at a global scale to disentangle the main demographic drivers of population dynamics in a long-lived species. We used 28 yr of Audouin's Gull demographic data encompassing 69 local patches (comprising 90% of the world population). Importantly, we took into account the unobservable fraction of non-breeders and also assessed the strength of density dependence for this fraction of the population. As predicted by life histories of long-lived organisms, temporal random variation in survival was highest for immature individuals (1.326, 95% credible interval [CRI] 1.290-1.940) and lowest for adults (0.499, 95% CRI 0.487-0.720). Large temporal fluctuations in the probability of taking a reproductive sabbatical would partly explain the consistency in adult survival, with individuals most likely refraining from breeding when environmental conditions were harsh. Immature survival and fertility were the main drivers of population dynamics during the study period (r2  = 0.83, 0.77-0.87 and 0.73, 0.63-0.79, respectively). We found strong evidence of density dependence, not only due to the number of breeders (r2  = -0.34, -0.43 to -0.24) but also due to individuals on sabbatical (r2  = -0.18, -0.33 to -0.01). From a conservation point of view, the species shows a 5% annual global decrease during the last 10 years, and we propose an update of its conservation status. Even though population dynamics of long-lived organisms are very sensitive to changes in adult survival, we show here that, in the absence of strong environmental perturbations affecting this vital rate, fluctuations in population density are mainly driven by variations in survival of immature individuals and fertility. Integrated models based on long-term monitoring at a global scale may enhance our ecological and evolutionary understanding of how demographic drivers influence population dynamics.

Spatial and temporal drivers of avian population dynamics across the annual cycle.

Untangling the spatial and temporal processes that influence population dynamics of migratory species is challenging, because changes in abundance are shaped by variation in vital rates across heterogeneous habitats and throughout the annual cycle. We developed a full-annual-cycle, integrated, population model and used demographic data collected between 2011 and 2014 in southern Indiana and Belize to estimate stage-specific vital rates of a declining migratory songbird, the Wood Thrush (Hylocichla mustelina). Our primary objective was to understand how spatial and temporal variation in demography contributes to local and regional population growth. Our full-annual-cycle model allowed us to estimate (1) age-specific, seasonal survival probabilities, including latent survival during both spring and autumn migration, and (2) how the relative contribution of vital rates to population growth differed among habitats. Wood Thrushes in our study populations experienced the lowest apparent survival rates during migration and apparent survival was lower during spring migration than during fall migration. Both mortality and high dispersal likely contributed to low apparent survival during spring migration. Population growth in high-quality habitat was most sensitive to variation in fecundity and apparent survival of juveniles during spring migration, whereas population growth in low-quality sites was most sensitive to adult apparent breeding-season survival. These results elucidate how full-annual-cycle vital rates, particularly apparent survival during migration, interact with spatial variation in habitat quality to influence population dynamics in migratory species.

Integrating count and detection-nondetection data to model population dynamics.

There is increasing need for methods that integrate multiple data types into a single analytical framework as the spatial and temporal scale of ecological research expands. Current work on this topic primarily focuses on combining capture-recapture data from marked individuals with other data types into integrated population models. Yet, studies of species distributions and trends often rely on data from unmarked individuals across broad scales where local abundance and environmental variables may vary. We present a modeling framework for integrating detection-nondetection and count data into a single analysis to estimate population dynamics, abundance, and individual detection probabilities during sampling. Our dynamic population model assumes that site-specific abundance can change over time according to survival of individuals and gains through reproduction and immigration. The observation process for each data type is modeled by assuming that every individual present at a site has an equal probability of being detected during sampling processes. We examine our modeling approach through a series of simulations illustrating the relative value of count vs. detection-nondetection data under a variety of parameter values and survey configurations. We also provide an empirical example of the model by combining long-term detection-nondetection data (1995-2014) with newly collected count data (2015-2016) from a growing population of Barred Owl (Strix varia) in the Pacific Northwest to examine the factors influencing population abundance over time. Our model provides a foundation for incorporating unmarked data within a single framework, even in cases where sampling processes yield different detection probabilities. This approach will be useful for survey design and to researchers interested in incorporating historical or citizen science data into analyses focused on understanding how demographic rates drive population abundance.

An integrated data model to estimate spatiotemporal occupancy, abundance, and colonization dynamics.

Ecological invasions and colonizations occur dynamically through space and time. Estimating the distribution and abundance of colonizing species is critical for efficient management or conservation. We describe a statistical framework for simultaneously estimating spatiotemporal occupancy and abundance dynamics of a colonizing species. Our method accounts for several issues that are common when modeling spatiotemporal ecological data including multiple levels of detection probability, multiple data sources, and computational limitations that occur when making fine-scale inference over a large spatiotemporal domain. We apply the model to estimate the colonization dynamics of sea otters (Enhydra lutris) in Glacier Bay, in southeastern Alaska.

Integrated population modeling reveals the impact of climate on the survival of juvenile emperor penguins.

Early-life demographic traits are poorly known, impeding our understanding of population processes and sensitivity to climate change. Survival of immature individuals is a critical component of population dynamics and recruitment in particular. However, obtaining reliable estimates of juvenile survival (i.e., from independence to first year) remains challenging, as immatures are often difficult to observe and to monitor individually in the field. This is particularly acute for seabirds, in which juveniles stay at sea and remain undetectable for several years. In this work, we developed a Bayesian integrated population model to estimate the juvenile survival of emperor penguins (Aptenodytes forsteri), and other demographic parameters including adult survival and fecundity of the species. Using this statistical method, we simultaneously analyzed capture-recapture data of adults, the annual number of breeding females, and the number of fledglings of emperor penguins collected at Dumont d'Urville, Antarctica, for the period 1971-1998. We also assessed how climate covariates known to affect the species foraging habitats and prey [southern annular mode (SAM), sea ice concentration (SIC)] affect juvenile survival. Our analyses revealed that there was a strong evidence for the positive effect of SAM during the rearing period (SAMR) on juvenile survival. Our findings suggest that this large-scale climate index affects juvenile emperor penguins body condition and survival through its influence on wind patterns, fast ice extent, and distance to open water. Estimating the influence of environmental covariates on juvenile survival is of major importance to understand the impacts of climate variability and change on the population dynamics of emperor penguins and seabirds in general and to make robust predictions on the impact of climate change on marine predators.