Search

Optimal design procedures provide a framework to leverage the learning generated by ecological models to flexibly and efficiently deploy future monitoring efforts. At the same time, Bayesian hierarchical models have become widespread in ecology and offer a rich set of tools for ecological learning and inference. However, coupling these methods with an optimal design framework can become computationally intractable. Recursive Bayesian computation offers a way to substantially reduce this computational burden, making optimal design accessible for modern Bayesian ecological models. We demonstrate the application of so-called prior-proposal recursive Bayes to optimal design using a simulated data binary regression and the real-world example of monitoring and modeling sea otters in Glacier Bay, Alaska. These examples highlight the computational gains offered by recursive Bayesian methods and the tighter fusion of monitoring and science that those computational gains enable., (© 2021 The Ecological Society of America. This article has been contributed to by US Government employees and their work is in the public domain in the USA.)

Accurate assessment of abundance forms a central challenge in population ecology and wildlife management. Many statistical techniques have been developed to estimate population sizes because populations change over time and space and to correct for the bias resulting from animals that are present in a study area but not observed. The mobility of individuals makes it difficult to design sampling procedures that account for movement into and out of areas with fixed jurisdictional boundaries. Aerial surveys are the gold standard used to obtain data of large mobile species in geographic regions with harsh terrain, but these surveys can be prohibitively expensive and dangerous. Estimating abundance with ground-based census methods have practical advantages, but it can be difficult to simultaneously account for temporary emigration and observer error to avoid biased results. Contemporary research in population ecology increasingly relies on telemetry observations of the states and locations of individuals to gain insight on vital rates, animal movements, and population abundance. Analytical models that use observations of movements to improve estimates of abundance have not been developed. Here we build upon existing multi-state mark-recapture methods using a hierarchical N-mixture model with multiple sources of data, including telemetry data on locations of individuals, to improve estimates of population sizes. We used a state-space approach to model animal movements to approximate the number of marked animals present within the study area at any observation period, thereby accounting for a frequently changing number of marked individuals. We illustrate the approach using data on a population of elk (Cervus elaphus nelsoni) in Northern Colorado, USA. We demonstrate substantial improvement compared to existing abundance estimation methods and corroborate our results from the ground based surveys with estimates from aerial surveys during the same seasons. We develop a hierarchical Bayesian N-mixture model using multiple sources of data on abundance, movement and survival to estimate the population size of a mobile species that uses remote conservation areas. The model improves accuracy of inference relative to previous methods for estimating abundance of open populations., (© 2018 by the Ecological Society of America.)

Population dynamics vary in space and time. Survey designs that ignore these dynamics may be inefficient and fail to capture essential spatio-temporal variability of a process. Alternatively, dynamic survey designs explicitly incorporate knowledge of ecological processes, the associated uncertainty in those processes, and can be optimized with respect to monitoring objectives. We describe a cohesive framework for monitoring a spreading population that explicitly links animal movement models with survey design and monitoring objectives. We apply the framework to develop an optimal survey design for sea otters in Glacier Bay. Sea otters were first detected in Glacier Bay in 1988 and have since increased in both abundance and distribution; abundance estimates increased from 5 otters to >5,000 otters, and they have spread faster than 2.7 km/yr. By explicitly linking animal movement models and survey design, we are able to reduce uncertainty associated with forecasting occupancy, abundance, and distribution compared to other potential random designs. The framework we describe is general, and we outline steps to applying it to novel systems and taxa., (© 2018 by the Ecological Society of America.)

Analyzing ecological data often requires modeling the autocorrelation created by spatial and temporal processes. Many seemingly disparate statistical methods used to account for autocorrelation can be expressed as regression models that include basis functions. Basis functions also enable ecologists to modify a wide range of existing ecological models in order to account for autocorrelation, which can improve inference and predictive accuracy. Furthermore, understanding the properties of basis functions is essential for evaluating the fit of spatial or time-series models, detecting a hidden form of collinearity, and analyzing large data sets. We present important concepts and properties related to basis functions and illustrate several tools and techniques ecologists can use when modeling autocorrelation in ecological data., (© 2016 by the Ecological Society of America.)

Ecological invasions and colonizations occur dynamically through space and time. Estimating the distribution and abundance of colonizing species is critical for efficient management or conservation. We describe a statistical framework for simultaneously estimating spatiotemporal occupancy and abundance dynamics of a colonizing species. Our method accounts for several issues that are common when modeling spatiotemporal ecological data including multiple levels of detection probability, multiple data sources, and computational limitations that occur when making fine-scale inference over a large spatiotemporal domain. We apply the model to estimate the colonization dynamics of sea otters (Enhydra lutris) in Glacier Bay, in southeastern Alaska., (© 2016 by the Ecological Society of America.)

Satellite telemetry devices collect valuable information concerning the sites visited by animals, including the location of central places like dens, nests, rookeries, or haul-outs. Existing methods for estimating the location of central places from telemetry data require user-specified thresholds and ignore common nuances like measurement error. We present a fully model-based approach for locating central places from telemetry data that accounts for multiple sources of uncertainty and uses all of the available locational data. Our general framework consists of an observation model to account for large telemetry measurement error and animal movement, and a highly flexible mixture model specified using a Dirichlet process to identify the location of central places. We also quantify temporal patterns in central place use by incorporating ancillary behavioral data into the model; however, our framework is also suitable when no such behavioral data exist. We apply the model to a simulated data set as proof of concept. We then illustrate our framework by analyzing an Argos satellite telemetry data set on harbor seals (Phoca vitulina) in the Gulf of Alaska, a species that exhibits fidelity to terrestrial haul-out sites., (© 2016 by the Ecological Society of America.)

Ecological processes operate across temporal and spatial scales. Anthropogenic disturbances impact these processes, but examinations of scale dependence in impacts are infrequent. Such examinations can provide important insight to wildlife-human interactions and guide management efforts to reduce impacts. We assessed spatiotemporal scale dependence in habitat selection of mule deer (Odocoileus hemionus) in the Piceance Basin of Colorado, USA, an area of ongoing natural gas development. We employed a newly developed animal movement method to assess habitat selection across scales defined using animal-centric spatiotemporal definitions ranging from the local (defined from five hour movements) to the broad (defined from weekly movements). We extended our analysis to examine variation in scale dependence between night and day and assess functional responses in habitat selection patterns relative to the density of anthropogenic features. Mule deer displayed scale invariance in the direction of their response to energy development features, avoiding well pads and the areas closest to roads at all scales, though with increasing strength of avoidance at coarser scales. Deer displayed scale-dependent responses to most other habitat features, including land cover type and habitat edges. Selection differed between night and day at the finest scales, but homogenized as scale increased. Deer displayed functional responses to development, with deer inhabiting the least developed ranges more strongly avoiding development relative to those with more development in their ranges. Energy development was a primary driver of habitat selection patterns in mule deer, structuring their behaviors across all scales examined. Stronger avoidance at coarser scales suggests that deer behaviorally mediated their interaction with development, but only to a degree. At higher development densities than seen in this area, such mediation may not be possible and thus maintenance of sufficient habitat with lower development densities will be a critical best management practice as development expands globally., (© 2016 by the Ecological Society of America.)

Evaluation of the progress of management programs for invasive species is crucial for demonstrating impacts to stakeholders and strategic planning of resource allocation. Estimates of abundance before and after management activities can serve as a useful metric of population management programs. However, many methods of estimating population size are too labor intensive and costly to implement, posing restrictive levels of burden on operational programs. Removal models are a reliable method for estimating abundance before and after management using data from the removal activities exclusively, thus requiring no work in addition to management. We developed a Bayesian hierarchical model to estimate abundance from removal data accounting for varying levels of effort, and used simulations to assess the conditions under which reliable population estimates are obtained. We applied this model to estimate site-specific abundance of an invasive species, feral swine (Sus scrofa), using removal data from aerial gunning in 59 site/time-frame combinations (480-19,600 acres) throughout Oklahoma and Texas, USA. Simulations showed that abundance estimates were generally accurate when effective removal rates (removal rate accounting for total effort) were above 0.40. However, when abundances were small (<50) the effective removal rate needed to accurately estimates abundances was considerably higher (0.70). Based on our post-validation method, 78% of our site/time frame estimates were accurate. To use this modeling framework it is important to have multiple removals (more than three) within a time frame during which demographic changes are minimized (i.e., a closed population; ≤3 months for feral swine). Our results show that the probability of accurately estimating abundance from this model improves with increased sampling effort (8+ flight hours across the 3-month window is best) and increased removal rate. Based on the inverse relationship between inaccurate abundances and inaccurate removal rates, we suggest auxiliary information that could be collected and included in the model as covariates (e.g., habitat effects, differences between pilots) to improve accuracy of removal rates and hence abundance estimates., (© 2016 by the Ecological Society of America.)

Statistical decision theory (SDT) is a sub-field of decision theory that formally incorporates statistical investigation into a decision-theoretic framework to account for uncertainties in a decision problem. SDT provides a unifying analysis of three types of information: statistical results from a data set, knowledge of the consequences of potential choices (i.e., loss), and prior beliefs about a system. SDT links the theoretical development of a large body of statistical methods, including point estimation, hypothesis testing, and confidence interval estimation. The theory and application of SDT have mainly been developed and published in the fields of mathematics, statistics, operations research, and other decision sciences, but have had limited exposure in ecology. Thus, we provide an introduction to SDT for ecologists and describe its utility for linking the conventionally separate tasks of statistical investigation and decision making in a single framework. We describe the basic framework of both Bayesian and frequentist SDT, its traditional use in statistics, and discuss its application to decision problems that occur in ecology. We demonstrate SDT with two types of decisions: Bayesian point estimation and an applied management problem of selecting a prescribed fire rotation for managing a grassland bird species. Central to SDT, and decision theory in general, are loss functions. Thus, we also provide basic guidance and references for constructing loss functions for an SDT problem., (© 2016 by the Ecological Society of America.)

While multi-species occupancy models (MSOMs) are emerging as a popular method for analyzing biodiversity data, formal checking and validation approaches for this class of models have lagged behind. Concurrent with the rise in application of MSOMs among ecologists, a quiet regime shift is occurring in Bayesian statistics where predictive model comparison approaches are experiencing a resurgence. Unlike single-species occupancy models that use integrated likelihoods, MSOMs are usually couched in a Bayesian framework and contain multiple levels. Standard model checking and selection methods are often unreliable in this setting and there is only limited guidance in the ecological literature for this class of models. We examined several different contemporary Bayesian hierarchical approaches for checking and validating MSOMs and applied these methods to a freshwater aquatic study system in Colorado, USA, to better understand the diversity and distributions of plains fishes. Our findings indicated distinct differences among model selection approaches, with cross-validation techniques performing the best in terms of prediction., (© 2016 by the Ecological Society of America.)

The dynamic, multi-season occupancy model framework has become a popular tool for modeling open populations with occupancies that change over time through local colonizations and extinctions. However, few versions of the model relate these probabilities to the occupancies of neighboring sites or patches. We present a modeling framework that incorporates this information and is capable of describing a wide variety of spatiotemporal colonization and extinction processes. A key feature of the model is that it is based on a simple set of small-scale rules describing how the process evolves. The result is a dynamic process that can account for complicated large-scale features. In our model, a site is more likely to be colonized if more of its neighbors were previously occupied and if it provides more appealing environmental characteristics than its neighboring sites. Additionally, a site without occupied neighbors may also become colonized through the inclusion of a long-distance dispersal process. Although similar model specifications have been developed for epidemiological applications, ours formally accounts for detectability using the well-known occupancy modeling framework. After demonstrating the viability and potential of this new form of dynamic occupancy model in a simulation study, we use it to obtain inference for the ongoing Common Myna (Acridotheres tristis) invasion in South Africa. Our results suggest that the Common Myna continues to enlarge its distribution and its spread via short distance movement, rather than long-distance dispersal. Overall, this new modeling framework provides a powerful tool for managers examining the drivers of colonization including short- vs. long-distance dispersal, habitat quality, and distance from source populations.

Multiple factors complicate the analysis of animal telemetry location data. Recent advancements address issues such as temporal autocorrelation and telemetry measurement error, but additional challenges remain. Difficulties introduced by complicated error structures or barriers to animal movement can weaken inference. We propose an approach for obtaining resource selection inference from animal location data that accounts for complicated error structures, movement constraints, and temporally autocorrelated observations. We specify a model for telemetry data observed with error conditional on unobserved true locations that reflects prior knowledge about constraints in the animal movement process. The observed telemetry data are modeled using a flexible distribution that accommodates extreme errors and complicated error structures. Although constraints to movement are often viewed as a nuisance, we use constraints to simultaneously estimate and account for telemetry error. We apply the model to simulated data, showing that it outperforms common ad hoc approaches used when confronted with measurement error and movement constraints. We then apply our framework to an Argos satellite telemetry data set on harbor seals (Phoca vitulina) in the Gulf of Alaska, a species that is constrained to move within the marine environment and adjacent coastlines.

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.

Habitat selection is a fundamental aspect of animal ecology, the understanding of which is critical to management and conservation. Global positioning system data from animals allow fine-scale assessments of habitat selection and typically are analyzed in a use availability framework, whereby animal locations are contrasted with random locations (the availability sample). Although most use-availability methods are in fact spatial point process models, they often are fit using logistic regression. This framework offers numerous methodological challenges, for which the literature provides little guidance. Specifically, the size and spatial extent of the availability sample influences coefficient estimates potentially causing interpretational bias. We examined the influence of availability on statistical inference through simulations and analysis of serially correlated mule deer GPS data. Bias in estimates arose from incorrectly assessing and sampling the spatial extent of availability. Spatial autocorrelation in covariates, which is common for landscape characteristics, exacerbated the error in availability sampling leading to increased bias. These results have strong implications for habitat selection analyses using GPS data, which are increasingly prevalent in the literature. We recommend that researchers assess the sensitivity of their results to their availability sample and, where bias is likely, take care with interpretations and use cross validation to assess robustness.

Power-law relationships are among the most well-studied functional relationships in biology. Recently the common practice of fitting power laws using linear regression (LR) on log-transformed data has been criticized, calling into question the conclusions of hundreds of studies. It has been suggested that nonlinear regression (NLR) is preferable, but no rigorous comparison of these two methods has been conducted. Using Monte Carlo simulations, we demonstrate that the error distribution determines which method performs better, with NLR better characterizing data with additive, homoscedastic, normal error and LR better characterizing data with multiplicative, heteroscedastic, lognormal error. Analysis of 471 biological power laws shows that both forms of error occur in nature. While previous analyses based on log-transformation appear to be generally valid, future analyses should choose methods based on a combination of biological plausibility and analysis of the error distribution. We provide detailed guidelines and associated computer code for doing so, including a model averaging approach for cases where the error structure is uncertain.

Ecological spatial data often come from multiple sources, varying in extent and accuracy. We describe a general approach to reconciling such data sets through the use of the Bayesian hierarchical framework. This approach provides a way for the data sets to borrow strength from one another while allowing for inference on the underlying ecological process. We apply this approach to study the incidence of eastern spruce dwarf mistletoe (Arceuthobium pusillum) in Minnesota black spruce (Picea mariana). A Minnesota Department of Natural Resources operational inventory of black spruce stands in northern Minnesota found mistletoe in 11% of surveyed stands, while a small, specific-pest survey found mistletoe in 56% of the surveyed stands. We reconcile these two surveys within a Bayesian hierarchical framework and predict that 35-59% of black spruce stands in northern Minnesota are infested with dwarf mistletoe.

Climate change could alter the population growth of dominant species, leading to profound effects on community structure and ecosystem dynamics. Understanding the links between historical variation in climate and population vital rates (survival, growth, recruitment) is one way to predict the impact of future climate change. Using a unique, long-term data set from eastern Idaho, USA, we parameterized integral projection models (IPMs) for Pseudoroegneria spicata, Hesperostipa comata, and Artemisia tripartita to identify the demographic rates and climate variables most important for population growth. We described survival, growth, and recruitment as a function of genet size using mixed-effect regression models that incorporated climate variables. Elasticites for the survival + growth portion of the kernel were larger than the recruitment portion for all three species, with survival + growth accounting for 87-95% of the total elasticity. The genet sizes with the highest elasticity values in each species were very close to the genet size threshold where survival approached 100%. We found strong effects of climate on the population growth rate of two of our three species. In H. comata, a 1% decrease in previous year's precipitation would lead to a 0.6% decrease in population growth. In A. tripartita, a 1% increase in summer temperature would result in a 1.3% increase in population growth. In both H. comata and A. tripartita, climate influenced population growth by affecting genet growth more than survival or recruitment. Late-winter snow was the most important climate variable for P. spicata, but its effect on population growth was smaller than the climate effects we found in H. comata or A. tripartita. For all three species, demographic responses lagged climate by at least one year. Our analysis indicates that understanding climate effects on genet growth may be crucial for anticipating future changes in the structure and function of sagebrush steppe vegetation.