Dynamics of a low‐density tiger population in Southeast Asia in the context of improved law enforcement

Recovering small populations of threatened species is an important global conservation strategy. Monitoring the anticipated recovery, however, often relies on uncertain abundance indices rather than on rigorous demographic estimates. To counter the severe threat from poaching of wild tigers (Panthera tigris), the Government of Thailand established an intensive patrolling system in 2005 to protect and recover its largest source population in Huai Kha Khaeng Wildlife Sanctuary. Concurrently, we assessed the dynamics of this tiger population over the next 8 years with rigorous photographic capture‐recapture methods. From 2006 to 2012, we sampled across 624–1026 km² with 137–200 camera traps. Cameras deployed for 21,359 trap days yielded photographic records of 90 distinct individuals. We used closed model Bayesian spatial capture‐recapture methods to estimate tiger abundances annually. Abundance estimates were integrated with likelihood‐based open model analyses to estimate rates of annual and overall rates of survival, recruitment, and changes in abundance. Estimates of demographic parameters fluctuated widely: annual density ranged from 1.25 to 2.01 tigers/100 km², abundance from 35 to 58 tigers, survival from 79.6% to 95.5%, and annual recruitment from 0 to 25 tigers. The number of distinct individuals photographed demonstrates the value of photographic capture–recapture methods for assessments of population dynamics in rare and elusive species that are identifiable from natural markings. Possibly because of poaching pressure, overall tiger densities at Huai Kha Khaeng were 82–90% lower than in ecologically comparable sites in India. However, intensified patrolling after 2006 appeared to reduce poaching and was correlated with marginal improvement in tiger survival and recruitment. Our results suggest that population recovery of low‐density tiger populations may be slower than anticipated by current global strategies aimed at doubling the number of wild tigers in a decade.     

Population viability of the Siberian Tiger in a changing landscape: Going, going and gone?

The Amur tiger (Panthera tigris altaica) is a flagship species of the boreal forest ecosystem in northeastern China and Russia Far East. During the past century, the tiger population has declined sharply from more than 3000 to fewer than 600 individuals, and its habitat has become much smaller and greatly fragmented. Poaching, habitat degradation, habitat loss, and habitat fragmentation have been widely recognized as the primary causes for the observed population decline. Using a population viability analysis tool (RAMAS/GIS), we simulated the effects of poaching, habitat degradation, habitat loss, and habitat fragmentation on the population dynamics and extinction risk of the Amur tiger, and then explored the relative effectiveness of three conservation strategies involving improving habitat quality and establishing movement corridors in China and Russia. A series of controlled simulation experiments were performed based on the current spatial distribution of habitat and field-observed vital rates. Our results showed that the Amur tiger population could be viable for the next 100 years if the current habitat area and quality were well-maintained, with poaching strictly prohibited of the tigers and their main prey species. Poaching and habitat degradation (mainly prey scarcity) had the largest negative impacts on the tiger population persistence. While the effect of habitat loss was also substantial, habitat fragmentation per se had less influence on the long-term fate of the tiger population. However, to sustain the subpopulations in both Russia and China would take much greater conservation efforts. The viability of the Chinese population of tigers would rely heavily on its connectivity with the largest patch on the other side of the border. Improving the habitat quality of small patches only or increasing habitat connectivity through movement corridors alone would not be enough to guarantee the long-term population persistence of the Amur tiger in both Russia and China. The only conservation strategy that allowed for long-term persistence of tigers in both countries required both the improvement of habitat quality and the establishment of a transnational reserve network. Our study provides new insights into the metapopulation dynamics and persistence of the Amur tiger, which should be useful in landscape and conservation planning for protecting the biggest cat species in the world.     

Killing Tigers to Save Them: Fallacies of the Farming Argument

Abstract: The lucrative, illegal trade in tigers (Panthera tigris) remains a major conservation problem. Tiger farming has been proposed as a potential solution, with farmed tigers substituting for wild tigers. At first glance, this argument's logic seems simple: farming will increase the supply of tigers, prices will fall, and poaching will no longer be profitable. We contend, however, that this supply‐side argument relies on mistaken assumptions. First, tiger markets are imperfect, meaning they are dominated by a few producers who control price. Second, consumers prefer wild tigers to farmed tigers and therefore the two are not pure substitutes. In economic terms, products from wild tigers are luxury goods, commanding a price premium. Third, there is no evidence that farmed tigers can be produced or sold more cheaply than wild tigers. In sum, it is unlikely that farming will drive down the price of wild‐caught tigers or decrease profitability for tiger poachers. Rather, tiger farming is more likely to increase aggregate demand for tiger products and stimulate higher levels of poaching.     

Genetic and demographic consequences of importing animals into a small population: a simulation model of the Texas State Bison Herd (USA)

The extant 40 bison (Bison bison) constituting the Texas State Bison Herd (TSBH; USA) are directly and exclusively descended from a bison herd assembled by Charles Goodnight in the 1880s, representing a historically and genetically valuable resource. The population currently suffers from low genetic variation, low heterozygosity, high calf mortality, and low natality rates compared with other closed bison populations. Population viability analysis using the VORTEX program previously indicated a 99% chance of population extinction within the next 41 years [J. Mamm. 85 (2004) in press]. We developed a stochastic simulation model to evaluate the genetic and demographic consequences of various management scenarios for the TSBH using genotypic data from 51 microsatellite loci and demographic information recorded over a 6-year period. Our results reveal that without the introduction of new genetic variation, approximately 37% of the representative microsatellite loci will become fixed as the TSBH continues to lose genetic variation at a staggering rate of 30–40% within the next 50 years. Furthermore, if the current trends in natality and mortality rates continue, our model indicates the TSBH faces a 99% chance of extinction in the next 51 years. With the importation of unrelated male bison into the TSBH, and under the assumption of increased fitness, the probability of population survival in the next 100 years increases to 100%, and the population will reach the approximate carrying capacity of 200 bison in 15–16 years. Furthermore, our model predicts increases in genetic diversity and heterozygosity of 24.7–48.4% and 17.5–36.5%, respectively, in the next 100 years following the addition of new genetic variation. We conclude that the importation of bison into the TSBH is necessary to prevent extinction and ensure long-term population survival.     

Implications of dispersal and life history strategies for the persistence of Linyphiid spider populations

Linyphiid spiders have evolved the ability to disperse long distances by a process known as ballooning. It has been hypothesized that ballooning may allow populations to persist in the highly disturbed agricultural areas that the spiders prefer. In this study, I develop a stochastic population model to explore how the propensity for this type of long distance dispersal influences long term population persistence in a heterogeneous landscape where catastrophic mortality events are common. Analysis of this model indicates that although some dispersal does indeed decrease the probability of extinction of the population, the frequency of dispersal is only important in certain extremes. Instead, both the mean population birth and death rates, and the landscape composition, are much more important in determining the probability of extinction than the dispersal process. Thus, in order to develop effective conservation strategies for these spiders, better understanding of life history processes should be prioritized over an understanding of dispersal strategies.     

Population Viability Analysis for an Endangered Plant

Abstract: Demographic modeling is used to understand the population viability of Furbish's lousewort, Pedicularis furbishiae , a perennial plant species endemic to the St. John River Valley in northern Maine. Environment-specific summaries of demographic parameters (survivorship, growth, and fecundity) over four years, organized into stage-based projection matrices, provide predictions of future population dynamics given a deterministic extension of past conditions. Stochastic modeling, using (I) empirically observed variation in demographic parameters, and (2) estimated rates of natural catastrophes, leads to predictions of extinction probability.
P. furbishiae viability has varied widely over the study period Viable populations with finite rates of increase > 1 are found where cover is low, woody plants do not dominate, and disturbance does not occur. Rates of increase vary over time, suggesting that stochastic analyses would be realistic. Stochastic measures of population viability incorporating environmental variation suggest that early successional environments, especially wetter sites, can support viable populations in the absence of disturbance. However; observed rates of natural catastrophe dominate viability estimates of individual populations. Metapopulation dynamics feature extinction rates that are greater than recolonization rates, and may be affected by land use in the watershed Species management needs to consider a large-scale view of the riverine corridor.     

Erosion of Heterozygosity in Fluctuating Populations

Abstract: Demographic, environmental, and genetic stochasticity threaten the persistence of isolated populations. The relative importance of these intertwining factors remains unresolved, but a common view is that random demographic and environmental events will usually drive small populations to the brink of extinction before genetic deterioration poses a serious threat. To evaluate the potential importance of genetic factors, we analyzed a model linking demographic and environmental conditions to the loss of genetic diversity in isolated populations undergoing natural levels of fluctuation. Nongenetic processes—environmental stochasticity and population demography—were modeled according to a bounded diffusion process. Genetic processes were modeled by quantifying the rate of drift according to the effective population size, which was predicted from the same parameters used to describe the nongenetic processes. We combined these models to predict the heterozygosity remaining at the time of extinction, as predicted by the nongenetic portion of the model. Our model predicts that many populations will lose most or all of their neutral genetic diversity before nongenetic random events lead to extinction. Given the abundant evidence for inbreeding depression and recent evidence for elevated extinction rates of inbred populations, our findings suggest that inbreeding may be a greater general threat to population persistence than is generally recognized. Therefore, conservation biologists should not ignore the genetic component of extinction risk when assessing species endangerment and developing recovery plans.     

Cover Caption

Cover : Farming of tigers (Panthera tigris) has been proposed as a disincentive to illegal trade that is decreasing the abundance and distribution of tigers in the wild. In this issue, Kirkpatrick and Emerton (pp. 655–659) examine the economic rationale for tiger farming. In theory, the pro. tability of poaching will decrease as availability of tigers increases. Given the consumer preference for wild tigers and price controls by traders, however, farming actually may increase demand for wild tigers and lead to a greater incidence of poaching.     

Population viability analysis for several populations using multivariate state-space models

The International Union for the Conservation of Nature and Natural Resources (IUCN), the world's largest and most important global conservation network, has listed approximately 16,000 species worldwide as threatened. The most important tool for recognizing and listing species as threatened is population viability analysis (PVA), which estimates the probability of extinction of a population or species over a specified time horizon. The most common PVA approach is to apply it to single time series of population abundance. This approach to population viability analysis ignores covariability of local populations. Covariability can be important because high synchrony of local populations reduces the effective number of local populations and leads to greater extinction risk. Needed is a way of extending PVA to model correlation structure among multiple local populations. Multivariate state-space modeling is applied to this problem and alternative estimation methods are compared. The multivariate state-space technique is applied to endangered populations of pacific salmon, USA. Simulations demonstrated that the correlation structure can strongly influence population viability and is best estimated using restricted maximum likelihood instead of maximum likelihood.     

Effects of Social Structure and Prey Dynamics on Extinction Risk in Gray Wolves

Extinction models based on diffusion theory generally fail to incorporate two important aspects of population biology—social structure and prey dynamics. We include these aspects in an individual-based extinction model for small, isolated populations of the gray wolf (Canis lupus). Our model predicts mean times to extinction significantly longer than those predicted by more general (diffusion) models. According to our model, an isolated population of 50 wolves has a 95% chance of surviving just 9 years and only a 30% chance of surviving beyond 100 years. Reflecting the influence of social structure, a wolf population initially comprising 50 individuals is expected to persist only a few years longer, on average (71 years), than is a population initially comprising just a single reproductive pair (62 years). In contrast, substantially greater average prey abundance leads to dramatically longer expected persistence times. Autocorrelated prey dynamics result in a more complex distribution of extinction times than predicted by many extinction models. We contend that demographic stochasticity may pose the greatest threat to small, isolated wolf populations, although environmental stochasticity and genetic effects may compound this threat. Our work highlights the importance of considering social structure and resource dynamics in the development of population viability analyses.     

Simulation of the Effect of Spatial and Temporal Variation in Fire Regimes on the Population Viability of a Banksia Species

Populations of plants that rely on seeds for recovery from disturbance by fire (obligate seeders) are sensitive to regimes of frequent fire. Obligate seeders are prominent in fire-prone heathlands of southern Australia and South Africa. Population extinction may occur if there are successive fires during a plant's juvenile period. Research on the population biology of obligate seeders has influenced the management of fire in these heath and shrublands, but work on the effects of the spatial variability of fires is lacking. We hypothesize that extinction maybe avoided under an adverse fire frequency if fires are patchy. We present a model that simulates the effects of spatial and temporal variations in fire regimes on the viability of a plant population in a grid landscape. Seedling establishment, maturation, senescence, and seed dispersal determine the presence or absence of plants in each cell. We used values typical of serotinous Banksia species to estimate probability of extinction in relation to fire frequency and size. We examined the sensitivity of predictions to dispersal, senescence, fire frequency, spatial burning pattern and size variance, and the size of the grid. Simulations 200 years in length indicated that extinction probability was lowest when mean fire frequency was intermediate and mean fire size was large. When fire frequency was high, extinction probability was high irrespective of fire size. Senescence was more important than high-frequency fire as a cause of extinction in cells. Interactions between dispersal, fire frequency, and size were complex, indicating that extinction is governed by intercell connectivity. The model indicates that fire patchiness cannot be assumed to ensure avoidance of extinction of populations. Conservation of populations is most likely when fire patchiness is relatively low—when the size of fires is moderate to large and when burned patches are contiguous.     

The Role of Translocation in Recovery of Woodland Caribou Populations

Abstract: Maintenance of viable populations of many endangered species will require conservation action in perpetuity. Efforts to conserve these species are more likely to be successful if their reliance on conservation actions is assessed at the population level. Woodland caribou (Rangifer tarandus caribou) were extirpated recently from Banff National Park, Canada, and translocations of caribou to Banff and neighboring Jasper National Park are being considered. We used population viability analysis to assess the relative need for and benefits from translocation of individuals among caribou populations. We measured stochastic growth rates and the probability of quasi extinction of four populations of woodland caribou with and without translocation. We used two vital rates in our analysis: mean adult female survival and mean number of calves per breeding‐age female as estimates of mean fecundity. We isolated process variance for each vital rate. Our results suggested the Tonquin caribou population in Jasper is likely to remain viable without translocation, but that translocation is probably insufficient to prevent eventual extirpation of the two other populations in Jasper. Simulated reintroductions of caribou into Banff resulted in a 53–98% probability of >8 females remaining after 20 years, which suggests translocation may be an effective recovery tool for some caribou populations.     

Habitat Change and Plant Demography: Assessing the Extinction Risk of a Formerly Common Grassland Perennial

Abstract:  An important aim of conservation biology is to understand how habitat change affects the dynamics and extinction risk of populations. We used matrix models to analyze the effect of habitat degradation on the demography of the declining perennial plant Trifolium montanum in 9 calcareous grasslands in Germany over 4 years and experimentally tested the effect of grassland management. Finite population growth rates (λ) decreased with light competition, measured as leaf-area index above T. montanum plants. At unmanaged sites λ was <1 due to lower recruitment and lower survival and flowering probability of large plants. Nevertheless, in stochastic simulations, extinction of unmanaged populations of 100 flowering plants was delayed for several decades. Clipping as a management technique rapidly increased population growth because of higher survival and flowering probability of large plants in managed than in unmanaged plots. Transition-matrix simulations from these plots indicated grazing or mowing every second year would be sufficient to ensure a growth rate ≥1 if conditions stayed the same. At frequently grazed sites, the finite growth rate was approximately 1 in most populations of T. montanum . In stochastic simulations, the extinction risk of even relatively small grazed populations was low, but about half the extant populations of T. montanum in central Germany are smaller than would be sufficient for a probability of survival of >95% over 100 years. We conclude that habitat change after cessation of management strongly reduces recruitment and survival of established individuals of this perennial plant. Nevertheless, our results suggest extinction processes may take a long time in perennial plants, resulting in an extinction debt. Even if management is frequent, many remnant populations of T. montanum may be at risk because of their small size, but even a slight increase in size could considerably reduce their extinction risk.     

Geometrical Factors and Metapopulation Dynamics of the Bush Cricket, Metrioptera bicolor Philippi (Orthoptera: Tettigoniidae)

This paper presents a metapopulation study of the bush cricket, Metrioptera bicolor , living in a recently fragmented landscape. The species inhabits grass and heathland patches of varying area and isolation. Analyses are made of how these geometrical factors affect local population size and density, distribution pattern, and the probability of local extinction and colonization. The proportion of available patches occupied varied between 72 and 79% during 1985–1990. Unoccupied patches were smaller and more isolated than those that were occupied. Patches where populations became extinct during this period were smaller than those with persisting populations. Since local population size was well correlated with patch area, it was clear that stochastic extinctions only occurred in small populations. Critical patch size for population extinction was approximately half a hectare. Colonized patches were less isolated than those that had not been colonized. Critical inter-patch distance for colonization was about 100 meters. The turnover was restricted to an identifiable share of the available patches. Only 33% of the patches were so small that extinction due to stochastic causes could be considered highly probable. This metapopulation will therefore most likely persist over a considerable period in its present spatial structure. There are apparent threats of further fragmentation, however, and nothing is known about the likelihood of large-scale extinctions resulting from extremely unfavorable weather conditions. Nevertheless, our results show that it is appropriate to include geometrical factors in metapopulation models.     

Genetic and Demographic Threats to the Black Rhinoceros Population in the Ngorongoro Crater

A resident population of 13 black rhinoceros ( Diceros bicornis ) persist in Ngorongoro Crater, Tanzania. The effective population size ( N _e ) may be as few as 5 animals. Projected growth for this population suggests that the effective population size will remain small for the near future, threatening this Iocal population with extinction due to the stochastic factors associated with small population size. A summary of historic and recent demographic data for this population reveals a population crash during the period of heavy poaching that affected this species throughout its range. Although poaching of this species has been brought under control the population remains small. These data and models of projected population growth argue for consideration of more-intensive management within the framework of the small population paradigm. This case is an example of applied conservation resulting from this paradigm used in conjunction with rather than competing with the declining population paradigm. We identify additional monitoring, particularly of density-dependent behaviors, that will be necessary for designing a successful management program. Finally, the use of molecular markers for developing an accurate pedigree for this population is suggested in order to maintain a genetically healthy population. These strategies have broad applicability to black rhinoceros conservation throughout Africa.     

Predicting the Deleterious Effects of Mutation Load in Fragmented Populations

Abstract:  Human-induced habitat fragmentation constitutes a major threat to biodiversity. Both genetic and demographic factors combine to drive small and isolated populations into extinction vortices. Nevertheless, the deleterious effects of inbreeding and drift load may depend on population structure, migration patterns, and mating systems and are difficult to predict in the absence of crossing experiments. We performed stochastic individual-based simulations aimed at predicting the effects of deleterious mutations on population fitness (offspring viability and median time to extinction) under a variety of settings (landscape configurations, migration models, and mating systems) on the basis of easy-to-collect demographic and genetic information. Pooling all simulations, a large part (70%) of variance in offspring viability was explained by a combination of genetic structure ( F_ST ) and within-deme heterozygosity ( H_S ). A similar part of variance in median time to extinction was explained by a combination of local population size ( N ) and heterozygosity ( H_S ). In both cases the predictive power increased above 80% when information on mating systems was available. These results provide robust predictive models to evaluate the viability prospects of fragmented populations.     

Rates of Movement of Threatened Bird Species between IUCN Red List Categories and toward Extinction

Abstract:  In recent centuries bird species have been deteriorating in status and becoming extinct at a rate that may be 2–3 orders of magnitude higher than in prehuman times. We examined extinction rates of bird species designated critically endangered in 1994 and the rate at which species have moved through the IUCN (World Conservation Union) Red List categories of extinction risk globally for the period 1988–2004 and regionally in Australia from 1750 to 2000. For Australia we drew on historical accounts of the extent and condition of species habitats, spread of invasive species, and changes in sighting frequencies. These data sets permitted comparison of observed rates of movement through the IUCN Red List categories with novel predictions based on the IUCN Red List criterion E, which relates to explicit extinction probabilities determined, for example, by population viability analysis. The comparison also tested whether species listed on the basis of other criteria face a similar probability of moving to a higher threat category as those listed under criterion E. For the rate at which species moved from vulnerable to endangered, there was a good match between observations and predictions, both worldwide and in Australia. Nevertheless, species have become extinct at a rate that, although historically high, is 2 (Australia) to 10 (globally) times lower than predicted. Although the extinction probability associated with the critically endangered category may be too high, the shortfall in realized extinctions can also be attributed to the beneficial impact of conservation intervention. These efforts may have reduced the number of global extinctions from 19 to 3 and substantially slowed the extinction trajectory of 33 additional critically endangered species. Our results suggest that current conservation action benefits species on the brink of extinction, but is less targeted at or has less effect on moderately threatened species.     

Effects of learning and adaptation on population viability

Cultural adaptation is one means by which conservationists may help populations adapt to threats. A learned behavior may protect an individual from a threat, and the behavior can be transmitted horizontally (within generations) and vertically (between generations), rapidly conferring population-level protection. Although possible in theory, it remains unclear whether such manipulations work in a conservation setting; what conditions are required for them to work; and how they might affect the evolutionary process. We examined models in which a population can adapt through both genetic and cultural mechanisms. Our work was motivated by the invasion of highly toxic cane toads (Rhinella marina) across northern Australia and the resultant declines of endangered northern quolls (Dasyurus hallucatus), which attack and are fatally poisoned by the toxic toads. We examined whether a novel management strategy in which wild quolls are trained to avoid toads can reduce extinction probability. We used a simulation model tailored to quoll life history. Within simulations, individuals were trained and a continuous evolving trait determined innate tendency to attack toads. We applied this model in a population viability setting. The strategy reduced extinction probability only when heritability of innate aversion was low (<20%) and when trained mothers trained >70% of their young to avoid toads. When these conditions were met, genetic adaptation was slower, but rapid cultural adaptation kept the population extant while genetic adaptation was completed. To gain insight into the evolutionary dynamics (in which we saw a transitory peak in cultural adaptation over time), we also developed a simple analytical model of evolutionary dynamics. This model showed that the strength of natural selection declined as the cultural transmission rate increased and that adaptation proceeded only when the rate of cultural transmission was below a critical value determined by the relative levels of protection conferred by genetic versus cultural mechanisms. Together, our models showed that cultural adaptation can play a powerful role in preventing extinction, but that rates of cultural transmission need to be high for this to occur.     

The Interplay between Climate Variability and Density Dependence in the Population Viability of Chinook Salmon

Abstract:  The viability of populations is influenced by driving forces such as density dependence and climate variability, but most population viability analyses (PVAs) ignore these factors because of data limitations. Additionally, simplified PVAs produce limited measures of population viability such as annual population growth rate (λ) or extinction risk. Here we developed a "mechanistic" PVA of threatened Chinook salmon ( Oncorhynchus tshawytscha ) in which, based on 40 years of detailed data, we related freshwater recruitment of juveniles to density of spawners, and third-year survival in the ocean to monthly indices of broad-scale ocean and climate conditions. Including climate variability in the model produced important effects: estimated population viability was very sensitive to assumptions of future climate conditions and the autocorrelation contained in the climate signal increased mean population abundance while increasing probability of quasi extinction. Because of the presence of density dependence in the model, however, we could not distinguish among alternative climate scenarios through mean λ values, emphasizing the importance of considering multiple measures to elucidate population viability. Our sensitivity analyses demonstrated that the importance of particular parameters varied across models and depended on which viability measure was the response variable. The density-dependent parameter associated with freshwater recruitment was consistently the most important, regardless of viability measure, suggesting that increasing juvenile carrying capacity is important for recovery.     

Regional Consequences of Local Population Demography and Genetics in Relation to Habitat Management in Gentiana pneumonanthe

Abstract:  A joint demographic and population genetics stage-based model for a subdivided population was applied to Gentiana pneumonanthe , an early successional perennial herb, at a regional (metapopulation) scale. We used numerical simulations to determine the optimal frequency of habitat disturbance (sod cutting) and the intensity of gene flow among populations of G. pneumonanthe to manage both population viability and genetic diversity in this species. The simulations showed that even small populations that initially had near-equal allele frequencies could, if managed properly through sod cutting every 6 to 7 years, sustain their high genetic variation over the long run without gene flow. The more the allele frequencies in the small populations are skewed, however, the higher the probability that in the absence of gene flow, some alleles will be lost and within-population genetic variation will decrease even under proper management. This implies that although local population dynamics should be the major target for management, regional dynamics become important when habitat fragmentation and decreased population size lead to the loss of local genetic diversity. The recommended strategy to improve genetic composition of small populations is the introduction of seeds or seedlings of nonlocal origin.