Linking irreplaceable landforms in a self‐organizing landscape to sensitivity of population vital rates for an ecological specialist

Irreplaceable, self‐organizing landforms and the endemic and ecologically specialized biodiversity they support are threatened globally by anthropogenic disturbances. Although the outcome of disrupting landforms is somewhat understood, little information exists that documents population consequences of landform disturbance on endemic biodiversity. Conservation strategies for species dependent upon landforms have been difficult to devise because they require understanding complex feedbacks that create and maintain landforms and the consequences of landform configuration on demography of species. We characterized and quantified links between landform configuration and demography of an ecological specialist, the dunes sagebrush lizard (Sceloporus arenicolus), which occurs only in blowouts (i.e., wind‐blown sandy depressions) of Shinnery oak (Quercus havardii) sand‐dune landforms. We used matrix models to estimate vital rates from a multisite mark‐recapture study of 6 populations occupying landforms with different spatial configurations. Sensitivity and elasticity analyses demonstrated demographic rates among populations varied in sensitivity to different landform configurations. Specifically, significant relationships between blowout shape complexity and vital rate elasticities suggested direct links between S. arenicolus demography and amount of edge in Shinnery oak sand‐dune landforms. These landforms are irreplaceable, based on permanent transition of disturbed areas to alternative grassland ecosystem states. Additionally, complex feedbacks between wind, sand, and Shinnery oak maintain this landform, indicating restoration through land management practices is unlikely. Our findings that S. arenicolus population dynamics depended on landform configuration suggest that failure to consider processes of landform organization and their effects on species’ population dynamics may lead to incorrect inferences about threats to endemic species and ineffective habitat management for threatened or endangered species. As such, successful conservation of these systems and the biodiversity they support must be informed by research linking process‐oriented studies of self‐organized landforms with studies of movement, behavior, and demography of species that dwell in them.     

Translating effects of inbreeding depression on component vital rates to overall population growth in endangered bighorn sheep

Evidence of inbreeding depression is commonly detected from the fitness traits of animals, yet its effects on population growth rates of endangered species are rarely assessed. We examined whether inbreeding depression was affecting Sierra Nevada bighorn sheep (Ovis canadensis sierrae), a subspecies listed as endangered under the U.S. Endangered Species Act. Our objectives were to characterize genetic variation in this subspecies; test whether inbreeding depression affects bighorn sheep vital rates (adult survival and female fecundity); evaluate whether inbreeding depression may limit subspecies recovery; and examine the potential for genetic management to increase population growth rates. Genetic variation in 4 populations of Sierra Nevada bighorn sheep was among the lowest reported for any wild bighorn sheep population, and our results suggest that inbreeding depression has reduced adult female fecundity. Despite this population sizes and growth rates predicted from matrix-based projection models demonstrated that inbreeding depression would not substantially inhibit the recovery of Sierra Nevada bighorn sheep populations in the next approximately 8 bighorn sheep generations (48 years). Furthermore, simulations of genetic rescue within the subspecies did not suggest that such activities would appreciably increase population sizes or growth rates during the period we modeled (10 bighorn sheep generations, 60 years). Only simulations that augmented the Mono Basin population with genetic variation from other subspecies, which is not currently a management option, predicted significant increases in population size. Although we recommend that recovery activities should minimize future losses of genetic variation, genetic effects within these endangered populations-either negative (inbreeding depression) or positive (within subspecies genetic rescue)-appear unlikely to dramatically compromise or stimulate short-term conservation efforts. The distinction between detecting the effects of inbreeding depression on a component vital rate (e.g., fecundity) and the effects of inbreeding depression on population growth underscores the importance of quantifying inbreeding costs relative to population dynamics to effectively manage endangered populations.     

Reliability of Conservation Actions Based on Elasticity Analysis of Matrix Models

Abstract: Matrix population models have entered the mainstream of conservation biology, with analysis of proportional sensitivities (elasticity analysis) of demographic rates becoming important components of conservation decision making. We identify areas where management applications using elasticity analysis potentially conflict with the mathematical basis of the technique, and we use a hypothetical example and three real data sets (Prairie Chicken [   Tympanuchus cupido ], desert tortoise [ Gopherus agassizii ], and killer whale [ Orcinus orca ]) to evaluate the extent to which conservation recommendations based on elasticities might be misleading. First, changes in one demographic rate can change the qualitative ranking of the elasticity values calculated from a population matrix, a result that dampens enthusiasm for ranking conservation actions based solely on which rates have the highest elasticity values. Second, although elasticities often provide accurate predictions of future changes in population growth rate under management perturbations that are large or that affect more than one rate concurrently, concordance frequently fails when different rates vary by different amounts. In particular, when vital rates change to their high or low values observed in nature, predictions of future growth rate based on elasticities of a mean matrix can be misleading, even predicting population increase when the population growth rate actually declines following a perturbation. Elasticity measures will continue to be useful tools for applied ecologists, but they should be interpreted with considerable care. We suggest that studies using analytical elasticity analysis explicitly consider the range of variation possible for different rates and that simulation methods are a useful tool to this end.     

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

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

Developing Metapopulation Connectivity Criteria from Genetic and Habitat Data to Recover the Endangered Mexican Wolf

Restoring connectivity between fragmented populations is an important tool for alleviating genetic threats to endangered species. Yet recovery plans typically lack quantitative criteria for ensuring such population connectivity. We demonstrate how models that integrate habitat, genetic, and demographic data can be used to develop connectivity criteria for the endangered Mexican wolf (Canis lupus baileyi), which is currently being restored to the wild from a captive population descended from 7 founders. We used population viability analysis that incorporated pedigree data to evaluate the relation between connectivity and persistence for a restored Mexican wolf metapopulation of 3 populations of equal size. Decreasing dispersal rates greatly increased extinction risk for small populations (<150–200), especially as dispersal rates dropped below 0.5 genetically effective migrants per generation. We compared observed migration rates in the Northern Rocky Mountains (NRM) wolf metapopulation to 2 habitat‐based effective distance metrics, least‐cost and resistance distance. We then used effective distance between potential primary core populations in a restored Mexican wolf metapopulation to evaluate potential dispersal rates. Although potential connectivity was lower in the Mexican wolf versus the NRM wolf metapopulation, a connectivity rate of >0.5 genetically effective migrants per generation may be achievable via natural dispersal under current landscape conditions. When sufficient data are available, these methods allow planners to move beyond general aspirational connectivity goals or rules of thumb to develop objective and measurable connectivity criteria that more effectively support species recovery. The shift from simple connectivity rules of thumb to species‐specific analyses parallels the previous shift from general minimum‐viable‐population thresholds to detailed viability modeling in endangered species recovery planning. Desarrollo de Criterios de Conectividad Metapoblacional a Partir de Datos Genéticos y de Hábitat para Recuperar al Lobo Mexicano en Peligro de Extinción     

Population size affects vital rates but not population growth rate of a perennial plant

Negative effects of habitat fragmentation on individual performance have been widely documented, but relatively little is known about how simultaneous effects on multiple vital rates translate into effects on population viability in long-lived species. In this study, we examined relationships between population size, individual growth, survival and reproduction, and population growth rate in the perennial plant Phyteuma spicatum. Population size positively affected the growth of seedlings, the survival of juveniles, the proportion of adults flowering, and potential seed production. Analyses with integral projection models, however, showed no relationship between population size and population growth rate. This was due to the fact that herbivores and pathogens eliminated the relationship between population size and seed production, and that population growth rate was not sensitive to changes in the vital rates that varied with population size. We conclude that effects of population size on vital rates must not translate into effects on population growth rate, and that populations of long-lived organisms may partly be able to buffer negative effects of small population size on vital rates that have a relatively small influence on population growth rate. Our study illustrates that we need to be cautious when assessing the consequences of habitat fragmentation for population viability based on effects on only one or a few vital rates.     

Population inertia and its sensitivity to changes in vital rates and population structure

Because the (st)age structure of a population may rarely be stable, studies of transient population dynamics and population momentum are becoming ever more popular. Yet, studies of "population momentum" are restricted in the sense that they describe the inertia of population size resulting from a demographic transition to the stationary population growth rate. Although rarely mentioned, inertia in population size is a general phenomenon and can be produced by any demographic transition or perturbation. Because population size is of central importance in demography, conservation, and management, formulas relating the sensitivity of population inertia to changes in underlying vital rates and population structure could provide much-needed insight into the dynamics of populations with unstable (st)age structure. Here, we derive such formulas, which are readily computable, and provide examples of their potential use in studies of life history and applied arenas of population study.     

Possible Extinction Vortex for a Population of Iberian Lynx on the Verge of Extirpation

Abstract: Theory suggests that demographic and genetic traits deteriorate (i.e., fitness and genetic diversity decrease) when populations become small, and that such deterioration could precipitate positive feedback loops called extinction vortices. We examined whether demographic attributes and genetic traits have changed over time in one of the 2 remaining small populations of the highly endangered Iberian lynx (Lynx pardinus) in Doñana, Spain. From 1983 to 2008, we recorded nontraumatic mortality rates, litter size, offspring survival, age at territory acquisition, and sex ratio. We combined these demographic attributes with measures of inbreeding and genetic diversity at neutral loci (microsatellites) and genes subjected to selection (major histocompatibility complex). Data on demographic traits were obtained through capture and radio tracking, checking dens during breeding, track surveys, and camera trapping. For genetic analyses, we obtained blood or tissue samples from captured or necropsied individuals or from museum specimens. Over time a female‐biased sex ratio developed, age of territory acquisition decreased, mean litter size decreased, and rates of nontraumatic mortality increased, but there were no significant changes in overall mortality rates, standardized individual heterozygosity declined steadily, and allelic diversity of exon 2 of class II major histocompatibility complex DRB genes remained constant (2 allelic variants present in all individuals analyzed). Changes in sex ratio and age of territory acquisition may have resulted from demographic stochasticity, whereas changes in litter size and nontraumatic mortality may be related to observed increases in inbreeding. Concomitant deterioration of both demographic attributes and genetic traits is consistent with an extinction vortex. The co‐occurrence, with or without interaction, of demographic and genetic deterioration may explain the lack of success of conservation efforts with the Doñana population of Iberian lynx.     

Reexamining the Minimum Viable Population Concept for Long‐Lived Species

For decades conservation biologists have proposed general rules of thumb for minimum viable population size (MVP); typically, they range from hundreds to thousands of individuals. These rules have shifted conservation resources away from small and fragmented populations. We examined whether iteroparous, long‐lived species might constitute an exception to general MVP guidelines. On the basis of results from a 10‐year capture‐recapture study in eastern New York (U.S.A.), we developed a comprehensive demographic model for the globally threatened bog turtle (Glyptemys muhlenbergii), which is designated as endangered by the IUCN in 2011. We assessed population viability across a wide range of initial abundances and carrying capacities. Not accounting for inbreeding, our results suggest that bog turtle colonies with as few as 15 breeding females have >90% probability of persisting for >100 years, provided vital rates and environmental variance remain at currently estimated levels. On the basis of our results, we suggest that MVP thresholds may be 1–2 orders of magnitude too high for many long‐lived organisms. Consequently, protection of small and fragmented populations may constitute a viable conservation option for such species, especially in a regional or metapopulation context. Reexaminando el Concepto de Población Mínima Viable para Especies Longevas Resumen     

A case study of bats and white‐nose syndrome demonstrating how to model population viability with evolutionary effects

Ecological factors generally affect population viability on rapid time scales. Traditional population viability analyses (PVA) therefore focus on alleviating ecological pressures, discounting potential evolutionary impacts on individual phenotypes. Recent studies of evolutionary rescue (ER) focus on cases in which severe, environmentally induced population bottlenecks trigger a rapid evolutionary response that can potentially reverse demographic threats. ER models have focused on shifting genetics and resulting population recovery, but no one has explored how to incorporate those findings into PVA. We integrated ER into PVA to identify the critical decision interval for evolutionary rescue (DIER) under which targeted conservation action should be applied to buffer populations undergoing ER against extinction from stochastic events and to determine the most appropriate vital rate to target to promote population recovery. We applied this model to little brown bats (Myotis lucifugus) affected by white‐nose syndrome (WNS), a fungal disease causing massive declines in several North American bat populations. Under the ER scenario, the model predicted that the DIER period for little brown bats was within 11 years of initial WNS emergence, after which they stabilized at a positive growth rate (λ = 1.05). By comparing our model results with population trajectories of multiple infected hibernacula across the WNS range, we concluded that ER is a potential explanation of observed little brown bat population trajectories across multiple hibernacula within the affected range. Our approach provides a tool that can be used by all managers to provide testable hypotheses regarding the occurrence of ER in declining populations, suggest empirical studies to better parameterize the population genetics and conservation‐relevant vital rates, and identify the DIER period during which management strategies will be most effective for species conservation.     

Bridging gaps in demographic analysis with phylogenetic imputation

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

Demographic Effects of Harvesting Epiphytic Bromeliads and an Alternative Approach to Collection

Abstract: Hundreds of epiphytic bromeliads species are harvested from the wild for trade and for cultural uses, but little is known about the effects of this harvest. We assessed the potential demographic effects of harvesting from the wild on 2 epiphytic bromeliads: Tillandsia macdougallii, an atmospheric bromeliad (adsorbs water and nutrients directly from the atmosphere), and T. violaceae, a tank bromeliad (accumulates water and organic material between its leaves). We also examined an alternative to harvesting bromeliads from trees—the collection of fallen bromeliads from the forest floor. We censused populations of T. macdougallii each year from 2005 to 2010 and of T. violaceae from 2005 to 2008, in Oaxaca, Mexico. We also measured monthly fall rates of bromeliads over 1 year and monitored the survival of fallen bromeliads on the forest floor. The tank bromeliad had significantly higher rates of survival, reproduction, and stochastic population growth rates (λ_s) than the atmospheric bromeliad, but λ_s for both species were <1, which suggests that the populations will decline even without harvest. Elasticity patterns differed between species, but in both, survival of large individuals had high elasticity values. No fallen bromeliads survived more than 1.5 years on the forest floor and the rate of bromeliad fall was comparable to current harvest rates. Low rates of population growth recorded for the species we studied and other epiphytic bromeliads and high elasticity values for the vital rates that were most affected by harvest suggest that commercial harvesting in the wild of these species is not sustainable. We propose the collection of fallen bromeliads as an ecologically and, potentially, economically viable alternative.     

Importance of Habitat Quality and Landscape Connectivity for the Persistence of Endangered Natterjack Toads

Abstract: The natterjack toad (Bufo calamita) is endangered in several parts of its distribution, including Belgium, where it occurs mainly in artificial habitats. We parameterized a general model for natterjack population viability analysis (PVA) and tested its sensitivity to changes in the values of basic parameters. Then we assessed the relative efficiency of various conservation measures in 2 situations: a small isolated population and a system of 4 populations connected by rare dispersal movements. We based the population viability analysis on a stage‐structured model of natterjack population dynamics. We parameterized the model in the RAMAS GIS platform with vital rates obtained from our own field experience and from published studies. Simulated natterjack populations were highly sensitive to habitat quality (particularly pond drying), to dispersal from surrounding local populations, and to a lesser extent to values of fecundity and survival of terrestrial stages. Population trajectories were nearly insensitive to initial abundances, carrying capacities, and the frequency of extreme climatic conditions. The simulations showed that in habitats with highly ephemeral ponds, where premetamorphosis mortality was high, natterjack populations nearly always had a very high extinction risk. We also illustrated how low dispersal rates (<1 dispersing individual/generation) efficiently rescued declining local populations. Such source‐sink dynamics demonstrate that the identification and management of source populations should be a high priority.     

Effects of life history and reproduction on recruitment time lags in reintroductions of rare plants

Reintroductions are important components of conservation and recovery programs for rare plant species, but their long-term success rates are poorly understood. Previous reviews of plant reintroductions focused on short-term (e.g., ≤3 years) survival and flowering of founder individuals rather than on benchmarks of intergenerational persistence, such as seedling recruitment. However, short-term metrics may obscure outcomes because the unique demographic properties of reintroductions, including small size and unstable stage structure, could create lags in population growth. We used time-to-event analysis on a database of unusually well-monitored and long-term (4–28 years) reintroductions of 27 rare plant species to test whether life-history traits and population characteristics of reintroductions create time-lagged responses in seedling recruitment (i.e., recruitment time lags [RTLs]), an important benchmark of success and indicator of persistence in reintroduced populations. Recruitment time lags were highly variable among reintroductions, ranging from <1 to 17 years after installation. Recruitment patterns matched predictions from life-history theory with short-lived species (fast species) exhibiting consistently shorter and less variable RTLs than long-lived species (slow species). Long RTLs occurred in long-lived herbs, especially in grasslands, whereas short RTLs occurred in short-lived subtropical woody plants and annual herbs. Across plant life histories, as reproductive adult abundance increased, RTLs decreased. Highly variable RTLs were observed in species with multiple reintroduction events, suggesting local processes are just as important as life-history strategy in determining reintroduction outcomes. Time lags in restoration outcomes highlight the need to scale success benchmarks in reintroduction monitoring programs with plant life-history strategies and the unique demographic properties of restored populations. Drawing conclusions on the long-term success of plant reintroduction programs is premature given that demographic processes in species with slow life-histories take decades to unfold.     

Ability of Matrix Models to Explain the Past and Predict the Future of Plant Populations

Uncertainty associated with ecological forecasts has long been recognized, but forecast accuracy is rarely quantified. We evaluated how well data on 82 populations of 20 species of plants spanning 3 continents explained and predicted plant population dynamics. We parameterized stage‐based matrix models with demographic data from individually marked plants and determined how well these models forecast population sizes observed at least 5 years into the future. Simple demographic models forecasted population dynamics poorly; only 40% of observed population sizes fell within our forecasts’ 95% confidence limits. However, these models explained population dynamics during the years in which data were collected; observed changes in population size during the data‐collection period were strongly positively correlated with population growth rate. Thus, these models are at least a sound way to quantify population status. Poor forecasts were not associated with the number of individual plants or years of data. We tested whether vital rates were density dependent and found both positive and negative density dependence. However, density dependence was not associated with forecast error. Forecast error was significantly associated with environmental differences between the data collection and forecast periods. To forecast population fates, more detailed models, such as those that project how environments are likely to change and how these changes will affect population dynamics, may be needed. Such detailed models are not always feasible. Thus, it may be wiser to make risk‐averse decisions than to expect precise forecasts from models. Habilidad de los Modelos Matriciales para Explicar el Pasado y Predecir el Futuro de las Poblaciones de Plantas     

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.     

Effects of Harvest of Nontimber Forest Products and Ecological Differences between Sites on the Demography of African Mahogany

Abstract: The demographic impacts of harvesting nontimber forest products (NTFP) have been increasingly studied because of reports of potentially unsustainable harvest. Nevertheless, our understanding of how plant demographic response to harvest is altered by variation in ecological conditions, which is critical for developing realistic sustainable‐use plans, is limited. We built matrix population models to test whether and how variation in ecological conditions affects population responses to harvest. In particular, we examined the effect of bark and foliage harvest on the demography of populations of African mahogany (Khaya senegalensis) in two contrasting ecological regions of Benin, West Africa. K. senegalensis bark and foliage harvest significantly reduced its stochastic population growth rates, but ecological differences between regions had a greater effect on population growth rates than did harvest. The effect of harvest on population growth rates (Δλ) was slightly stronger in the moist than in the drier region. Life‐table response experiments revealed that the mechanism by which harvesting reduced λ differed between ecological regions. Lowered stasis (persistence) of larger life stages lead to a reduction in λ in the drier region, whereas lowered growth of all life stages lowered λ in moist region. Potential strategies to increase population growth rates should include decreasing the proportion of individuals harvested, promoting harvester‐owned plantations of African mahogany, and increasing survival and growth by promoting no‐fire zones in gallery forests. Our results show how population responses to harvest of NTFP may be altered by ecological differences across sites and emphasize the importance of monitoring populations over the climatic range in which they occur to develop more realistic recommendations for conservation.     

Testing sustainability by prospective and retrospective demographic analyses: evaluation for palm leaf harvest.

Harvesting nontimber forest products (NTFPs) is a major economic activity in tropical forests. As many NTFPs are overexploited, sustainability analyses are required to set harvest guidelines. Here we introduce and apply a new approach to evaluating sustainability, which combines prospective (elasticity) and retrospective (Life Table Response Experiments [LTRE]) demographic analyses of matrix population models. We relate the elasticity of vital rates (representing their importance for population growth rate, lamda) to their contribution to harvest-induced change in lamda ("LTRE contribution"). When high-elasticity vital rates have a low LTRE contribution, exploitation is potentially sustainable as negative effects for population growth are buffered. If the reverse is found, there is little scope for sustainability because crucial vital rates are affected. Our approach is less sensitive to chance fluctuations than the commonly used sustainability criterion of lamda = 1.0, as it does not depend on the absolute value of lamda. We applied this analysis to Geonoma deversa, a clustered forest understory palm. We studied three experimentally defoliated and control populations in a Bolivian rainforest during two years. Cutting all leaves of large ramets did not change mortality but strongly affected growth and reproduction. In spite of severe changes in some vital rates, population growth rate was not significantly reduced after defoliation. A literature review revealed that six other understory palms species responded very similarly to defoliation. The combination of LTRE contributions and elasticity analyses showed that low-elasticity vital rates were mainly responsible for the defoliation-induced change in lamda for Geonoma deversa. For two other understory palms (Astrocaryum mexicanum and Chamaedorea radicalis) new demographic analyses yielded very similar results. For Geonoma, the LTRE contribution-elasticity relation strongly changed when we mimicked harvest damage. Adding 5% mortality to defoliated palms caused stronger change in lamda, mainly due to changes in a high-elasticity vital rate (survival). Therefore, harvest practices that involve stem killing are clearly unsustainable. Our results show that commercial leaf cutting in Geonoma deversa is potentially sustainable, and that this is likely the case for understory palms in general. Our approach to evaluating harvest sustainability can be applied to other NTFPs.     

Conserving Slow-Growing, Long-Lived Tree Species: Input from the Demography of a Rare Understory Conifer, Taxus floridana

Abstract:  Although land preservation and promotion of successful regeneration are important conservation actions, their ability to increase population growth rates of slow-growing, long-lived trees is limited. We investigated the demography of Taxus floridana Nutt., a rare understory conifer, in three populations in different ravine forests spanning its entire geographic range along the Apalachicola River Bluffs in northern Florida (U.S.A.). We examined spatial and temporal patterns in demographic parameters and projected population growth rates by using four years of data on the recruitment and survival of seedlings and established stems, and on diameter growth from cross-sections of dead stems. All populations experienced a roughly 10-fold increase in seedling recruitment in 1996 compared with other years. The fates of seedlings and stems between 8 and 16 mm differed among populations. The fates of stems in two other size classes (the 2- to 4-mm class and the 4- to 8-mm class) differed among both populations and years. Individual stems in all populations exhibited similarly slow growth rates. Stochastic matrix models projected declines in all populations. Stochastic matrix analysis revealed the high elasticity of a measure of stochastic population growth rate to perturbations in the stasis of large reproductive stems for all populations. Additional analyses also indicated that occasional episodes of high recruitment do not greatly affect population growth rates. Conservation efforts directed at long-lived, slow-growing rare plants like Taxus floridana should both protect established reproductive individuals and further enhance survival of individuals in other life-history stages, such as juveniles, that often do not appear to contribute greatly to population growth rates.     

Use of Integrated Modeling to Enhance Estimates of Population Dynamics Obtained from Limited Data

Abstract:  Demographic data of rare and endangered species are often too sparse to estimate vital rates and population size with sufficient precision for understanding population growth and decline. Yet, the combination of different sources of demographic data into one statistical model holds promise. We applied Bayesian integrated population modeling to demographic data from a colony of the endangered greater horseshoe bats (Rhinolophus ferrumequinum) . Available data were the number of subadults and adults emerging from the colony roost at dusk, the number of newborns from 1991 to 2005, and recapture data of subadults and adults from 2004 and 2005. Survival rates did not differ between sexes, and demographic rates remained constant across time. The greater horseshoe bat is a long-lived species with high survival rates (first year: 0.49 [SD 0.06]; adults: 0.91 [SD 0.02]) and low fecundity (0.74 [SD 0.12]). The yearly average population growth was 4.4% (SD 0.1%) and there were 92 (SD 10) adults in the colony in year 2005. Had we analyzed each data set separately, we would not have been able to estimate fecundity, the estimates of survival would have been less precise, and the estimate of population growth biased. Our results demonstrate that integrated models are suitable for obtaining crucial demographic information from limited data.