No abstract

An understanding of species relationships is critical in the management and conservation of populations facing climate change, yet few studies address how climate alters species interactions and other population drivers. We use a long-term, broad-scale data set of relative abundance to examine the influence of climate, predators, and density dependence on the population dynamics of declining scaup (Aythya) species within the core of their breeding range. The state-space modeling approach we use applies to a wide range of wildlife species, especially populations monitored over broad spatiotemporal extents. Using this approach, we found that immediate snow cover extent in the preceding winter and spring had the strongest effects, with increases in mean snow cover extent having a positive effect on the local surveyed abundance of scaup. The direct effects of mesopredator abundance on scaup population dynamics were weaker, but the results still indicated a potentil interactive process between climate and food web dynamics (mesopredators, alternative prey, and scaup). By considering climate variables and other potential effects on population dynamics, and using a rigorous estimation framework, we provide insight into complex ecological processes for guiding. conservation and policy actions aimed at mitigating and reversing the decline of scaup.

A common goal in ecology and wildlife management is to determine the causes of variation in population dynamics over long periods of time and across large spatial scales. Many assumptions must nevertheless be overcome to make appropriate inference about spatio-temporal variation in population dynamics, such as autocorrelation among data points, excess zeros, and observation error in count data. To address these issues, many scientists and statisticians have recommended the use of Bayesian hierarchical models. Unfortunately, hierarchical statistical models remain somewhat difficult to use because of the necessary quantitative background needed to implement them, or because of the computational demands of using Markov Chain Monte Carlo algorithms to estimate parameters. Fortunately, new tools have recently been developed that make it more feasible for wildlife biologists to fit sophisticated hierarchical Bayesian models (i.e., Integrated Nested Laplace Approximation, ‘INLA’). We present a case study using two important game species in North America, the lesser and greater scaup, to demonstrate how INLA can be used to estimate the parameters in a hierarchical model that decouples observation error from process variation, and accounts for unknown sources of excess zeros as well as spatial and temporal dependence in the data. Ultimately, our goal was to make unbiased inference about spatial variation in population trends over time.

No abstract

Managing for species using current weather patterns fails to incorporate the uncertainty associated with future climatic conditions; without incorporating potential changes in climate into conservation strategies, management and conservation efforts may fall short or waste valuable resources. Understanding the effects of climate change on species in the Great Plains of North America is especially important, as this region is projected to experience an increased magnitude of climate change. Of particular ecological and conservation interest is the lesser prairie‐chicken (Tympanuchus pallidicinctus), which was listed as “threatened” under the U.S. Endangered Species Act in May 2014. We used Bayesian hierarchical models to quantify the effects of extreme climatic events (extreme values of the Palmer Drought Severity Index [PDSI]) relative to intermediate (changes in El Niño Southern Oscillation) and long‐term climate variability (changes in the Pacific Decadal Oscillation) on trends in lesser prairie‐chicken abundance from 1981 to 2014. Our results indicate that lesser prairie‐chicken abundance on leks responded to environmental conditions of the year previous by positively responding to wet springs (high PDSI) and negatively to years with hot, dry summers (low PDSI), but had little response to variation in the El Niño Southern Oscillation and the Pacific Decadal Oscillation. Additionally, greater variation in abundance on leks was explained by variation in site relative to broad‐scale climatic indices. Consequently, lesser prairie‐chicken abundance on leks in Kansas is more strongly influenced by extreme drought events during summer than other climatic conditions, which may have negative consequences for the population as drought conditions intensify throughout the Great Plains.

Covariance matrix estimation arises in multivariate problems including multivariate normal sampling models and regression models where random effects are jointly modeled, e.g. random-intercept, random-slope models. A Bayesian analysis of these problems requires a prior on the covariance matrix. Here we compare an inverse Wishart, scaled inverse Wishart, hierarchical inverse Wishart, and a separation strategy as possible priors for the covariance matrix. We evaluate these priors through a simulation study and application to a real data set. Generally all priors work well with the exception of the inverse Wishart when the true variance is small relative to prior mean. In this case, the posterior for the variance is biased toward larger values and the correlation is biased toward zero. This bias persists even for large sample sizes and therefore caution should be used when using the inverse Wishart prior.

Hellbenders Cryptobranchus alleganiensis are critically imperiled amphibians throughout the eastern USA. Rock-lifting is widely used to monitor hellbenders but can severely disturb habitat. We asked whether artificial shelter occupancy (the proportion of occupied shelters in an array) would function as a proxy for hellbender abundance and thereby serve as a viable alternative to rock-lifting. We hypothesized that shelter occupancy would vary spatially in response to hellbender density, natural shelter density, or both, and would vary temporally with hellbender seasonal activity patterns and time since shelter deployment. We established shelter arrays (n = 30 shelters each) in 6 stream reaches and monitored them monthly for up to 2 yr. We used Bayesian mixed logistic regression and model ranking criteria to assess support for hypotheses concerning drivers of shelter occupancy. In all reaches, shelter occupancy was highest from June-August each year and was higher in Year 2 relative to Year 1. Our best-supported model indicated that the extent of boulder and bedrock (hereafter, natural shelter) in a reach mediated the relationship between hellbender abundance and shelter occupancy. More explicitly, shelter occupancy was positively correlated with abundance when natural shelter covered <20% of a reach, but uncorrelated with abundance when natural shelter was more abundant. While shelter occupancy should not be used to infer variation in hellbender relative abundance when substrate composition varies among reaches, we showed that artificial shelters can function as valuable monitoring tools when reaches meet certain criteria, though regular shelter maintenance is critical.

No abstract