Bridging gaps in demographic analysis with phylogenetic imputation

Phylogenetically informed imputation methods have rarely been applied to estimate missing values in demographic data but may be a powerful tool for reconstructing vital rates of survival, maturation, and fecundity for species of conservation concern. Imputed vital rates could be used to parameterize demographic models to explore how populations respond when vital rates are perturbed. We used standardized vital rate estimates for 50 bird species to assess the use of phylogenetic imputation to fill gaps in demographic data. We calculated imputation accuracy for vital rates of focal species excluded from the data set either singly or in combination and with and without phylogeny, body mass, and life-history trait data. We used imputed vital rates to calculate demographic metrics, including generation time, to validate the use of imputation in demographic analyses. Covariance among vital rates and other trait data provided a strong basis to guide imputation of missing vital rates in birds, even in the absence of phylogenetic information. Mean NRMSE for null and phylogenetic models differed by <0.01 except when no vital rates were available or for vital rates with high phylogenetic signal (Pagel's λ > 0.8). In these cases, including body mass and life-history trait data compensated for lack of phylogenetic information: mean normalized root mean square error (NRMSE) for null and phylogenetic models differed by <0.01 for adult survival and <0.04 for maturation rate. Estimates of demographic metrics were sensitive to the accuracy of imputed vital rates. For example, mean error in generation time doubled in response to inaccurate estimates of maturation time. Accurate demographic data and metrics, such as generation time, are needed to inform conservation planning processes, for example through International Union for Conservation of Nature Red List assessments and population viability analysis. Imputed vital rates could be useful in this context but, as for any estimated model parameters, awareness of the sensitivities of demographic model outputs to the imputed vital rates is essential.     

Accuracy of Short‐Term Demographic Data in Projecting Long‐Term Fate of Populations

Short‐term surveys are useful in conservation of species if they can be used to reliably predict the long‐term fate of populations. However, statistical evaluations of reliability are rare. We studied how well short‐term demographic data (1999–2002) of tartar catchfly (Silene tatarica), a perennial riparian plant, projected the fate and growth of 23 populations of this species up to the year 2010. Surveyed populations occurred along a river with natural flood dynamics and along a regulated river. Riparian plant populations are affected by flooding, which maintains unvegetated shores, while forest succession proceeds in areas with little flooding. Flooding is less severe along the regulated river, and vegetation overgrowth reduces abundance of tartar catchfly on unvegetated shores. We built matrix models to calculate population growth rates and estimated times to population extinction in natural and in regulated rivers, 13 and 10 populations, respectively. Models predicted population survival well (model predictions matched observed survival in 91% of populations) and accurately predicted abundance increases and decreases in 65% of populations. The observed and projected population growth rates differed significantly in all but 3 populations. In most cases, the model overestimated population growth. Model predictions did not improve when data from more years were used (1999–2006). In the regulated river, the poorest model predictions occurred in areas where cover of other plant species changed the fastest. Although vegetation cover increased in most populations, it decreased in 4 populations along the natural river. Our results highlight the need to combine disturbance and succession dynamics in demographic models and the importance of habitat management for species survival along regulated rivers. Precisión de Datos Demográficos de Corto Plazo en la Proyección del Destino de Poblaciones a Largo Plazo     

Estimating the functional form for the density dependence from life history data

Two contrasting approaches to the analysis of population dynamics are currently popular: demographic approaches where the associations between demographic rates and statistics summarizing the population dynamics are identified; and time series approaches where the associations between population dynamics, population density, and environmental covariates are investigated. In this paper, we develop an approach to combine these methods and apply it to detailed data from Soay sheep (Ovis aries). We examine how density dependence and climate contribute to fluctuations in population size via age- and sex-specific demographic rates, and how fluctuations in demographic structure influence population dynamics. Density dependence contributes most, followed by climatic variation, age structure fluctuations and interactions between density and climate. We then simplify the density-dependent, stochastic, age-structured demographic model and derive a new phenomenological time series which captures the dynamics better than previously selected functions. The simple method we develop has potential to provide substantial insight into the relative contributions of population and individual-level processes to the dynamics of populations in stochastic environments.     

Joint modelling of breeding and survival in the kittiwake using frailty models

Assessment of population dynamics is central to population dynamics and conservation. In structured populations, matrix population models based on demographic data have been widely used to assess such dynamics. Although highlighted in several studies, the influence of heterogeneity among individuals in demographic parameters and of the possible correlation among these parameters has usually been ignored, mostly because of difficulties in estimating such individual-specific parameters. In the kittiwake (Rissa tridactyla), a long-lived seabird species, differences in survival and breeding probabilities among individual birds are well documented. Several approaches have been used in the animal ecology literature to establish the association between survival and breeding rates. However, most are based on observed heterogeneity between groups of individuals, an approach that seldom accounts for individual heterogeneity. Few attempts have been made to build models permitting estimation of the correlation between vital rates. For example, survival and breeding probability of individual birds were jointly modelled using logistic random effects models by [Cam, E., Link, W.A., Cooch, E.G., Monnat, J., Danchin, E., 2002. Individual covariation in life-history traits: seeing the trees despite the forest. Am. Naturalist, 159, in press]. This is the only example in wildlife animal populations we are aware of. Here we adopt the survival analysis approaches from epidemiology. We model the survival and the breeding probability jointly using a normally distributed random effect (frailty). Conditionally on this random effect, the survival time is modelled assuming a lognormal distribution, and breeding is modelled with a logistic model. Since the deaths are observed in year-intervals, we also take into account that the data are interval censored. The joint model is estimated using classic frequentist methods and also MCMC techniques in Winbugs. The association between survival and breeding attempt is quantified using the standard deviation of the random frailty parameters. We apply our joint model on a large data set of 862 birds, that was followed from 1984 to 1995 in Brittany (France). Survival is positively correlated with breeding indicating that birds with greater inclination to breed also had higher survival.     

Using Age Structure to Detect Impacts on Threatened Populations: a Case Study with Steller Sea Lions

Abstract:   A delayed response to change is often a characteristic of long-lived species and presents a major challenge to monitoring their status. However, rapid shifts in age structure can occur even while population size remains relatively static. We used time-varying matrix models to study age-structure information as a tool for improving detection of survivorship and fecundity change and status. We applied the methods to Steller sea lions (  Eumetopias jubatus ), a long-lived endangered marine mammal found throughout the North Pacific Rim. Population and newborn counts were supplemented with information on the fraction of the population that was juvenile, obtained by measuring animals in aerial photographs taken during range-wide censuses. By fitting the model to 1976–1998 data, we obtained maximum-likelihood estimates and 95% confidence intervals for juvenile survivorship, adult survivorship, and adult fecundity in the mid-1980s, late 1980s, and 1990s. We used a series of nested models to test whether the data were best fit by a model with one, two, or three temporal changes in demographic rates, and we fit the models to different lengths of data to test the number of years of data needed to detect a demographic change. The declines in the early 1980s were associated with severely low juvenile survivorship, whereas declines in the 1990s were associated with disproportionately low fecundity. We repeated these analyses, fitting only to the count data without the juvenile-fraction information, to determine whether the age-structure information changed the conclusions and/or changed the certainty and speed with which demographic-rate changes could be detected. The juvenile-fraction data substantially improved the degree to which estimates from the model were consistent with field data and significantly improved the speed and certainty with which changes in demographic rates were detected.     

Reliability of Absolute and Relative Predictions of Population Persistence Based on Time Series

Abstract:  Conventional population viability analysis (PVA) is often impractical because data are scarce for many threatened species. For this reason, simple count-based models are being advocated. The simplest of these models requires nothing more than a time series of abundance estimates, from which variance and autocorrelation in growth rate are estimated and predictions of population persistence are generated. What remains unclear, however, is how many years of data are needed to generate reliable estimates of these parameters and hence reliable predictions of persistence. By analyzing published and simulated time series, we show that several decades of data are needed. Predictions based on short time series were very unreliable mainly because limited data yielded biased, unreliable estimates of variance in growth rate, especially when growth rate was strongly autocorrelated. More optimistically, our results suggest that count-based PVA is sometimes useful for relative risk assessment (i.e., for ranking populations by extinction risk), even when time series are only a decade long. However, some conditions consistently lead to backward rankings. We explored the limited conditions under which simple count-based PVA may be useful for relative risk assessment.     

Relating chronic toxicity responses to population-level effects: A comparison of population-level parameters for three salmon species as a function of low-level toxicity

Standard laboratory toxicity tests assess the physiological responses of individual organisms to exposure to toxic substances under controlled conditions. Time and space restrictions often prevent the assessment of population-level responses to a toxic substance. Contaminants can affect various biological functions (e.g. growth, fecundity or behavior), which may alter different demographic traits, leading to population-level impacts. In this study, immune suppression, reproductive dysfunction and somatic growth impairment were examined using life history matrix models for coho salmon (Oncorhynchus kisutch), sockeye salmon (Oncorhynchus nerka) and chinook salmon (Oncorhynchus tshawytscha). Our intent was to gauge the relative magnitude of response to toxic effects among species and between life history stages, not provide a specific estimate of population growth rate or abundance. Effects due to immune suppression were modeled as reductions in age-specific survival. Toxic impacts on reproductive function were modeled as a 10% reduction in reproductive contribution for all reproductively mature age groups. Model runs that examined the effect of somatic growth reduction on population parameters incorporated both survival and reproductive impacts. All impacts were modeled as 10% reductions in the affected population demographic parameters. First-year survival and reproductive impacts produced similar population growth rates (λ), but resulted in different sensitivity and stable age distributions. Modeled somatic growth reduction produced additive effects on survival and reproduction. Toxic stressors producing similar changes in λ did not necessarily produce similar changes in the age distributions. Sensitivity and elasticity analyses demonstrated that changes to the first-year survival rate produced the greatest per-unit effect on λ for each species. Alteration in abundance of mature females also strongly influenced λ. Differences observed between species showed that the number of reproductive ages and time to reproductive maturity were important components for population-level responses. These results emphasize the importance of linking toxicity responses at low concentrations to the demographic traits they affect, and help to highlight the toxicity tests that are more suitable for assessing impacts on the focal species. Additionally, life history modeling is a useful tool for developing testable hypotheses regarding impacts on specific populations as well as for conducting comparisons between populations.     

Spatiotemporal variation in survival rates: implications for population dynamics of yellow-bellied marmots

Spatiotemporal variation in age-specific survival rates can profoundly influence population dynamics, but few studies of vertebrates have thoroughly investigated both spatial and temporal variability in age-specific survival rates. We used 28 years (1976-2003) of capture-mark-recapture (CMR) data from 17 locations to parameterize an age-structured Cormack-Jolly-Seber model, and investigated spatial and temporal variation in age-specific annual survival rates of yellow-bellied marmots (Marmota flaviventris). Survival rates varied both spatially and temporally, with survival of younger animals exhibiting the highest degree of variation. Juvenile survival rates varied from 0.52 +/- 0.05 to 0.78 +/- 0.10 among sites and from 0.15 +/- 0.14 to 0.89 +/- 0.06 over time. Adult survival rates varied from 0.62 +/- 0.09 to 0.80 +/- 0.03 among sites, but did not vary significantly over time. We used reverse-time CMR models to estimate the realized population growth rate (lamda), and to investigate the influence of the observed variation in age-specific survival rates on lamda. The realized growth rate of the population closely covaried with, and was significantly influenced by, spatiotemporal variation in juvenile survival rate. High variability in juvenile survival rates over space and time clearly influenced the dynamics of our study population and is also likely to be an important determinant of the spatiotemporal variation in the population dynamics of other mammals with similar life history characteristics.     

Demographic Forecasting in Koala Conservation

Abstract: The koala currently needs conservation intervention. There is clear evidence of decline in many populations, but the existence of other stable or expanding populations offers the possibility of a variety of creative solutions to their conservation problems. The 1998 National Koala Conservation Strategy emphasizes the need to obtain demographic information and to use this information to assess management options for koalas. We need accurate diagnoses of the status of koala populations and forecasts of their demographic future with and without particular management actions. In a qualitative fashion, this process has been undertaken many times on a local and national scale. Quantitative demographic forecasting tools are increasingly available, and koala management could benefit from their application both at the scale of individual populations and more broadly. There is already a considerable body of suitable data on the dispersal, effects of normal and catastrophic environmental variation on reproduction and survival, and on the effects of habitat change. Demographic forecasting, however, is hampered because the full suite of information is rarely available from a single population. In two Queensland populations, retrospective population viability analyses provided forecasts that were in agreement with observed population trends. Work is needed to determine whether data from one population can be applied to other populations. Models can then be developed to make projections at a multipopulation level on the basis of local population dynamics and dispersal. Certain koala populations, because of their long history of study, offer the opportunity to test demographic models retrospectively. These tests will not only aid in fine-tuning the models for koala biology and data but will also assist with the more general process of validating the models.     

Using demography and movement behavior to predict range expansion of the southern sea otter

In addition to forecasting population growth, basic demographic data combined with movement data provide a means for predicting rates of range expansion. Quantitative models of range expansion have rarely been applied to large vertebrates, although such tools could be useful for restoration and management of many threatened but recovering populations. Using the southern sea otter (Enhydra lutris nereis) as a case study, we utilized integro-difference equations in combination with a stage-structured projection matrix that incorporated spatial variation in dispersal and demography to make forecasts of population recovery and range recolonization. In addition to these basic predictions, we emphasize how to make these modeling predictions useful in a management context through the inclusion of parameter uncertainty and sensitivity analysis. Our models resulted in hind-cast (1989-2003) predictions of net population growth and range expansion that closely matched observed patterns. We next made projections of future range expansion and population growth, incorporating uncertainty in all model parameters, and explored the sensitivity of model predictions to variation in spatially explicit survival and dispersal rates. The predicted rate of southward range expansion (median = 5.2 km/yr) was sensitive to both dispersal and survival rates; elasticity analysis indicated that changes in adult survival would have the greatest potential effect on the rate of range expansion, while perturbation analysis showed that variation in subadult dispersal contributed most to variance in model predictions. Variation in survival and dispersal of females at the south end of the range contributed most of the variance in predicted southward range expansion. Our approach provides guidance for the acquisition of further data and a means of forecasting the consequence of specific management actions. Similar methods could aid in the management of other recovering populations.     

Reconstructing the historic demography of an endangered seabird

Reducing extinction risk for threatened species requires determining which demographic parameters are depressed and causing population declines. Museum collections may constitute a unique, underutilized resource for measuring demographic changes over long time periods using age-ratio analysis. We reconstruct the historic demography of a U.S. federally endangered seabird, the Marbled Murrelet (Brachyramphus marmoratus), from specimens collected approximately 100 years ago for comparison with predictions from comparative analyses and with results from contemporary field studies using both age-ratio analysis and conventional demographic estimators. Reproduction in the late 1800s and early 1900s matched predictions from comparative analysis, but was 8-9 times greater than contemporary estimates, whereas adult survival was unchanged. Historic reproductive rates would support stable populations, but contemporary levels should result in population declines. Contemporary demographic estimates derived from age-ratio analysis were similar to estimates from conventional estimators. Using museum specimens to reconstruct historic demography provides a unique approach to identify causes of decline and to set demographic benchmarks for recovery of endangered species that meet most assumptions of age-ratio analysis.     

An integrated approach for predicting fates of reintroductions with demographic data from multiple populations

We devised a novel approach to model reintroduced populations whereby demographic data collected from multiple sites are integrated into a Bayesian hierarchical model. Integrating data from multiple reintroductions allows more precise population-growth projections to be made, especially for populations for which data are sparse, and allows projections that account for random site-to-site variation to be made before new reintroductions are attempted. We used data from reintroductions of the North Island Robin (Petroica longipes), an endemic New Zealand passerine, to 10 sites where non-native mammalian predators are controlled. A comparison of candidate models that we based on deviance information criterion showed that rat-tracking rate (an index of rat density) was a useful predictor of robin fecundity and adult female survival, that landscape connectivity and a binary measure of whether sites were on a peninsula were useful predictors of apparent juvenile survival (probably due to differential dispersal away from reintroduction sites), and that there was unexplained random variation among sites in all demographic rates. We used the two best supported models to estimate the finite rate of increase (λ) for populations at each of the 10 sites, and for a proposed reintroduction site, under different levels of rat control. Only three of the reintroduction sites had λ distributions completely >1 for either model. At two sites, λ was expected to be >1 if rat-tracking rates were <5%. At the other five reintroduction sites, λ was predicted to be close to 1, and it was unclear whether growth was expected. Predictions of λ for the proposed reintroduction site were less precise than for other sites because distributions incorporated the full range of site-to-site random variation in vital rates. Our methods can be applied to any species for which postrelease data on demographic rates are available and potentially can be extended to model multiple species simultaneously.     

Effect of stage-specific vital rates on population growth rates and effective population sizes in an endangered iteroparous plant

Effective population size (N(e)) determines the strength of genetic drift and can influence the level of genetic diversity a population can maintain. Assessing how changes in demographic rates associated with environmental variables and management actions affect N(e) thus can be crucial to the conservation of endangered species. Calculation of N(e) through demographic models makes it possible to use elasticity analyses to study this issue. The elasticity of N(e) to a given vital rate is the proportional change in N(e) associated with a proportional increase in that vital rate. In addition, demographic models can be used to study N(e) and population growth rate (λ) simultaneously. Simultaneous examination is important because some vital rates differ diametrically in their associations with λ and N(e). For example, in some cases increasing these vital rates increases λ and decreases N(e). We used elasticity analysis to study the effect of stage-specific survival and flowering rates on N(e), annual effective population size (N(a)), and λ in seven populations of the endangered plant Austrian dragonhead (Dracocephalum austriacum). In populations with λ ≥ 1, the elasticities of N(e) and N(a) were similar to those of λ. Survival rates of adults were associated with greater elasticities than survival rates of juveniles, flowering rates, or fecundity. In populations with λ < 1, N(e) and N(a) exhibited greater elasticities to juvenile than to adult vital rates. These patterns are similar to those observed in other species with similar life histories. We did not observe contrasting effects of any vital rate on λ and N(e); thus, management actions that increase the λ of populations of Austrian dragonhead will not increase genetic drift. Our results show that elasticity analyses of N(e) and N(a) can complement elasticity analysis of λ. Moreover, such analyses do not require more data than standard matrix models of population dynamics.     

Combining demographic and count-based approaches to identify source-sink dynamics of a threatened seabird.

Identifying source-sink dynamics is of fundamental importance for conservation but is often limited by an inability to determine how immigration and emigration influence population processes. We demonstrate two ways to assess the role of immigration on population processes without directly observing individuals dispersing from one population to another and apply these methods to a population of Marbled Murrelets (Brachyramphus marmoratus) in California (USA). In the first method, the rate of immigration (i) is estimated by subtracting local recruitment (recruitment from within the population due to reproduction) estimated with demographic data from total recruitment (f; recruitment from within the population plus recruitment from other populations) estimated using temporal symmetry mark-recapture models developed by R. Pradel. The second method compares population growth rates estimated with temporal symmetry models (lambdaTS) and/or population growth rates estimated from counts of individuals over multiple sampling periods (lambdaC) with growth estimates from a stage-structured projection matrix model (lambdaM). Both lambdaTS and lambdaC incorporate all demographic processes affecting population change (birth, death, immigration, and emigration), whereas matrix models are usually constructed without incorporating immigration. Thus, if lambdaTS and lambdaC are > or = 1 and lambdaM < 1, the population is sustained by immigration and is considered to be a sink. Using the first method, recruitment estimated with temporal symmetry models was high (f= 0.182, SE = 0.058), the mean adult birth rate, as estimated using the ratio of juveniles to > or = 1 year old individuals (observed during ship-based surveys) was low (bA = 0.039, SE = 0.014), and immigration was 0.160 (SE = 0.057). Using the second method, murrelet numbers in central California were stable (lambdaC = 1.058, SE = 0.047; lambdaTS = 1.064, SE = 0.033), but were projected to decline 9.5% annually in the absence of immigration (lambdaM = 0.905, SE = 0.053). Our results suggest that Marbled Murrelets in central California represent a sink population that is stable but would decline in the absence of immigration from larger populations to the north. However, the extent to which modeled immigration is due to permanent recruitment or temporarily dispersing individuals that simply mask population declines is uncertain.     

Effective Population Size in Red-Cockaded Woodpeckers: Population and Model Differences

Loss of genetic variability in isolated populations is an important issue for conservation biology. Most studies involve only a single population of a given species and a single method of estimating rate of loss. Here we present analyses for three different Red-cockaded Woodpecker ( Picoides borealis ) populations from different geographic regions. We compare two different models for estimating the expected rate of loss of genetic variability, and test their sensitivity to model parameters. We found that the simpler model (Reed et al. 1988) consistently estimated a greater rate of loss of genetic variability from a population than did the Emigh and Pollak (1979) model. The ratio of effective population size (which describes the expected rate of loss of genetic variability) to breeder population size varied widely among Red-cockaded Woodpecker populations due to geographic variation in demography. For this species, estimates of effective size were extremely sensitive to survival parameters, but not to the probability of breeding or reproductive success. Sensitivity was sufficient that error in estimating survival rates in the field could easily mask true population differences in effective size. Our results indicate that accurate and precise demographic data are prerequisites to determining effective population size for this species using genetic models, and that a single estimate of rate of loss of genetic variability is not valid across populations.     

Demography of central and marginal populations of monkeyflowers (Mimulus cardinalis and M. lewisii)

Every species occupies a limited geographic area, but how spatiotemporal environmental variation affects individual and population fitness to create range limits is not well understood. Because range boundaries arise where, on average, populations are more likely to go extinct than to persist, range limits are an inherently population-level problem for which a demographic framework is useful. In this study, I compare demographic parameters and population dynamics between central and marginal populations of monkeyflowers, Mimulus cardinalis and M. lewisii, along an elevation gradient spanning both species' ranges. Central and marginal populations of both species differed in survival and fecundity. For M. lewisii, these components of fitness were higher in central than in marginal populations, but for M. cardinalis the converse was true. To assess spatiotemporal variation in population dynamics, I used transition matrix models to estimate asymptotic population growth rates (lambda) and found that population growth rates of M. lewisii were highest at the range center and reduced at the range margin. Population growth rates of M. cardinalis were highest at the range margin and greatly reduced at the range center. Life table response analysis decomposed spatiotemporal variation in lambda into contributions from each transition between life stages, finding that transitions from large nonreproductive and reproductive plants to the seed class and stasis in the reproductive class made the largest contributions to spatial differences in lambda. These transitions had only low to moderate sensitivities, indicating that differences in projected population growth rates resulted mainly from observed differences in transition matrix parameters and their underlying vital rates.     

Life history and viability of a long-lived marine invertebrate: the octocoral Paramuricea clavata

The red gorgonian Paramuricea clavata is a long-lived, slow-growing sessile invertebrate of ecological and conservation importance in the northwestern Mediterranean Sea. We develop a series of size-based matrix models for two Paramuricea clavata populations. These models were used to estimate basic life history traits for this species and to evaluate the viability of the red gorgonian populations we studied. As for many other slow-growing species, sensitivity and elasticity analysis demonstrate that gorgonian population growth is far more sensitive to changes in survival rates than to growth, shrinkage, or reproductive rates. The slow growth and low mortality of red gorgonians results in low damping ratios, indicating slow convergence to stable size structures (at least 50 years). The stable distributions predicted by the model did not differ from the observed ones. However, our simulations point out the fragility of this species, showing both populations in decline and high risk of extinction over moderate time horizons. These declines appear to be related to a recent increase in anthropogenic disturbances. Relative to their life span, the values of recruitment elasticity for Paramuricea clavata are lower than those reported for other marine organisms but are similar to those reported for some long-lived plants. These values and the delayed age of sexual maturity, in combination with the longevity of the species, show a clear fecundity/mortality trade-off. Full demographic studies of sessile marine species are quite scarce but can provide insight into population dynamics and life history patterns for these difficult and under-studied species. While our work shows clear results for the red gorgonian, the variability in some of our estimates suggest that future work should include data collection over longer temporal and spatial scales to better understand the long-term effects of natural and anthropogenic disturbances on red gorgonian populations.     

Evaluating mortality rates with a novel integrated framework for nonmonogamous species

The conservation of wildlife requires management based on quantitative evidence, and especially for large carnivores, unraveling cause‐specific mortalities and understanding their impact on population dynamics is crucial. Acquiring this knowledge is challenging because it is difficult to obtain robust long‐term data sets on endangered populations and, usually, data are collected through diverse sampling strategies. Integrated population models (IPMs) offer a way to integrate data generated through different processes. However, IPMs are female‐based models that cannot account for mate availability, and this feature limits their applicability to monogamous species only. We extended classical IPMs to a two‐sex framework that allows investigation of population dynamics and quantification of cause‐specific mortality rates in nonmonogamous species. We illustrated our approach by simultaneously modeling different types of data from a reintroduced, unhunted brown bear (Ursus arctos) population living in an area with a dense human population. In a population mainly driven by adult survival, we estimated that on average 11% of cubs and 61% of adults died from human‐related causes. Although the population is currently not at risk, adult survival and thus population dynamics are driven by anthropogenic mortality. Given the recent increase of human‐bear conflicts in the area, removal of individuals for management purposes and through poaching may increase, reversing the positive population growth rate. Our approach can be generalized to other species affected by cause‐specific mortality and will be useful to inform conservation decisions for other nonmonogamous species, such as most large carnivores, for which data are scarce and diverse and thus data integration is highly desirable.     

Age-structured modeling reveals long-term declines in the natality of western Steller sea lions.

Since the mid-1970s, the western Steller sea lion (Eumetopias jubatus), inhabiting Alaskan waters from Prince William Sound west through the Aleutian Islands, has declined by over 80%. Changing oceanographic conditions, competition from fishing operations, direct human-related mortality, and predators have been suggested as factors driving the decline, but the indirect and interactive nature of their effects on sea lions have made it difficult to attribute changes in abundance to specific factors. In part, this is because only changes in abundance, not changes in vital rates, are known. To determine how vital rates of the western Steller sea lion have changed during its 28-year decline, we first estimated the changes in Steller sea lion age structure using measurements of animals in aerial photographs taken during population surveys since 1985 in the central Gulf of Alaska (CGOA). We then fit an age-structured model with temporally varying vital rates to the age-structure data and to total population and pup counts. The model fits indicate that birth rate in the CGOA steadily declined from 1976 to 2004. Over the same period, survivorship first dropped severely in the early 1980s, when the population collapsed, and then survivorship steadily recovered. The best-fitting model indicates that in 2004, the birth rate in the central Gulf of Alaska was 36% lower than in the 1970s, while adult and juvenile survivorship were close to or slightly above 1970s levels. These predictions and other model predictions concerning population structure match independent field data from mark-recapture studies and photometric analyses. The dominant eigenvalue for the estimated 2004 Leslie matrix is 1.0014, indicating a stable population. The stability, however, depends on very high adult survival, and the shift in vital rates results in a population that is more sensitive to changes in adult survivorship. Although our modeling analysis focused exclusively on the central Gulf of Alaska, the western Gulf of Alaska and eastern Aleutians show a similar pattern of declining pup fraction with no increase in the juvenile, or pre-breeding, fraction. This suggests that declining birth rate may be a problem for western Steller sea lions across the Gulf of Alaska and into the Aleutian Islands.