Geology topics

Kathryn Irvine

Publications and source records attributed to Kathryn Irvine.

At least 19 recordsLinked to original sources

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

Estimating disease prevalence from preferentially sampled, pooled data

After the onset of the COVID-19 pandemic, scientific interest in coronaviruses endemic in animal populations has increased dramatically. However, investigating the prevalence of disease in animal populations across the landscape, which requires finding and capturing animals can be difficult. Spatial random sampling over a grid could be extremely inefficient because animals can be hard to locate, and the total number of samples may be small. Alternatively, preferential sampling, using existing knowledge to inform sample location, can guarantee larger numbers of samples, but estimates derived from this sampling scheme may exhibit bias if there is a relationship between higher probability sampling locations and the disease prevalence. Sample specimens are commonly grouped and tested in pools which can also be an added challenge when combined with preferential sampling. Here we present a Bayesian method for estimating disease prevalence with preferential sampling in pooled presence-absence data motivated by estimating factors related to coronavirus infection among Mexican free-tailed bats ( Tadarida brasiliensis ) in California. We demonstrate the efficacy of our approach in a simulation study, where a naive model, not accounting for preferential sampling, returns biased estimates of parameter values; however, our model returns unbiased results regardless of the degree of preferential sampling. Our model framework is then applied to data from California to estimate factors related to coronavirus prevalence. After accounting for preferential sampling impacts, our model suggests small prevalence differences between male and female bats.

California↗

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

Bayesian model selection to investigate meaningful spatial scales

Ecologists and other statistical practitioners with access to high-resolution spatial data 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 concentric circles 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 MCMC techniques, enables efficient estimation of posterior model probabilities to assist with the choice of meaningful spatial scales for environmental covariates.

Authorea↗

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↗

USGS2025

Statistical review of systematic reconnaissance flight (SRF) surveys (2008-2023) for monitoring American alligator (Alligator mississippiensis) nests in Everglades National Park

We reviewed current (2008-2023) protocols for monitoring American alligator nests in Everglades National Park (ENP) using a qualitative statistical review framework that emphasizes connecting measurable objectives with evaluation of survey design elements. Our review outlined the statistical assumptions that, if severely violated, would prevent defensible inferences from being drawn regarding the three alligator nesting metrics monitored in ENP’s systematic reconnaissance flight (SRF) surveys: nest effort (number of nests), nest distribution, and nest fate. After evaluating the current survey and inference methods, we offer potential adjustments or augmentations to the design and models currently used to better account for potential error sources and align the stated objectives with statistical inferences.

Florida↗

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↗

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

Using mobile acoustic monitoring and false-positive N-mixture models to estimate bat abundance and population trends

Estimating the abundance of unmarked animal populations from acoustic data is challenging due to the inability to identify individuals and the need to adjust for observation biases including detectability (false negatives), species misclassification (false positives), and sampling exposure. Acoustic surveys conducted along mobile transects were designed to avoid counting individuals more than once, where raw counts are commonly treated as an index of abundance. More recently, false-positive abundance models have been developed to estimate abundance while accounting for imperfect detection and misclassification. We adapted these methods to model summertime abundance and trends of three species of bats at multiple spatial scales using acoustic recordings collected along mobile transects by partners of the North American Bat Monitoring Program (NABat) from 2012 to 2020. This multiscale modeling spanned individual transect routes, larger NABat grid cells (10 km × 10 km), and across the entire extent of modeled species ranges. We estimated relationships between species abundances and a suite of abiotic and biotic predictors (landcover types, climatological variables, physiographic diversity, building density, and the impacts of white-nose syndrome [WNS]) and found varying levels of support between species. We present clear evidence of substantial declines in populations of tricolored bats ( Perimyotis subflavus ) and little brown bats ( Myotis lucifugus ), declines that corresponded in space and time with the progression of WNS, a devastating disease of hibernating bats. In contrast, our analysis revealed that similar population-wide declines probably have not occurred in big brown bats ( Eptesicus fuscus ), a species known to be less affected by WNS. This study provides the first abundance-based species distribution predictions and population trends for bats in their summer ranges in North America. These models will probably be applicable to assessing wildlife populations in other monitoring programs where acoustic data are used or where false-negative and false-positive detections are present. Finally, our abundance framework (as a spatial point pattern process) can serve as a foundation from which more sophisticated integrated species distribution models that incorporate additional streams of monitoring data (e.g., stationary acoustics, captures) can be developed for North American bats.

Ecological Monographs↗

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

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

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

NABat: A top-down, bottom-up solution to collaborative continental-scale monitoring

Collaborative monitoring over broad scales and levels of ecological organization can inform conservation efforts necessary to address the contemporary biodiversity crisis. An important challenge to collaborative monitoring is motivating local engagement with enough buy-in from stakeholders while providing adequate top-down direction for scientific rigor, quality control, and coordination. Collaborative monitoring must reconcile this inherent tension between top-down control and bottom-up engagement. Highly mobile and cryptic taxa, such as bats, present a particularly acute challenge. Given their scale of movement, complex life histories, and rapidly expanding threats, understanding population trends of bats requires coordinated broad-scale collaborative monitoring. The North American Bat Monitoring Program (NABat) reconciles top-down, bottom-up tension with a hierarchical master sample survey design, integrated data analysis, dynamic data curation, regional monitoring hubs, and knowledge delivery through web-based infrastructure. NABat supports collaborative monitoring across spatial and organizational scales and the full annual lifecycle of bats.

Ambio↗

USGS2020

Scientist’s guide to developing explanatory statistical models using causal analysis principles

Recent discussions of model selection and multimodel inference highlight a general challenge for researchers, which is how to clearly convey the explanatory content of a hypothesized model or set of competing models. The advice from statisticians for scientists employing multimodel inference is to develop a well‐thought‐out set of candidate models for comparison, though precise instructions for how to do that are typically not given. A coherent body of knowledge, which falls under the general term causal analysis, now exists for examining the explanatory scientific content of candidate models. Much of the literature on causal analysis has been recently developed and we suspect may not be familiar to many ecologists. This body of knowledge comprises a set of graphical tools and axiomatic principles to support scientists in their endeavors to create “well‐formed hypotheses”, as statisticians are asking them to do. Causal analysis is complementary to methods such as structural equation modeling, which provides the means for evaluation of proposed hypotheses against data. In this paper, we summarize and illustrate a set of principles that can guide scientists in their quest to develop explanatory hypotheses for evaluation. The principles presented in this paper have the capacity to close the communication gap between statisticians, who urge scientists to develop well‐thought‐out coherent models, and scientists, who would like some practical advice for exactly how to do that.

Ecology↗

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

Climatic correlates of white pine blister rust infection in whitebark pine in the Greater Yellowstone Ecosystem

Whitebark pine, a foundation species at tree line in the Western U.S. and Canada, has declined due to native mountain pine beetle epidemics, wildfire, and white pine blister rust. These declines are concerning for the multitude of ecosystem and human benefits provided by this species. Understanding climatic correlates associated with spread is needed to successfully manage impacts from forest pathogens. In the Greater Yellowstone Ecosystem since 2000 mountain pine beetles have killed 75 percent of the mature cone-bearing trees, and 40.9 percent of monitored trees have been infected with white pine blister rust. We identified models of white pine blister rust infection that indicate an August and September interaction between relative humidity and temperature were better predictors of white pine blister rust infection in whitebark pine than location and site characteristics in the Greater Yellowstone Ecosystem. The climate conditions conducive to white pine blister rust occur throughout the ecosystem, but larger trees in relatively warm and humid conditions were more likely to be infected between 2000 and 2018. We mapped the infection probability over the past two decades to identify coarse-scale patterns of climate conditions conducive to white pine blister rust infection in whitebark pine.

Wyoming↗