The founding of a southern elephant seal colony

The only large mainland colony of southern elephant seals (Mirounga leonina) is on Península Valdés, at 42°S, in Argentine Patagonia. Censuses of pups have been carried out regularly there since 1970, and the population grew five‐fold by 2010. Here we use Bayesian modeling tools to make rigorous estimates of the rate of population growth, r, and to estimate survival and recruitment parameters that could account for the growth, incorporating observation error across different census methods. In the 1970s, r= 8%/yr, but has slowed to <1%/yr over the past decade. Using explicit demographic models, we established that the high growth of the 1970s was consistent with adult and juvenile survival at the upper end of published values (0.87/yr adult female survival; 0.40 juvenile survivorship to age four); the decline in the rate of population growth from 1970 to 2010 can be described by density‐dependent reductions in adult and juvenile survival that fall well within published variation. Extrapolating empirical models of population growth rate backwards illustrates that the population could have been an established colony, with 100 pups born per year, between 1915 and 1945, consistent with qualitative observations prior to 1950. We conclude that the Valdés colony was founded by a few immigrants early in the 20th century and has been growing mostly by internal recruitment, with unknown density‐dependent processes causing a reduction in growth and stabilization at 15,000–16,000 pups born.     

The demography of a newly established Osprey Pandion haliaetus population in France

Using long‐term mark–resighting data acquired over 27 years in continental France, we estimated demographic parameters and modelled the dynamics of a newly established population of Ospreys Pandion haliaetus using a life‐history model. We then performed prospective and retrospective analyses to estimate the sensitivity of the population growth rate to demographic parameters, and to quantify their contribution to the observed variation in abundance. The observed population growth rate was estimated at 1.150 (from one to 38 pairs in the period 1985–2011), and the stochastic population growth rate was estimated at 1.156. The number of fledglings per nest made the largest contribution to the variance of the observed population growth rate. Breeding productivity was stable across years. In contrast, the prospective analysis indicated that the sensitivity of the population growth rate was greatest for immigration and adult survival. Our results suggest that the increase of a new and recently established breeding population of Ospreys was mainly driven by local dynamics (high productivity and high proportion of breeding individuals), with no sign of density‐dependence except for juvenile survival. This probably reflects highly favourable conditions for breeding. Our results show that productivity can be a major driver in recovering raptor populations, and conservation work should aim to protect occupied nest‐sites and their surrounding habitat and to maintain highly favourable foraging areas in the vicinity of breeding sites.     

Viability analysis of a threatened amphibian population: modelling the past,present and future

Conservation actions that have maintained populations in the past may not necessarily do so in the future. Population viability analysis provides one tool for exploring the impact of management actions on large temporal scales. However, there are relatively few long‐term data sets that provide the demographic and environmental data demanded by such models. Using a 37‐yr data set, we used RAMAS Metapop to model the persistence of natterjack toads Epidalea (Bufo) calamita on a heathland in southern Britain. A retrospective analysis showed that the best fit between the predicted population trajectories and the real population was when the management carried out was modelled as an increase in K of 150 toads yr^?1. However, even if ongoing management continues to improve K by a further 40–60 toads yr^?1 over the next 50 yr, the population still has an extinction risk of at least 60% if other factors remain unchanged. Sensitivity analyses and simulated management scenarios indicated that the population was most sensitive to changes in the survival of juvenile (i.e. 1–2 yr old) toads. In addition, if the frequency and severity of pond desiccation increases, the risk of extinction was predicted to increase as a result of reduced recruitment. Low levels of extinction risk occurred irrespective of K when juvenile survival was enhanced in combination with low frequency and severity of pond desiccation. The models suggest that populations that are responding to management against a background of natural fluctuations may remain vulnerable to extinction for several decades. These extinction risks may increase if habitat management fails to offset reductions in recruitment and juvenile survival caused by environmental change.     

Rapid recovery of Dutch gray seal colonies fueled by immigration

Gray seals were first observed breeding in the Dutch Wadden Sea in 1985, after centuries of absence. The breeding colony there is now the largest on the European continent. We describe the changes in gray seal numbers and their geographical expansion, and estimate how these processes were influenced by immigration from other colonies. Counts of hauled out animals were carried out between 1985 and 2013, monitoring three different periods of the seals' annual cycle. Using priors determined for the UK population, a Bayesian demographic model was fitted to pup numbers to estimate the population parameters driving the growth. This included immigration of subadults into the breeding population, which contributed to an average growth rate in the pup counts of 19%/yr, much higher than expected in a closed population. This immigration may account for approximately 35% of the total annual growth. In addition, at least 200 gray seals from the UK visit the area temporarily. Recovery of the population in the Netherlands occurred more than 50 yr after gray seals were protected in the UK. These time scales should be taken into account when studying long living marine mammals, e.g., in impact and conservation studies.     

A combined inverse method and multivariate approach for exploring population trends of Florida manatees

Given the expense and time required to monitor marine mammal populations effectively, approaches that fully exploit the resulting data certainly are warranted. We employed a two‐step modeling approach to estimate key demographic parameters, including immigration, from aerial surveys of manatees (Trichechus manatus latirostris) in the Northwest management unit of Florida. Abundances of adults and calves were predicted by multivariate adaptive regression spline (MARS) models, after accounting for heterogeneous detection rates caused by variable environmental conditions. The resulting predictions were incorporated into a stage‐structured, deterministic model that used an inverse method to estimate parameters with and without immigration. The model without immigration estimated mean survival probabilities of 0.966, 0.923, and 0.794 for adults, subadults and calves, respectively, with a per capita reproductive rate of 0.135. These parameter estimates yielded an overall mean population growth rate of approximately 1.037, which is comparable to rates from mark‐recapture studies. When we added an immigration term that accounted for the greater slope in adult counts since 1999, as identified by the MARS model, the estimated per capita reproductive rate was 0.122, with survival probabilities for adults, subadults and calves of 0.926, 0.920, and 0.833, respectively. These rates were coupled with an estimated mean winter immigration rate corresponding to roughly 5.2% of the adult and subadult population. In this latter scenario, the number of manatees in the core population of the Northwest management unit was predicted to remain constant, with a population growth rate near one, and additional manatees counted during aerial surveys were deemed to be immigrants. While further studies could certainly expound on the potential effects of migrants on population indices, we present this first published immigration estimate for wintering manatees in northwest Florida.     

The demographic drivers of local population dynamics in two rare migratory birds

The exchange of individuals among populations can have strong effects on the dynamics and persistence of a given population. Yet, estimation of immigration rates remains one of the greatest challenges for animal demographers. Little empirical knowledge exists about the effects of immigration on population dynamics. New integrated population models fitted using Bayesian methods enable simultaneous estimation of fecundity, survival and immigration, as well as the growth rate of a population of interest. We applied this novel analytical framework to the demography of two populations of long-distance migratory birds, hoopoe Upupa epops and wryneck Jynx torquilla, in a study area in south-western Switzerland. During 2002–2010, the hoopoe population increased annually by 11%, while the wryneck population remained fairly stable. Apparent juvenile and adult survival probability was nearly identical in both species, but fecundity and immigration were slightly higher in the hoopoe. Hoopoe population growth rate was strongly correlated with juvenile survival, fecundity and immigration, while that of wrynecks strongly correlated only with immigration. This indicates that demographic components impacting the arrival of new individuals into the populations were more important for their dynamics than demographic components affecting the loss of individuals. The finding that immigration plays a crucial role in the population growth rates of these two rare species emphasizes the need for a broad rather than local perspective for population studies, and the development of wide-scale conservation actions.     

Estimate of population extinction risk and its application to ecological risk management

Environmental threats, such as habitat size reduction or environmental pollution, may not cause immediate extinction of a population but may shorten the expected time to extinction. We developed a method to estimate the mean time to extinction for a density-dependent population with environmental fluctuation and to compare the impacts of different risk factors. We first derived a formula of the mean extinction time for a population with logistic growth and environmental and demographic stochasticities expressed as a stochastic differential equation model (canonical model). The relative importance of different risk factors is evaluated by the decrease in the mean extinction time. We studied an approximated formula for the reduction in habitat size that enhances extinction risk by the same magnitude as a given decrease in survivorship caused by toxic chemical exposure. In a large population (large K) or in a slowly growing population (small r), a small decrease in survivorship can cause the extinction risk to increase, corresponding to a significant reduction in the habitat size. Finally, we studied an approximate maximum likelihood estimate of three parameters (intrinsic growth rate r, carrying capacity K, and environmental stochasticity σ ² _e ) from time series data. By Monte Carlo sampling, we can remove the bias very effectively and determine the confidence interval. We discuss here how the reliability of the estimate changes with the length of time series. If we know the intrinsic rate of population growth r, the mean extinction time is estimated quite accurately even when only a short time series is available for parameter estimation. Received: March 31, 1999 / Accepted: November 9, 1999     

Demography of Agile Gibbons (<Emphasis Type="Italic">Hylobates agilis</Emphasis>) in a Lowland Tropical Rain Forest of Southern Sumatra,Indonesia: Problems in Paradise

Gibbons are among the best-studied Asian primates, but few studies address their demography and life history strategies. We used annual censuses to study the demography of agile gibbons (Hylobates agilis) between 1998 and 2009 in rain forests of Bukit Barisan Selatan National Park, Indonesia. The population declined from 22 individuals (9 groups) to 14 individuals (5 groups) over the 12 yr of study. Infant survival to the juvenile age class was 33.3%, and 16.7% of infants survived to the subadult age class. The interbirth interval was 3.83 ± 1.15 yr and birth rate was 0.22–0.28 infants female^–1 yr^–1. Two groups colonized the study area but subsequently disappeared. We documented 7 immigrations, 17 disappearances, and ≥10 transients in the population. Compared to lar gibbons (Hylobates lar) and Bornean white-bearded gibbons (Hylobates albibarbis), Way Canguk’s agile gibbon population is characterized by slow reproduction, low survival, and high group turnover. We hypothesize that, although the habitat is high in fruit resources, agile gibbons may be displaced or excluded from the best fruit resources by larger and more numerous competitors, incurring costs of decreased opportunities to forage and increased travel, and leading to higher mortality for young agile gibbons. The reproductive potential of this agile gibbon population is insufficient to compensate for high mortality, and the population is unlikely to persist without immigration from outside the area. Given the agile gibbons’ endangered status and limited capacity to respond demographically to change, it is likely that intensive management interventions will be required to conserve this species.     

Demography and dynamics of three wild horse populations in the Australian Alps

Wild horses (Equus caballus) are a non‐native species occupying over 2800 km² of the nationally significant Australian Alps National Parks. We estimated key demographic parameters (fecundity, adult and juvenile survival and annual finite population growth rate) over 3 years and related these to horse body condition and available food for three populations under natural conditions, and found a trend consistent with food limitation. The populations were independent, with different site characteristics and occupied areas, identified by land managers, as areas of concern about possible conservation impacts. Annual fecundity and juvenile survival varied across sites averaging between 0.21 and 0.31 female young per adult female, and 0.83 and 0.90 per annum, respectively, and annual adult survival was consistent across sites averaging 0.91 per annum. One population was increasing (λ = 1.09 year^?1; 95% CI 1.04–1.14) and two populations were stable (λ ～ 1.0 year^?1). Mean body condition of horses was positively correlated with mean pasture biomass rank. Across the three populations, fecundity, recruitment, body condition and annual finite population growth rate were lowest when mean pasture biomass rank was lowest and conversely highest when pasture rank was highest. We conclude that food limitation appears to be operating across these three sites. We used our results to assess the sensitivity of annual finite rate of increase (λ) to changes in key demographic parameters and found that λ was most sensitive to a change in adult survival, with the second most sensitive parameter being fecundity. Thus, if the aim of management is to reduce the size of the wild horse population then targeting adult survival is most important, followed by fecundity. Finally, we estimated the linear, negative, numerical response for wild horses between annual λ and horses per unit pasture biomass.     

Inferring causal factors for a declining population of bottlenose dolphins via temporal symmetry capture–recapture modeling

We applied temporal symmetry capture–recapture (TSCR) models to assess the strength of evidence for factors potentially responsible for population decline in bottlenose dolphins (Tursiops truncatus) in Doubtful Sound, New Zealand from 1995 to 2008. Model selection was conducted to estimate recruitment and population growth rates. There were similar levels of support for three different models, each reflecting distinct trends in recruitment. Modeling yielded low overall estimates of recruitment (0.0249, 95% CI: 0.0174–0.0324) and population growth rate (0.9642, 95% CI: 0.9546–0.9737). The TSCR rate of population decline was consistent with an estimate derived from trends in abundance (lambda = 0.9632, 95% CI: 0.9599–0.9665). The TSCR model selection confirmed the influence of a decline in the survival of calves (<1 yr old) since 2002 for population trends. However, TSCR population growth rates did not exceed 1 in any year between 1995 and 2008, indicating the population was declining prior to 2002. A separate reduction in juvenile survival (1–3 yr old) prior to 2002 was identified as a likely contributing factor in the population decline. Thus, TSCR modeling indicated the potential cause of the population decline in Doubtful Sound: cumulative impacts on individuals <3 yr old resulting in a reduced recruitment.     

Survival and abundance of short‐finned pilot whales in the archipelago of Madeira,NE Atlantic

Estimates of population parameters for the short‐finned pilot whale, Globicephala macrorhynchus, are scarce in literature, contributing to an International Union for Conservation of Nature (IUCN) status of Data Deficient. In this study, photo‐identification data collected over 7 yr from Madeira were used to estimate for the first time survivorship, capture probability, and abundance in this species using mark‐recapture methodology. The Cormack‐Jolly‐Seber model estimated that the adult island‐associated (i.e., resident and regular visitor) whales had a constant survival rate of 0.960 (95% CI: 0.853–0.990) and an annual capture probability varying between 0.372 (CI: 0.178–0.619) and 0.843 (CI: 0.619–0.947). A parameterization of the Jolly‐Seber model estimated that 140 island‐associated whales (CI: 131–151) used the area throughout the course of the study. Based on a closed population model, the most precise (lower CV) annual estimate of the total number of pilot whales using the southern and eastern waters of Madeira (~900 km²) in a 3 mo period covering summer/autumn was 334 animals (CI: 260–437). No trend was observed. Despite including biases, the approach used in this study provided plausible estimates of population parameters, which can contribute to the regional conservation strategies.     

Contemporary effective population size and predicted maintenance of genetic diversity in the endangered kea (Nestor notabilis)

Population size and the potential for maintenance of genetic diversity are critical information for the monitoring of species of conservation concern. However, direct estimates of population size are not always feasible, making indirect genetic approaches a valuable alternative. We estimated contemporary effective population size (N_e) in the endangered kea (Nestor notabilis) using three different methods. We then inferred the census size (N_C) using published N_e/N_C ratios and modelled the future maintenance of genetic diversity assuming a number of demographic parameters. Short-term N_e was small with a range-wide N_e?N_C was within the range of the current estimate (c. 1000–5000). Forward simulations showed low probability of retaining 90% of rare alleles without immigration. However, the probability of maintaining genetic diversity was high with immigration, juvenile survival of?≥?30%, and an initial sex ratio of c. 0.5–0.6. Despite the low N_e in kea, predator control and/or artificial immigration might be sufficient to maintain the present genetic diversity.     

Growth and survival of peritrich ciliates in an urban stream

Summary A computerized system was devised to trace the attachment, growth and disappearance of peritrich ciliates in an urban river. By tracing the development of each Carchesium polypinum colony in the river, it was possible to estimate the actual growth rate of this species without considering the effect of immigration. The survival and colonization rates of the colonies and individuals of solitary species could also be estimated. C. polypinum showed high growth rates (r=1.370 day^-1, doubling time 12.14 h). The number of daily colonizers also increased at a high rate, and the combination of growth and colonization caused very high population increase rates at the area level. Because of its low survivorship, the rate of Vorticella microstoma increase was much lower than previously reported. The survival of solitary species and small colonies was lower than with large colonies, and fewer new colonizers survived with the development of the attached microbial community.     

Estimating Survival for Elusive Juvenile Pond-Breeding Salamanders

Juvenile vital rates have important effects on population dynamics for many species, but this demographic is often difficult to locate and track. As such, we frequently lack reliable estimates of juvenile survival, which are necessary for accurately assessing population stability and potential management approaches to conserve biodiversity. We estimated survival rates for elusive juveniles of 3 species, the ringed salamander (Ambystoma annulatum), spotted salamander (A. maculatum), and small-mouthed salamander (A. texanum), using 2 approaches. First, we conducted an 11-month (2016–2017) mark-recapture study within semi-natural enclosures and used Bayesian Cormack-Jolly-Seber models to estimate survival and recapture probabilities. Second, we inferred the expected annual juvenile survival rate given published vital rates for pre-metamorphic and adult ambystomatids assuming stable population growth. For all 3 species, juvenile survival probabilities were constant across recapture occasions, whereas recapture probability estimates were time-dependent. Further, survival and recapture probabilities among study species were similar. Post-study sampling revealed that the initial study period median estimate of annual survival probability (0.39) underestimated the number of salamanders known alive at 11 months. We therefore appended approximately 1 year of opportunistic data, which produced a median annual survival probability of 0.50, encompassing salamanders that we knew to have been alive. Calculation from literature values suggested a mean annual terrestrial juvenile ambystomatid survival probability of 0.49. Similar results among our approaches indicated that juvenile survival estimates for the study species were robust and likely comparable to rates in nature. These estimates can now be confidently applied to research, monitoring, and management efforts for the study species and ecologically similar taxa. Our findings indicated that similarly robust vital rate estimates for subsets of ecologically and phylogenetically similar species can provide reasonable surrogate demographic information that can be used to reveal key factors influencing population viability for data-deficient species. © 2020 The Wildlife Society.     

Population size and survival of the Brazilian Torrent Frog Hylodes heyeri (Anura,Hylodidae)

Population dynamics are influenced by environmental variability and understanding the abundance and persistence of individuals and populations is a fundamental goal of population ecology. Thus, estimating demographic parameters to identify the factors important for population variability is required to understand temporal and spatial dynamics. The stream-living diurnal frog Hylodes heyeri is endemic to the Atlantic Forest of Brazil in the states of Paraná, São Paulo and Santa Catarina. Here we use capture-mark-recapture methods to estimate survival rates and population size of this Brazilian Torrent Frog in Pico do Marumbi State Park, Paraná. We used CJS models for an open population to estimate apparent survival, capturability and population size in two streams. The number of captures during each session was positively correlated with the minimum weekly temperature. Despite that correlation, the most parsimonious model of survival and capturability was the constant model for both parameters, resulting in a monthly survival rate of 0.38 (95% CI = 0.30–0.46). Thus, only the abundance of the frog differed in the two streams (79 vs. 36), with the population size estimate of 187 individuals. Reproduction is seasonal in the Brazilian Torrent Frog and so the low monthly survival rate suggests that animals move over time rather than die, because 38% month⁻¹ survival should result in <1% of the population remaining after 5 months. Thus, researchers must recognize that populations are labile and individuals often move or are washed downstream during heavy rainfall, generating apparently rapid local turnover that is unlikely to reflect true mortality.     

Survival of Adult Martens in Northern Wisconsin

Abstract: Low adult marten (Martes americana) survival may be one factor limiting their population growth >30 yr after their reintroduction in Wisconsin, USA. We estimated annual adult marten survival at 0.81 in northern Wisconsin, with lower survival during winter (0.87) than summer-fall (1.00). Fisher (Martes pennanti) and raptor kills were infrequent, and each reduced marten adult annual survival <10%. Annual adult survival was similar to or higher than survival in other areas, suggesting that it was not unusually low and therefore did not limit recovery of marten populations in northern Wisconsin. We captured few juvenile martens, suggesting low reproduction or reduced juvenile survival.     

Source populations in carnivore management: cougar demography and emigration in a lightly hunted population

Wildlife agencies typically attempt to manage carnivore numbers in localized game management units through hunting, and do not always consider the potential influences of immigration and emigration on the outcome of those hunting practices. However, such a closed population structure may not be an appropriate model for management of carnivore populations where immigration and emigration are important population parameters. The closed population hypothesis predicts that high hunting mortality will reduce numbers and densities of carnivores and that low hunting mortality will increase numbers and densities. By contrast, the open population hypothesis predicts that high hunting mortality may not reduce carnivore densities because of compensatory immigration, and low hunting mortality may not result in more carnivores because of compensatory emigration. Previous research supported the open population hypothesis with high immigration rates in a heavily hunted (hunting mortality rate=0.24) cougar population in northern Washington. We test the open population hypothesis and high emigration rates in a lightly hunted (hunting mortality rate=0.11) cougar population in central Washington by monitoring demography from 2002 to 2007. We used a dual sex survival/fecundity Leslie matrix to estimate closed population growth and annual census counts to estimate open population growth. The observed open population growth rate of 0.98 was lower than the closed survival/fecundity growth rates of 1.13 (deterministic) and 1.10 (stochastic), and suggests a 12–15% annual emigration rate. Our data support the open population hypothesis for lightly hunted populations of carnivores. Low hunting mortality did not result in increased numbers and densities of cougars, as commonly believed because of compensatory emigration.     

Quantifying the contribution of immigration to population dynamics: a review of methods,evidence and perspectives in birds and mammals

The demography of a population is often reduced to the apparent (or local) survival of individuals and their realised fecundity within a study area defined according to logistical constraints rather than landscape features. Such demographics are then used to infer whether a local population contributes positively to population dynamics across a wider landscape context. Such a simplistic approach ignores a fundamental process underpinning population dynamics: dispersal. Indeed, it has long been accepted that immigration contributed by dispersers that emigrated from neighbouring populations may strongly influence the net growth of a local population. To date however, we lack a clear picture of how widely immigration rate varies both among and within populations, in relation to extrinsic and intrinsic ecological conditions, even for the best‐studied avian and mammalian populations. This empirical knowledge gap precludes the emergence of a sound conceptual framework that ought to inform conservation and population ecology. This review, conducted on both birds and mammals, has thus three complementary objectives. First, we describe and evaluate the relative merits of methods used to quantify immigration and how they relate to widely applicable metrics. We identify two simple and unifying metrics to measure immigration: the immigration rate i_t defined as the ratio of the number of immigrants present in the population at time t + 1 and the total breeding population in year t, and π_t , the proportion of immigrants among new recruits (i.e. new breeders). Two recently developed methods are likely to provide the most valuable data on immigration in the near future: individual parentage (rather than population) assignments based on genetic sampling, and spatially explicit integrated population models combining multiple sources of demographic data (survival, fecundity and population counts). Second, we report on a systematic literature review of studies providing a quantitative measure of immigration. Although the diversity of methods employed precludes detailed analyses, it appears that the number of immigrants exceeds locally born individuals in recruitment for most avian populations (median π_t  = 0.57, N = 45 estimates from 37 studies), a figure twofold higher than estimated for mammalian populations (median π_t  = 0.26, N = 33 estimates from 11 studies). Third, recent quantitative studies reveal that immigration can be the main driver of temporal variation in population growth rates, across a wide array of demographic and spatial contexts. To what extent immigration acts as a regulatory process has however been considered only rarely to date and deserves more attention. Overall, it is likely that most populations benefit from immigrants without necessarily being sink populations. Furthermore, we suggest that quantitative estimates of immigration should be core to future demographic studies and plead for more empirical evidence about the ways in which immigration interacts with local demographic processes to shape population dynamics. Finally, we discuss how to tackle spatial population dynamics by exploring, beyond the classical source–sink framework, the extent to which populations exchange individuals according to spatial scale and type of population distribution throughout the landscape.     

Inferring the temperature dependence of population parameters: the effects of experimental design and inference algorithm

Understanding and quantifying the temperature dependence of population parameters, such as intrinsic growth rate and carrying capacity, is critical for predicting the ecological responses to environmental change. Many studies provide empirical estimates of such temperature dependencies, but a thorough investigation of the methods used to infer them has not been performed yet. We created artificial population time series using a stochastic logistic model parameterized with the Arrhenius equation, so that activation energy drives the temperature dependence of population parameters. We simulated different experimental designs and used different inference methods, varying the likelihood functions and other aspects of the parameter estimation methods. Finally, we applied the best performing inference methods to real data for the species Paramecium caudatum. The relative error of the estimates of activation energy varied between 5% and 30%. The fraction of habitat sampled played the most important role in determining the relative error; sampling at least 1% of the habitat kept it below 50%. We found that methods that simultaneously use all time series data (direct methods) and methods that estimate population parameters separately for each temperature (indirect methods) are complementary. Indirect methods provide a clearer insight into the shape of the functional form describing the temperature dependence of population parameters; direct methods enable a more accurate estimation of the parameters of such functional forms. Using both methods, we found that growth rate and carrying capacity of Paramecium caudatum scale with temperature according to different activation energies. Our study shows how careful choice of experimental design and inference methods can increase the accuracy of the inferred relationships between temperature and population parameters. The comparison of estimation methods provided here can increase the accuracy of model predictions, with important implications in understanding and predicting the effects of temperature on the dynamics of populations.     

Effects of Coyote Population Reduction on Swift Fox Demographics in Southeastern Colorado

ABSTRACT The distribution and abundance of swift foxes (Vulpes velox) has declined from historic levels. Causes for the decline include habitat loss and fragmentation, incidental poisoning, changing land use practices, trapping, and predation by other carnivores. Coyotes (Canis latrans) overlap the geographical distribution of swift foxes, compete for similar resources, and are a significant source of mortality amongst many swift fox populations. Current swift fox conservation and management plans to bolster declining or recovering fox populations may include coyote population reduction to decrease predation. However, the role of coyote predation in swift fox population dynamics is not well-understood. To better understand the interactions of swift foxes and coyotes, we compared swift fox population demographics (survival rates, dispersal rates, reproduction, density) between areas with and without coyote population reduction. On the Piñon Canyon Maneuver Site, Colorado, USA, we monitored 141 swift foxes for 65,226 radio-days from 15 December 1998 to 14 December 2000 with 18,035 total telemetry locations collected. Juvenile swift fox survival rate was increased and survival was temporarily prolonged in the coyote removal area. Adult fox survival patterns were also altered by coyote removal, but only following late-summer coyote removals and, again, only temporarily. Coyote predation remained the main cause of juvenile and adult fox mortality in both areas. The increase in juvenile fox survival in the coyote removal area resulted in a compensatory increase in the juvenile dispersal rate and an earlier pulse in dispersal movements. Adult fox dispersal rate was more consistent throughout the year in the coyote removal area. Coyote removal did not influence the reproductive parameters of the swift foxes. Even though juvenile survival increased, swift fox density remained similar between the areas due to the compensatory dispersal rate among juvenile foxes. We concluded that the swift fox population in the area was saturated. Although coyote predation appeared additive in the juvenile cohort, it was compensatory with dispersal.