Modeling the impacts of hunting on the population dynamics of red howler monkeys (Alouatta seniculus)

Overexploitation of wildlife populations occurs across the humid tropics and is a significant threat to the long-term survival of large-bodied primates. To investigate the impacts of hunting on primates and ways to mitigate them, we developed a spatially explicit, individual-based model for a landscape that included hunted and un-hunted areas. We used the large-bodied neotropical red howler monkey (Alouatta seniculus) as our case study species because its life history characteristics make it vulnerable to hunting. We modeled the influence of different rates of harvest and proportions of landscape dedicated to un-hunted reserves on population persistence, population size, social dynamics, and hunting yields of red howler monkeys. In most scenarios, the un-hunted populations maintained a constant density regardless of hunting pressure elsewhere, and allowed the overall population to persist. Therefore, the overall population was quite resilient to extinction; only in scenarios without any un-hunted areas did the population go extinct. However, the total and hunted populations did experience large declines over 100 years under moderate and high hunting pressure. In addition, when reserve area decreased, population losses and losses per unit area increased disproportionately. Furthermore, hunting disrupted the social structure of troops. The number of male turnovers and infanticides increased in hunted populations, while birth rates decreased and exacerbated population losses due to hunting. Finally, our results indicated that when more than 55% of the landscape was harvested at high (30%) rates, hunting yields, as measured by kilograms of biomass, were less than those obtained from moderate harvest rates. Additionally, hunting yields, expressed as the number of individuals hunted/year/km², increased in proximity to un-hunted areas, and suggested that dispersal from un-hunted areas may have contributed to hunting sustainability. These results indicate that un-hunted areas serve to enhance hunting yields, population size, and population persistence in hunted landscapes. Therefore, spatial regulation of hunting via a reserve system may be an effective management strategy for sustainable hunting, and we recommend it because it may also be more feasible to implement than harvest quotas or restrictions on season length.     

Interactive Effects of Harvest and Deer Herbivory on the Population Dynamics of American Ginseng

Abstract:  Few demographic models for any species consider the role of multiple, interacting ecological threats. Many forest herbs are heavily browsed by white-tailed deer ( Odocoileus virginianus ) and a number of these are also harvested for the medicinal, floral, or horticultural trades. Previous studies of the viability of American ginseng ( Panax quinquefolius ) have separately examined the effects of harvesting and deer herbivory. We followed individually marked ginseng plants in 6 populations for 8 years and documented deer browse levels, conducted helicopter surveys to estimate the deer herd size, and documented 2 ginseng harvests. We used this long-term data set to develop a stochastic demographic model that quantified the separate and interactive role of these threats to ginseng viability. Although harvesting and deer herbivory negatively affected ginseng population growth, their effects were not additive. Deer herbivory negatively affected population growth in the absence but not in the presence of harvesting. Life table response experiments revealed that in the presence of harvesting, deer herbivory had some positive effects on vital rates because browsed plants were less apparent to harvesters. Ginseng populations that were harvested responsibly (i.e., planting seeds from harvested individuals) had higher growth rates than those that were harvested irresponsibly. We concluded that both deer populations and harvesting must be managed to ensure sustainable populations of American ginseng. Our findings underscore the importance of long-term monitoring to assess threats to viability and the need for a broad ecological understanding of the complexity of ecosystem management.     

Adjusting age and stage distributions for misclassification errors

Ecologists often use samples from the age or stage structure of a population to make inferences about population-level processes and to parameterize matrix models. Typically, researchers make a simplifying assumption that age and stage classes are determined without error, when in fact some level of misclassification often can be expected. If unaccounted for, misclassification will lead to overly optimistic levels of precision and can cause biased estimates of age or stage structure. Although several studies have used information from known-age individuals to quantify errors in age or stage distribution, the problem of estimating the age or stage structure in face of such errors has received comparably little attention. In this paper, we describe a general statistical framework for estimating the true stage distribution of a sample when misclassification rates can be estimated. The estimation process requires auxiliary information on misclassification rates, such as data from individuals of known age. We analyze age-structured harvest records from black bears in Pennsylvania to illustrate how incorporating misclassification errors leads to changes in point estimates and provides a measure of precision.     

Habitat and population modelling of roe deer using an interactive geographic information system

Management of German roe deer (Capreolus capreolus) populations is a challenge for wildlife managers and foresters because population densities are difficult to estimate in forests and forest regeneration can be negatively affected when roe deer density is high. We describe a model to determine deer population densities compatible with forest management goals, and to assess harvest rates necessary to maintain desired deer densities. A geographic information system (GIS) was used to model wildlife habitat and population dynamics over time. Our model interactively incorporates knowledge of field biologists and foresters via a graphical user interface (GUI). Calibration of the model with deer damage maps allowed us to evaluate density dependence of a roe deer population. Incorporation of local knowledge into temporally dynamic and spatial models increases understanding of population dynamics and improves wildlife management.     

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.     

Calculating Ne and Ne/N in age-structured populations: a hybrid Felsenstein-Hill approach

The concept of effective population size (Ne) was developed under a discrete-generation model, but most species have overlapping generations. In the early 1970s, J. Felsenstein and W. G. Hill independently developed methods for calculating Ne in age-structured populations; the two approaches produce the same answer under certain conditions and have contrasting advantages and disadvantages. Here, we describe a hybrid approach that combines useful features of both. Like Felsenstein's model, the new method is based on age-specific survival and fertility rates and therefore can be directly applied to any species for which life table data are available. Like Hill, we relax the restrictive assumption in Felsenstein's model regarding random variance in reproductive success, which allows more general application. The basic principle underlying the new method is that age structure stratifies a population into winners and losers in the game of life: individuals that live longer have more opportunities to reproduce and therefore have a higher mean lifetime reproductive success. This creates different classes of individuals within the population, and grouping individuals by age at death provides a simple means of calculating lifetime variance in reproductive success of a newborn cohort. The new method has the following features: (1) it can accommodate unequal sex ratio and sex-specific vital rates and overdispersed variance in reproductive success; (2) it can calculate effective size in species that change sex during their lifetime; (3) it can calculate Ne and the ratio Ne/N based on various ways of defining N; (4) it allows one to explore the relationship between Ne and the effective number of breeders per year (Nb), which is a quantity that genetic estimators of contemporary Ne commonly provide information about; and (5) it is implemented in freely available software (AgeNe).     

A spatially explicit empirical model of structural development processes in natural forests based on climate and topography

Stand structure develops with stand age. Old-growth forests with well-developed stand structure support many species. However, development rates of stand structure likely vary with climate and topography. We modeled structural development of 4 key stand variables and a composite old-growth index as functions of climatic and topographic covariates. We used a hierarchical Bayesian method for analysis of extensive snap-shot National Forest Inventory (NFI) data in Japan (n = 9244) to account for differences in stand age. Development rates of structural variables and the old-growth index exhibited curvilinear responses to environmental covariates. Flat sites were characterized by high rates of structural development. Approximately 150 years were generally required to attain high values (approximately 0.8) of the old-growth index. However, the predicted age to achieve specific values varied depending on environmental conditions. Spatial predictions highlighted regional variation in potential structural development rates. For example, sometimes there were differences of >100 years among sites, even in the same catchment, in attainment of a medium index value (0.5) after timber harvesting. The NFI data suggested that natural forests, especially old natural forests (>150 years), remain generally on unproductive ridges, steep slopes, or areas with low temperature and deep snow, where many structural variables show slow development rates. We suggest that maintenance and restoration of old natural forests on flat sites should be prioritized for conservation due to the likely rapid development of stand structure, although remaining natural forests on low-productivity sites are still important and should be protected.     

Ungulate population dynamics and optimization models

Optimal harvesting strategies for an ungulate population are estimated using stochastic dynamic programming. Data on the Llano Basin white-tailed deer (Odocoileus virginianus) population were used to construct a 2-variable population dynamics model. The model provided the basis for estimating optimal harvesting strategies as a feedback function of the current values of the state variables (prefawning older deer and juveniles). Optimal harvest strategies were insensitive to assumptions about the probability distributions of the stochastic variable (rainfall). The response of the population components to harvesting and the returns obtained from applying optimal strategies were explored through simulation. Mean annual harvest is about 15% of the population. Simplified harvesting strategies based on age-ratios as well as a simplified version based on optimal strategies—but assuming persisting equilibrium juvenile deer density—were compared to optimal strategies through examining values of information. Simplified harvesting strategies lead to a lower harvest over a 50-year simulation period.     

The Influence of Model Structure on Conclusions about the Viability and Harvesting of Serengeti Wildebeest

We investigate how the viability and harvestability predicted by population models are affected by details of model construction. Based on this analysis we discuss some of the pitfalls associated with the use of classical statistical techniques for resolving the uncertainties associated with modeling population dynamics. The management of the Serengeti wildebeest (Connochaetes taurinus) is used as a case study. We fitted a collection of age-structured and unstructured models to a common set of available data and compared model predictions in terms of wildebeest viability and harvest. Models that depicted demographic processes in strikingly different ways fitted the data equally well. However, upon further analysis it became clear that models that fit the data equally well could nonetheless have very different management implications. In general, model structure had a much larger effect on viability analysis (e.g., time to collapse) than on optimal harvest analysis (e.g., harvest rate that maximizes harvest). Some modeling decisions, such as including age-dependent fertility rates, did not affect management predictions, but others had a strong effect (e.g., choice of model structure). Because several suitable models of comparable complexity fitted the data equally well, traditional model selection methods based on the parsimony principle were not practical for judging the value of alternative models. Our results stress the need to implement analytical frameworks for population management that explicitly consider the uncertainty about the behavior of natural systems.     

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.     

Emerging prion disease drives host selection in a wildlife population

Infectious diseases are increasingly recognized as an important force driving population dynamics, conservation biology, and natural selection in wildlife populations. Infectious agents have been implicated in the decline of small or endangered populations and may act to constrain population size, distribution, growth rates, or migration patterns. Further, diseases may provide selective pressures that shape the genetic diversity of populations or species. Thus, understanding disease dynamics and selective pressures from pathogens is crucial to understanding population processes, managing wildlife diseases, and conserving biological diversity. There is ample evidence that variation in the prion protein gene (PRNP) impacts host susceptibility to prion diseases. Still, little is known about how genetic differences might influence natural selection within wildlife populations. Here we link genetic variation with differential susceptibility of white-tailed deer to chronic wasting disease (CWD), with implications for fitness and disease-driven genetic selection. We developed a single nucleotide polymorphism (SNP) assay to efficiently genotype deer at the locus of interest (in the 96th codon of the PRNP gene). Then, using a Bayesian modeling approach, we found that the more susceptible genotype had over four times greater risk of CWD infection; and, once infected, deer with the resistant genotype survived 49% longer (8.25 more months). We used these epidemiological parameters in a multi-stage population matrix model to evaluate relative fitness based on genotype-specific population growth rates. The differences in disease infection and mortality rates allowed genetically resistant deer to achieve higher population growth and obtain a long-term fitness advantage, which translated into a selection coefficient of over 1% favoring the CWD-resistant genotype. This selective pressure suggests that the resistant allele could become dominant in the population within an evolutionarily short time frame. Our work provides a rare example of a quantifiable disease-driven selection process in a wildlife population, demonstrating the potential for infectious diseases to alter host populations. This will have direct bearing on the epidemiology, dynamics, and future trends in CWD transmission and spread. Understanding genotype-specific epidemiology will improve predictive models and inform management strategies for CWD-affected cervid populations.     

Growth and population dynamic model for the non-zooxanthellate temperate solitary coral Leptopsammia pruvoti (Scleractinia, Dendrophylliidae)

In corals where complex life history processes decoupling age from size (e.g., fragmentation, fusion, partial colony mortality) are rare or clearly detectable, individual age may be determined from size, and age-based growth and population dynamic models may be applied. One example is the solitary Mediterranean coral Leptopsammia pruvoti Lacaze-Duthiers 1897, whose population size and structure, and growth rates were determined at Calafuria (43°28′N and 10°20′E), Ligurian Sea, from December 2007 to June 2009. Growth rate decreased with increasing size. The growth curve derived from field measurements confirmed the one obtained by growth band analysis. The frequency of individuals decreased exponentially with age, indicating a steady state population. Turnover time was 2.3 years. Maximum life span was 13 years. Most reproductive output was from intermediate age classes (6 years), while older individuals (>7 years), although having higher fecundity, were rare and accounted for a minority of population reproductive output. In comparison with other solitary dendrophylliids, L. pruvoti life strategy was characterized by a reproduction with r-strategy correlates (high fecundity, short embryo incubation, small planula size, fast achievement of sexual maturity), and a rate of demographic renewal occurring halfway along the r–K continuum (intermediate turnover time and maximum longevity).     

Optimal control of Atlantic population Canada geese

Management of Canada geese (Branta canadensis) can be a balance between providing sustained harvest opportunity while not allowing populations to become overabundant and cause damage. In this paper, we focus on the Atlantic population of Canada geese and use stochastic dynamic programming to determine the optimal harvest strategy over a range of plausible models for population dynamics. There is evidence to suggest that the population exhibits significant age structure, and it is possible to reconstruct age structure from surveys. Consequently the harvest strategy is a function of the age composition, as well as the abundance, of the population. The objective is to maximize harvest while maintaining the number of breeding adults in the population between specified upper and lower limits. In addition, the total harvest capacity is limited and there is uncertainty about the strength of density-dependence. We find that under a density-independent model, harvest is maximized by maintaining the breeding population at the highest acceptable abundance. However if harvest capacity is limited, then the optimal long-term breeding population size is lower than the highest acceptable level, to reduce the risk of the population growing to an unacceptably large size. Under the proposed density-dependent model, harvest is maximized by maintaining the breeding population at an intermediate level between the bounds on acceptable population size; limits to harvest capacity have little effect on the optimal long-term population size. It is clear that the strength of density-dependence and constraints on harvest significantly affect the optimal harvest strategy for this population. Model discrimination might be achieved in the long term, while continuing to meet management goals, by adopting an adaptive management strategy.     

Conserving spawning stocks through harvest slot limits and no-take protected areas

The key to the conservation of harvested species is the maintenance of reproductive success. Yet for many marine species large, old, individuals are targeted despite their disproportionate contribution to reproduction. We hypothesized that a combination of no-take marine protected areas (MPAs) and harvest slot limits (maximum and minimum size limits) would result in the conservation of large spawning individuals under heavy harvest. We tested this approach under different harvest intensities with a 2-sex, stage-structured metapopulation model for the Caribbean spiny lobster (Panulirus argus). P. argus is intensively harvested in the Caribbean, and in many localities large, mature individuals no longer exist. No-take MPAs and harvest slot limits combined, rebuilt and maintained large mature individuals even under high harvest pressure. The most conservative model (a 30% MPA and harvest slot limit of 75–105 mm) increased spawner abundance by 5.53E¹² compared with the fishing status quo at the end of 30 years. Spawning stock abundance also increased by 2.76–9.56E¹² individuals at a high harvest intensity over 30 years with MPAs alone. Our results demonstrate the potential of MPAs and harvest slot limits for the conservation of large breeding individuals in some marine and freshwater environments. Decisions on which management strategy best suits a fishery, however, requires balancing what is ecologically desirable with what is economically and socially feasible.     

Demographic Side Effects of Selective Hunting in Ungulates and Carnivores

Abstract:  Selective harvesting regimes are often implemented because age and sex classes contribute differently to population dynamics and hunters show preferences associated with body size and trophy value. We reviewed the literature on how such cropping regimes affect the demography of the remaining population (here termed demographic side effects ). First, we examined the implications of removing a large proportion of a specific age or sex class. Such harvesting strategies often bias the population sex ratio toward females and reduce the mean age of males, which may consequently delay birth dates, reduce birth synchrony, delay body mass development, and alter offspring sex ratios. Second, we reviewed the side effects associated with the selective removal of relatively few specific individuals, often large trophy males. Such selective harvesting can destabilize social structures and the dominance hierarchy and may cause loss of social knowledge, sexually selected infanticide, habitat changes among reproductive females, and changes in offspring sex ratio. A common feature of many of the reported mechanisms is that they ultimately depress recruitment and in some extreme cases even cause total reproductive collapse. These effects could act additively and destabilize the dynamics of populations, thus having a stronger effect on population growth rate than first anticipated. Although more experimental than observational studies reported demographic side effects, we argue that this may reflect the quite subtle mechanisms involved, which are unlikely to be detected in observational studies without rigorous monitoring regimes. We call for more detailed studies of hunted populations with marked individuals that address how the expression of these effects varies across mating systems, habitats, and with population density. Theoretical models investigating how strongly these effects influence population growth rates are also required.     

Patterns of body mass senescence and selective disappearance differ among three species of free-living ungulates

Declines in survival and reproduction with age are prevalent in wild vertebrates, but we know little about longitudinal changes in behavioral, morphological, or physiological variables that may explain these demographic declines. We compared age-related variation in body mass of adult females in three free-living ungulate populations that have been the focus of long-term, individual-based research: bighorn sheep (Ovis canadensis) at Ram Mountain, Canada; roe deer (Capreolus capreolus) at Trois Fontaines, France; and Soay sheep (Ovis aries) on St. Kilda, Scotland. We use two recently proposed approaches to separate contributions to age-dependent variation at the population level from within-individual changes and between-individual selective disappearance. Selective disappearance of light individuals in all three populations was most evident at the youngest and oldest ages. In later adulthood, bighorn sheep and roe deer showed a continuous decline in body mass that accelerated with age while Soay sheep showed a precipitous decrease in mass in the two years preceding death. Our results highlight the importance of mass loss in explaining within-individual demographic declines in later adulthood in natural populations. They also reveal that the pattern of senescence, and potentially also the processes underlying demographic declines in late life, can differ markedly across related species with similar life histories.     

Selective harvesting and habitat loss produce long-term life history changes in a mouflon population.

We examined the long-term effects (28 years) of habitat loss and phenotype-based selective harvest on body mass, horn size, and horn shape of mouflon (Ovis gmelini musimon) in southern France. This population has experienced habitat deterioration (loss of 50.8% of open area) since its introduction in 1956 and unrestricted selective hunting of the largest horned males since 1973. Both processes are predicted to lead to a decrease in phenotype quality by decreasing habitat quality and by reducing the reproductive contribution of individuals carrying traits that are targeted by hunters. Body mass and body size of both sexes and horn measurements of males markedly decreased (by 3.4-38.3%) in all age classes from the 1970s. Lamb body mass varied in relation to the spatiotemporal variation of habitat closure within the hunting-free reserve, suggesting that habitat closure explains part of these changes. However, the fact that there was no significant spatial variation in body mass in the early part of the study, when a decline in phenotypic quality already had occurred, provided support for the influence of selective harvesting. We also found that the allometric relationship between horn breadth and horn length changed over the study period. For a given horn length, horn breadth was lower during the second part of the study. This result, as well as changes in horn curve diameter, supports the interpretation that selective harvesting of males based on their horn configuration had evolutionary consequences for horn shape, since this phenotypic trait is less likely to be affected by changes in habitat characteristics. Moreover, males required more time (approximately four years) to develop a desirable trophy, suggesting that trophy hunting favors the reproductive contribution of animals with slow-growing horns. Managers should exploit hunters' desire for trophy males to finance management strategies which ensure a balance between the population and its environment. However, for long-term sustainable exploitation, harvest strategy should also ensure that selectively targeted males are allowed to contribute genetically to the next generations.     

Estimating survival rates with time series of standing age-structure data

It has long been recognized that age-structure data contain useful information for assessing the status and dynamics of wildlife populations. For example, age-specific survival rates can be estimated with just a single sample from the age distribution of a stable, stationary population. For a population that is not stable, age-specific survival rates can be estimated using techniques such as inverse methods that combine time series of age-structure data with other demographic data. However, estimation of survival rates using these methods typically requires numerical optimization, a relatively long time series of data, and smoothing or other constraints to provide useful estimates. We developed general models for possibly unstable populations that combine time series of age-structure data with other demographic data to provide explicit maximum likelihood estimators of age-specific survival rates with as few as two years of data. As an example, we applied these methods to estimate survival rates for female bison (Bison bison) in Yellowstone National Park, USA. This approach provides a simple tool for monitoring survival rates based on age-structure data.     

A survival-analysis-based simulation model for Russian wheat aphid population dynamics

A simulation model for Russian wheat aphid (RWA), Diuraphis noxia (Mordvilko), populations is built by integrating survival-analysis-based development and survivor functions and the same-shape reproduction distribution model in the framework of Leslie [Leslie, P.H., 1945. On the use of matrices in certain population mathematics. Biometrika 33, 183–212] matrix structure. Survival analysis is utilized to model both the development and survival of RWA populations, and the Cox (1972) proportional hazards model is fitted with the data sets from our laboratory observation of 1800 RWA individuals under 25 factorial combinations of five temperature regimes and five barley plant-growth stages. Rather than using simple age-specific survivor rates as in the traditional Leslie matrix, the survivor functions based on survival analysis describe age-specific, temperature and plant stage-dependent RWA survival probabilities. Similarly, a probability model from survival analysis to estimate the probability that an individual will reach mature adult stage is utilized to describe the development process; this makes the transition from nymphal stage to mature adult stage dependent on RWA age as well as temperature and plant-growth stage.Inspired by the same-shape distribution and rate-summation approach for modeling insect development, a similar approach for modeling insect reproduction under variable temperature is developed. This new same-shape reproduction distribution model incorporates individual variation in reproduction capability, as well as the effects of RWA age, temperature and plant-growth stage. Consequently, the same-shape reproduction distribution model replaces the simple age-specific fecundities in Leslie matrix model. To the best of our knowledge, this work is the first to introduce survival analysis to simulation modeling in entomology and ecology and also the first to integrate our newly developed same-shape reproduction distribution model into application.