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.

The effect of prey abundance and fisheries on the survival, reproduction, and social structure of killer whales (Orcinus orca) at subantarctic Marion Island.

Most marine apex predators are keystone species that fundamentally influence their ecosystems through cascading top-down processes. Reductions in worldwide predator abundances, attributed to environmental- and anthropogenic-induced changes to prey availability and negative interactions with fisheries, can have far-reaching ecosystem impacts. We tested whether the survival of killer whales (Orcinus orca) observed at Marion Island in the Southern Indian Ocean correlated with social structure and prey variables (direct measures of prey abundance, Patagonian toothfish fishery effort, and environmental proxies) using multistate models of capture-recapture data spanning 12 years (2006-2018). We also tested the effect of these same variables on killer whale social structure and reproduction measured over the same period. Indices of social structure had the strongest correlation with survival, with higher sociality associated with increased survival probability. Survival was also positively correlated with Patagonian toothfish fishing effort during the previous year, suggesting that fishery-linked resource availability is an important determinant of survival. No correlation between survival and environmental proxies of prey abundance was found. At-island prey availability influenced the social structure of Marion Island killer whales, but none of the variables explained variability in reproduction. Future increases in legal fishing activity may benefit this population of killer whales through the artificial provisioning of resources they provide.

Determinants of moult haulout phenology and duration in southern elephant seals.

Phenological shifts are among the most obvious biological responses to environmental change, yet documented responses for Southern Ocean marine mammals are extremely rare. Marine mammals can respond to environmental changes through phenological flexibility of their life-history events such as breeding and moulting. Southern elephant seals (Mirounga leonina) undergo an obligatory annual moult which involves the rapid shedding of epidermal skin and hair while seals fast ashore. We quantified the timing (phenology) and duration (the time from arrival ashore to departure) of the moult haulout of 4612 female elephant seals at Marion Island over 32 years. Using linear mixed-effects models, we investigated age, breeding state and environmental drivers of moult timing and haulout duration. We found no clear evidence for a temporal shift in moult phenology or its duration. Annual variation in moult arrival date and haulout duration was small relative to age and breeding effects, which explained more than 90% of the variance in moult arrival date and 25% in moult haulout duration. All environmental covariates we tested explained minimal variation in the data. Female elephant seals moulted progressively later as juveniles, but adults age 4 and older had similar moult start dates that depended on the breeding state of the female. In contrast, moult haulout duration was not constant with age among adults, but instead became shorter with increasing age. Moulting is energetically expensive and differences in the moult haulout duration are possibly due to individual variation in body mass and associated metabolizable energy reserves, although other drivers (e.g. hormones) may also be present. Individual-based data on moult arrival dates and haulout duration can be used as auxiliary data in demographic modelling and may be useful proxies of other important biological parameters such as body condition and breeding history.

Positive early-late life-history trait correlations in elephant seals.

Correlations between early- and late-life performance are a major prediction of life-history theory. Negative early-late correlations can emerge because biological processes are optimized for early but not late life (e.g., rapid development may accelerate the onset of senescence; "developmental theory of aging") or because allocation to early-life performance comes at a cost in terms of late-life performance (as in the disposable soma theory). But variation in genetic and environmental challenges that each individual has to cope with during early life may also lead to positive early-late life-history trait correlations (the "fixed heterogeneity" or "individual quality" hypothesis). We analyzed individual life-history trajectories of 7,420 known-age female southern elephant seals (Mirounga leonina) monitored over 36 yr to determine how actuarial senescence (a proxy for late-life performance) correlate with age at first reproduction (a proxy for early-life performance). As some breeding events may not be detected in this field study, we used a custom "multievent" hierarchical model to estimate the age at first reproduction and correlate it to other life-history traits. The probability of first reproduction was 0.34 at age 3, with most females breeding for the first time at age 4, and comparatively few at older ages. Females with an early age of first reproduction outperformed delayed breeders in all aspects we considered (survival, rate of senescence, net reproductive output) but one: early breeders appeared to have an onset of actuarial senescence 1 yr earlier compared to late breeders. Genetics and environmental conditions during early life likely explain the positive correlation between early- and late-life performance. Our results provide the first evidence of actuarial senescence in female southern elephant seals.

Abundance of Indo-Pacific bottlenose dolphins (Tursiops aduncus) along the south coast of South Africa.

Coastally distributed dolphin species are vulnerable to a variety of anthropogenic pressures, yet a lack of abundance data often prevents data-driven conservation management strategies from being implemented. We investigated the abundance of Indo-Pacific bottlenose dolphins (Tursiops aduncus) along the south coast of South Africa, from the Goukamma Marine Protected Area (MPA) to the Tsitsikamma MPA, between 2014 and 2016. During this period, 662.3h of boat-based photo-identification survey effort was carried out during 189 surveys. The sighting histories of 817 identified individuals were used to estimate abundance using capture-recapture modelling. Using open population (POPAN) models, we estimated that 2,155 individuals (95% CI: 1,873-2,479) occurred in the study area, although many individuals appeared to be transients. We recorded smaller group sizes and an apparent decline in abundance in a subset of the study area (Plettenberg Bay) compared to estimates obtained in 2002-2003 at this location. We recorded declines of more than 70% in both abundance and group size for a subset of the study area (Plettenberg Bay), in relation to estimates obtained in 2002-2003 at this location. We discuss plausible hypotheses for causes of the declines, including anthropogenic pressure, ecosystem change, and methodological inconsistencies. Our study highlights the importance of assessing trends in abundance at other locations to inform data-driven conservation management strategies of T. aduncus in South Africa.

Trade-offs between age-related breeding improvement and survival senescence in highly polygynous elephant seals: Dominant males always do better.

Life history trade-off theory predicts that current reproduction can negatively affect survival and future reproduction. Few studies have assessed breeding costs for males of polygynous species compared to females, despite substantial variation in breeding success among individual males (e.g. subordinate cf. dominant breeders). Specifically, differentiating between the cost of attending breeding seasons, and the additional cost of successfully securing and mating females is lacking. We investigated whether trade-offs are present in the highly polygynous male southern elephant seal (Mirounga leonina) using 34-years of individual-level data. We compare age-specific survival, recruitment and future breeding success probabilities of pre-breeders (males yet to recruit) and breeders (subordinate and dominant social ranks) using multievent models. Pre-breeders and breeders of overlapping ages had similar survival probabilities, suggesting that there was no attendance cost for early recruits. In addition, the probability of recruiting as a dominant breeder never exceeded recruitment probability as a subordinate breeder of the same age. Therefore, older pre-breeders that delayed attendance costs generally did not improve their breeding success (probability of being dominant) at recruitment more than younger recruits. Rather, recruitment age may be a function of individual quality, with lower quality individuals requiring more time to socially mature. When comparing subordinate and dominant breeders, we found clear evidence for survival senescence, with subordinate breeders having a higher baseline mortality. In contrast, age-specific future breeding success (probability of being dominant at t + 1) increased with age, with dominant breeders maintaining higher subsequent breeding success than subordinate breeders. The opposite trends in survival and future breeding success for both subordinate and dominant breeders may indicate a lifetime, population-level trade-off. However, we found no evidence to suggest that being a dominant breeder consecutively (and having a higher accumulated breeding cost) accelerated the rate of senescence when compared to individuals that were previously subordinate. Thus, males experienced actuarial senescence regardless of social rank, with dominant (and possibly high quality) breeders showing a reduced trade-off between survival and future breeding success. We make several novel contributions to understanding polygynous male life histories and southern elephant seal demography.

Making use of multiple surveys: Estimating breeding probability using a multievent-robust design capture-recapture model.

Increased environmental stochasticity due to climate change will intensify temporal variance in the life-history traits, and especially breeding probabilities, of long-lived iteroparous species. These changes may decrease individual fitness and population viability and is therefore important to monitor. In wild animal populations with imperfect individual detection, breeding probabilities are best estimated using capture-recapture methods. However, in many vertebrate species (e.g., amphibians, turtles, seabirds), nonbreeders are unobservable because they are not tied to a territory or breeding location. Although unobservable states can be used to model temporary emigration of nonbreeders, there are disadvantages to having unobservable states in capture-recapture models. The best solution to deal with unobservable life-history states is therefore to eliminate them altogether. Here, we achieve this objective by fitting novel multievent-robust design models which utilize information obtained from multiple surveys conducted throughout the year. We use this approach to estimate annual breeding probabilities of capital breeding female elephant seals (Mirounga leonina). Conceptually, our approach parallels a multistate version of the Barker/robust design in that it combines robust design capture data collected during discrete breeding seasons with observations made at other times of the year. A substantial advantage of our approach is that the nonbreeder state became "observable" when multiple data sources were analyzed together. This allowed us to test for the existence of state-dependent survival (with some support found for lower survival in breeders compared to nonbreeders), and to estimate annual breeding transitions to and from the nonbreeder state with greater precision (where current breeders tended to have higher future breeding probabilities than nonbreeders). We used program E-SURGE (2.1.2) to fit the multievent-robust design models, with uncertainty in breeding state assignment (breeder, nonbreeder) being incorporated via a hidden Markov process. This flexible modeling approach can easily be adapted to suit sampling designs from numerous species which may be encountered during and outside of discrete breeding seasons.

Correction: Return Customers: Foraging Site Fidelity and the Effect of Environmental Variability in Wide-Ranging Antarctic Fur Seals.

[This corrects the article DOI: 10.1371/journal.pone.0120888.].

Return customers: foraging site fidelity and the effect of environmental variability in wide-ranging antarctic fur seals.

Strategies employed by wide-ranging foraging animals involve consideration of habitat quality and predictability and should maximise net energy gain. Fidelity to foraging sites is common in areas of high resource availability or where predictable changes in resource availability occur. However, if resource availability is heterogeneous or unpredictable, as it often is in marine environments, then habitat familiarity may also present ecological benefits to individuals. We examined the winter foraging distribution of female Antarctic fur seals, Arctocephalus gazelle, over four years to assess the degree of foraging site fidelity at two scales; within and between years. On average, between-year fidelity was strong, with most individuals utilising more than half of their annual foraging home range over multiple years. However, fidelity was a bimodal strategy among individuals, with five out of eight animals recording between-year overlap values of greater than 50%, while three animals recorded values of less than 5%. High long-term variance in sea surface temperature, a potential proxy for elevated long-term productivity and prey availability, typified areas of overlap. Within-year foraging site fidelity was weak, indicating that successive trips over the winter target different geographic areas. We suggest that over a season, changes in prey availability are predictable enough for individuals to shift foraging area in response, with limited associated energetic costs. Conversely, over multiple years, the availability of prey resources is less spatially and temporally predictable, increasing the potential costs of shifting foraging area and favouring long-term site fidelity. In a dynamic and patchy environment, multi-year foraging site fidelity may confer a long-term energetic advantage to the individual. Such behaviours that operate at the individual level have evolutionary and ecological implications and are potential drivers of niche specialization and modifiers of intra-specific competition.

Satellite tagging and biopsy sampling of killer whales at subantarctic Marion Island: effectiveness, immediate reactions and long-term responses.

Remote tissue biopsy sampling and satellite tagging are becoming widely used in large marine vertebrate studies because they allow the collection of a diverse suite of otherwise difficult-to-obtain data which are critical in understanding the ecology of these species and to their conservation and management. Researchers must carefully consider their methods not only from an animal welfare perspective, but also to ensure the scientific rigour and validity of their results. We report methods for shore-based, remote biopsy sampling and satellite tagging of killer whales Orcinus orca at Subantarctic Marion Island. The performance of these methods is critically assessed using 1) the attachment duration of low-impact minimally percutaneous satellite tags; 2) the immediate behavioural reactions of animals to biopsy sampling and satellite tagging; 3) the effect of researcher experience on biopsy sampling and satellite tagging; and 4) the mid- (1 month) and long- (24 month) term behavioural consequences. To study mid- and long-term behavioural changes we used multievent capture-recapture models that accommodate imperfect detection and individual heterogeneity. We made 72 biopsy sampling attempts (resulting in 32 tissue samples) and 37 satellite tagging attempts (deploying 19 tags). Biopsy sampling success rates were low (43%), but tagging rates were high with improved tag designs (86%). The improved tags remained attached for 26±14 days (mean ± SD). Individuals most often showed no reaction when attempts missed (66%) and a slight reaction-defined as a slight flinch, slight shake, short acceleration, or immediate dive-when hit (54%). Severe immediate reactions were never observed. Hit or miss and age-sex class were important predictors of the reaction, but the method (tag or biopsy) was unimportant. Multievent trap-dependence modelling revealed considerable variation in individual sighting patterns; however, there were no significant mid- or long-term changes following biopsy sampling or tagging.