Evaluating recovery strategies for an ocelot (Leopardus pardalis) population in the united states

The ocelot Leopardus pardalis population in the United States was listed as endangered in 1982, with only two known isolated breeding populations occurring in southern Texas. Conservation concerns for ocelots include loss of dense thornshrub habitat, mortality from ocelot-vehicle collisions, and genetic erosion. In this study, we used a population viability analysis (PVA) to evaluate four recovery strategies (i.e., supplementation of additional ocelots, reduced road mortality, habitat protection and restoration, and linkage of two breeding populations) for ocelot conservation management. We used the VORTEX (Version 9.42) program to conduct our PVA for an ocelot population located in Cameron County, Texas. Each scenario was simulated 500 times over 100 years. We compared the effectiveness of recovery strategies and combinations thereof with estimates of extinction probability and final population size. Model scenarios with no recovery strategies predicted an extinction probability of 0.65 for the Cameron population of ocelots over 100 years. The protection and restoration of thornshrub habitat was the most effective recovery strategy, followed by population linkage and reduced road mortality, with the supplementation of ocelots being the least effective strategy. Protection and restoration of ocelot habitat cannot be accomplished without the participation of private landowners. Using an adaptive management approach, future actions need to be taken to monitor ocelot populations and habitats and help to reduce the high probability of extinction predicted in our PVA for the ocelot population in Cameron County.     

A habitat-based population viability analysis for ocelots (Leopardus pardalis) in the United States

Strategies are needed to recover the ocelot Leopardus pardalis from the endangered species list. Recently, a population viability analysis (PVA) was developed which concluded that combinations of different recovery strategies were needed to effectively reduce ocelot extinction probability in the United States (US), with habitat protection and restoration identified as the most effective recovery scenario. We expanded this PVA model by incorporating landscape data to develop a more realistic habitat-based PVA for ocelots in southern Texas. We used RAMAS/gis software to conduct a habitat-based PVA by linking landscape data with a demographic metapopulation model. The primary goal of this study was to provide a model for evaluating ocelot recovery strategies in the US. Each model scenario was simulated 1000 times over 50 years and we defined extinction as one individual remaining. Using the RAMAS/gis program we identified 11 possible ocelot habitat patches (i.e., subpopulations) occurring in southern Texas. In addition, based on the habitat-based PVA model we found that combinations of different recovery strategies were needed to effectively reduce ocelot extinction probability in the US, with reducing road mortality the single most effective strategy. Short-term recovery strategies should include reducing ocelot road mortality, and translocation of ocelots into the US from northern Mexico. Long-term recovery strategies should include the restoration of habitat between and around existing ocelot habitat patches and the establishment of a dispersal corridor between ocelot breeding populations.     

A quantitative framework to evaluate incidental take and endangered species population viability

Incidental take is the permitted unintended or collateral killing, harassment, or habitat destruction of a protected species under endangered species law and is permissible as long as the take does not jeopardize the species’ persistence in the wild. However, take is seldom addressed in a quantitative or population modeling context and the criterion of “jeopardy” has no universal, quantitative definition. To model the effect of incidental take on population viability, we modified an existing population model and simulated abundance, population growth, and quasi-extinction probability (the probability of falling below a predetermined abundance threshold) for Piping Plovers (Charadrius melodus) in the Great Plains, USA. The model incorporated environmental stochasticity and variation due to sampling variance. Eggs and chicks were taken out of the population as a “harvest” to simulate incidental take that currently occurs in the Missouri River system. We used least-squares regression and an AIC model selection approach to evaluate the population’s elasticity to incidental take covariates. Even in the absence of take the population declined by 7.5% annually. Population growth and final abundance were reduced and the probability of quasi-extinction was increased in simulations where egg take and chick take were applied. The model selection analysis indicated that incidental take of eggs and chicks depresses population viability and the probability of recovery of Piping Plovers in the Great Plains. Though the model was useful in putting take and jeopardy in a quantitative setting, the question remains as to whether permitted levels of take causes jeopardy for Piping Plovers since there are no decision standards defined by the USFWS. However, evaluating take in a quantitative framework, as we have, will make jeopardy decisions more explicit in terms of viability and recovery metrics.     

Minimum viable population and conservation status of the Atlantic Forest spiny rat Trinomys eliasi

A population viability analysis (PVA) was conducted to assess the minimum viable population (MVP) of the Atlantic Forest spiny rat Trinomys eliasi, a species threatened by habitat loss and restricted geographical distribution. Objectives were to suggest quasi-extinction thresholds, estimate minimum areas of suitable habitat (MASH) and MVPs, and compare results with the species’ current status. The computer package VORTEX was used. The model predicted sizes of 200 animals to achieve demographic stability, but buffering declines in genetic variability required populations of 2000 animals. Estimated MASHs were approximately 250 and 2500 ha for demographic and genetic stability, respectively. Mortality rate and mean litter size were the most sensitive parameters to changes in model assumptions. The protection of known populations and the search for extant populations are the first steps in conservation. T. eliasi's issue could help protecting the coastal shrubland ecosystem of Rio de Janeiro state. Observing IUCN's criteria for listing threatened species, it is suggested that T. eliasi should be ranked as vulnerable in red lists.     

Population viability of Cape mountain zebra in Gamka Mountain Nature Reserve, South Africa: the influence of habitat and fire

The small Cape mountain zebra population in Gamka Mountain Nature Reserve represents a third of the entire gene pool of this endangered species and is thus vital for it’s conservation. Presently, management of this population is largely hands off, with the belief that it will grow to levels which will allow it to form a source for the mixing of mountain zebra stocks in the future. The growth of this population however, has been slow and we investigated the influence of habitat and fire on this growth. Firstly, we used a diffusion model to perform a population viability analysis. This analysis indicated that the population had a low probability of attaining quasi-extinction in the next 50 years (G = 0.0032). However, our findings indicated that less than 30% of the reserve was suitable for mountain zebra and that the preferred habitat would have to be burnt at unnaturally short intervals to sustain the present growth. We therefore argue that the risk of quasi-extinction to this population is greater than predicted and suggest that management options need to be implemented to reduce this risk. These options include; translocation to another protected area; acquisition of adjacent land; burning preferred habitat at unnaturally short intervals; forming a conservancy with adjacent landowners; leasing cultivated land for pasture. We suggest that only the latter two options are likely to stimulate mountain zebra population growth in the short term and that these should receive immediate attention.     

Potential factors controlling the population viability of newly introduced endangered marble trout populations

While several population viability analyses (PVAs) have been performed on anadromous salmonids, less attention has been given to stream-living salmonids. In this work, we explore the role of PVA as a tool in the recovery of threatened stream-living salmonid species. The analysis has been performed with reference to marble trout Salmo marmoratus, a salmonid with a limited geographic distribution and at risk of extinction due to hybridization with the non-native introduced brown trout. Demographic parameters, such as survival, fecundity and density-dependent patterns were estimated from an eight year on-going monitoring program of two translocated marble trout populations in pristine, previously fishless streams (Zakojska and Gorska) in the Soca and Idrijca river basins (Slovenia). To explore the importance of disturbance events such as floods on marble trout population dynamics, we performed a PVA under three scenarios: (1) occurrence of both severe and moderate floods; (2) occurrence of only moderate floods; (3) no flood events. Our analysis shows that population viability is threatened only by severe flood events, otherwise the two populations prove to be fairly stable with population abundance fluctuating around stream carrying capacity. A sensitivity analysis performed on model parameters highlighted that density-dependence in first-year survival and the magnitude of reduction in population size after a severe flood are the two most crucial parameters affecting population abundance and quasi-extinction probability, respectively. While only extreme floods can drive the population to extinction, the increase in juvenile survival when population abundance collapses after a major flood may allow the populations to quickly recover from few reproductive individuals back to stream carrying capacity.     

Viability of mountain caribou in British Columbia, Canada: Effects of habitat change and population density

Population viability analyses (PVA) are frequently employed to develop recovery plans and inform management of endangered species. Translating results from PVA into meaningful management recommendations often depends on an understanding of how population parameters change with environmental conditions as well as population density. The decline of mountain caribou (Rangifer tarandus caribou) in British Columbia, Canada, is believed to be caused by apparent competition with alternative prey species following changes to the forest age structure from timber harvest and wildfire. In addition, populations have been shown to decline at faster rates at low population density. To evaluate the potential effects of habitat change and population density on population persistence, we used stochastic projection models for 10 distinct populations varying in initial size from <10 to approximately 150 females. In an initial model, we used estimates of vital rates based on information sampled from >350 radiocollared caribou between 1984 and 2004. We then compared the results of the initial model to a set of models that evaluated the effects of habitat conditions and population density via their expected relationships to female adult survival. Assuming that vital rates remain constant over a 200-year time frame, only three populations have high probabilities (>0.95) of extinction. When models incorporate the declines in adult female survival know to occur with increasing proportions of young forest and declining population densities, all 10 populations are predicted to decline to extinction within <200 years. Based on our results, we suggest that PVA models that fail to incorporate the effects of changes in vital rates with habitat and population density may lead to overly optimistic assessments of the probability of population persistence in endangered species.     

Preliminary modelling of sediment production and delivery in the Xihanshui River basin, Gansu, China

This paper outlines an analysis of the spatial distribution of sediment production, delivery and yield in the Xihanshui River basin, South Gansu, China, using the modelling tools of SedNet (Prosser et al., 2001). This model can assess the delivery efficiency to downstream locations, as well as identifying locations with high rates of sediment production. Preliminary model experiments assist understanding of the spatial dynamics of these sediment processes and evaluation of the effectiveness of soil conservation practices since the mid-1980s. Three scenario years (dry, average and wet) from the 1983–2005 record are identified and modelled, and land use and management are represented in the model to reflect known changes since the 1980s. Results show hillslope erosion to be a dominant source of sediment supply, causing the latter to decrease ten-fold between 1984 and 1997/2000. Estimated bank erosion and floodplain deposition rates are sensitive to parameter values, but bank erosion appears less sensitive than hillslope supply to rainfall. The model can be used to assess net changes in floodplain storage; for default parameters, floodplain deposition rates are 25–200 times the rates of bank erosion depending on the climate scenario. Comparing simulation results with measured sediment yields at the three gauging stations indicates encouraging agreement in 2000. In 1984 (the wet year), the model under-predicts, suggesting that additional unmodelled sediment production processes, especially mass movement and gully erosion, may be important in wet years. Mass movement inventory data could close the gap between the high yields measured in the wet scenario year and the estimated yield due to hillslope erosion alone. In 1997 (the dry year), the model over-predicts; this suggests that the land use change parameters required to reflect the effects of conservation may not have been sufficient, implying that conservation has been generally effective, and that evidence of declining sediment yield is not simply a reflection of drier conditions.     

Preliminary modelling of sediment production and delivery in the Xihanshui River basin,Gansu, China

This paper outlines an analysis of the spatial distribution of sediment production, delivery and yield in the Xihanshui River basin, South Gansu, China, using the modelling tools of SedNet (Prosser et al., 2001). This model can assess the delivery efficiency to downstream locations, as well as identifying locations with high rates of sediment production. Preliminary model experiments assist understanding of the spatial dynamics of these sediment processes and evaluation of the effectiveness of soil conservation practices since the mid-1980s. Three scenario years (dry, average and wet) from the 1983–2005 record are identified and modelled, and land use and management are represented in the model to reflect known changes since the 1980s. Results show hillslope erosion to be a dominant source of sediment supply, causing the latter to decrease ten-fold between 1984 and 1997/2000. Estimated bank erosion and floodplain deposition rates are sensitive to parameter values, but bank erosion appears less sensitive than hillslope supply to rainfall. The model can be used to assess net changes in floodplain storage; for default parameters, floodplain deposition rates are 25–200 times the rates of bank erosion depending on the climate scenario. Comparing simulation results with measured sediment yields at the three gauging stations indicates encouraging agreement in 2000. In 1984 (the wet year), the model under-predicts, suggesting that additional unmodelled sediment production processes, especially mass movement and gully erosion, may be important in wet years. Mass movement inventory data could close the gap between the high yields measured in the wet scenario year and the estimated yield due to hillslope erosion alone. In 1997 (the dry year), the model over-predicts; this suggests that the land use change parameters required to reflect the effects of conservation may not have been sufficient, implying that conservation has been generally effective, and that evidence of declining sediment yield is not simply a reflection of drier conditions.     

Estimating the rate of quasi-extinction of the Australian grey nurse shark (Carcharias taurus) population using deterministic age- and stage-classified models

Grey nurse sharks off the east coast of Australia are listed nationally as “critically endangered” under Schedule 1 of the Environmental Protection and Biodiversity Conservation Act (1999) and may number no more than 300 in New South Wales and southern Queensland waters. They are an inshore, coastal dwelling species and were severely depleted by spearfishing in the 1960s. The population has continued to decline despite protection since 1984. Their life history (long-lived to 25+ years), late maturation (6-8 years), low fecundity (maximum 2 live young biennially), specific habitat requirements, limited inshore distribution, and small population size render them particularly vulnerable to extinction. We estimated the time to quasi-extinction (years elapsed for the population to consist of ?50 females) for the grey nurse shark population off the east coast of Australia based on current estimates of abundance and known anthropogenic rates of mortality. Estimated minimum population size was 300 as of 2002, and minimum anthropogenic mortality assessed from recovered carcasses was 12/year of which 75% were females. We modelled time to quasi-extinction using deterministic age- and stage-classified models for worst-, likely and best-case scenarios. Population size was estimated at 300 (worst), 1000 (likely) and 3000 (best). Anthropogenic mortality was added to the model assuming either all carcasses are being recovered (best), or conservatively, that only 50% are reported (realistic). Depending on model structure, if all carcasses are being reported, quasi-extinction times for worst-, likely and best-case scenarios range from 13 to 16 years, 84-98 years and 289-324 years, respectively. If under-reporting is occurring, time to quasi-extinction ranges from 6 to 8 years, 45-53 years and 173-200 years, respectively. In all scenarios modelled the grey nurse shark population will decline if no further steps are taken to remove anthropogenic sources of mortality. Because estimates of quasi-extinction rate depend on initial population size, and sensitivity analysis revealed that population rate of change was most sensitive to changes in the survival probability of the smallest length classes, obtaining precise estimates of abundance and annual survival of young females is critical.     

Probable consequences of high female mortality for speckled warblers living in habitat remnants

The effects of habitat fragmentation on the Australian avifauna have been widespread with species richness and abundance declining with reduced remnant size and habitat quality and increased habitat isolation. The speckled warbler, Chthonicola sagittata is one species from the highly fragmented temperate woodlands of eastern Australia that has declined across its range and populations that remain appear to be patchily distributed in habitat remnants. Specific causes of decline are unknown but several aspects of its biology make the species particularly vulnerable to decline in fragmented landscapes. Here, we analyse survey data (presence/absence) of speckled warblers in a large sample of habitat remnants from three regions to identify patterns of occupancy. We explore the effects of patch size on extinction risk using population viability analyses (PVA) and detailed demographic data from a behavioural study of individuals in the Australian Capital Territory, south-eastern Australia. Patch size was a strong predictor of the persistence of speckled warblers in habitat remnants. High density populations had higher probabilities of persisting, and inclusion of an Allee effect during drought decreased the probability of persistence. In the absence of an Allee effect, only high density populations in patches greater than 300 ha and low density populations in patches greater than700 ha had more than an 80% probability of persisting over 100 years. The accelerating decrease in population persistence below approximately 200-400 ha suggests that small populations were particularly vulnerable to stochastic demographic and environmental events. Adult female mortality was the single most important factor in driving population extinction. Our PVA model predictions matched the survey data for the Australian Capital Territory region remarkably well, but failed to predict occupancy rates in remnants in other regions. Differences in occupancy patterns between regions may, however, have resulted from inbreeding depression. This study demonstrates both the strengths and limitations of PVA analysis. PVA can predict occupancy patterns with reasonable accuracy, given good demographic data, but data for one region cannot be used universally for all regions. We highlight the need for studies of demography in different regions to interpret regional patterns of occupancy and to identify mechanisms of decline in remnant habitat.     

Integrating demographic and genetic approaches in plant conservation

We summarize the problems that populations of formerly common plants may encounter when habitat fragmentation isolates them and reduces population size. Genetic erosion, inbreeding depression, Allee-effects on reproductive success, catastrophes and environmental stochasticity are illustrated with studies on species that have recently become rare in The Netherlands due to habitat fragmentation. These clearly indicate that population viability is negatively affected. We also show that in the recent literature (since 1980), most studies on the conservation of rare plants have addressed population genetic structure and relationships between genetic variation and population size. Though important, these studies are not suitable for assessing the importance of genetics for population viability. In turn, demographic studies can detect changes in vital rates in small populations, but cannot reveal underlying genetic causes. Fitness and demographic studies are also well-represented in the literature, but remarkably few studies have attempted to integrate empirical demographic and genetic studies. We discuss two approaches to fill this very important lacuna in our knowledge. One of these constructs matrix-projection models on the basis of demographic censuses of—if possible—large and viable populations, and combines these with the results of experiments to determine inbreeding effects on demographic transitions and, subsequently, population growth and extinction. The other approach is to demographically monitor experimentally created small and large populations with low and high genetic variation and measure their actual growth rates and probabilities of extinction. We conclude that demography and demographic-genetic experiments should play a central role in plant conservation genetics.     

Viability and patch occupancy of a swamp rabbit metapopulation at the northern edge of its distribution

Swamp rabbits (Sylvilagus aquaticus) are state-endangered in Indiana, USA, and population decline has been attributed to habitat loss. We conducted pellet surveys as part of a long-term survey effort that has been conducted at approximate 10-year intervals over the last 40 years. We modeled patch occupancy and conducted a spatially-explicit population viability analysis (PVA). Although occupancy of individual patches varied over time, occupancy rate has been constant for the last 30 years, and Indiana swamp rabbits exist as a metapopulation that appears to be stable. Metapopulation dynamics were best characterized as being stationary, but area was an important factor in extinction rates; occupied patches (142 ± 37 ha) were significantly larger (P = 0.01) than unoccupied patches (79 ± 20 ha). We did not find strong support for models with colonization rates as a function of distance to neighboring patches, nor was distance to contiguous patches of habitat significantly different (P = 0.12) for occupied and unoccupied sites. Population viability analysis corroborated our findings based on occupancy modeling, and evaluation of the PVA model using occupancy data for the period 1985–2006 resulted in predictions that nearly matched our field observations (33% observed patch occupancy vs. 25% predicted patch occupancy). Population viability was most sensitive to reductions in survival and fecundity rates, but was otherwise robust to changes in parameters such as initial abundance and carrying capacity. Our findings provide novel insights into a poorly studied member of Sylvilagus and into species metapopulation dynamics at the edge of the range.     

Metapopulation dynamics and conservation of the marsh fritillary butterfly: Population viability analysis and management options for a critically endangered species in Western Europe

This study investigates the dynamics and viability of a marsh fritillary butterfly Euphydryas aurinia metapopulation in a Belgian successional landscape. Based on capture-mark-recapture and winter nest census data, we first estimated demography (survival and recruitment rates, population size, density dependence) and dispersal parameters (emigration rate, effect of patch connectivity on dispersal, mortality during dispersal). Then using RAMAS/GIS platform, we parameterised a population viability analysis (PVA) model with these parameters to simulate the future of this metapopulation under different scenarios.The metapopulation does not seem viable even if natural reforestation is controlled by adequate management. In its present state, the patch system is not able to sustain enough individuals: due to the large temporal fluctuations in demographic parameters, a carrying capacity far higher than currently would be necessary to limit extinction risk to 1%, suggesting the existence of an extinction debt for the species in Belgium. The situation of E. aurinia appears much worse compared to two other fritillary species threatened in Belgium, for which similar PVA are available. It is therefore urgent to increase the carrying capacity of the patch system. How and where it is achieved are of secondary importance for the gain in viability: improvement of habitat quality through restoration, or increase of habitat quantity via enlargement of existing patches and/or creation of new habitat in the matrix. A regime of management based on regular re-opening and maintenance of habitat patches may be the only guarantee of long-term persistence for this critically endangered species in Belgium.     

Population viability and species interactions: Life outside the single-species vacuum

Population viability analysis (PVA) has become a widely used set of tools for evaluating relative extinction risk and prioritizing management options among imperiled populations. While PVA is a widely sanctioned tool in conservation biology, the field of population viability is in its infancy with respect to species interactions. In this paper, I review available methods for evaluating extinction risk when species interactions contribute significantly to population viability. This review includes an evaluation of six broad categories of species interactions (predation, disease, competition, mutualism, parasitism and host-parasitoid interactions) in population viability analysis, with a particular focus on predation as a case study. I first evaluated how often species interactions are considered when PVA is applied to population data from imperiled species. I identified 378 articles in commonly cited conservation journals, of which 24 attempted a viability analysis for populations threatened by interactions with other species. Most of these PVA’s treat a putative species interaction as a constant source of mortality rather than a coupled, dynamic population process. Second, I reviewed the literature to identify the availability of time-series of abundance data for two interacting species in which at least one species was threatened or endangered. Adequate time-series data were available for both species comprising an interacting pair in only 9 out of 407 papers reviewed. Third, I used a stochastic, fully stage-structured predator prey model to create time-series data (vital rates and projection matrices) in order to quantify the efficacy of two matrix-based, single-species PVA approaches. Simple single-species PVAs confound stochastic variation with population cycles induced by species interactions (in this case predation). As a result these models provide conservatively biased forecasts of viability. Unfortunately, the data needed to construct more complex PVA’s with feedback and multi-species stochasticity are rarely collected. I close with a discussion of key advances needed to “escape the population vacuum” in a move toward more realistic estimates of extinction risk.     

The role of demographic compensation theory in incidental take assessments for endangered species

Many endangered species laws provide exceptions to legislated prohibitions through incidental take provisions as long as take is the result of unintended consequences of an otherwise legal activity. These allowances presumably invoke the theory of demographic compensation, commonly applied to harvested species, by allowing limited harm as long as the probability of the species’ survival or recovery is not reduced appreciably. Demographic compensation requires some density-dependent limits on survival or reproduction in a species’ annual cycle that can be alleviated through incidental take. Using a population model for piping plovers in the Great Plains, we found that when the population is in rapid decline or when there is no density dependence, the probability of quasi-extinction increased linearly with increasing take. However, when the population is near stability and subject to density-dependent survival, there was no relationship between quasi-extinction probability and take rates. We note however, that a brief examination of piping plover demography and annual cycles suggests little room for compensatory capacity. We argue that a population’s capacity for demographic compensation of incidental take should be evaluated when considering incidental allowances because compensation is the only mechanism whereby a population can absorb the negative effects of take without incurring a reduction in the probability of survival in the wild. With many endangered species there is probably little known about density dependence and compensatory capacity. Under these circumstances, using multiple system models (with and without compensation) to predict the population’s response to incidental take and implementing follow-up monitoring to assess species response may be valuable in increasing knowledge and improving future decision making.     

Lynx reintroductions in fragmented landscapes of Germany: Projects with a future or misunderstood wildlife conservation?

Eurasian lynx are slowly recovering in Germany after an absence of about 100 years, and additional reintroduction programs have been launched. However, suitable habitat is patchily distributed in Germany, and whether patches could host a viable population or contribute to the potential spread of lynx is uncertain. We combined demographic scenarios with a spatially explicit population simulation model to evaluate the viability and colonization success of lynx in the different patches, the aim being to conclude guidelines for reintroductions. The spatial basis of our model is a validated habitat model for the lynx in Germany. The dispersal module stems from a calibrated dispersal model, while the demographic module uses plausible published information on the lynx’ life history. The results indicate that (1) a viable population is possible, but that (2) source patches are not interconnected except along the German-Czech border, and that (3) from a demographic viewpoint at least 10 females and 5 males are required for a start that will develop into a viable population with an extinction probability of less than 5% in 50 years. The survival rate of resident adults was the most sensitive parameter, and the best management strategy for the success of reintroduction would be to reduce the mortality of residents in the source patches. Nevertheless, the extremely low probability of connectivity between suitable patches makes most of the reintroduction plans isolated efforts, and they are therefore questionable in the long run. With such a model, the suitability of the single habitat patches can be assessed and the most appropriate management scheme applied. This study shows that simulation models are useful tools for establishing the comparative effectiveness of reintroduction plans aimed at increasing the viability of the species.     

Flawed population viability analysis can result in misleading population assessment: A case study for wolves in Algonquin park, Canada

For many populations and species, population viability analysis (PVA) plays a critical role in developing defensible conservation strategies and recovery plans. Although technical aspects of PVA have been well scrutinized, misapplication of PVA and misinterpretation of its results have received less attention. To illustrate potential hazards of improper use of PVA we reanalyzed data from a recent study on viability of wolves in Algonquin Provincial Park (APP), Ontario, Canada. The original PVA predicted extirpation of wolves from APP and prompted both a ban on wolf harvesting in a 10-16 km buffer zone surrounding the 7571 km² park and an intensive research program to evaluate efficacy of the ban. Our reanalysis showed that limited and imprecise wolf population density and demographic rate estimates, as well as flawed population assessment and reconstruction methods, led to overly pessimistic evaluation of wolf population status in APP. In fact, our analyses suggest that wolves in APP are unlikely to decline significantly over the next 20 years. Further, contrary to earlier conclusions we suggest that rapid wolf population recovery following protection from human exploitation would be likely and readily detectable. These findings highlight the need for adequate data, appropriate methodology, and proper analytical context when conducting PVA. Because the original PVA prompted substantial redirection of staff and financial resources from other significant conservation initiatives, we conclude that improper PVA may undermine execution of effective wildlife management and ultimately provide a disservice to conservation biology.     

Small mammals, habitat patches and PVA models: a field test of model predictive ability

The results are described of comparisons between actual values for patch occupancy for two species of Australian small mammals (Bush Rat Rattus fuscipes and Agile Antechinus Antechinus agilis) determined from field sampling and predictions of patch occupancy made using VORTEX, a generic simulation model for Population Viability Analysis (PVA). The work focussed on a fragmented forest in south-eastern Australia comprised of a network of 39 patches of native eucalypt forest surrounded by extensive stands of exotic softwood Radiata Pine (Pinus radiata) plantation. A range of modelling scenarios were completed in which four broad factors were varied: (1) inter-patch variation in habitat quality; (2) the pattern of inter-patch dispersal; (3) the rate of inter-patch dispersal; and (4) the population sink effects of the Radiata Pine matrix that surrounded the eucalypt patches. Model predictions were made for the total number of animals, the distribution of animal density among patches, the total number of occupied patches, and the probability of patch occupancy. Predictions were then compared with observed values for these same measures based on extensive field surveys of small mammals in the patch system. For most models for the Bush Rat, the predicted relative density of animals per patch correlated well with the values estimated from field surveys. Predictions of patch occupancy were not significantly different from the actual value for the number of occupied patches in half the models tested. The better models explained 10-16% of the log-likelihood of the probability of patch occupancy. While some of the models gave reasonable forecasts of the number of occupied patches, even in these cases, they had only moderate ability to predict which patches were occupied. Field surveys revealed there was no relationship between patch area and population density for the Agile Antechinus—an outcome correctly predicted by only a few models. Five of the 18 scenarios completed for the Agile Antechinus gave predicted numbers of occupied patches not significantly different from the observed number. In each of these five cases, large standard deviations around the mean predicted value meant uncertainty generated by the simulation model limited the predictive power of the PVA. Some of the models gave reasonable predictions for the number of occupied patches, but those models were unable to predict which ones were actually occupied. The results of our study suggest that key processes influencing which specific patches would be occupied were not modelled appropriately. High levels of variability and fecundity drive the population dynamics of the Bush Rat and Agile Antechinus, making the patch system unpredictable and difficult to model accurately. Despite the fact that both the Bush Rat and the Agile Antechinus are two of the most studied mammals in Australia, there are attributes of their biology that are presently poorly understood (which were not included in the VORTEX model), but which could strongly influence patch occupancy. For example, local landscape features may be important determinants of inter-patch movement and habitat utilisation in the patch system. Further empirical studies are needed to explore this aspect of small mammal biology.     

Two-sex population projection of the endemic and dioecious rainforest shrub, Gardenia actinocarpa (Rubiaceae)

Gardenia actinocarpa (Rubiaceae)—an endemic and dioecious shrub-small tree of the lowland rainforest of Cape Tribulation, NE Australia—is evaluated for its long-term population persistence using five-year demographic data from three permanent plots of the species. Two-sex, size-specific Lefkovitch matrices were constructed to evaluate the species’ vital rates and extinction thresholds. One site indicated a sign of significant population growth (annual λ=1.056), while the remaining two gave λs close to unity (λ=1.0102 and 1.007). Overall, λ was greater than 1 (λ=1.0281) indicating that G. actinocarpa will increase in population size with time, though male plants were projected to grow faster in number than female plants. Elasticity analyses indicated that the largest size class (>22 cm diameter) made the greatest contribution to population growth (28%), especially the male plants (23%). Also λ was equally sensitive to survival/stasis (45%) and growth (44%) and least sensitive to fecundity (11%). Stochastic simulations gave 0.505 probability of quasi-extinction to 10% of its initial size in 100 years of population growth projection- suggesting that the endemic G. actinocarpa, despite a positive rate of population increase, is indeed threatened. Sensitivity analyses and numerical simulations showed that to maintain population persistence of G. actinocarpa, efforts should concentrate, following cost-benefit analyses, on increasing seed production coupled with enhancing establishment of individuals in the seedling, juvenile and vegetative phases.