Geology topics

Kathryn M. Irvine

Publications and source records attributed to Kathryn M. Irvine.

At least 19 recordsLinked to original sources

USGS2026

Bayesian model selection to investigate meaningful spatial scales

Ecologists and other statistical practitioners with access to high-resolution spatial data often lack guidance on best approaches for discerning meaningful spatial scales for environmental covariates which is necessary when spatial factors influence environmental processes. Recently developed methods have attempted to automate investigating spatial scales for covariates by evaluating models for which potential explanatory variables are derived from circular extents of increasing size centered at survey locations. However, these methods make a strong assumption on the inclusion of the covariate and do not help discern whether a covariate should be included in the model. We present an approach that utilizes researcher guidance to create informative priors on the model space that, along with parallelizable Reversible Jump Markov chain Monte Carlo techniques, enables efficient estimation of posterior model probabilities to assist with the choice of meaningful spatial scales for environmental covariates.

Ecology and Evolution↗

USGS2025

Missing data in ecology: Syntheses, clarifications, and considerations

In ecology and related sciences, missing data are common and occur in a variety of different contexts. When missing data are not handled properly, subsequent statistical estimates tend to be biased, inefficient, and lack proper confidence interval coverage. Missing data are often grouped into three categories: missing completely at random (MCAR), missing at random (MAR), and missing not at random (MNAR). We review each category and compare their benefits and drawbacks. We review several approaches to handling missing data including complete case analysis, imputation, inverse probability weighting, and data augmentation. We clarify what types of variables should accompany imputation methods and how those variables are influenced by the analysis methods. Additionally, we discuss missing data that lack a formal basis for measurement and hence are fundamentally different from MCAR, MAR, and MNAR missing data. Throughout, we introduce concepts and numeric examples using both simulated data and data from the United States Environmental Protection Agency's 2016 National Wetland Condition Assessment. We conclude by providing five considerations for ecologists and other scientists handling missing data.

Ecological Monographs↗

USGS2025

Predicting bat roosts in bridges using Bayesian Additive Regression Trees

Human-built structures can provide important habitat for wildlife, but predicting which structures are most likely to be used remains challenging. To evaluate the predictive capabilities of data-driven ensemble modeling approaches, we conducted surveys for bats and signs of bat use, such as urine and guano staining, at bridges across the southwestern United States. We developed a bat roost discovery tool using Bayesian Additive Regression Trees (BART) and evaluated the predictive ability of this model against other commonly used approaches. We found that the lack of nearby water resources was associated with a lower predicted probability of bat presence or signs of bat use at bridges. While the presence of nearby water resources was associated with higher average predicted probability of bat presence or signs of bat use, high uncertainty surrounding these estimates indicates that other factors also play a role in determining which bridge roosts bats are more likely to use. As such, our model could be particularly useful for predicting which bridges can be excluded from survey efforts due to low probability of bat presence or signs of bat use. We extrapolated our model to unsurveyed bridges across the study region and provide an interactive dashboard application interface for the exploration of these results. Overall, this study demonstrates the application of BART as a predictive tool for prioritizing future bridge surveys for bats roosting in transportation structures.

Arizona, California, Nevada, New Mexico, Texas↗

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

Trends in plant cover derived from vegetation plot data using ordinal zero-augmented beta regression

Questions Plant cover values in vegetation plot data are bounded between 0 and 1, and cover is typically recorded in discrete classes with non-equal intervals. Consequently, cover data are skewed and heteroskedastic, which hampers the application of conventional regression methods. Recently developed ordinal beta regression models consider these statistical difficulties. Our primary question is whether we can detect species trends in vegetation plot time series data with this modelling approach. A second question is whether trends in cover have additional value compared to trends in occurrence, which are easier to assess for practitioners. Location The Netherlands, Western Europe. Methods We used vegetation plot data collected from 10,000 fixed plots which were surveyed once every four years during 1999–2022. We used the ordinal zero-augmented beta regression (OZAB) model, a hierarchical model consisting of a logistic regression for presence and an ordinal beta regression for cover. We adapted the OZAB model for longitudinal data and produced estimates of cover and occurrence for each four-year period. Thereafter we assessed trends in cover and in occurrence across all periods. Results We found evidence of a trend in cover in 318 out of the 721 species (44%) with sufficient data. Most species showed similar directional trends in occurrence and percent cover. No trend in occurrence was detected for 64 species that had evidence of a trend in cover. Declining species had stronger relative changes in cover than in occurrence. Conclusions Our model enables researchers to detect trends in cover using longitudinal vegetation plot data. Cover trends often corroborated trends in occurrence, but we also regularly found trends in cover even in the absence of evidence for trends in occurrence. Our approach thus contributes to a more complete picture of (changes in) vegetation composition based on large monitoring data sets.

Journal of Vegetation Science↗

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

Spatial Gaussian processes improve multi-species occupancy models when range boundaries are uncertain and nonoverlapping

Species distribution models enable practitioners to analyze large datasets of encounter records and make predictions about species occurrence at unsurveyed locations. In omnibus surveys that record data on multiple species simultaneously, species ranges are often nonoverlapping and misaligned with the administrative unit defining the spatial domain of interest (e.g., a state or province). Consequently, some species display differentially restricted extents within a study area. Assuming hard boundaries based on expert opinion or published range maps to restrict species occurrence predictions implies a false sense of certainty in model-based inferences. We propose a multi-species occupancy model with a spatial Gaussian process on site-specific effects for each species as a model-based solution. Specifying informative Bayesian hyperpriors on the spatial hyperparameters encapsulates broad-scale correlation among site occupancy probabilities for each species. We fit this model to acoustic detection/nondetection data collected with autonomous recording units during summer of 2016–2019 throughout Oregon and Washington, USA, on 15 bat species. We found vast improvements in spatial predictions of spotted bat ( Euderma maculatum ), canyon bat ( Parastrellus hesperus ), and Brazilian free-tailed bat ( Tadarida brasiliensis ) when the available environmental predictors were insufficient for characterizing their restricted ranges within the region. In contrast, widespread species ( Lasionycteris noctivagans , Myotis californicus , Myotis evotis , Myotis volans ) were appropriately modeled using only environmental predictors, such as percentage forest cover and cliff and canyon cover. Utilizing spatial Gaussian processes within a community or multi-species model incorporates uncertainty in range boundaries and allows for simultaneous predictions for the entire faunal assemblage even if species have nonoverlapping or restricted ranges within a spatial domain of interest. Such modeling improvements are essential if species distribution models are to accurately inform monitoring, species recovery plans, and other conservation efforts.

Ecology and Evolution↗

USGS2021

Adaptive monitoring in action: Reconsidering design-based estimators reveals underestimation of whitebark pine disease prevalence in the Greater Yellowstone Ecosystem

Identifying and understanding status and trends in ecological indicators motivates continual monitoring over decades. Many programs rely on probability surveys and their companion design-based estimators for status assessments (e.g. Horvitz–Thompson). Design-based estimators do not easily extend to trend estimation nor situations with observation errors. Field-based monitoring efforts inevitably have turnover of field crew members which may affect consistency and accuracy of data collection over time. Additionally, design-based estimators ignore the complexities of spatial and temporal heterogeneity in an ecological indicator and how this variability may be linked to environmental or biological dynamics. We propose monitoring programs should re-evaluate their prescribed statistical methods, consider model-based approaches and adapt their sampling designs as needed to improve inferences. The Greater Yellowstone Ecosystem, home to two of the most iconic U.S. National Parks, has experienced significant declines in whitebark pine Pinus albicaulis communities due to forest pathogens, insect outbreaks, wildland fires and drought. Whitebark pine is a keystone species found in mountainous environments throughout the Western U.S. and Canada. We assessed the design-based ratio estimator originally recommended for estimating prevalence of white pine blister rust Cronartium ribicola . We compared the design-based estimator to a model-based approach that accounts for the sampling design, imperfect detection and allows for infection probabilities to vary over space and time. Ignoring observation errors led to lower estimated prevalence of white pine blister rust in the general population. Using model-based approaches, we found that the probability of infection has increased since 2004. However, overall prevalence likely has not changed because of the mountain pine beetle Dendroctonus ponderosae -induced shift towards smaller diameter trees that have a lower probability of infection compared to their larger cohorts. Synthesis and Applications . Using a design-based approach to detect change in ecological indicators falls short because of the inability to account for observation errors or to explore environmental or biological factors explaining temporal dynamics. Inherently understanding the mechanisms leading to changes in an ecological indicator over time informs potential management actions. Our assessment underscores the need for continued evaluation and updating of a monitoring program's sampling design and analytical procedures to maintain relevancy.

Wyoming↗

USGS2021

Resilience to fire and resistance to annual grass invasion in sagebrush ecosystems of US National Parks

Western North American sagebrush shrublands and steppe face accelerating risks from fire-driven feedback loops that transition these ecosystems into self-reinforcing states dominated by invasive annual grasses. In response, sagebrush conservation decision-making is increasingly done through the lens of resilience to fire and annual grass invasion resistance. Operationalizing resilience and resistance concepts requires place-based understanding of resilience and resistance variation among landscapes over time. Place-based insights allow for landscape prioritization in targeted areas of significance such as protected-area sagebrush ecosystems that exhibit inherently low resilience and are therefore at high risk of loss. We used a multi-scale approach to evaluate sagebrush resiliency and strategic planning across 1) the US National Park system, 2) a regional suite of five parks, and 3) for two specific park case studies. First, we summarized broad patterns of relative resilience to fire and resistance to annual grass invasion across all parks with sagebrush ecosystems. We found that national parks represented ~11% of US protected-area sagebrush ecosystems and reflected a similar low-resilience bias that occurs across the biome, broadly. Climate change is likely to shift both low- and high-resilience park sagebrush ecosystems towards moderate resiliency, creating new opportunities and constraints for park conservation. Approximately seventy park units include at least some sagebrush shrublands or steppe, but we identified 40 parks with substantial amounts (>20% of park area) that can be included in an agency-wide conservation strategy. Second, we examined detailed patterns of resilience and resistance, fire history and fire risk, cheatgrass ( Bromus tectorum ) invasion, and sagebrush shrub ( Artemisia spp.) persistence in five national park units in Columbia Basin and Snake River Plain sagebrush steppe, contextualized by the broader summary. In these five parks, fire frequency and size increased in recent decades. Cheatgrass invasion and sagebrush persistence correlated strongly with resilience, burn frequency (0–3 fires since ~1940), and burn probability, but with important variation, in part mediated by local-scale topography. Third, we used these insights to assemble strategic sagebrush ecosystem fire protection mapping scenarios in two additional parks – Lava Beds National Monument and Great Basin National Park. Readily available and periodically updated geospatial data including soil surveys, fire histories, vegetation inventories, and long-term monitoring support resiliency-based adaptive management through tactical planning of pre-fire protection, post-fire restoration, and triage. Our assessment establishes the precarious importance of the US national park system to sagebrush ecosystem conservation and an operational strategy for place-based and science-supported conservation.

Arizona, California, Colorado, Idaho, Montana, New↗

USGS2021

Advancements in analytical approaches improve whitebark pine monitoring results

Long-term monitoring programs track the status and trends of species in increasingly vulnerable environments. These monitoring results provide critical information for evaluating, understanding, and managing natural resources. To accurately interpret if and how conditions may be changing for select ecological indicators, it is essential that monitoring programs adopt methods to ensure exceptional data quality. To do this and remain relevant, practitioners need to be flexible and willing to embrace a degree of adaptivity in their protocols. They must periodically re-evaluate their statistical methods and field data collection techniques to provide contemporary, significant, and applicable inferences.

Idaho, Montana, Wyoming↗

USGS2020

Post-fire vegetation response in a repeatedly burned low-elevation sagebrush steppe protected area provides insights about resilience and invasion resistance

Sagebrush steppe ecosystems are threatened by human land-use legacies, biological invasions, and altered fire and climate dynamics. Steppe protected areas are therefore of heightened conservation importance but are few and vulnerable to the same impacts broadly affecting sagebrush steppe. To address this problem, sagebrush steppe conservation science is increasingly emphasizing a focus on resilience to fire and resistance to non-native annual grass invasion as a decision framework. It is well-established that the positive feedback loop between fire and annual grass invasion is the driving process of most contemporary steppe degradation. We use a newly developed ordinal zero-augmented beta regression model fit to large-sample vegetation monitoring data from John Day Fossil Beds National Monument, USA, spanning 7 years to evaluate fire responses of two native perennial foundation bunchgrasses and two non-native invasive annual grasses in a repeatedly burned, historically grazed, and inherently low-resilient protected area. We structured our model hierarchically to support inferences about variation among ecological site types and over time after also accounting for growing-season water deficit, fine-scale topographic variation, and burn severity. We use a state-and-transition conceptual diagram and abundances of plants listed in ecological site reference conditions to formalize our hypothesis of fire-accelerated transition to ecologically novel annual grassland. Notably, big sagebrush ( Artemisia tridentata ) and other woody species were entirely removed by fire. The two perennial grasses, bluebunch wheatgrass ( Pseudoroegneria spicata ) and Thurber's needlegrass ( Achnatherum thurberianum ) exhibited fire resiliency, with no apparent trend after fire. The two annual grasses, cheatgrass ( Bromus tectorum ) and medusahead ( Taeniatherum caput-medusae ), increased in response to burn severity, most notably medusahead. Surprisingly, we found no variation in grass cover among ecological sites, suggesting fire-driven homogenization as shrubs were removed and annual grasses became dominant. We found contrasting responses among all four grass species along gradients of topography and water deficit, informative to protected-area conservation strategies. The fine-grained influence of topography was particularly important to variation in cover among species and provides a foothold for conservation in low-resilient, aridic steppe. Broadly, our study demonstrates how to operationalize resilience and resistance concepts for protected areas by integrating empirical data with conceptual and statistical models.

Oregon↗

USGS2020

The use of Bayesian priors in Ecology: The good, the bad and the not great

Bayesian data analysis (BDA) is a powerful tool for making inference from ecological data, but its full potential has yet to be realized. Despite a generally positive trajectory in research surrounding model development and assessment, far too little attention has been given to prior specification. Default priors, a sub‐class of non‐informative prior distributions that are often chosen without critical thought or evaluation, are commonly used in practice. We believe the fear of being too ‘subjective’ has prevented many researchers from using any prior information in their analyses despite the fact that defending prior choice (informative or not) promotes good statistical practice. In this commentary, we provide an overview of how BDA is currently being used in a random sample of articles, discuss implications for inference if current bad practices continue, and highlight sub‐fields where knowledge about the system has improved inference and promoted good statistical practices through the careful and justified use of informative priors. We hope to inspire a renewed discussion about the use of Bayesian priors in Ecology with particular attention paid to specification and justification. We also emphasize that all priors are the result of a subjective choice, and should be discussed in that way.

Methods in Ecology and Evolution↗

USGS2019

Evidence of region‐wide bat population decline from long‐term monitoring and Bayesian occupancy models with empirically informed priors

Strategic conservation efforts for cryptic species, especially bats, are hindered by limited understanding of distribution and population trends. Integrating long‐term encounter surveys with multi‐season occupancy models provides a solution whereby inferences about changing occupancy probabilities and latent changes in abundance can be supported. When harnessed to a Bayesian inferential paradigm, this modeling framework offers flexibility for conservation programs that need to update prior model‐based understanding about at‐risk species with new data. This scenario is exemplified by a bat monitoring program in the Pacific Northwestern United States in which results from 8 years of surveys from 2003 to 2010 require updating with new data from 2016 to 2018. The new data were collected after the arrival of bat white‐nose syndrome and expansion of wind power generation, stressors expected to cause population declines in at least two vulnerable species, little brown bat ( Myotis lucifugus ) and the hoary bat ( Lasiurus cinereus ). We used multi‐season occupancy models with empirically informed prior distributions drawn from previous occupancy results (2003–2010) to assess evidence of contemporary decline in these two species. Empirically informed priors provided the bridge across the two monitoring periods and increased precision of parameter posterior distributions, but did not alter inferences relative to use of vague priors. We found evidence of region‐wide summertime decline for the hoary bat ( = 0.86 ± 0.10) since 2010, but no evidence of decline for the little brown bat ( = 1.1 ± 0.10). White‐nose syndrome was documented in the region in 2016 and may not yet have caused regional impact to the little brown bat. However, our discovery of hoary bat decline is consistent with the hypothesis that the longer duration and greater geographic extent of the wind energy stressor (collision and barotrauma) have impacted the species. These hypotheses can be evaluated and updated over time within our framework of pre–post impact monitoring and modeling. Our approach provides the foundation for a strategic evidence‐based conservation system and contributes to a growing preponderance of evidence from multiple lines of inquiry that bat species are declining.

Oregon, Washington↗

USGS2019

Cohesive framework for modeling plant cover class data

The study of plant distribution and abundance is a fundamental pursuit in ecology and conservation biology. Measuring plant abundance by visually assessing percent cover and recording a cover class is a common field method that yields ordinal data. Statistical models for ordinal data exist but entail cumbersome interpretations and sometimes restrictive assumptions. We propose a Bayesian hierarchical framework for analysing cover class data that allows for linking ordinal observations to a latent beta distribution and accounts for zero inflation. Harnessing a latent beta distribution supports interpreting changes in abundance in terms of mean percent cover rather than odds ratios of cumulative cover classes as for cumulative link models. The zero augmentation allows for simultaneous inferences on both occurrence (distribution) and abundance. We show how our model can account for true and false zeros, misclassification of cover classes, multiple species and hierarchical sampling designs, using empirical examples and simulations. Simulated observation errors, when ignored, led to models overestimating abundance and underestimating occurrence. Based on simulations, we found no substantial difference between mean percent cover estimates when analyzing ordinal cover classes versus continuous percent cover as the response. Our empirical datasets displayed high probability of detection (>0.85 on average for all species), likely due to the sampling design used and training of observers. Probability of occurrence was slightly underestimated for bare ground, Artemisia tridentata , Elycap medusae , and Poa secunda using a model that ignored imperfect detection. Estimated mean percent cover was not substantially impacted by ignoring measurement error for five plant species and bare ground. Our modelling framework for cover class data allows for an explicit separation of distribution from abundance and, importantly, allows for interpreting species–environment relationships in terms of variation in mean percent cover as compared to cumulative odds ratios. The beta distribution inherently accommodates heteroscedasticity and skewness, statistical properties that are a consequence of spatially aggregated patterns common to plant survey data. Recording cover classes provides a reliable, efficient way to measure plants and our simulations suggest little loss of information compared to assuming continuous percent cover. We provide JAGS and Stan model code for implementation.

Methods in Ecology and Evolution↗

USGS2019

Using the beta distribution to analyze plant cover data

Most plant species are spatially aggregated. Local demographic and ecological processes (e.g. vegetative growth and limited seed dispersal) result in a clustered spatial pattern within an environmentally homogenous area. Spatial aggregation should be considered when modelling plant abundance data. Commonly, plant abundance is quantified by measuring cover within multiple areal plots or along multiple lines randomly placed within a study area. A common practice for analyzing plant cover is to use statistical methods that rely on the normal distribution for quantifying uncertainty. This is problematic because plant cover data tend to be left‐skewed (J‐shaped), right skewed (L‐shaped) or U‐shaped and, therefore, commonly violate classic statistical assumptions, such as normality. We outline statistical analyses that explicitly account for spatial aggregation by assuming that plant cover is beta‐distributed. The beta distribution is a flexible choice because within the open unit interval it can take on a wide range of shapes (L, J, U, or a bell‐shaped). We discuss and introduce extensions to the beta distribution that address common analysis issues encountered in plant cover datasets, such as i) the treatment of zero and one cover values, ii) hierarchical data structures, and iii) observations errors. For heuristic purposes, we focus on single species analyses, but we demonstrate how the outlined methods can be generalized to more species. The assumption that plant cover is beta‐distributed allows us to estimate the degree of spatial aggregation, and the ecological significance of this new knowledge is discussed. We provide a summary of available software for analyses (emphasizing standard R packages) and include worked examples and a simulation study comparing analysis options as supplemental information. Synthesis . Previously, the state of the statistical software made it practically difficult for empirical plant ecologists to analyze their cover data correctly, but new theory and R‐packages have been developed, and this difficulty no longer exists. We recommend that empirical plant ecologists embrace the new statistical possibilities for exploring the exciting ecological features in spatial variation of plant cover.

Journal of Ecology↗

USGS2019

Identifying occupancy model inadequacies: Can residuals separately assess detection and presence?

Occupancy models are widely applied to estimate species distributions, but few methods exist for model checking. Thorough model assessments can uncover inadequacies and allow for deeper ecological insight by exploring structure in the observed data not accounted for by a model. We introduce occupancy model residual definitions that utilize the posterior distribution of the partially latent occupancy states. Residual‐based assessments are valuable because they can target specific assumptions and identify ways to improve a model, such as adding spatial correlation or meaningful covariates. Our approach defines separate residuals for occupancy and detection, and we use simulation to examine whether missing structure for modeling detection probabilities can be distinguished from that for occupancy probabilities. In many scenarios, our residual diagnostics were able to successfully separate inadequacies at the different model levels, but we describe other situations when this may not be the case. Applying Moran's I residual diagnostics to assess models for silver‐haired (Lasionycteris noctivagans) and little brown (Myotis lucifugus) bats only provided evidence of residual spatial correlation among detections. Targeting specific model assumptions using carefully chosen residual diagnostics is valuable for any analysis, and we remove previous barriers for occupancy analyses — lack of examples and practical advice.

Ecology↗