Geology topics

Devin S. Johnson

Publications and source records attributed to Devin S. Johnson.

At least 19 recordsLinked to original sources

USGS2023

Examining the effect of environmental variability on the viability of endangered Steller sea lions using an integrated population model

Understanding spatio-temporal variability in demography and the influence of environmental conditions offers insight into the factors underlying population dynamics. This is particularly true for species with divergent demographic patterns across large geographic areas. The contrasting abundance trends observed across the range of Steller sea lions (Eumetopias jubatus) have been studied extensively, with research suggesting that the primary drivers of localized population dynamics vary over time and space. We developed a Bayesian integrated population model for the endangered western distinct population segment of Steller sea lions that combines mark-recapture and count data from 2000-2021 to estimate demographic rates, abundance trends, and the effects of environmental variability on population growth. Our results highlight subregional demographic differences, including reduced pup survival in the central Aleutian Islands and reduced yearling survival west of Samalga Pass. Range-wide abundance increased by 1.7% yr-1 (95% credible interval: 0.14; 3.4%) over the study period, with a positive annual growth rate of 3.0% (1.1; 5.1%) yr-1 east of Samalga Pass, a negative growth rate of -2.1% (-4.6; 0.5%) yr-1 west of Samalga Pass, and an overall low probability of local extirpation (<2%) in 100 years even in subregions experiencing continued decline. The effect of environmental variability on population growth varied depending on subpopulation size and vital rates and was strongest in the area of greatest decline. Our model improves upon existing approaches for estimating abundance, accounts for environmental variability within the viability analysis, and can facilitate evaluating the efficacy of conservation actions and progress toward recovery goals.

Alaska↗

USGS2022

Environmental drivers of demography and potential factors limiting the recovery of an endangered marine top predator

Understanding what drives changes in wildlife demography is fundamental to the conservation and management of depleted or declining populations, though making inference about the intrinsic and extrinsic factors that influence survival and reproduction remains challenging. Here we use mark–resight data from 2000 to 2018 to examine the effects of environmental variability on age-specific survival and natality for the endangered western distinct population segment (wDPS) of Steller sea lions ( Eumetopias jubatus ) in Alaska, USA. Though this population has been studied extensively over the last four decades, the causes of divergent abundance trends that have been observed across the wDPS range remain unknown. We developed a Bayesian multievent mark–resight model that accounts for female reproductive state uncertainty. Annual survival probabilities for male pups (0.44; 0.36–0.53), female yearlings (0.63; 0.49–0.73), and male yearlings (0.62; 0.51–0.71) born in the western portion of the wDPS range, estimated here for the first time, were lower than those in the eastern portion of the wDPS range, estimated as: male pups (0.69; 0.65–0.74), female yearlings (0.76; 0.71–0.81), and male yearlings (0.71; 0.65–0.78). There was a higher proportion of young female breeders in the western portion of the range, but overall natality was lower (0.69; 0.47–0.96) than in the eastern portion of the range (0.80; 0.74–0.84). Additionally, pup mass had a positive effect on pup survival in the eastern portion of the range and a negative effect in the western portion of the range, potentially due to earlier weaning of heavier pups. Local- and basin-scale oceanographic features such as the Aleutian Low, the Arctic Oscillation Index, the North Pacific Gyre Oscillation, chlorophyll concentration, upwelling, and wind in certain seasons were correlated with vital rates. However, drawing strong inferences from these correlations is challenging given that relationships between ocean conditions and an adaptive top predator in a dynamic ecosystem are exceedingly complex. This study provides the first demographic rate estimates for the western portion of the range where abundance estimates continue to decline. These results will advance efforts to identify factors driving regionally divergent abundance trends, with implications for population-level responses to future climate variability.

Alaska↗

USGS2022

Greater than the sum of its parts: Computationally flexible Bayesian hierarchical modeling

We propose a multistage method for making inference at all levels of a Bayesian hierarchical model (BHM) using natural data partitions to increase efficiency by allowing computations to take place in parallel form using software that is most appropriate for each data partition. The full hierarchical model is then approximated by the product of independent normal distributions for the data component of the model. In the second stage, the Bayesian maximum a posteriori (MAP) estimator is found by maximizing the approximated posterior density with respect to the parameters. If the parameters of the model can be represented as normally distributed random effects, then the second-stage optimization is equivalent to fitting a multivariate normal linear mixed model. We consider a third stage that updates the estimates of distinct parameters for each data partition based on the results of the second stage. The method is demonstrated with two ecological data sets and models, a generalized linear mixed effects model (GLMM) and an integrated population model (IPM). The multistage results were compared to estimates from models fit in single stages to the entire data set. In both cases, multistage results were very similar to a full MCMC analysis. Supplementary materials accompanying this paper appear online.

Journal of Agricultural, Biological and Environmen↗

USGS2021

An integrated path for spatial capture–recapture and animal movement modeling

Ecologists and conservation biologists increasingly rely on spatial capture–recapture (SCR) and movement modeling to study animal populations. Historically, SCR has focused on population-level processes (e.g., vital rates, abundance, density, and distribution), whereas animal movement modeling has focused on the behavior of individuals (e.g., activity budgets, resource selection, migration). Even though animal movement is clearly a driver of population-level patterns and dynamics, technical and conceptual developments to date have not forged a firm link between the two fields. Instead, movement modeling has typically focused on the individual level without providing a coherent scaling from individual- to population-level processes, whereas SCR has typically focused on the population level while greatly simplifying the movement processes that give rise to the observations underlying these models. In our view, the integration of SCR and animal movement modeling has tremendous potential for allowing ecologists to scale up from individuals to populations and advancing the types of inferences that can be made at the intersection of population, movement, and landscape ecology. Properly accounting for complex animal movement processes can also potentially reduce bias in estimators of population-level parameters, thereby improving inferences that are critical for species conservation and management. This introductory article to the Special Feature reviews recent advances in SCR and animal movement modeling, establishes a common notation, highlights potential advantages of linking individual-level (Lagrangian) movements to population-level (Eulerian) processes, and outlines a general conceptual framework for the integration of movement and SCR models. We then identify important avenues for future research, including key challenges and potential pitfalls in the developments and applications that lie ahead.

Ecology↗

USGS2020

Estimating population-specific predation effects on Chinook salmon via data integration

Recent success in the conservation of many marine mammals has resulted in new management challenges due to increasing conflict with fisheries. Increasing predation by pinnipeds on threatened salmon is of particular concern. Seemingly, pinniped conservation is now in conflict with the recovery of threatened salmon, creating a dilemma for managers. We use the Lower Columbia River as a case study for examining the relationship between seasonal California sea lion Zalophus californianus abundance and survival of threatened salmon. To quantify mortality associated with increasing sea lion abundance, we examined the effect of seasonal sea lion abundance on adult Chinook salmon Oncorhynchus tshawytscha survival during migrations through the Lower Columbia River. We integrated data on survival with data on population-specific migration timing, allowing quantification of the relationship between sea lion abundance and survival in 18 populations of spring–summer Chinook salmon listed as Threatened or Endangered under the U.S. Endangered Species Act. Of the 18 populations examined, earlier migrating populations experienced lower survival in association with increased exposure to higher sea lion abundance. We estimated that in years with high sea lion abundance, the nine earliest-migrating populations experienced an additional 21.1% (95% CI = 16.3–26.1) mortality compared to years with baseline sea lion abundance, while the nine latest migrating populations experienced an additional 10.1% (7.5–13.0). Synthesis and applications . Integrating datasets on seasonal survival and migration timing made it possible for us to estimate population-specific mortality associated with increased sea lion abundance in the Lower Columbia River. This information could not be produced from any one dataset, highlighting the utility of data integration approaches. The mortality experienced by early migrating Chinook salmon suggests the potential for demographic and evolutionary consequences. Management actions such as hazing, relocating, or removing individuals that are frequent predators on salmon have been proposed. Identifying the management actions that will allow for socially and legally acceptable trade-offs between multiple conservation and other social values will be facilitated by development of explicit multi-species management frameworks. Continued monitoring will help to reduce the substantial uncertainty about the effect of pinnipeds on salmon and the predicted outcomes of alternative management actions.

Oregon, Washington↗

USGS2019

Making Recursive Bayesian inference accessible

Bayesian models provide recursive inference naturally because they can formally reconcile new data and existing scientific information. However, popular use of Bayesian methods often avoids priors that are based on exact posterior distributions resulting from former studies. Two existing Recursive Bayesian methods are: Prior- and Proposal-Recursive Bayes. Prior-Recursive Bayes uses Bayesian updating, fitting models to partitions of data sequentially, and provides a way to accommodate new data as they become available using the posterior from the previous stage as the prior in the new stage based on the latest data. ProposalRecursive Bayes is intended for use with hierarchical Bayesian models and uses a set of transient priors in first stage independent analyses of the data partitions. The second stage of Proposal-Recursive Bayes uses the posteriors from the first stage as proposals in an MCMC algorithm to fit the full model. We combine Prior- and Proposal-Recursive concepts to fit any Bayesian model, and often with computational improvements. We demonstrate our method with two case studies. Our approach has implications for big data, streaming data, and optimal adaptive design situations.

American Statistician↗

USGS2018

A guide to Bayesian model checking for ecologists

Checking that models adequately represent data is an essential component of applied statistical inference. Ecologists increasingly use hierarchical Bayesian statistical models in their research. The appeal of this modeling paradigm is undeniable, as researchers can build and fit models that embody complex ecological processes while simultaneously accounting for observation error. However, ecologists tend to be less focused on checking model assumptions and assessing potential lack of fit when applying Bayesian methods than when applying more traditional modes of inference such as maximum likelihood. There are also multiple ways of assessing the fit of Bayesian models, each of which has strengths and weaknesses. For instance, Bayesian P values are relatively easy to compute, but are well known to be conservative, producing P values biased toward 0.5. Alternatively, lesser known approaches to model checking, such as prior predictive checks, cross‐validation probability integral transforms, and pivot discrepancy measures may produce more accurate characterizations of goodness‐of‐fit but are not as well known to ecologists. In addition, a suite of visual and targeted diagnostics can be used to examine violations of different model assumptions and lack of fit at different levels of the modeling hierarchy, and to check for residual temporal or spatial autocorrelation. In this review, we synthesize existing literature to guide ecologists through the many available options for Bayesian model checking. We illustrate methods and procedures with several ecological case studies including (1) analysis of simulated spatiotemporal count data, (2) N‐mixture models for estimating abundance of sea otters from an aircraft, and (3) hidden Markov modeling to describe attendance patterns of California sea lion mothers on a rookery. We find that commonly used procedures based on posterior predictive P values detect extreme model inadequacy, but often do not detect more subtle cases of lack of fit. Tests based on cross‐validation and pivot discrepancy measures (including the “sampled predictive P value”) appear to be better suited to model checking and to have better overall statistical performance. We conclude that model checking is necessary to ensure that scientific inference is well founded. As an essential component of scientific discovery, it should accompany most Bayesian analyses presented in the literature.

Ecological Monographs↗

USGS2018

Process convolution approaches for modeling interacting trajectories

Gaussian processes are a fundamental statistical tool used in a wide range of applications. In the spatiotemporal setting, several families of covariance functions exist to accommodate a wide variety of dependence structures arising in different applications. These parametric families can be restrictive and are insufficient in some situations. In contrast, process convolutions represent a flexible, interpretable approach to defining the covariance of a Gaussian process and have modest requirements to ensure validity. We introduce a generalization of the process convolution approach that employs multiple convolutions sequentially to form a “process convolution chain”. In our proposed multistage framework, complex dependencies that arise from a combination of different interacting mechanisms are decomposed into a series of interpretable kernel smoothers. We demonstrate an application of process convolution chains to model killer whale movement, in which the paths taken by multiple individuals are not independent but reflect dynamic social interactions within the population. Our proposed model for dependent movement provides inference for the latent dynamic social structure in the study population. Additionally, by leveraging the positive dependence among individual paths, we achieve a reduction in uncertainty for the estimated locations of the killer whales, compared with a model that treats paths as independent.

Environmetrics↗

USGS2017

Basis function models for animal movement

Advances in satellite-based data collection techniques have served as a catalyst for new statistical methodology to analyze these data. In wildlife ecological studies, satellite-based data and methodology have provided a wealth of information about animal space use and the investigation of individual-based animal–environment relationships. With the technology for data collection improving dramatically over time, we are left with massive archives of historical animal telemetry data of varying quality. While many contemporary statistical approaches for inferring movement behavior are specified in discrete time, we develop a flexible continuous-time stochastic integral equation framework that is amenable to reduced-rank second-order covariance parameterizations. We demonstrate how the associated first-order basis functions can be constructed to mimic behavioral characteristics in realistic trajectory processes using telemetry data from mule deer and mountain lion individuals in western North America. Our approach is parallelizable and provides inference for heterogenous trajectories using nonstationary spatial modeling techniques that are feasible for large telemetry datasets. Supplementary materials for this article are available online.

Journal of the American Statistical Association↗

USGS2017

Imputation approaches for animal movement modeling

The analysis of telemetry data is common in animal ecological studies. While the collection of telemetry data for individual animals has improved dramatically, the methods to properly account for inherent uncertainties (e.g., measurement error, dependence, barriers to movement) have lagged behind. Still, many new statistical approaches have been developed to infer unknown quantities affecting animal movement or predict movement based on telemetry data. Hierarchical statistical models are useful to account for some of the aforementioned uncertainties, as well as provide population-level inference, but they often come with an increased computational burden. For certain types of statistical models, it is straightforward to provide inference if the latent true animal trajectory is known, but challenging otherwise. In these cases, approaches related to multiple imputation have been employed to account for the uncertainty associated with our knowledge of the latent trajectory. Despite the increasing use of imputation approaches for modeling animal movement, the general sensitivity and accuracy of these methods have not been explored in detail. We provide an introduction to animal movement modeling and describe how imputation approaches may be helpful for certain types of models. We also assess the performance of imputation approaches in two simulation studies. Our simulation studies suggests that inference for model parameters directly related to the location of an individual may be more accurate than inference for parameters associated with higher-order processes such as velocity or acceleration. Finally, we apply these methods to analyze a telemetry data set involving northern fur seals ( Callorhinus ursinus ) in the Bering Sea. Supplementary materials accompanying this paper appear online.

Journal of Agricultural, Biological, and Environme↗

USGS2017

Animal movement: Statistical models for telemetry data

The study of animal movement has always been a key element in ecological science, because it is inherently linked to critical processes that scale from individuals to populations and communities to ecosystems. Rapid improvements in biotelemetry data collection and processing technology have given rise to a variety of statistical methods for characterizing animal movement. The book serves as a comprehensive reference for the types of statistical models used to study individual-based animal movement.

Book↗

USGS2016

Dynamic occupancy models for explicit colonization processes

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.

Ecology↗

USGS2016

Dynamic social networks based on movement

Network modeling techniques provide a means for quantifying social structure in populations of individuals. Data used to define social connectivity are often expensive to collect and based on case-specific, ad hoc criteria. Moreover, in applications involving animal social networks, collection of these data is often opportunistic and can be invasive. Frequently, the social network of interest for a given population is closely related to the way individuals move. Thus, telemetry data, which are minimally invasive and relatively inexpensive to collect, present an alternative source of information. We develop a framework for using telemetry data to infer social relationships among animals. To achieve this, we propose a Bayesian hierarchical model with an underlying dynamic social network controlling movement of individuals via two mechanisms: an attractive effect and an aligning effect. We demonstrate the model and its ability to accurately identify complex social behavior in simulation, and apply our model to telemetry data arising from killer whales. Using auxiliary information about the study population, we investigate model validity and find the inferred dynamic social network is consistent with killer whale ecology and expert knowledge.

The Annals of Applied Statistics↗

USGS2015

Using spatiotemporal statistical models to estimate animal abundance and infer ecological dynamics from survey counts

Ecologists often fit models to survey data to estimate and explain variation in animal abundance. Such models typically require that animal density remains constant across the landscape where sampling is being conducted, a potentially problematic assumption for animals inhabiting dynamic landscapes or otherwise exhibiting considerable spatiotemporal variation in density. We review several concepts from the burgeoning literature on spatiotemporal statistical models, including the nature of the temporal structure (i.e., descriptive or dynamical) and strategies for dimension reduction to promote computational tractability. We also review several features as they specifically relate to abundance estimation, including boundary conditions, population closure, choice of link function, and extrapolation of predicted relationships to unsampled areas. We then compare a suite of novel and existing spatiotemporal hierarchical models for animal count data that permit animal density to vary over space and time, including formulations motivated by resource selection and allowing for closed populations. We gauge the relative performance (bias, precision, computational demands) of alternative spatiotemporal models when confronted with simulated and real data sets from dynamic animal populations. For the latter, we analyze spotted seal ( Phoca largha ) counts from an aerial survey of the Bering Sea where the quantity and quality of suitable habitat (sea ice) changed dramatically while surveys were being conducted. Simulation analyses suggested that multiple types of spatiotemporal models provide reasonable inference (low positive bias, high precision) about animal abundance, but have potential for overestimating precision. Analysis of spotted seal data indicated that several model formulations, including those based on a log-Gaussian Cox process, had a tendency to overestimate abundance. By contrast, a model that included a population closure assumption and a scale prior on total abundance produced estimates that largely conformed to our a priori expectation. Although care must be taken to tailor models to match the study population and survey data available, we argue that hierarchical spatiotemporal statistical models represent a powerful way forward for estimating abundance and explaining variation in the distribution of dynamical populations.

Ecological Monographs↗

USGS2015

Estimating demographic parameters using a combination of known-fate and open N -mixture models

Accurate estimates of demographic parameters are required to infer appropriate ecological relationships and inform management actions. Known-fate data from marked individuals are commonly used to estimate survival rates, whereas N-mixture models use count data from unmarked individuals to estimate multiple demographic parameters. However, a joint approach combining the strengths of both analytical tools has not been developed. Here we develop an integrated model combining known-fate and open N-mixture models, allowing the estimation of detection probability, recruitment, and the joint estimation of survival. We demonstrate our approach through both simulations and an applied example using four years of known-fate and pack count data for wolves (Canis lupus). Simulation results indicated that the integrated model reliably recovered parameters with no evidence of bias, and survival estimates were more precise under the joint model. Results from the applied example indicated that the marked sample of wolves was biased toward individuals with higher apparent survival rates than the unmarked pack mates, suggesting that joint estimates may be more representative of the overall population. Our integrated model is a practical approach for reducing bias while increasing precision and the amount of information gained from mark–resight data sets. We provide implementations in both the BUGS language and an R package.

Alaska↗

USGS2014

Hidden Markov model for dependent mark loss and survival estimation

Mark-recapture estimators assume no loss of marks to provide unbiased estimates of population parameters. We describe a hidden Markov model (HMM) framework that integrates a mark loss model with a Cormack–Jolly–Seber model for survival estimation. Mark loss can be estimated with single-marked animals as long as a sub-sample of animals has a permanent mark. Double-marking provides an estimate of mark loss assuming independence but dependence can be modeled with a permanently marked sub-sample. We use a log-linear approach to include covariates for mark loss and dependence which is more flexible than existing published methods for integrated models. The HMM approach is demonstrated with a dataset of black bears ( Ursus americanus ) with two ear tags and a subset of which were permanently marked with tattoos. The data were analyzed with and without the tattoo. Dropping the tattoos resulted in estimates of survival that were reduced by 0.005–0.035 due to tag loss dependence that could not be modeled. We also analyzed the data with and without the tattoo using a single tag. By not using.

Journal of Agricultural, Biological, and Environme↗

USGS2014

When to be discrete: The importance of time formulation in understanding animal movement

Animal movement is essential to our understanding of population dynamics, animal behavior, and the impacts of global change. Coupled with high-resolution biotelemetry data, exciting new inferences about animal movement have been facilitated by various specifications of contemporary models. These approaches differ, but most share common themes. One key distinction is whether the underlying movement process is conceptualized in discrete or continuous time. This is perhaps the greatest source of confusion among practitioners, both in terms of implementation and biological interpretation. In general, animal movement occurs in continuous time but we observe it at fixed discrete-time intervals. Thus, continuous time is conceptually and theoretically appealing, but in practice it is perhaps more intuitive to interpret movement in discrete intervals. With an emphasis on state-space models, we explore the differences and similarities between continuous and discrete versions of mechanistic movement models, establish some common terminology, and indicate under which circumstances one form might be preferred over another. Counter to the overly simplistic view that discrete- and continuous-time conceptualizations are merely different means to the same end, we present novel mathematical results revealing hitherto unappreciated consequences of model formulation on inferences about animal movement. Notably, the speed and direction of movement are intrinsically linked in current continuous-time random walk formulations, and this can have important implications when interpreting animal behavior. We illustrate these concepts in the context of state-space models with multiple movement behavior states using northern fur seal ( Callorhinus ursinus ) biotelemetry data.

Movement Ecology↗

USGS2013

Estimating animal resource selection from telemetry data using point process models

Analyses of animal resource selection functions (RSF) using data collected from relocations of individuals via remote telemetry devices have become commonplace. Increasing technological advances, however, have produced statistical challenges in analysing such highly autocorrelated data. Weighted distribution methods have been proposed for analysing RSFs with telemetry data. However, they can be computationally challenging due to an intractable normalizing constant and cannot be aggregated (i.e. collapsed) over time to make space-only inference. In this study, we take a conceptually different approach to modelling animal telemetry data for making RSF inference. We consider the telemetry data to be a realization of a space–time point process. Under the point process paradigm, the times of the relocations are also considered to be random rather than fixed. We show the point process models we propose are a generalization of the weighted distribution telemetry models. By generalizing the weighted model, we can access several numerical techniques for evaluating point process likelihoods that make use of common statistical software. Thus, the analysis methods can be readily implemented by animal ecologists. In addition to ease of computation, the point process models can be aggregated over time by marginalizing over the temporal component of the model. This allows a full range of models to be constructed for RSF analysis at the individual movement level up to the study area level. To demonstrate the analysis of telemetry data with the point process approach, we analysed a data set of telemetry locations from northern fur seals (Callorhinus ursinus) in the Pribilof Islands, Alaska. Both a space–time and an aggregated space-only model were fitted. At the individual level, the space–time analysis showed little selection relative to the habitat covariates. However, at the study area level, the space-only model showed strong selection relative to the covariates.

Alaska↗