A cloudy forecast for species distribution models: Predictive uncertainties abound for California birds after a century of climate and land-use change.

Correlative species distribution models are widely used to quantify past shifts in ranges or communities, and to predict future outcomes under ongoing global change. Practitioners confront a wide range of potentially plausible models for ecological dynamics, but most specific applications only consider a narrow set. Here, we clarify that certain model structures can embed restrictive assumptions about key sources of forecast uncertainty into an analysis.

Estimating spatially variable and density-dependent survival using open-population spatial capture-recapture models.

Open-population spatial capture-recapture (OPSCR) models use the spatial information contained in individual detections collected over multiple consecutive occasions to estimate not only occasion-specific density, but also demographic parameters. OPSCR models can also estimate spatial variation in vital rates, but such models are neither widely used nor thoroughly tested. We developed a Bayesian OPSCR model that not only accounts for spatial variation in survival using spatial covariates but also estimates local density-dependent effects on survival within a unified framework.

A flexible and efficient Bayesian implementation of point process models for spatial capture-recapture data.

Spatial capture-recapture (SCR) is now routinely used for estimating abundance and density of wildlife populations. A standard SCR model includes sub-models for the distribution of individual activity centers (ACs) and for individual detections conditional on the locations of these ACs. Both sub-models can be expressed as point processes taking place in continuous space, but there is a lack of accessible and efficient tools to fit such models in a Bayesian paradigm.

Comparing N-mixture models and GLMMs for relative abundance estimation in a citizen science dataset.

To analyze species count data when detection is imperfect, ecologists need models to estimate relative abundance in the presence of unknown sources of heterogeneity. Two candidate models are generalized linear mixed models (GLMMs) and hierarchical N-mixture models. GLMMs are computationally robust but do not explicitly separate detection from abundance patterns.

Identifying engaging bird species and traits with community science observations.

Identifying rates at which birders engage with different species can inform the impact and efficacy of conservation outreach and the scientific use of community-collected biodiversity data. Species that are thought to be “charismatic” are often prioritized in conservation, and previous researchers have used sociological experiments and digital records to estimate charisma indirectly. In this study, we take advantage of community science efforts as another record of human engagement with animals that can reveal observer biases directly, which are in part driven by observer preference.

Building integral projection models with nonindependent vital rates.

Population dynamics are functions of several demographic processes including survival, reproduction, somatic growth, and maturation. The rates or probabilities for these processes can vary by time, by location, and by individual. These processes can co-vary and interact to varying degrees, e.

To rarefy or not to rarefy: robustness and efficiency trade-offs of rarefying microbiome data.

Motivation: Microbiome datasets provide rich information about microbial communities. However, vast library size variations across samples present great challenges for proper statistical comparisons. To deal with these challenges, rarefaction is often used in practice as a normalization technique, although there has been debate whether rarefaction should ever be used.

Multi-scale trend analysis of water quality using error propagation of generalized additive models.

Effective stewardship of ecosystems to sustain current ecological status or mitigate impacts requires nuanced understanding of how conditions have changed over time in response to anthropogenic pressures and natural variability. Detecting and appropriately characterizing changes requires accurate and flexible trend assessment methods that can be readily applied to environmental monitoring datasets. A key requirement is complete propagation of uncertainty through the analysis.

Estimating and forecasting spatial population dynamics of apex predators using transnational genetic monitoring.

The ongoing recovery of terrestrial large carnivores in North America and Europe is accompanied by intense controversy. On the one hand, reestablishment of large carnivores entails a recovery of their most important ecological role, predation. On the other hand, societies are struggling to relearn how to live with apex predators that kill livestock, compete for game species, and occasionally injure or kill people.

One size does not fit all: Customizing MCMC methods for hierarchical models using NIMBLE.

Improved efficiency of Markov chain Monte Carlo facilitates all aspects of statistical analysis with Bayesian hierarchical models. Identifying strategies to improve MCMC performance is becoming increasingly crucial as the complexity of models, and the run times to fit them, increases. We evaluate different strategies for improving MCMC efficiency using the open-source software NIMBLE (R package nimble) using common ecological models of species occurrence and abundance as examples.

A pathway for multivariate analysis of ecological communities using copulas.

We describe a new pathway for multivariate analysis of data consisting of counts of species abundances that includes two key components: copulas, to provide a flexible joint model of individual species, and dissimilarity-based methods, to integrate information across species and provide a holistic view of the community. Individual species are characterized using suitable (marginal) statistical distributions, with the mean, the degree of over-dispersion, and/or zero-inflation being allowed to vary among a priori groups of sampling units. Associations among species are then modeled using copulas, which allow any pair of disparate types of variables to be coupled through their cumulative distribution function, while maintaining entirely the separate individual marginal distributions appropriate for each species.

Proximity of restored hedgerows interacts with local floral diversity and species' traits to shape long-term pollinator metacommunity dynamics.

Disconnected habitat fragments are poor at supporting population and community persistence; restoration ecologists, therefore, advocate for the establishment of habitat networks across landscapes. Few empirical studies, however, have considered how networks of restored habitat patches affect metacommunity dynamics. Here, using a 10-year study on restored hedgerows and unrestored field margins within an intensive agricultural landscape, we integrate occupancy modelling with network theory to examine the interaction between local and landscape characteristics, habitat selection and dispersal in shaping pollinator metacommunity dynamics.

A century of climate and land-use change cause species turnover without loss of beta diversity in California's Central Valley.

Climate and land-use changes are thought to be the greatest threats to biodiversity, but few studies have directly measured their simultaneous impacts on species distributions. We used a unique historic resource-early 20th-century bird surveys conducted by Joseph Grinnell and colleagues-paired with contemporary resurveys a century later to examine changes in bird distributions in California's Central Valley, one of the most intensively modified agricultural zones in the world and a region of heterogeneous climate change. We analyzed species- and community-level occupancy using multispecies occupancy models that explicitly accounted for imperfect detection probability, and developed a novel, simulation-based method to compare the relative influences of climate and land-use covariates on site-level species richness and beta diversity (measured by Jaccard similarity).

Commentary to: a cross-validation-based approach for delimiting reliable home range estimates.

Background: Continued exploration of the performance of the recently proposed cross-validation-based approach for delimiting home ranges using the Time Local Convex Hull (T-LoCoH) method has revealed a number of issues with the original formulation.

Main Text: Here we replace the ad hoc cross-validation score with a new formulation based on the total log probability of out-of-sample predictions. To obtain these probabilities, we interpret the normalized LoCoH hulls as a probability density.

Origin and maintenance of chemical diversity in a species-rich tropical tree lineage.

Plant secondary metabolites play important ecological and evolutionary roles, most notably in the deterrence of natural enemies. The classical theory explaining the evolution of plant chemical diversity is that new defences arise through a pairwise co-evolutionary arms race between plants and their specialized natural enemies. However, plant species are bombarded by dozens of different herbivore taxa from disparate phylogenetic lineages that span a wide range of feeding strategies and have distinctive physiological constraints that interact differently with particular plant metabolites.

Phylogenetic occupancy models integrate imperfect detection and phylogenetic signal to analyze community structure.

Biological communities are structured phylogenetically-closely related species are typically more likely to be found at the same sites. This may be, in part, because they respond similarly to environmental gradients. Accurately surveying biological communities is, however, made difficult by the fact that detection of species is not perfect.

Warming and provenance limit tree recruitment across and beyond the elevation range of subalpine forest.

Climate niche models project that subalpine forest ranges will extend upslope with climate warming. These projections assume that the climate suitable for adult trees will be adequate for forest regeneration, ignoring climate requirements for seedling recruitment, a potential demographic bottleneck. Moreover, local genetic adaptation is expected to facilitate range expansion, with tree populations at the upper forest edge providing the seed best adapted to the alpine.

Monte Carlo estimation of stage structured development from cohort data.

Cohort data are frequently collected to study stage-structured development and mortalities of many organisms, particularly arthropods. Such data can provide information on mean stage durations, among-individual variation in stage durations, and on mortality rates. Current statistical methods for cohort data lack flexibility in the specification of stage duration distributions and mortality rates.

Diversification practices reduce organic to conventional yield gap.

Agriculture today places great strains on biodiversity, soils, water and the atmosphere, and these strains will be exacerbated if current trends in population growth, meat and energy consumption, and food waste continue. Thus, farming systems that are both highly productive and minimize environmental harms are critically needed. How organic agriculture may contribute to world food production has been subject to vigorous debate over the past decade.

Statistical ecology comes of age.

The desire to predict the consequences of global environmental change has been the driver towards more realistic models embracing the variability and uncertainties inherent in ecology. Statistical ecology has gelled over the past decade as a discipline that moves away from describing patterns towards modelling the ecological processes that generate these patterns. Following the fourth International Statistical Ecology Conference (1-4 July 2014) in Montpellier, France, we analyse current trends in statistical ecology.

An approximate Bayesian computation approach to parameter estimation in a stochastic stage-structured population model.

Complex population processes may require equally complex models, which can lead to analytically intractable estimation problems. Approximate Bayesian computation (ABC) is a computational tool for parameter estimation in situations where likelihoods cannot be computed. Instead of using likelihoods, ABC methods quantify the similarities between an observed data set and repeated simulations from a model.

The importance of individual developmental variation in stage-structured population models.

Population stage structure is fundamental to ecology, and models of this structure have proven useful in many different systems. Many ecological variables other than stage, such as habitat type, site occupancy and metapopulation status are also modelled using transitions among discrete states. Transitions among life stages can be characterised by the distribution of time spent in each stage, including the mean and variance of each stage duration and within-individual correlations among multiple stage durations.