Geology topics

Christian Stratton

Publications and source records attributed to Christian Stratton.

12 recordsLinked to original sources

USGS2025

Integrated distribution modeling resolves asynchrony between bat population impacts and occupancy trends through latent abundance

Monitoring populations is challenging for cryptic species with seasonal life cycles, where data from multiple field techniques are commonly collected and analyzed as multiple lines of evidence. Data integration can provide comprehensive inferences while improving accuracy, precision, and scope but faces challenges in modeling misaligned resolutions and observational uncertainties. We developed a multi-scale, integrated species distribution model (MS-iSDM) for North American bats to combine data across monitoring types and seasons using joint likelihood methods, observational models with false-negatives and false-positives, and seasonal migratory connectivity. We applied this model to 11 years of data for an imperiled bat species (tricolored bat, Perimyotis subflavus ). Relative abundance and occupancy were linked with multi-scale predictors, revealing clear patterns of population declines, but with important differences in spatial trends (abundance: corresponded with white-nose syndrome impacts, occupancy: at the range periphery) and overall severity (abundance: -74.8%, 95% CRI: -79.7% to -69.3%; occupancy: -35.5%, 95% CRI: -41.1% to -30.2%). The asynchrony between occupancy trends and population impacts was explained as an emergent pattern of spatiotemporal variation in abundance in the integrated distribution model. Compared to multiple lines of evidence, the integrated model provided consensus-estimates, increased precision and spatiotemporal scope, and strengthened evidence of population declines.

Communications Biology↗

USGS2025

Leveraging an observed-data likelihood improves the use of machine learning labels in a Bayesian hierarchical model for bioacoustic data

Classification of massive datasets by machine learning (ML) algorithms is promising for many scientific domains, especially wildlife monitoring programs that rely on passive acoustic surveys for detecting species. However, treating ML-predicted class labels (e.g., species identity) as truth biases inferences of focal parameters within common modeling frameworks. One solution is to model the misclassification process explicitly using human-validated true-class labels for a subset of observations. Validation by experts can present a substantial bottleneck in otherwise efficient workflows that use ML predictions. Bioacoustics practitioners seek guidance on both the quantity and process for selecting ML-labeled data to validate by an expert. We derive an alternative model formulation that jointly models human-validated and ML-predicted class labels with an observed-data likelihood (ODL) and use empirically informed simulations motivated by a real-data application to explore different probability designs for selecting class labels for validation. Simulation results suggest that with smaller validation sets the ODL formulation increases computational speed and reduces estimation error compared to a default MCMC data augmentation routine. Our methodology is transferable to applications that treat predictions from classification algorithms as the response variable of interest.

Annals of Applied Statistics↗

USGS2024

Clustering and unconstrained ordination with Dirichlet process mixture models

Assessment of similarity in species composition or abundance across sampled locations is a common goal in multi-species monitoring programs. Existing ordination techniques provide a framework for clustering sample locations based on species composition by projecting high-dimensional community data into a low-dimensional, latent ecological gradient representing species composition. However, these techniques require specification of the number of distinct ecological communities present in the latent space, which can be difficult to determine in advance. We develop an ordination model capable of simultaneous clustering and ordination that allows for estimation of the number of clusters present in the latent ecological gradient. This model draws latent coordinates for each sample location from a Dirichlet process mixture model, affording researchers with probabilistic statements about the number of clusters present in the latent ecological gradient. The model is compared to existing methods for simultaneous clustering and ordination via simulation and applied to two empirical datasets; JAGS code to fit the proposed model is provided in an appendix. The first dataset concerns presence-absence records of fish in the Doubs river in eastern France and the second dataset describes presence-absence records of plant species in Craters of the Moon National Monument and Preserve (CRMO) in Idaho, USA. Results from both analyses align with existing ecological gradients at each location. Development of the Dirichlet process ordination model provides wildlife managers with data-driven inferences about the number of distinct communities present across monitored locations, allowing for more cost-effective monitoring and reliable decision-making for conservation management.

Methods in Ecology and Evolution↗

USGS2024

Joint spatial modeling bridges the gap between disparate disease surveillance and population monitoring efforts informing conservation of at-risk bat species

White-Nose Syndrome (WNS) is a wildlife disease that has decimated hibernating bats since its introduction in North America in 2006. As the disease spreads westward, assessing the potentially differential impact of the disease on western bat species is an urgent conservation need. The statistical challenge is that the disease surveillance and species response monitoring data are not co-located, available at different spatial resolutions, non-Gaussian, and subject to observation error requiring a novel extension to spatially misaligned regression models for analysis. Previous work motivated by epidemiology applications has proposed two-step approaches that overcome the spatial misalignment while intentionally preventing the human health outcome from informing estimation of exposure. In our application, the impacted animals contribute to spreading the fungus that causes WNS, motivating development of a joint framework that exploits the known biological relationship. We introduce a Bayesian, joint spatial modeling framework that provides inferences about the impact of WNS on measures of relative bat activity and accounts for the uncertainty in estimation of WNS presence at non-surveyed locations. Our simulations demonstrate that the joint model produced more precise estimates of disease occurrence and unbiased estimates of the association between disease presence and the count response relative to competing two-step approaches. Our statistical framework provides a solution that leverages disparate monitoring activities and informs species conservation across large landscapes. Stan code and documentation are provided to facilitate access and adaptation for other wildlife disease applications.

Journal of Agricultural, Biological and Environmen↗

USGS2022

NABat ML: Utilizing deep learning to enable crowdsourced development of automated, scalable solutions for documenting North American bat populations

Bats play crucial ecological roles and provide valuable ecosystem services, yet many populations face serious threats from various ecological disturbances. The North American Bat Monitoring Program (NABat) aims to use its technology infrastructure to assess status and trends of bat populations, while developing innovative and community-driven conservation solutions. Here, we present NABat ML , an automated machine-learning algorithm that improves the scalability and scientific transparency of NABat acoustic monitoring. This model combines signal processing techniques and convolutional neural networks (CNNs) to detect and classify recorded bat echolocation calls. We developed our CNN model with internet-based computing resources (‘cloud environment’), and trained it on >600,000 spectrogram images. We also incorporated species range maps to improve the robustness and accuracy of the model for future ‘unseen’ data. We evaluated model performance using a comprehensive, independent, holdout dataset. NABat ML successfully distinguished 31 classes (30 species and a noise class) with overall weighted-average accuracy and precision rates of 92%, and ≥90% classification accuracy for 19 of the bat species. Using a single cloud-environment computing instance, the entire model training process took <16 h. Synthesis and applications . Our convolutional neural network (CNN)-based model, NABat ML , classifies 30 North American bat species using their recorded echolocation calls with an overall accuracy of 92%. In addition to providing highly accurate species-level classification, NABat ML and its outputs are compatible with Bayesian and other statistical techniques for measuring uncertainty in classification. Our model is open-source and reproducible, enabling future implementations as software on end-user devices and cloud-based web applications. These qualities make NABat ML highly suitable for applications ranging from grassroots community science initiatives to big-data methods developed and implemented by researchers and professional practitioners. We believe the transparency and accessibility of NABat ML will encourage broad-scale participation in bat monitoring, and enable development of innovative solutions needed to conserve North American bat species.

Journal of Applied Ecology↗

USGS2022

Statistical assessment on determining local presence of rare bat species

Surveying cryptic, sparsely distributed taxa using autonomous recording units, although cost-effective, provides imperfect knowledge about species presence. Summertime bat acoustic surveys in North America exemplify the challenges with characterizing sources of uncertainty: observation error, inability to census populations, and natural stochastic variation. Statistical uncertainty, if not considered thoroughly, hampers determining rare species presence accurately and/or estimating rangewide status and trends with suitable precision. Bat acoustic data are processed using an automated workflow in which proprietary or open-source algorithms assign a species label to each recorded high-frequency echolocation sequence. A false-negative occurs, if a species is actually present but not recorded and/or all recordings from the species are of such poor quality that a correct species identity cannot be assigned to any observation. False positives for a focal species are a direct result of the presence and incorrect identification of a recording from another species. We compare four analytical approaches in terms of parameter estimation and their resulting (in)correct decisions regarding species presence or absence using realistic data-generating scenarios for bat acoustic data within a simulation study. The current standard for deciding species presence or absence uses a multinomial likelihood-ratio test p value (maximum likelihood estimate [MLE]-metric) that accounts for known species misidentifications, but not imperfect detection and only returns a binary outcome (evidence of presence or not). We found that the MLE-metric had estimated median correct decisions less than 60% for presence and greater than 85% for absence. Alternatively, a multispecies count detection model was equivalent to or better than the MLE-metric for correct claims of rare species presence or absence using the posterior probability a species was present at a site and, importantly, provided unbiased estimates of relative activity and probability of occurrence, creating opportunities for reducing posterior uncertainty through the inclusion of meaningful covariates. Single-species occupancy models with and without false-positive detections removed were insufficient for determining local presence because of substantially biased occurrence and detection probabilities. We propose solutions to potential barriers for integrating local, short-term and rangewide, long-term acoustic surveys within a cohesive statistical framework that facilitates determining local species presence with uncertainty concurrent with estimating species–environment relationships.

Ecosphere↗

USGS2022

Coupling validation effort with in situ bioacoustic data improves estimating relative activity and occupancy for multiple species with cross-species misclassifications

The increasing complexity and pace of ecological change requires natural resource managers to consider entire species assemblages. Acoustic recording units (ARUs) require minimal cost and effort to deploy and inform relative activity, or encounter rates, for multiple species simultaneously. ARU-based surveys require post-processing of the recordings via software algorithms that assign a species label to each recording. The automated classification process can result in cross-species misidentifications that should be accounted for when employing statistical modelling for conservation decision-making. Using simulation and ARU-based detection counts from 17 bat species in British Columbia, Canada, we investigate three strategies for adjusting statistical inference for species misclassification: (a) ‘coupling’ ambiguous and unambiguous detections by validating a subset of survey events post-hoc, (b) using a calibration dataset on the software algorithm's (in)accuracy for species identification or (c) specifying informative Bayesian priors on classification probabilities. We explore the impact of different Bayesian prior specifications for the classification probabilities on posterior estimation. We then consider how the quantity of data validated post-hoc impacts model convergence and resulting inferences for bat species relative activity as related to nightly conditions and yearly site occupancy after accounting for site-level environmental variables. Coupled methods resulted in less bias and uncertainty when estimating relative activity and species classification probabilities relative to calibration approaches. We found that species that were difficult-to-detect and those that were often inaccurately identified by the software required more validation effort than more easily detected and/or identified species. Our results suggest that, when possible, acoustic surveys should rely on coupled validated detection information to account for false-positive detections, rather than uncoupled calibration datasets. However, if the assemblage of interest contains a large number of rarely detected or less prevalent species, an intractable amount of effort may be required, suggesting there are benefits to curating a calibration dataset that is representative of the observation process. Our findings provide insights into the practical challenges associated with statistical analyses of ARU data and possible analytical solutions to support reliable and cost-effective decision-making for wildlife conservation/management in the face of known sources of observation errors.

British Columbia↗

USGS2022

Status and trends of North American bats: Summer occupancy analysis 2010-2019

• We developed an analytical pipeline supported by web-based infrastructure for integrating continental scale bat monitoring data (stationary acoustic, mobile acoustic, and capture records) to estimate summer (May 1–Aug 31) occupancy probabilities and changes in occupancy over time for 12 North American bat species. This serves as one of multiple lines of evidence that inform the status and trends of bat populations. • We analyzed data from a total of 12 bat species (Table 1), 11 of which have tested positive for Pseudogymnoascus destructans (Pd), a fungal pathogen that causes white-nose syndrome (WNS)—a disease that has led to significant rates of mortality for subterranean hibernating bat species in North America. A twelfth species was also selected because of high rates of mortality at wind energy facilities. Additional species were considered but not selected due to data limitations. • We estimated occupancy probabilities for 2010 through 2019 for three species ( Myotis lucifugus , MYLU; Myotis septentrionalis , MYSE; and P erimyotis subflavus , PESU). For an additional nine species, we estimated occupancy probabilities for 2016 through 2019 ( Myotis evotis , MYEV; Myotis grisescens , MYGR; Myotis leibii , MYLE; Myotis thysanodes , MYTH; Myotis volans , MYVO; Myotis yumanensis , MYYU; Eptesicus fuscus , EPFU; Lasionycteris noctivagans , LANO; and Lasiurus cinereus , LACI). • For each species, we provide range-wide occupancy probability predictions (e.g., predicted summer occupancy distribution maps) each year at a spatial resolution of 100 km2 and provide regional estimates of mean occupancy probability aggregated at larger spatial scales (state/province/territory, range-wide). • For each species, we also provide trends over time (average annual change rate and total change rate) in mean occupancy probabilities at multiple spatial scales (state/province/territory, range-wide) and when possible, over multiple timescales (short, medium, long). • Results suggest that over the short-term (2016-2019), two ( Myotis lucifugus and Perimyotis subflavus ) of 12 species have experienced declines in range-wide average occupancy probability with at least 95% certainty. Seven species showed either minor increases or decreases in range-wide average occupancy probability but with less than 95% certainty in both trend indicators. Results over the longer term (eight years and 10 years of sampling) suggest that three hibernating species known to be highly affected by white-nose syndrome ( Myotis lucifugus , Myotis septentrionalis , and Perimyotis subflavus ) have experienced marked declines in range-wide average occupancy probabilities, with severity varying by species and region. Finally, the results for three species ( Eptesicus fuscus , Lasiurus cinereus , Lasionycteris noctivagans ) were inconclusive due to 1) borderline convergence issues in the model fitting procedure which suggests potentially unreliable estimates, 2) failure to reliably distinguish between false positives and true positive detections for ambiguous detections, and 3) largely uninformative covariates for occupancy and detection. • For Myotis lucifugus , Myotis septentrionalis , and Perimyotis subflavus we found meaningful associations in space and time between declining winter populations (likely a result of WNS) and summer occupancy distributions. • The representativeness of sampling data for each species’ status and trend estimates (e.g., state/province/territory) were also evaluated based on the percent of grid cells sampled each year with a goal of understanding the reliability of regional estimates and improving future monitoring efforts. • This work represents the most comprehensive effort to date to model North American bat distributions across their continental ranges. Despite current limitations highlighted in the discussion, the analytical methods and resulting status and trends estimates provide the best available science on summer bat populations across North America and will continue to improve over time as monitoring data sets and analytical methods improve. • Moving forward, our occupancy analyses will continue to improve with submission of more 1) data from currently underrepresented areas (i.e., improved geographic representation), 2) manually-vetted acoustic recordings, 3) capture records, and 4) roost location and count data (summer and winter).

Report↗

USGS2021

Integrating environmental DNA results with diverse data sets to improve biosurveillance of river health

Autonomous, robotic environmental (e)DNA samplers now make it possible for biological observations to match the scale and quality of abiotic measurements collected by automated sensor networks. Merging these automated data streams may allow for improved insight into biotic responses to environmental change and stressors. Here, we merged eDNA data collected by robotic samplers installed at three U.S. Geological Survey (USGS) streamgages with gridded daily weather data, and daily water quality and quantity data into a cloud-hosted database. The eDNA targets were a rare fish parasite and a more common salmonid fish. We then used computationally expedient Bayesian hierarchical occupancy models to evaluate associations between abiotic conditions and eDNA detections and to simulate how uncertainty in result interpretation changes with the frequency of autonomous robotic eDNA sample collection. We developed scripts to automate data merging, cleaning and analysis steps into a chained-step, workflow. We found that inclusion of abiotic covariates only provided improved insight for the more common salmonid fish since its DNA was more frequently detected. Rare fish parasite DNA was infrequently detected, which caused occupancy parameter estimates and covariate associations to have high uncertainty. Our simulations found that collecting samples at least once per day resulted in more detections and less parameter uncertainty than less frequent sampling. Our occupancy and simulation results together demonstrate the advantages of robotic eDNA samplers and how these samples can be combined with easy to acquire, publicly available data to foster real-time biosurveillance and forecasting.

Idaho, Wyoming, Montana↗

USGS2020

msocc: Fit and analyse computationally efficient multi‐scale occupancy models in R

Environmental DNA (eDNA) sampling is a promising tool for the detection of rare and cryptic taxa, such as aquatic pathogens, parasites and invasive species. Environmental DNA sampling workflows commonly rely on multi‐stage hierarchical sampling designs that induce complicated dependencies within the data. This complex dependence structure can be intuitively modelled with Bayesian multi‐scale occupancy models. However, current software for such models are computationally demanding, impeding their use. We present an r package, msocc , that implements a data augmentation strategy to fit fully Bayesian, computationally efficient multi‐scale occupancy models. The msocc package allows users to fit multi‐scale occupancy models, to estimate and visualize posterior summaries of site, sample and replicate‐level occupancy, and to compare different models using Bayesian information criterion. Additionally, we provide a supplemental web application that allows users to investigate study design for multi‐scale occupancy models and acts as a graphical user interface to the msocc package. The utility of the msocc package is illustrated on a published dataset and the functions in msocc are compared to the primary Bayesian toolkit for multi‐scale occupancy modelling, eDNAoccupancy , using various computational benchmarks. These benchmarks indicate that msocc is capable of fitting models 50 times faster than eDNAoccupancy . We hope that access to software that efficiently fits, analyses and conducts study design investigations for multi‐scale occupancy models facilitates their implementation by the research and wildlife management communities.

Methods in Ecology and Evolution↗

USGS2019

Adding invasive species bio-surveillance to the U.S. Geological Survey streamgage network

The costs of invasive species in the United States alone are estimated to exceed US$100 billion per year so a critical tactic in minimizing the costs of invasive species is the development of effective, early-detection systems. To this end, we evaluated the efficacy of adding environmental (e)DNA surveillance to the U.S. Geological Survey (USGS) streamgage network, which consists of > 8,200 streamgages nationwide systemically visited by USGS hydrologic technicians. Incorporating strategic eDNA sample collection during routine streamgage visits could provide early detection surveillance of aquatic invasive species with minimal additional cost. For this evaluation, USGS hydrologic technicians collected monthly eDNA water samples, May – September 2018, from streamgages downstream of reservoirs in the Columbia River Basin thought to be vulnerable to invasive dreissenid mussel ( Dreissenidae spp.) establishment. We tested water samples for dreissenid mussel DNA and also for kokanee ( Oncorhynchus nerka ) and yellow perch ( Perca flavescens ) DNA; the two fishes were used to assess if streamgages are adequately located to provide early-detection eDNA surveillance of taxa known to be present in upstream reservoirs. No Columbia River Basin streamgage samples met our criteria for being scored as positive for dreissenid DNA. We did detect kokanee and yellow perch DNA at all streamgages downstream of reservoirs where these species are known to occur. Field collection, laboratory analyses, and personnel time required for collection of four eDNA samples at a streamgage site cost US$500 -US$600 (net). Given these results, incorporating eDNA biosurveillance into routine streamgage visits might decrease costs associated with an invasion since early detection maximizes the potential for eradication, containment, and mitigation.

Ecosphere↗

USGS2019

Assessing spatial and temporal patterns in sagebrush steppe vegetation communities 2012-2018: Grand Teton National Park

Visual cover class data were collected on over 80 species across 30 permanent sampling frames in sagebrush steppe vegetation communities in Grand Teton National Park from 2012 to 2018. In this report, temporal and spatial patterns in species composition were assessed and used to inform potential sampling strategies for future monitoring. Specifically, the viability of a reduction in sampling effort was evaluated based on the similarity in species composition within each frame over time and among frames within each year. Using distance-based ordination techniques, we found little to no evidence of differences in species composition within each frame over time. Furthermore, there was little evidence of heterogeneity in species composition among frames within each year, though there was some evidence of differences in composition between the two principle sagebrush community types (sagebrush dry shrubland and sagebrush-bitterbrush) aggregated across frames. Based on these results, we propose that a reduction in sampling effort is viable and suggest a new monitoring schedule.

Wyoming↗