Effects of Trophy Hunting on Lion and Leopard Populations in Tanzania

Abstract: Tanzania holds most of the remaining large populations of African lions (Panthera leo) and has extensive areas of leopard habitat (Panthera pardus), and both species are subjected to sizable harvests by sport hunters. As a first step toward establishing sustainable management strategies, we analyzed harvest trends for lions and leopards across Tanzania's 300,000 km² of hunting blocks. We summarize lion population trends in protected areas where lion abundance has been directly measured and data on the frequency of lion attacks on humans in high‐conflict agricultural areas. We place these findings in context of the rapidly growing human population in rural Tanzania and the concomitant effects of habitat loss, human‐wildlife conflict, and cultural practices. Lion harvests declined by 50% across Tanzania between 1996 and 2008, and hunting areas with the highest initial harvests suffered the steepest declines. Although each part of the country is subject to some form of anthropogenic impact from local people, the intensity of trophy hunting was the only significant factor in a statistical analysis of lion harvest trends. Although leopard harvests were more stable, regions outside the Selous Game Reserve with the highest initial leopard harvests again showed the steepest declines. Our quantitative analyses suggest that annual hunting quotas be limited to 0.5 lions and 1.0 leopard/1000 km² of hunting area, except hunting blocks in the Selous Game Reserve, where harvests should be limited to 1.0 lion and 3.0 leopards/1000 km².     

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

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

Modeling the Effects of Trophy Selection and Environmental Disturbance on a Simulated Population of African Lions

Abstract:  Tanzania is a premier destination for trophy hunting of African lions (Panthera leo) and is home to the most extensive long-term study of unhunted lions. Thus, it provides a unique opportunity to apply data from a long-term field study to a conservation dilemma: How can a trophy-hunted species whose reproductive success is closely tied to social stability be harvested sustainably? We used an individually based, spatially explicit, stochastic model, parameterized with nearly 40 years of behavioral and demographic data on lions in the Serengeti, to examine the separate effects of trophy selection and environmental disturbance on the viability of a simulated lion population in response to annual harvesting. Female population size was sensitive to the harvesting of young males (≥3 years), whereas hunting represented a relatively trivial threat to population viability when the harvest was restricted to mature males (≥6 years). Overall model performance was robust to environmental disturbance and to errors in age assessment based on nose coloration as an index used to age potential trophies. Introducing an environmental disturbance did not eliminate the capacity to maintain a viable breeding population when harvesting only older males, and initially depleted populations recovered within 15–25 years after the disturbance to levels comparable to hunted populations that did not experience a catastrophic event. These results are consistent with empirical observations of lion resilience to environmental stochasticity .     

Conservation and monitoring of a persecuted African lion population by Maasai warriors

Although Africa has many threatened species and biological hot spots, there are few citizen science schemes, particularly in rural communities, and there has been limited evaluation of existing programs. We engaged traditional Maasai warriors (pastoralist men aged 15 to 35) in community‐based conservation and demographic monitoring of a persecuted African lion (Panthera leo) population. Through direct engagement, we investigated whether a citizen science approach employing local warriors, who had no formal education, could produce reliable data on the demographics, predation, and movements of a species with which their communities have been in conflict for generations. Warriors were given benefits such as literacy training and skill enhancement and engaged in the monitoring of the lions. The trained warriors reported on lion sign across an area nearly 4000 km². Scientists worked together with the warriors to verify their reports and gather observations on the lion population. Using the verified reports and collected observations, we examined our scientific knowledge relative to the lion population preceding and during the citizen science program. Our observations showed that data quality and quantity improved with the involvement and training of the participants. Furthermore, because they engaged in conservation and gained personal benefits, the participants came to appreciate a species that was traditionally their foe. We believe engaging other local communities in biodiversity conservation and monitoring may be an effective conservation approach in rural Africa.     

Use of Spatial Capture‐Recapture Modeling and DNA Data to Estimate Densities of Elusive Animals

Abstract: Assessment of abundance, survival, recruitment rates, and density (i.e., population assessment) is especially challenging for elusive species most in need of protection (e.g., rare carnivores). Individual identification methods, such as DNA sampling, provide ways of studying such species efficiently and noninvasively. Additionally, statistical methods that correct for undetected animals and account for locations where animals are captured are available to efficiently estimate density and other demographic parameters. We collected hair samples of European wildcat (Felis silvestris) from cheek‐rub lure sticks, extracted DNA from the samples, and identified each animals’ genotype. To estimate the density of wildcats, we used Bayesian inference in a spatial capture‐recapture model. We used WinBUGS to fit a model that accounted for differences in detection probability among individuals and seasons and between two lure arrays. We detected 21 individual wildcats (including possible hybrids) 47 times. Wildcat density was estimated at 0.29/km² (SE 0.06), and 95% of the activity of wildcats was estimated to occur within 1.83 km from their home‐range center. Lures located systematically were associated with a greater number of detections than lures placed in a cell on the basis of expert opinion. Detection probability of individual cats was greatest in late March. Our model is a generalized linear mixed model; hence, it can be easily extended, for instance, to incorporate trap‐ and individual‐level covariates. We believe that the combined use of noninvasive sampling techniques and spatial capture‐recapture models will improve population assessments, especially for rare and elusive animals.     

Importance of Accounting for Detection Heterogeneity When Estimating Abundance: the Case of French Wolves

Abstract: Assessing conservation strategies requires reliable estimates of abundance. Because detecting all individuals is most often impossible in free‐ranging populations, estimation procedures have to account for a <1 detection probability. Capture–recapture methods allow biologists to cope with this issue of detectability. Nevertheless, capture–recapture models for open populations are built on the assumption that all individuals share the same detection probability, although detection heterogeneity among individuals has led to underestimating abundance of closed populations. We developed multievent capture–recapture models for an open population and proposed an associated estimator of population size that both account for individual detection heterogeneity (IDH). We considered a two‐class mixture model with weakly and highly detectable individuals to account for IDH. In a noninvasive capture–recapture study of wolves we based on genotypes identified in feces and hairs, we found a large underestimation of population size (27% on average) occurred when IDH was ignored.     

Estimating Sonoran pronghorn abundance and survival with fecal DNA and capture–recapture methods

Population abundance estimates are important for management but can be challenging to determine in low‐density, wide‐ranging, and endangered species, such as Sonoran pronghorn (Antilocapra americana sonoriensis). The Sonoran pronghorn population has been increasing; however, population estimates are currently derived from a biennial aerial count that does not provide survival or recruitment estimates. We identified individuals through noninvasively collected fecal DNA and used robust‐design capture–recapture to estimate abundance and survival for Sonoran pronghorn in the United States from 2013 to 2014. In 2014 we generated separate population estimates for pronghorn gathered near 13 different artificial water holes and for pronghorn not near water holes. The population using artificial water holes had 116 (95% CI 102–131) and 121 individuals (95% CI 112–132) in 2013 and 2014, respectively. For all locations, we estimated there were 144 individuals (95% CI 132–157). Adults had higher annual survival probabilities (0.83, 95% CI 0.69–0.92) than fawns (0.41, 95% CI 0.21–0.65). Our use of targeted noninvasive genetic sampling and capture–recapture with Sonoran pronghorn fecal DNA was an effective method for monitoring a large proportion of the population. Our results provided the first survival estimates for this population in over 2 decades and precise estimates of the population using artificial water holes. Our method could be used for targeted sampling of broadly distributed species in other systems, such as in African savanna ecosystems, where many species congregate at watering sites.     

Efficacy of Two Lion Conservation Programs in Maasailand,Kenya

Lion (Panthera leo) populations are in decline throughout most of Africa. The problem is particularly acute in southern Kenya, where Maasai pastoralists have been spearing and poisoning lions at a rate that will ensure near term local extinction. We investigated 2 approaches for improving local tolerance of lions: compensation payments for livestock lost to predators and Lion Guardians, which draws on local cultural values and knowledge to mitigate livestock‐carnivore conflict and monitor carnivores. To gauge the overall influence of conservation intervention, we combined both programs into a single conservation treatment variable. Using 8 years of lion killing data, we applied Manski's partial identification approach with bounded assumptions to investigate the effect of conservation treatment on lion killing in 4 contiguous areas. In 3 of the areas, conservation treatment was positively associated with a reduction in lion killing. We then applied a generalized linear model to assess the relative efficacy of the 2 interventions. The model estimated that compensation resulted in an 87–91% drop in the number of lions killed, whereas Lion Guardians (operating in combination with compensation and alone) resulted in a 99% drop in lion killing. Eficacia de Dos Programas de Conservación de Leones en Maasailand, Kenia     

Optimization of capture–recapture monitoring of elusive species illustrated with a threatened grasshopper

Information on population sizes and trends of threatened species is essential for their conservation, but obtaining reliable estimates can be challenging. We devised a method to improve the precision of estimates of population size obtained from capture–recapture studies for species with low capture and recapture probabilities and short seasonal activity, illustrated with population data of an elusive grasshopper (Prionotropis rhodanica). We used data from 5 capture–recapture studies to identify methodological and environmental factors affecting capture and recapture probabilities and estimates of population size. In a simulation, we used the population size and capture and recapture probability estimates obtained from the field studies to identify the minimum number of sampling occasions needed to obtain unbiased and robust estimates of population size. Based on these results we optimized the capture–recapture design, implemented it in 2 additional studies, and compared their precision with those of the nonoptimized studies. Additionally, we simulated scenarios based on thresholds of population size in criteria C and D of the International Union for Conservation of Nature (IUCN) Red List to investigate whether estimates of population size for elusive species can reliably inform red-list assessments. Identifying parameters that affect capture and recapture probabilities (for the grasshopper time since emergence of first adults) and optimizing field protocols based on this information reduced study effort (−6% to −27% sampling occasions) and provided more precise estimates of population size (reduced coefficient of variation) compared with nonoptimized studies. Estimates of population size from the scenarios based on the IUCN thresholds were mostly unbiased and robust (only the combination of very small populations and little study effort produced unreliable estimates), suggesting capture–recapture can be considered reliable for informing red-list assessments. Although capture–recapture remains difficult and costly for elusive species, our optimization procedure can help determine efficient protocols to increase data quality and minimize monitoring effort.     

Limitation of African Wild Dogs by Competition with Larger Carnivores

African wild dogs ( Lycaon pictus ) are endangered largely because their population-density is low under all conditions. Interspecific competition with larger carnivores may be a factor limiting wild dog density. The density of wild dogs on a 2600-km² area of the Selous Game Reserve (Tanzania) was 0.04 adults/km². Spotted hyaena ( Crocuta crocuta ) density for the same area was estimated by audio playbacks as 0.32 hyaenas/km². Lion ( Panthera leo ) density, determined from the ratio of hyaenas to lions, was 0.11 lions/km². Across six ecosystems including Selous, there were strong negative correlations between wild dog and hyaena densities (r = −0.92; p = 0.01) and between wild dog and lion densities (r = −0.91; p = 0.03). Hyaenas out-numbered wild dogs by ratios ranging from 8:1 to 122:1. Ratios of lions to wild dogs ranged from 3:1 to 21:1. The diets of hyaenas and wild dogs overlap extensively; those of wild dogs and lions show less overlap. Where hyaenas are common and visibility is good, interference competition from hyaenas at wild dog kills is common and reduces wild dogs' feeding time. Where hyaena density is lower and visibility is poor, interference competition at wild dog kills is rare. Wild dogs are commonly killed by lions and occasionally by hyaenas. These data suggest that competition with spotted hyaenas may limit or exclude wild dogs when hyaena density is high. Competition with lions appears less intense, but direct predation by lions on wild dogs is important. Competition and predation by larger carnivores may be of broad importance to the conservation of wild dogs and other medium-sized carnivores.     

Determining Optimal Population Monitoring for Rare Butterflies

Abstract: Determining population viability of rare insects depends on precise, unbiased estimates of population size and other demographic parameters. We used data on the endangered St. Francis' satyr butterfly (Neonympha mitchellii francisci) to evaluate 2 approaches (mark–recapture and transect counts) for population analysis of rare butterflies. Mark–recapture analysis provided by far the greatest amount of demographic information, including estimates (and standard errors) of population size, detection, survival, and recruitment probabilities. Mark–recapture analysis can also be used to estimate dispersal and temporal variation in rates, although we did not do this here. Models of seasonal flight phenologies derived from transect counts (Insect Count Analyzer) provided an index of population size and estimates of survival and statistical uncertainty. Pollard–Yates population indices derived from transect counts did not provide estimates of demographic parameters. This index may be highly biased if detection and survival probabilities vary spatially and temporally. In terms of statistical performance, mark–recapture and Pollard–Yates indices were least variable. Mark–recapture estimates were less likely to fail than Insect Count Analyzer, but mark–recapture estimates became less precise as sampling intensity decreased. In general, count‐based approaches are less costly and less likely to cause harm to rare insects than mark–recapture. The optimal monitoring approach must reconcile these trade‐offs. Thus, mark–recapture should be favored when demographic estimates are needed, when financial resources enable frequent sampling, and when marking does not harm the insect populations. The optimal sampling strategy may use 2 sampling methods together in 1 overall sampling plan: limited mark–recapture sampling to estimate survival and detection probabilities and frequent but less expensive transect counts.     

Pheromone Monitoring of Rare and Threatened Insects: Exploiting a Pheromone–Kairomone System to Estimate Prey and Predator Abundance

Abstract: Pheromone‐based monitoring is a promising new method for assessing the conservation status of many threatened insect species. We examined the versatility and usefulness of pheromone‐based monitoring by integrating a pheromone–kairomone trapping system and pitfall trapping system in the monitoring of two saproxylic beetles, the hermit beetle Osmoderma eremita (Coleoptera: Scarabaeidae) and its predator Elater ferrugineus (Coleoptera: Elateridae), which live inside hollow trees. We performed mark–recapture studies of both species with unbaited pitfall traps in oak hollows combined with pheromone‐baited funnel traps suspended from oak branches to intercept dispersing individuals. For O. eremita, the integrated trapping system showed that the population in the study sites may be considerably higher than estimates based on extrapolation from pitfall trapping alone (approximately 3400 vs. 1100 or 1800 individuals, respectively). Recaptures between odor‐baited funnel traps showed that males and females had similar dispersal rates, but estimating the number of dispersing individuals was problematic due to declining recapture probability between subsequent capture events. Our conservative estimate, assuming a linear decrease in capture probability, suggested that around 1900 individuals, or at least half of the O. eremita population, may perform flights from their natal host trees, representing higher dispersal rates than previous estimates. E. ferrugineus was rarely caught in pitfall traps. One hundred thirty‐nine individuals, likely almost exclusively females, were caught in odor‐baited funnel traps with approximately 4% recapture probability. If recapture probability over consecutive capture events follows that of O. eremita, this would correspond to a total population size of 2500–3000 individuals of the predator; similar to its supposed prey O. eremita. Our results demonstrate that pheromone‐based monitoring is a valuable tool in the study of species or life‐history stages that would otherwise be inaccessible.     

Multiple data sources improve DNA-based mark-recapture population estimates of grizzly bears.

A fundamental challenge to estimating population size with mark-recapture methods is heterogeneous capture probabilities and subsequent bias of population estimates. Confronting this problem usually requires substantial sampling effort that can be difficult to achieve for some species, such as carnivores. We developed a methodology that uses two data sources to deal with heterogeneity and applied this to DNA mark-recapture data from grizzly bears (Ursus arctos). We improved population estimates by incorporating additional DNA "captures" of grizzly bears obtained by collecting hair from unbaited bear rub trees concurrently with baited, grid-based, hair snag sampling. We consider a Lincoln-Petersen estimator with hair snag captures as the initial session and rub tree captures as the recapture session and develop an estimator in program MARK that treats hair snag and rub tree samples as successive sessions. Using empirical data from a large-scale project in the greater Glacier National Park, Montana, USA, area and simulation modeling we evaluate these methods and compare the results to hair-snag-only estimates. Empirical results indicate that, compared with hair-snag-only data, the joint hair-snag-rub-tree methods produce similar but more precise estimates if capture and recapture rates are reasonably high for both methods. Simulation results suggest that estimators are potentially affected by correlation of capture probabilities between sample types in the presence of heterogeneity. Overall, closed population Huggins-Pledger estimators showed the highest precision and were most robust to sparse data, heterogeneity, and capture probability correlation among sampling types. Results also indicate that these estimators can be used when a segment of the population has zero capture probability for one of the methods. We propose that this general methodology may be useful for other species in which mark-recapture data are available from multiple sources.     

Estimating abundance without recaptures of marked pallid sturgeon in the Mississippi River

Abundance estimates are essential for assessing the viability of populations and the risks posed by alternative management actions. An effort to estimate abundance via a repeated mark‐recapture experiment may fail to recapture marked individuals. We devised a method for obtaining lower bounds on abundance in the absence of recaptures for both panmictic and spatially structured populations. The method assumes few enough recaptures were expected to be missed by random chance. The upper Bayesian credible limit on expected recaptures allows probabilistic statements about the minimum number of individuals present in the population. We applied this method to data from a 12‐year survey of pallid sturgeon (Scaphirhynchus albus) in the lower and middle Mississippi River (U.S.A.). None of the 241 individuals marked was recaptured in the survey. After accounting for survival and movement, our model‐averaged estimate of the total abundance of pallid sturgeon ≥3 years old in the study area had a 1%, 5%, or 25% chance of being <4,600, 7,000, or 15,000, respectively. When we assumed fish were distributed in proportion to survey catch per unit effort, the farthest downstream reach in the survey hosted at least 4.5–15 fish per river kilometer (rkm), whereas the remainder of the reaches in the lower and middle Mississippi River hosted at least 2.6–8.5 fish/rkm for all model variations examined. The lower Mississippi River had an average density of pallid sturgeon ≥3 years old of at least 3.0–9.8 fish/rkm. The choice of Bayesian prior was the largest source of uncertainty we considered but did not alter the order of magnitude of lower bounds. Nil‐recapture estimates of abundance are highly uncertain and require careful communication but can deliver insights from experiments that might otherwise be considered a failure.     

Multimethod,multistate Bayesian hierarchical modeling approach for use in regional monitoring of wolves

In many cases, the first step in large‐carnivore management is to obtain objective, reliable, and cost‐effective estimates of population parameters through procedures that are reproducible over time. However, monitoring predators over large areas is difficult, and the data have a high level of uncertainty. We devised a practical multimethod and multistate modeling approach based on Bayesian hierarchical‐site‐occupancy models that combined multiple survey methods to estimate different population states for use in monitoring large predators at a regional scale. We used wolves (Canis lupus) as our model species and generated reliable estimates of the number of sites with wolf reproduction (presence of pups). We used 2 wolf data sets from Spain (Western Galicia in 2013 and Asturias in 2004) to test the approach. Based on howling surveys, the naïve estimation (i.e., estimate based only on observations) of the number of sites with reproduction was 9 and 25 sites in Western Galicia and Asturias, respectively. Our model showed 33.4 (SD 9.6) and 34.4 (3.9) sites with wolf reproduction, respectively. The number of occupied sites with wolf reproduction was 0.67 (SD 0.19) and 0.76 (0.11), respectively. This approach can be used to design more cost‐effective monitoring programs (i.e., to define the sampling effort needed per site). Our approach should inspire well‐coordinated surveys across multiple administrative borders and populations and lead to improved decision making for management of large carnivores on a landscape level. The use of this Bayesian framework provides a simple way to visualize the degree of uncertainty around population‐parameter estimates and thus provides managers and stakeholders an intuitive approach to interpreting monitoring results. Our approach can be widely applied to large spatial scales in wildlife monitoring where detection probabilities differ between population states and where several methods are being used to estimate different population parameters.     

A Behaviorally Explicit Demographic Model Integrating Habitat Selection and Population Dynamics in California Sea Lions

Abstract: Although there has been a call for the integration of behavioral ecology and conservation biology, there are few tools currently available to achieve this integration. Explicitly including information about behavioral strategies in population viability analyses may enhance the ability of conservation biologists to understand and estimate patterns of extinction risk. Nevertheless, most behavioral‐based PVA approaches require detailed individual‐based data that are rarely available for imperiled species. We present a mechanistic approach that incorporates spatial and demographic consequences of behavioral strategies into population models used for conservation. We developed a stage‐structured matrix model that includes the costs and benefits of movement associated with 2 habitat‐selection strategies (philopatry and direct assessment). Using a life table for California sea lions (Zalophus californianus), we explored the sensitivity of model predictions to the inclusion of these behavioral parameters. Including behavioral information dramatically changed predicted population sizes, model dynamics, and the expected distribution of individuals among sites. Estimated population sizes projected in 100 years diverged up to 1 order of magnitude among scenarios that assumed different movement behavior. Scenarios also exhibited different model dynamics that ranged from stable equilibria to cycles or extinction. These results suggest that inclusion of behavioral data in viability models may improve estimates of extinction risk for imperiled species. Our approach provides a simple method for incorporating spatial and demographic consequences of behavioral strategies into population models and may be easily extended to other species and behaviors to understand the mechanisms of population dynamics for imperiled populations.     

Evaluating mortality rates with a novel integrated framework for nonmonogamous species

The conservation of wildlife requires management based on quantitative evidence, and especially for large carnivores, unraveling cause‐specific mortalities and understanding their impact on population dynamics is crucial. Acquiring this knowledge is challenging because it is difficult to obtain robust long‐term data sets on endangered populations and, usually, data are collected through diverse sampling strategies. Integrated population models (IPMs) offer a way to integrate data generated through different processes. However, IPMs are female‐based models that cannot account for mate availability, and this feature limits their applicability to monogamous species only. We extended classical IPMs to a two‐sex framework that allows investigation of population dynamics and quantification of cause‐specific mortality rates in nonmonogamous species. We illustrated our approach by simultaneously modeling different types of data from a reintroduced, unhunted brown bear (Ursus arctos) population living in an area with a dense human population. In a population mainly driven by adult survival, we estimated that on average 11% of cubs and 61% of adults died from human‐related causes. Although the population is currently not at risk, adult survival and thus population dynamics are driven by anthropogenic mortality. Given the recent increase of human‐bear conflicts in the area, removal of individuals for management purposes and through poaching may increase, reversing the positive population growth rate. Our approach can be generalized to other species affected by cause‐specific mortality and will be useful to inform conservation decisions for other nonmonogamous species, such as most large carnivores, for which data are scarce and diverse and thus data integration is highly desirable.     

A generalized estimating equations approach for capture–recapture closed population models

The estimation of population density animal population parameters, such as capture probability, population size, or population density, is an important issue in many ecological applications. Capture–recapture data may be considered as repeated observations that are often correlated over time. If these correlations are not taken into account then parameter estimates may be biased, possibly producing misleading results. We propose a generalized estimating equations (GEE) approach to account for correlation over time instead of assuming independence as in the traditional closed population capture–recapture studies. We also account for heterogeneity among observed individuals and over-dispersion, modelling capture probabilities as a function of covariates. The GEE versions of all closed population capture–recapture models and their corresponding estimating equations are proposed. We evaluate the effect of accounting for correlation structures on capture–recapture model selection based on the quasi-likelihood information criterion (QIC). An example is used for an illustrative application and for comparison to currently used methodology. A Horvitz–Thompson-like estimator is used to obtain estimates of population size based on conditional arguments. A simulation study is conducted to evaluate the performance of the GEE approach in capture-recapture studies. The GEE approach performs well for estimating population parameters, particularly when capture probabilities are high. The simulation results also reveal that estimated population size varies on the nature of the existing correlation among capture occasions.     

Closed population estimation models and their extensions in Program MARK

Program MARK provides > 65 data types in a common configuration for the estimation of population parameters from mark-encounter data. Encounter information from live captures, live resightings, and dead recoveries can be incorporated to estimate demographic parameters. Available estimates include survival (S or ϕ), rate of population change (λ), transition rates between strata (Ψ), emigration and immigration rates, and population size (N). Although N is the parameter most often desired by biologists, N is one of the most difficult parameters to estimate precisely without bias for a geographically and demographically closed population. The set of closed population estimation models available in Program MARK incorporate time (t) and behavioral (b) variation, and individual heterogeneity (h) in the estimation of capture and recapture probabilities in a likelihood framework. The full range of models from M ₀ (null model with all capture and recapture probabilities equal) to M _tbh are possible, including the ability to include temporal, group, and individual covariates to model capture and recapture probabilities. Both the full likelihood formulation of Otis et al. (1978) and the conditional model formulation of Huggins (1989, 1991) and Alho (1990) are provided in Program MARK, and all of these models are incorporated into the robust design (Kendall et al. 1995, 1997; Kendall and Nichols 1995) and robust-design multistrata (Hestbeck et al. 1991, Brownie et al. 1993) data types. Model selection is performed with AICc (Burnham and Anderson 2002) and model averaging (Burnham and Anderson 2002) is available in Program MARK to provide estimates of N with standard error that reflect model selection uncertainty.     

Spatial and temporal variability modify density dependence in populations of large herbivores

A central challenge in ecology is to understand the interplay of internal and external controls on the growth of populations. We examined the effects of temporal variation in weather and spatial variation in vegetation on the strength of density dependence in populations of large herbivores. We fit three subsets of the model ln(Nt) = a + (1 + b) x ln(N(t-1)) + c x ln(N(t-2)) to five time series of estimates (Nt) of abundance of ungulates in the Rocky Mountains, USA. The strength of density dependence was estimated by the magnitude of the coefficient b. We regressed the estimates of b on indices of temporal heterogeneity in weather and spatial heterogeneity in resources. The 95% posterior intervals of the slopes of these regressions showed that temporal heterogeneity strengthened density-dependent feedbacks to population growth, whereas spatial heterogeneity weakened them. This finding offers the first empirical evidence that density dependence responds in different ways to spatial heterogeneity and temporal heterogeneity.