- Open Access
Unravelling the sex- and age-specific impact of poaching mortality with multievent modeling
Frontiers in Zoology volume 16, Article number: 20 (2019)
Poaching is a prominent source of ‘hidden hurdles’, cryptic impacts of human activities that may hinder the conservation of animal populations. Estimating poaching mortality is challenging, as the evidence for illegal killing is not outwardly obvious. Using resighting and recovery data collected on 141 marked red deer Cervus elaphus within the Stelvio National Park (central Italian Alps), we show how multievent models allow to assess the direct impacts of illegal harvesting on age- and sex-specific survival, accounting for uncertainty over mortality causes.
Mortality caused by poaching was consistently higher for males than for females in all age classes. In males, the probability of dying from poaching was higher for extreme age classes, while in females all age classes showed fairly similar values of poaching mortality. The strong bias in sex-specific poaching mortality was possibly due to trophy killing in adult males and ‘bushmeat-like’ killing for private or commercial gain in young males and in females.
A robust assessment of age- and sex-specific prevalence of poaching in wildlife populations is pivotal when illegal killing is of conservation concern. This provides timely information on what segment of the population is most likely to be affected. Besides obvious demographic consequences on small populations, age- and sex-biased poaching prevalence may contrast with the need to maintain ecosystem complexity and may alter behavioral responses to human presence. The information provided by multievent models, whose flexibility makes them adaptable to many systems where individual-based data is part of population monitoring, offers a support to design appropriate strategies for the conservation of wildlife populations.
Knowing the impacts of human activities on animal behavior , population structure  and, ultimately, population dynamics, is pivotal to design management strategies for the conservation of wildlife populations . Through this knowledge, conservation management aims to secure ecosystem complexity, while allowing for the sustainable use of wildlife – be it for consumptive or non-consumptive purposes.
Some impacts, however, are difficult to either predict, study or quantify. Recreational activities or trophy hunting, for example, may have negativeimpacts on animals’ vigilance activity and personality traits, e.g., in elk Cervus canadensis [4, 5], population genetic quality, e.g., in bighorn sheep Ovis canadensis  or, more generally, social structure in ungulates . Despite their potential negative effects on wild populations, such ‘hidden hurdles’ remain hard to integrate into decision-making protocols. One of the most prominent human-related sources of ‘hidden hurdles’ in wildlife conservation is poaching, an illegal and lethal activity whose direct and indirect impacts are, by their very nature, cryptic .
For critically endangered mammals, the relevance of illegal killing is obvious from a demographic standpoint (e.g., rhinos  or Amur tiger Panthera tigris altaica ). Poaching, however, may have less predictable consequences, including for example the disruption of social structures and of mating system . Poaching consequences may be relevant also in species of less conservation concern. For instance, in the European roe deer Capreolus capreolus illegal killing can undo the effects of planned management strategies, possibly jeopardizing the sustainable use of resources . Poaching on common species can also have negative effects on ecosystem complexity and food webs, as poachers may compete with large carnivores over food resources and threaten their survival . These effects may be expected to aggravate when legal and illegal killing co-occur, in absence of compensatory mortality . The assessment of age- and sex-specific mortalities associated with poaching provides timely information on the sensitivity of individuals to illegal harvesting. This, in turn, may help to design appropriate management strategies where poaching is a concern. Nonetheless, estimating poaching-related mortality is challenging .
Unbiased estimates of the prevalence of illegal hunting may be obtained through long term monitoring of marked animals and analysis with capture–mark–recapture (CMR) models to account for detection failure, e.g. [15, 16]. To date, however, the models proposed to estimate poaching mortality seldom included tag loss / GPS malfunction (but see ) and never accounted for uncertainty over the causes of death. Multievent CMR models are particularly suited for estimating mortality probabilities by different sources  while accounting for changes in mark type carried by the individuals [19, 20], different detection probability among marks, as well as for uncertainty in the causes of death of recovered carcasses.
We used a multievent CMR modeling approach  to model age- and sex-specific poaching prevalence in red deer Cervus elaphus within a protected area in the Italian Alps using data collected on differently marked individuals. The red deer is widely distributed across Europe, albeit patchily , but some taxa are of conservation concern (e.g., Corsican red deer Cervus elaphus corsicanus, Mesola deer Cervus elaphus italicus, barbary deer Cervus elaphus barbarus) . The assessment of the impact of poaching mortality may thus provide timely information to maintain viable populations for, e.g., ecosystem complexity, contemplation by tourists, education or research where the species is abundant, and to improve conservation status where poaching is a threat (cf.  for Corsican red deer).
Our analytical framework for the estimation of poaching mortality incorporates tag loss, detection failure, heterogeneity of recapture probability and uncertainty over mortality causes. We show how model flexibility can be used to provide robust estimates of the direct impact of poaching on age- and sex-specific mortality probabilities, relative to other causes of death. Given the cryptic nature of poaching, no a priori hypotheses could be formulated on the age- and sex- specific impacts of illegal killing in our study population.
Study site and data collection
The study was conducted from 2007 to 2017 in an area of ca. 27,900 ha between 1200 and 3850 m a.s.l. within the Stelvio National Park, central Italian Alps (10°25′N, 46°27′E). Here, red deer abundance increased sharply over the past decades, reaching current densities of about 27.4 ind./km2 (± 2.5 SD), cf. . In 2011 the Park Authority established a culling plan to limit deer impacts on forest regeneration and agricultural activities. Culling always occurred between late October and early February, with the support of professional hunters under the supervision of the Park Authority. As the aim was to reduce deer density, the sex-ratio in the culling plan was biased towards females (2:1); no restrictions were imposed on females in terms of age and lactating status, whereas only males < 6 years of age could be culled. Marked deer were not allowed for culling. During the study period 56 males and 85 females were captured by darting between October and April. After sedation, each animal was assigned to a given sex and age, based on inspection of tooth wear and eruption. All 141 individuals were marked with ear tags; 35 males and 83 females were additionally marked with colored belts with unique pattern of colors and reflectors (hereafter ‘optical collars’). Twenty-nine of these animals (13 males and 16 females) had a Global Positioning System (GPS, Vectronic Aerospace GmbH) mounted on the optical collar. In summary, animals were of three ‘mark types’: ear tags only, ear tags plus optical collars, ear tags plus optical collars equipped with a GPS device (Fig. 1). All GPS collars had a timer-controlled drop-off device set to trigger after 3 years. Therefore, GPS-collared individuals were allowed to become ‘plainly ear-tagged deer’ (in case of successful drop-off of the collar) or ‘optical-collared deer’ (in case of drop-off failure and cease of the GPS signal). Individuals with a working GPS device had perfect detection probabilities, as GPS locations were automatically collected on an hourly basis until the battery ran out, the GPS failed permanently, or until the collar dropped-off. Resightings of deer with inactive or dropped-off GPS collars, optical collars and/or ear tags only were conducted using spotlight counts between mid-April and mid-May (details in ).
Data coding, model definition and selection
Live encounter data were collected from mid-April to mid-May of each year between 2008 and 2017 and coded into individual encounter histories. Observations were stratified according to the mark-type. We assigned ‘1’ to individuals detected with working GPS devices; ‘2’ to individuals observed alive with optical collar; ‘3’ to individuals observed with ear tags; ‘4’ to individuals resighted with inactive GPS collars (see Additional file 1 for further details on observation coding). Animals alive with inactive GPS collar have a detection probability similar to those with optical collars, but this code had to be specified to account for drop-off failure. A ‘0’ was assigned to individuals not observed on a given year. In addition to live encounter data, dead recoveries (n = 50; 20 males and 30 females) were reported by the park personnel or by local people and tourists on a non-systematic basis. Deer carcasses found in the field were subsequently inspected by a veterinarian, who established the likely causes of death, thereby minimizing the uncertainty associated to death cause assignment. Causes of death were defined as: ‘poaching’ (n = 16 individuals found within the park boundaries in no-hunting periods, with gunshot wounds or in abnormal situations, e.g., buried under snow), ‘others’ (n = 25, e.g., starvation, collision with cars or legally shot outside of the Park boundaries) and ‘unknown’ (n = 9 deer for which causes of death could not be unequivocally assigned). Dead encounters were coded depending on mark type and cause of death: ‘5’ (poached, with working GPS); ‘6’ (poached, with inactive or dropped-off GPS, optical collar or ear tags); ‘7’ (dead by other causes, with working GPS); ‘8’ (dead by other causes, with inactive or dropped-off GPS, optical collar or ear tags); ‘9’ (unknown causes, with working GPS); ‘10’ (unknown causes, with inactive or dropped-off GPS, optical collar or ear tags; cf. Additional file 1). Recovery data for individuals with optical collars (n = 29) and ear tags only (n = 6) were pooled to avoid identifiability problems caused by the small sample size for the latter mark type.
The model for the encounter histories described above was specified using a multievent statistical framework . As no goodness-of-fit (GOF) test is available for multievent models, we first tested the adequacy of a model assuming survival and resight parameters as time-dependent (i.e., the Cormack-Jolly-Seber model), using live observations only (codes 2, 3, and 4) in program U-CARE 2.3.2 . Observations collected by GPS-transmitters and recoveries were not considered for GOF because the former have perfect probability of recapture and the latter were too sparse for full time-dependent models. In the following paragraph we provide a verbal description of our multievent modeling procedure, technical details are reported in Additional file 1.
We specified the multievent models with three sets of parameters: 1) the initial state probabilities, i.e., the initial proportions of individuals with different mark type; 2) the probabilities of transition between the states, decomposed in the probabilities of different types of mark loss and the probability of mortality; 3) the probabilities of resight and recovery (Additional file 2). We considered an additional parameter to accommodate the uncertainty over the cause of death, so that the unknown causes were automatically ‘redistributed’ between the known ones. Due to identifiability issues, this parameter was modelled to independent of the cause of death. The model included 9 different states representing the type and status of the mark (working GPS collar, inactive GPS collar, optical collar, ear tags), the state alive or dead, and the cause of mortality (poaching vs. others) and 11 possible events (see above). Parameters were estimated simultaneously by maximum likelihood with program E-SURGE . We applied the same model design in males and females, separately.
Model selection began by fitting a general model for either sex, assuming the following structure: two age classes for the parameter ‘GPS-signal loss’ (GPS collars were set for a 2 year duration), a constant probability of GPS-transmitter drop-off, an age effect in survival using different age structures for young, adult and old animals based on the existing literature, e.g. [2, 27, 28], time- and mark type-dependent resighting probabilities, time-dependent recovery probabilities and a constant probability to determine the cause of death. We first selected the age structure that minimized the value of Akaike’s Information Criterion adjusted for the effective sample size (AICc ). We used this structure to test for time dependence in the recovery probability, and subsequently assess the effects of time and type of mark on resighting probabilities. Finally, we modelled different sources of mortality across age classes, in either sex, by contrasting 3 alternative scenarios: 1) mortalities vary interactively with age; 2) mortalities vary in parallel with age (additive model); 3) mortalities vary independently with age (constant model). Model selection was based on AICc values. Akaike’s weight (wj) for each model j in the final candidate set was calculated as an index of relative model plausibility and used to obtain averaged estimates .
The goodness-of-fit test of the Cormack-Jolly-Seber time-dependent model – with only live observations – indicated an adequate description of the data (χ2 = 45.86, df = 47, p = 0.52). The modeling of age effects on female mortality probabilities (Models 1F–3F, Table 1) supported different survival probabilities between yearlings (1 year), adults (2–7 years) and old animals (8+ years) (Model 1F, Table 1). For males (Models 1 M–3 M, Table 1) a slightly different age structure proved a better fit to the data: young (1–3 years), adults (4–7 years) and old animals (8+ years) (Model 2 M, Table 1). Models with constant recovery probabilities were selected in both sexes (Models 4F and 4 M, Table 1) whereas resighting probabilities showed temporal variation in females (Model 6F, Table 1) and varied depending on mark type in males (Model 7 M, Table 1). In females, the model with the lowest AICc value indicated that poaching mortality was constant across all age classes (Model 10F, Table 1). This model had a similar AICc value to the model assuming parallel variation of mortality causes with age (Model 9F, Table 1). In males, the age-dependent model with additive effects of mortality causes returned the lowest AICc value (Models 9 M, Table 1).
Averaged estimates indicated that GPS-transmitters began to fail the second year after deployment: the probability of signal loss was slightly higher for females (0.39, 95%CI = 0.17–0.66) than for males (0.27, 95%CI = 0.07–0.64). Annual probability of GPS drop-off was 0 for females and 0.15 for males (95%CI = 0.05–0.39). Dead recovery probabilities were 0.47 (95%CI = 0.33–0.61) for females and 0.38 (95%CI = 0.24–0.55) for males. Female resighting probabilities varied over time between 0.26 and 0.96 (mean = 0.75) and did not depend on the type of mark. Males marked with optical collars showed higher resighting probabilities (0.47, 95%CI = 0.33–0.61) than males marked with ear tags only (0.23, 95%CI = 0.13–0.36). Averaged survival estimates in males and females showed age-dependent variations, with highest values for the intermediate classes (Table 2). However, while females had fairly consistent survival probabilities across ages, males showed much greater fluctuations (Table 2). The estimates of mortality caused by poaching were consistently higher for males than for females in all age classes (Table 2, Fig. 2). In males, the probability of dying from poaching was higher for young and – especially – old individuals while in females all age classes showed fairly similar values of poaching mortality (Table 2, Fig. 2). Notably, in males the point estimates of poaching mortality were consistently higher than those of other mortality events, accounting for 52% of the total mortality in young individuals, 67% in adults and 56% in old deer (Fig. 2). In females, poaching accounted for c. 22% of the total mortality in yearling, 33% in adult hinds, and 22% in old individuals (Fig. 2).
Assessing poaching prevalence in wildlife populations is difficult. The multievent capture-mark-recapture modeling approach used here allowed to efficiently estimate sex- and age-specific mortality due to illegal killing, while accounting for different types of marks, tag loss and uncertainty over the causes of death. This information is of paramount importance for appropriate conservation management strategies, as it allows to identify the segments of the population most affected by poaching. Furthermore, the flexibility of multievent modeling makes this approach suitable for any study system where the collection of individual-based data is part of population monitoring (e.g., wolf Canis lupus ; brown bear Ursus arctos ).
Poaching-related mortality is among the main issues for the conservation of many threatened species, from large mammals to reptiles (e.g., rhinos ; cheetah Acinonyx jubatus ; leatherback turtles Dermochelys coriacea ). Less intuitively, poaching may have negative effects also in non-threatened species, when it causes depletion that may affect natural food webs (e.g., in carnivore-ungulate systems ) or when it contrasts with management plans aimed at securing a sustainable use of resources for consumptive and non-consumptive purposes. Here we found that male deer suffered a much greater poaching mortality than females in all age classes, with peaks of 25% in old stags (Table 2). In females, poaching accounted for 3–4% of the mortality, regardless of age. Although poaching may not raise immediate conservation concerns in super-abundant species, consequences may be expected to occur at different levels.
Management plans typically aim to achieve multiple goals. In our study site, for example, deer management aimed to lower consequences of deer presence on the ecosystem by reducing negative impacts on forest regeneration and agricultural activities, while conserving a viable population and a fearless behavior in deer, for contemplation by tourists. In this respect, the consequences of poaching are not straightforward to predict. While, in principle, poaching may help reducing browsing impacts by lowering densities in super-abundant ungulates, it may contrast with the long-term interest of maintaining demographically viable populations. From a demographic standpoint, in fact, age- and sex-biased poaching mortality is expected to bear fitness consequences . Our data, for example, show evidence for the occurrence of trophy-killing on adult males. Although there is no evidence of a positive association between breeding values for antler size and fitness in red deer, breeding success correlates positively with antler size . Furthermore, the removal of adult individuals from the population may lead to a greater reproductive investment by young adults, with potential detrimental consequences on their overwinter survival . To ensure the maintenance of viable populations, which are crucial to ecosystem complexity – especially when large predators are present , investigations on long-term demographic consequences of the sex- and age-specific impact of illegal killing are warranted, for example through the integrated analysis of surveys and individual-based information on poaching prevalence [31, 36].
Poaching may also bear consequences at the behavioral level, and the investigation of poaching prevalence may help depicting negative effects on wildlife populations. The occurrence of trophy-killing, for example, suggests that poachers may concentrate their activity during the rutting season, when the opportunity for successful poaching on adult males is greatest. This may select for fearful behavioral responses (cf. ) and make animals avoid human presence, in contrast with conservation strategies that aim to guarantee the possibility of contemplation by tourists (e.g., during the mating season). Furthermore, the assessment of age- and sex-specific poaching prevalence may provide clues to optimize survey activity by wildlife wardens which, in our case, may concentrate during the rutting period, to avoid undesirable consequences of poaching on males.
The information provided by multievent models, whose flexibility makes them adaptable to many systems where individual-based data is part of population monitoring, offers a support to design appropriate strategies for the conservation of wildlife populations. It is worth noting that our approach to investigate poaching prevalence can be extended to different systems where individual-based data are available. Monitoring individual-based information in natural populations entails the need to account for detection probability. When multiple marks are considered, there is the additional problem of dealing with a detection probability that depends on mark type and possibly with mark loss.  proposed a multievent model based on conditional processes to incorporate tag loss and different type of detection probabilities. Here, we have expanded this model to include a larger number of mark types and, more importantly, the uncertainty over mortality causes. This expansion, which allowed to include recoveries from uncertain causes, has also entailed a new event and a new set of parameters. The trade-off between additional information and model complexity should be carefully evaluated. A similar trade-off exists in including additional age classes or mark types. Despite these limitations, the multievent modelling approach allowed to investigate the age- and sex-specific impacts of poaching relative to other causes of mortality, framing the heterogeneity due to mark type, multiple data sources and uncertainty over mortality causes.
In our study, several deer were equipped with radio-transmitters. Radio-tracking studies have two important limitations: due to the costs of transmitters, only a relatively small number of animals can be marked and the information only lasts the lifespan of the battery. We showed how to integrate this information with those collected by using only ear tags. The current model can be additionally used, for example, to explore through simulations the costs and benefits of different mark methods and provide indications for planning the optimal study design given the available funds (e.g. radio-transmitters vs. ear tags) for species of conservation concern.
Akaike Information Criterion corrected for small samples
Goodness Of Fit
Global Positioning System
Berger-Tal O, Saltz D. Introduction: the whys and hows of conservation behavior. In: Berger-Tal O, Saltz D, editors. Conservation behavior. Applying behavioral Ecology to Wildlife conservation and management. Cambridge: Cambridge University Press; 2016. p. 3–35.
Gaillard J-M, Festa-Bianchet M, Yoccoz NG, Loison A, Toïgo C. Temporal variation in fitness components and population dynamics of large herbivores. Ann Rev Ecol Evol Syst. 2000;31:367–93.
Fryxell JM, Sinclair ARE, Caughley G. Wildlife Ecology, Conservation, and Management. Chichester: Wiley Blackwell; 2006.
Ciuti S, Northrup JM, Muhly TB, Simi S, Musiani M, Pitt JA, Boyce MS. Effects of humans on behaviour of wildlife exceed those of natural predators in a landscape of fear. PLoS One. 2012;7:e50611.
Ciuti S, Muhly TB, Paton DG, McDevitt AD, Musiani M, Boyce MS. Human selection of elk behavioural traits in a landscape of fear. Proc R Soc Lond B. 2012;279:4407–16.
Coltman DW, O’Donoghue P, Jorgenson JT, Hogg JT, Strobeck C, Festa-Bianchet M. Undesirable evolutionary consequences of trophy hunting. Nature. 2003;426:655–8.
Putman R, Apollonio M. Behaviour and management of European ungulates. In: Behaviour and Management of European Ungulates (eds. Putman R, Apollonio M). Dunbeath: Whittles Publishing; 2014. p. 1–12.
Liberg O, Chapron G, Wabakken P, Pedersen HC, Thompson Hobbs N, Sand H. Shoot, shovel and shut up: cryptic poaching slows restoration of a large carnivore in Europe. Proc R Soc Lond B. 2012;279:910–5.
Koen H, de Villiers JP, Roodt H, de Waal A. An expert-driven causal model of the rhino poaching problem. Ecol Model. 2017;347:29–39.
Goodrich JM, Kerley LL, Smirnov EN, Miquelle DG, McDonald L, Quigley HB, Hornocker MG, McDonald T. Survival rates and causes of mortality of Amur tigers on and near the Sikhote-Alin biosphere Zapovednik. J Zool. 2008;276:323–9.
Ishengoma DRS, Shedlock AM, Foley CAH, Foley LJ, Wasser SK, Balthazary ST, Mutayoba BM. Effects of poaching on bull mating success in a free ranging African elephant (Loxodonta africana) population in Tarangire National Park, Tanzania. Cons Gen. 2008;9:247–55.
Sönnichsen L, Borowik T, Podgórski T, Plis K, Berger A, Jędrzejewska B. Survival rates and causes of mortality of roe deer Capreolus capreolus in a rural landscape, eastern Poland. Mamm Res. 2017;62:141–7.
Ghoddousi A, Soofi M, Hamidi AK, Lumetsberger T, Egli L, Ashayeri A, Khorozyan I, Kiabi BH, Waltert M. When pork is not on the menu: assessing trophic competition between large carnivores and poachers. Biol Conserv. 2017;209:223–9.
Sæther BE, Engen S, Odden J, Linnell JDC, Grøtan V, Andrén H. Sustainable harvest strategies for age-structured Eurasian lynx populations: the use of reproductive value. Biol Conserv. 2010;143:1970–9.
Persson J, Ericsson G, Segerström P. Human caused mortality in the endangered Scandinavian wolverine population. Biol Conserv. 2009;142:325–31.
Andrén H, Linnell JDC, Liberg O, Andersen R, Danella A, Karlsson J, Odden J, Moa PF, Ahlqvista P, Kvamb T, Franzénd R, Segerström P. Survival rates and causes of mortality in Eurasian lynx (Lynx lynx) in multi-use landscapes. Biol Conserv. 2006;131:23–32.
Tavecchia G, Adrover J, Muñoz Navarro A, Pradel R. Modelling mortality causes in longitudinal data in the presence of tag loss: application to raptor poisoning and electrocution. J Appl Ecol. 2012;49:297–305.
Pradel R. Multievent: an extension of multistate capture–recapture models to uncertain states. Biometrics. 2005;61:442–7.
Sanz-Aguilar A, De Pablo F, Donázar JA. Age-dependent survival of island vs. mainland populations of two avian scavengers: delving into migration costs. Oecologia. 2015;179:405–14.
Souchay G, Schaub M. Investigating rates of hunting and survival in declining European lapwing populations. PLoS One. 2016;11:e0163850.
Lovari S, Lorenzini R, Masseti M, Pereladova O, Carden RF, Brook SM. Cervus elaphus. The IUCN red list of threatened species 2016;e.T55997072A22155320.
Iacolina L, Corlatti L, Buzan E, Safner T, Šprem N. Hybridisation in European ungulates: an overview of the current status, causes, and consequences. Mamm Rev. 2019;49:45–59.
Kidjo N, Feracci G, Bideau E, Gonzalez G, Mattéi C, Marchand B, Aulagnier S. Extirpation and reintroduction of the Corsican red deer Cervus elaphus corsicanus in Corsica. Oryx. 2007;41:488–94.
Corlatti L, Gugiatti A, Pedrotti L. Spring spotlight counts provide reliable indices to track changes in population size of mountain-dwelling red deer Cervus elaphus. Wildl Biol. 2016;22:268–76.
Choquet R, Lebreton J, Gimenez O, Reboulet A, Pradel R. U-CARE: utilities for performing goodness of fit tests and manipulating CApture–REcapture data. Ecography. 2009;32:1071–4.
Choquet R, Rouan L, Pradel R. Program E-SURGE: a software application for fitting multievent models. In: Thomson DL, Cooch EG, Conroy MJ, editors. Modeling Demographic Processes in Marked Populations. New York: Springer; 2009. p. 845–65.
Loison A, Festa-Bianchet M, Gaillard J-M, Jorgenson JT, Jullien J-M. Age-specific survival in five populations of ungulates: evidence of senescence. Ecology. 1999;80:2539–5254.
Catchpole EA, Fan Y, Morgan BJT, Clutton-Brock TH, Coulson T. Sexual dimorphism, survival and dispersal in Red Deer. J Agric Biol Env Stat. 2004;9:1–26.
Burnham KP, Anderson DR. Model selection and multi-model inference: a practical information-theoretic approach. New York: Springer; 2002.
Mech LD. Where can wolves live and how can we live with them? Biol Conserv. 2017;210:310–7.
Tenan S, Iemma A, Bragalanti N, Pedrini P, Barba M, Randi E, Groff C, Genovart M. Evaluating mortality rates with a novel integrated framework for non-monogamous species. Cons Biol. 2016;30:1307–19.
Durant SM, Mitchell N, Groom R, Pettorelli N, Ipavec A, Jacobson AP, et al. The global decline of cheetah Acinonyx jubatus and what it means for conservation. PNAS. 2017;114:528–33.
Santidrián Tomillo P, Saba VS, Piedra R, Paladino FV, Spotila JR. Effects of illegal harvest of eggs on the population decline of leatherback turtles in las Baulas marine National Park, Costa Rica. Cons Biol. 2008;22:1216–24.
Ginsberg JR, Milner-Gulland EJ. Sex biased harvesting and population dynamics in ungulates: implications for conservation and sustainable use. Cons Biol. 1994;8:157–66.
Kruuk LEB, Slate J, Pemberton JM, Brotherstone S, Guinness F, Clutton-Brock T. Antler size in red deer: heritability and selection but no evolution. Evolution. 2002;56:1683–95.
Tavecchia G, Besbeas P, Coulson T, Morgan BJT, Clutton-Brock TH. Estimating population size and hidden demographic parameters with state-space modeling. Am Nat. 2009;173:722–33.
We are grateful to A. Zanoli for his invaluable help in capturing deer. We thank all the people that in various ways helped us collecting resighting and recovery data, and four anonymous reviewers for helpful comments on previous drafts of the manuscript.
To conduct the research leading to these results, LC has received funding from the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme (FP7/2007–2013) under REA grant agreement n° . ASA was supported by a Ramón y Cajal contract (RYC-2017- 22796) co-funded by the Ministerio de Ciencia, Innovación y Universidades, the Agencia Estatal de Investigación and the European Social Fund (“El Fondo Social Europeo invierte en tu futuro”). The article processing charge was funded by the German Research Foundation (DFG) and the University of Freiburg in the funding programme Open Access Publishing.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.