Amy J. Davis

An integrated modeling approach to estimating Gunnison sage-grouse population dynamics: combining index and demographic data

Evaluation of population dynamics for rare and declining species is often limited to data that are sparse and/or of poor quality. Frequently, the best data available for rare bird species are based on large-scale, population count data. These data are commonly based on sampling methods that lack consistent sampling effort, do not account for detectability, and are complicated by observer bias. For some species, short-term studies of demographic rates have been conducted as well, but the data from such studies are typically analyzed separately. To utilize the strengths and minimize the weaknesses of these two data types, we developed a novel Bayesian integrated model that links population count data and population demographic data through population growth rate (λ) for Gunnison sage-grouse (Centrocercus minimus). The long-term population index data available for Gunnison sage-grouse are annual (years 1953–2012) male lek counts. An intensive demographic study was also conducted from years 2005 to 2010. We were able to reduce the variability in expected population growth rates across time, while correcting for potential small sample size bias in the demographic data. We found the population of Gunnison sage-grouse to be variable and slightly declining over the past 16 years.

Potential effects of incorporating fertility control into typical culling regimes in wild pig populations

Effective management of widespread invasive species such as wild pigs (Sus scrofa) is limited by resources available to devote to the effort. Better insight of the effectiveness of different management strategies on population dynamics is important for guiding decisions of resource allocation over space and time. Using a dynamic population model, we quantified effects of culling intensities and time between culling events on population dynamics of wild pigs in the USA using empirical culling patterns and data-based demographic parameters. In simulated populations closed to immigration, substantial population declines (50–100%) occurred within 4 years when 20–60% of the population was culled annually, but when immigration from surrounding areas occurred, there was a maximum of 50% reduction, even with the maximum culling intensity of 60%. Incorporating hypothetical levels of fertility control with realistic culling intensities was most effective in reducing populations when they were closed to immigration and when intrinsic population growth rate was too high (> = 1.78) to be controlled by culling alone. However, substantial benefits from fertility control used in conjunction with culling may only occur over a narrow range of net population growth rates (i.e., where net is the result of intrinsic growth rates and culling) that varies depending on intrinsic population growth rate. The management implications are that the decision to use fertility control in conjunction with culling should rely on concurrent consideration of achievable culling intensity, underlying demographic parameters, and costs of culling and fertility control. The addition of fertility control reduced abundance substantially more than culling alone, however the effects of fertility control were weaker than in populations without immigration. Because these populations were not being reduced substantially by culling alone, fertility control could be an especially helpful enhancement to culling for reducing abundance to target levels in areas where immigration can’t be prevented.

Nest Success of Gunnison Sage-Grouse in Colorado, USA.

Gunnison Sage-Grouse (Centrocercus minimus) is a species of concern for which little demographic information exists. To help fill this information gap, we investigated factors affecting nest success in two populations of Gunnison Sage-Grouse. We assessed the relative effects of (1) vegetation characteristics (e.g., shrub height, shrub cover, grass cover, and grass height), (2) temporal factors (e.g., year, timing of incubation initiation, and nest age), (3) precipitation, and (4) age of the nesting female (yearling or adult) on nest success rates. We found 177 nests in the Gunnison Basin population (that contains 85–90% of the species) from 2005–2010 and 20 nests in the San Miguel population (that contains < 10% of the species) from 2007–2010. Temporal factors had the greatest impact on nest success compared to vegetation characteristics, precipitation, and female age. Nest success varied considerably among years ranging from 4.0%-60.2% in Gunnison Basin and from 12.9%- 51.9% in San Miguel. Nests that were initiated earlier in the breeding season had higher nest success (at least one egg hatches). Daily nest survival rates decreased during the course of incubation. None of the vegetation characteristics we examined were strongly related to nest success.

Detection and persistence of environmental DNA from an invasive, terrestrial mammal

Abstract Invasive Sus scrofa, a species commonly referred to as wild pig or feral swine, is a destructive invasive species with a rapidly expanding distribution across the United States. We used artificial wallows and small waterers to determine the minimum amount of time needed for pig eDNA to accumulate in the water source to a detectable level. We removed water from the artificial wallows and tested eDNA detection over the course of 2 weeks to understand eDNA persistence. We show that our method is sensitive enough to detect very low quantities of eDNA shed by a terrestrial mammal that has limited interaction with water. Our experiments suggest that the number of individuals shedding into a water system can affect persistence of eDNA. Use of an eDNA detection technique can benefit management efforts by providing a sensitive method for finding even small numbers of individuals that may be elusive using other methods.

Declining recruitment of Gunnison Sage-Grouse highlights the need to monitor juvenile survival

ABSTRACT Recruitment of juveniles is an important vital rate that influences population growth and is fundamental to understanding trends in population size. Estimates of recruitment are often focused on the period just after hatching (prefledgling stage), which is typically the lowest survival period and often the most variable. Few studies examine true recruitment—survival from hatching to entering the breeding population—although this information is more relevant to understanding population trends. We studied the recruitment of Gunnison Sage-Grouse (Centrocercus minimus), a federally threatened species in the U.S., to examine the relative importance of chick and juvenile survival to recruitment patterns. We evaluated recruitment from 2005 to 2010 by combining separate estimates of chick survival (hatching to 30 days of age) and juvenile survival (31 days of age to the start of the first breeding season). To explain variation in these survival rates, we examined the effects of population, individual (i.e. age), and temporal (within-year and among-year differences) factors associated with recruitment of Gunnison Sage-Grouse. The factors that most explained juvenile survival rates were temporal (among-year trends and within-year seasonal effects). Chick survival rates varied by population, and daily chick survival increased with chick age. We found a slight negative trend in chick survival and a strong negative trend in juvenile survival from 2005 to 2010. The overall recruitment rate declined from 0.32 (± 0.09 SE) in 2005 to 0.04 (± 0.03 SE) in 2010. This decline coincided with a decline observed in population index data, which was not reflected in other demographic data. If survival had not been monitored past 30 days of age, estimates of recruitment would have remained relatively stable. This work highlights the importance of monitoring juvenile survival, as it may influence population dynamics.

Accounting for location uncertainty in azimuthal telemetry data improves ecological inference

BackgroundCharacterizing animal space use is critical for understanding ecological relationships. Animal telemetry technology has revolutionized the fields of ecology and conservation biology by providing high quality spatial data on animal movement. Radio-telemetry with very high frequency (VHF) radio signals continues to be a useful technology because of its low cost, miniaturization, and low battery requirements. Despite a number of statistical developments synthetically integrating animal location estimation and uncertainty with spatial process models using satellite telemetry data, we are unaware of similar developments for azimuthal telemetry data. As such, there are few statistical options to handle these unique data and no synthetic framework for modeling animal location uncertainty and accounting for it in ecological models.We developed a hierarchical modeling framework to provide robust animal location estimates from one or more intersecting or non-intersecting azimuths. We used our azimuthal telemetry model (ATM) to account for azimuthal uncertainty with covariates and propagate location uncertainty into spatial ecological models. We evaluate the ATM with commonly used estimators (Lenth (1981) maximum likelihood and M-Estimators) using simulation. We also provide illustrative empirical examples, demonstrating the impact of ignoring location uncertainty within home range and resource selection analyses. We further use simulation to better understand the relationship among location uncertainty, spatial covariate autocorrelation, and resource selection inference.ResultsWe found the ATM to have good performance in estimating locations and the only model that has appropriate measures of coverage. Ignoring animal location uncertainty when estimating resource selection or home ranges can have pernicious effects on ecological inference. Home range estimates can be overly confident and conservative when ignoring location uncertainty and resource selection coefficients can lead to incorrect inference and over confidence in the magnitude of selection. Furthermore, our simulation study clarified that incorporating location uncertainty helps reduce bias in resource selection coefficients across all levels of covariate spatial autocorrelation.ConclusionThe ATM can accommodate one or more azimuths when estimating animal locations, regardless of how they intersect; this ensures that all data collected are used for ecological inference. Our findings and model development have important implications for interpreting historical analyses using this type of data and the future design of radio-telemetry studies.

Quantifying site-level usage and certainty of absence for an invasive species through occupancy analysis of camera-trap data

Efficient implementation of management programs for invasive species depends on accurate surveillance for guiding prioritization of surveillance and control resources in space and time. Occupancy probabilities can be used to determine where surveillance should occur. Conversely, knowledge of the certainty of site-level absence is of special interest in situations where the objective is to completely remove populations despite substantial risk of re-invasion. Indeed, the decision to shift from emphasizing control activities over the full range to emphasizing reinvasion prevention, surveillance, and response near the borders, depends on accurate knowledge of absence across space. We used a dynamic occupancy model to monitor changes in the distribution of an invasive species, feral swine (Sus scrofa), based on camera-trap data collected as part of a management program from June 2014 to January 2016 in San Diego County, California. Site usage of feral swine declined overall. The most informative predictors of site usage were spatial (latitude and longitude). Site-level non-usage rates increased over time and in response to management removal efforts; and site-level usage rates were heavily impacted by having neighboring sites that were used. Combining the detection probability estimated from the occupancy model and Bayes Theorem, we demonstrated how certainty of local (site-level) absence can be estimated iteratively in time in areas with negative surveillance (no detections) data. Our framework provides a means for using management-based surveillance data to quantify certainty of site-level absence of an invasive species, allowing for adaptive prioritization of surveillance and control resources. Our approach is flexible for application to other species and types of surveillance (e.g., track-plates, eDNA).

Predicting spatial spread of rabies in skunk populations using surveillance data reported by the public

Background Prevention and control of wildlife disease invasions relies on the ability to predict spatio-temporal dynamics and understand the role of factors driving spread rates, such as seasonality and transmission distance. Passive disease surveillance (i.e., case reports by public) is a common method of monitoring emergence of wildlife diseases, but can be challenging to interpret due to spatial biases and limitations in data quantity and quality. Methodology/Principal findings We obtained passive rabies surveillance data from dead striped skunks (Mephitis mephitis) in an epizootic in northern Colorado, USA. We developed a dynamic patch-occupancy model which predicts spatio-temporal spreading while accounting for heterogeneous sampling. We estimated the distance travelled per transmission event, direction of invasion, rate of spatial spread, and effects of infection density and season. We also estimated mean transmission distance and rates of spatial spread using a phylogeographic approach on a subsample of viral sequences from the same epizootic. Both the occupancy and phylogeographic approaches predicted similar rates of spatio-temporal spread. Estimated mean transmission distances were 2.3 km (95% Highest Posterior Density (HPD95): 0.02, 11.9; phylogeographic) and 3.9 km (95% credible intervals (CI95): 1.4, 11.3; occupancy). Estimated rates of spatial spread in km/year were: 29.8 (HPD95: 20.8, 39.8; phylogeographic, branch velocity, homogenous model), 22.6 (HPD95: 15.3, 29.7; phylogeographic, diffusion rate, homogenous model) and 21.1 (CI95: 16.7, 25.5; occupancy). Initial colonization probability was twice as high in spring relative to fall. Conclusions/Significance Skunk-to-skunk transmission was primarily local (< 4 km) suggesting that if interventions were needed, they could be applied at the wave front. Slower viral invasions of skunk rabies in western USA compared to a similar epizootic in raccoons in the eastern USA implies host species or landscape factors underlie the dynamics of rabies invasions. Our framework provides a straightforward method for estimating rates of spatial spread of wildlife diseases.