Geology Reports⌕ Search

Geology topics

Thomas J. Rodhouse

Publications and source records attributed to Thomas J. Rodhouse.

15 recordsLinked to original sources
USGS2020

Case Study 4: NABat acoustic monitoring allows inferences about bat populations at multiple scales

North American bats face unprecedented risks from continuing and emerging threats including white-nose syndrome, wind energy development, and habitat loss. Many species of bats are thought to be recently experiencing unparalleled population declines unlike any previously observed (O’Shea et al. 2016). The North American Bat Monitoring Program (NABat) was conceived to better understand the true ecological consequences of these large-scale population reductions (Loeb et al. 2015). NABat aims is to improve the state of conservation science for the 47 species of bats shared by Canada, United States, and Mexico. To meet this objective, NABat offers standardize protocols and a unifying sample design facilitating a multi-agency, multinational, collaborative monitoring effort. A key element of NABat is cross-boundary partner coordination and sharing of limited resources for the collection of bat echolocation data. Here we provide three compelling examples of how NABat provides a convenient framework for using acoustic data to assess the potential impacts of current and future threats to North American bats across multiple spatial scales.

Oregon↗
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

North American Bat Monitoring Program regional protocol for surveying with stationary deployments of echolocation recording devices: Narrative version 1.0, Pacific Northwestern US

The outbreak of white-nose syndrome (WNS) and the growing awareness of the risks to bats from wind power generating facilities have driven radical changes to North American bat conservation. Over the last decade, formerly common species such as the little brown myotis (Myotis lucifugus) and hoary bat (Lasiurus cinereus) have experienced unprecedented mortality rates and are now facing non-trivial extinction risk. In response to this change, federal land management agencies such as the US National Park Service, US Fish and Wildlife Service, US Forest Service, US Bureau of Land Management and state wildlife management agencies such as the Oregon Department of Fish and Wildlife and Idaho Fish and Game have invested in collaborative, interagency bat monitoring to close the gap in information about bat welfare and to inform bat conservation strategies. Bats are notoriously difficult to track and study and there remains a paucity of fundamental information about the seasonal patterns of bat activity and habitat use and population distributions and abundances. Moreover, because bats are so highly mobile and difficult to survey (e.g., nocturnal flight), this information needs to be contextualized at broad regional (e.g., 10,000 km2) and range-wide extents. Delimiting bat populations at local scales (e.g., 100 km2) is very difficult and it is not clear, for example, how a declining trend in local (e.g., a small park unit) patterns of bat activity or relative abundance should be interpreted without broader context. In recognition of these challenges, a plan for coordinated continental-scale monitoring of bats, the North American Bat Monitoring Program (NABat) was developed (Loeb et al. 2015). The centerpiece of the plan is the use of a spatially-balanced randomized master sample of grid-cell sample units from a grid-based sampling frame to provide the architecture for collaboration and the statistical foundation for making inferences about bat populations across broad regions and entire bat geographic ranges. The plan outlines general goals, survey design, and field methods for both summertime acoustic surveys of bats as well as winter and summer counts of bats in hibernacula and maternity colonies but it does not provide field-level protocol and standard operating procedures for consistent and efficient implementation. This regional protocol provides these details for one component of NABat, the deployment of stationary acoustic detectors to record bats during summer, as is called for by the NABat plan. This protocol was written specifically to provide guidance and consistency across the Pacific Northwestern US (N. California [California Department of Fish and Wildlife Northern Region], Idaho, Washington, and Oregon; US Fish and Wildlife Service Region 1 and portion of Region 8 [in Northern California and Klamath Basin]; US Forest Service Region 6 and portions of Regions 1 and 5 in Idaho; and the Upper Columbia Basin, North Coast Cascades, and Klamath Networks of the National Park Service). This region has internal cohesion, sharing a distinct bat faunal assemblage of 15 species (with several additional species occurring on the southern periphery of the region), and a long history of collaborative bat monitoring beginning with the interagency Bat Grid Program which operated from 2003-2010 across Oregon and Washington (US Forest Service Region 6). This protocol will be coordinated and implemented by the Northwestern Bat Hub, on behalf of the collective interagency partnership. The Northwestern Bat Hub is housed on the Oregon State University-Cascades campus and leverages pooled partner funds and resources to maintain a small staff that coordinates and conducts monitoring, provides training and oversight, ensures high-quality data quality and control, and analyzes data and reports on results.

California, Idaho, Oregon, Washington↗
USGS2018

Occupancy modeling species–environment relationships with non‐ignorable survey designs

Statistical models supporting inferences about species occurrence patterns in relation to environmental gradients are fundamental to ecology and conservation biology. A common implicit assumption is that the sampling design is ignorable and does not need to be formally accounted for in analyses. The analyst assumes data are representative of the desired population and statistical modeling proceeds. However, if data sets from probability and non‐probability surveys are combined or unequal selection probabilities are used, the design may be non‐ignorable. We outline the use of pseudo‐maximum likelihood estimation for site‐occupancy models to account for such non‐ignorable survey designs. This estimation method accounts for the survey design by properly weighting the pseudo‐likelihood equation. In our empirical example, legacy and newer randomly selected locations were surveyed for bats to bridge a historic statewide effort with an ongoing nationwide program. We provide a worked example using bat acoustic detection/non‐detection data and show how analysts can diagnose whether their design is ignorable. Using simulations we assessed whether our approach is viable for modeling data sets composed of sites contributed outside of a probability design. Pseudo‐maximum likelihood estimates differed from the usual maximum likelihood occupancy estimates for some bat species. Using simulations we show the maximum likelihood estimator of species–environment relationships with non‐ignorable sampling designs was biased, whereas the pseudo‐likelihood estimator was design unbiased. However, in our simulation study the designs composed of a large proportion of legacy or non‐probability sites resulted in estimation issues for standard errors. These issues were likely a result of highly variable weights confounded by small sample sizes (5% or 10% sampling intensity and four revisits). Aggregating data sets from multiple sources logically supports larger sample sizes and potentially increases spatial extents for statistical inferences. Our results suggest that ignoring the mechanism for how locations were selected for data collection (e.g., the sampling design) could result in erroneous model‐based conclusions. Therefore, in order to ensure robust and defensible recommendations for evidence‐based conservation decision‐making, the survey design information in addition to the data themselves must be available for analysts. Details for constructing the weights used in estimation and code for implementation are provided.

Ecological Applications↗
USGS2018

Improving geographically extensive acoustic survey designs for modeling species occurrence with imperfect detection and misidentification

Acoustic recording units (ARUs) enable geographically extensive surveys of sensitive and elusive species. However, a hidden cost of using ARU data for modeling species occupancy is that prohibitive amounts of human verification may be required to correct species identifications made from automated software. Bat acoustic studies exemplify this challenge because large volumes of echolocation calls could be recorded and automatically classified to species. The standard occupancy model requires aggregating verified recordings to construct confirmed detection/non‐detection datasets. The multistep data processing workflow is not necessarily transparent nor consistent among studies. We share a workflow diagramming strategy that could provide coherency among practitioners. A false‐positive occupancy model is explored that accounts for misclassification errors and enables potential reduction in the number of confirmed detections. Simulations informed by real data were used to evaluate how much confirmation effort could be reduced without sacrificing site occupancy and detection error estimator bias and precision. We found even under a 50% reduction in total confirmation effort, estimator properties were reasonable for our assumed survey design, species‐specific parameter values, and desired precision. For transferability, a fully documented r package, OCacoustic, for implementing a false‐positive occupancy model is provided. Practitioners can apply OCacoustic to optimize their own study design (required sample sizes, number of visits, and confirmation scenarios) for properly implementing a false‐positive occupancy model with bat or other wildlife acoustic data. Additionally, our work highlights the importance of clearly defining research objectives and data processing strategies at the outset to align the study design with desired statistical inferences.

Ecology and Evolution↗
USGS2016

A goodness-of-fit test for occupancy models with correlated within-season revisits

Occupancy modeling is important for exploring species distribution patterns and for conservation monitoring. Within this framework, explicit attention is given to species detection probabilities estimated from replicate surveys to sample units. A central assumption is that replicate surveys are independent Bernoulli trials, but this assumption becomes untenable when ecologists serially deploy remote cameras and acoustic recording devices over days and weeks to survey rare and elusive animals. Proposed solutions involve modifying the detection-level component of the model (e.g., first-order Markov covariate). Evaluating whether a model sufficiently accounts for correlation is imperative, but clear guidance for practitioners is lacking. Currently, an omnibus goodnessof- fit test using a chi-square discrepancy measure on unique detection histories is available for occupancy models (MacKenzie and Bailey, Journal of Agricultural, Biological, and Environmental Statistics, 9, 2004, 300; hereafter, MacKenzie– Bailey test). We propose a join count summary measure adapted from spatial statistics to directly assess correlation after fitting a model. We motivate our work with a dataset of multinight bat call recordings from a pilot study for the North American Bat Monitoring Program. We found in simulations that our join count test was more reliable than the MacKenzie–Bailey test for detecting inadequacy of a model that assumed independence, particularly when serial correlation was low to moderate. A model that included a Markov-structured detection-level covariate produced unbiased occupancy estimates except in the presence of strong serial correlation and a revisit design consisting only of temporal replicates. When applied to two common bat species, our approach illustrates that sophisticated models do not guarantee adequate fit to real data, underscoring the importance of model assessment. Our join count test provides a widely applicable goodness-of-fit test and specifically evaluates occupancy model lack of fit related to correlation among detections within a sample unit. Our diagnostic tool is available for practitioners that serially deploy survey equipment as a way to achieve cost savings.

Ecology and Evolution↗
USGS2016

Distribution of a climate-sensitive species at an interior range margin

Advances in understanding the factors that limit a species’ range, particularly in the context of climate change, have come disproportionately through investigations at range edges or margins. The margins of a species’ range might often correspond with anomalous microclimates that confer habitat suitability where the species would otherwise fail to persist. We addressed this hypothesis using data from an interior, climatic range margin of the American pika ( Ochotona princeps ), an indicator of relatively cool, mesic climates in rocky habitats of western North America. Pikas in Lava Beds National Monument, northeastern California, USA, occur at elevations much lower than predicted by latitude and longitude. We hypothesized that pika occurrence within Lava Beds would be associated primarily with features such as “ice caves” in which sub-surface ice persists outside the winter months. We used data loggers to monitor sub-surface temperatures at cave entrances and at non-cave sites, confirming that temperatures were cooler and more stable at cave entrances. We surveyed habitat characteristics and evidence of pika occupancy across a random sample of cave and non-cave sites over a 2-yr period. Pika detection probability was high (~0.97), and the combined occupancy of cave and non-cave sites varied across the 2 yr from 27% to 69%. Contrary to our hypothesis, occupancy was not higher at cave sites. Vegetation metrics were the best predictors of site use by pikas, followed by an edge effect and elevation. The importance of vegetation as a predictor of pika distribution at this interior range margin is congruent with recent studies from other portions of the species’ range. However, we caution that vegetation composition depends on microclimate, which might be the proximal driver of pika distribution. The microclimates available in non-cave crevices accessible to small animals have not been characterized adequately for lava landscapes. We advocate innovation in the acquisition and use of microclimatic data for understanding the distributions of many taxa. Appropriately scaled microclimatic data are increasingly available but rarely used in studies of range dynamics.

Ecosphere↗
USGS2015

Establishing conservation baselines with dynamic distribution models for bat populations facing imminent decline

Aim Bat mortality rates from white-nose syndrome and wind power development are unprecedented. Cryptic and wide-ranging behaviours of bats make them difficult to survey, and population estimation is often intractable. We advance a model-based framework for making spatially explicit predictions about summertime distributions of bats from capture and acoustic surveys. Motivated by species-energy and life-history theory, our models describe hypotheses about spatio-temporal variation in bat distributions along environmental gradients and life-history attributes, providing a statistical basis for conservation decision-making. Location Oregon and Washington, USA. Methods We developed Bayesian hierarchical models for 14 bat species from an 8-year monitoring dataset across a ~430,000 km 2 study area. Models accounted for imperfect detection and were temporally dynamic. We mapped predicted occurrence probabilities and prediction uncertainties as baselines for assessing future declines. Results Forest cover, snag abundance and cliffs were important predictors for most species. Species occurrence patterns varied along elevation and precipitation gradients, suggesting a potential hump-shaped diversity–productivity relationship. Annual turnover in occurrence was generally low, and occurrence probabilities were stable among most species. We found modest evidence that turnover covaried with the relative riskiness of bat roosting and migration. The fringed myotis ( Myotis thysanodes ), canyon bat ( Parastrellus hesperus ) and pallid bat ( Antrozous pallidus ) were rare; fringed myotis occurrence probabilities declined over the study period. We simulated anticipated declines to demonstrate that mapped occurrence probabilities, updated over time, provide an intuitive way to assess bat conservation status for a broad audience. Main conclusions Landscape keystone structures associated with roosting habitat emerged as regionally important predictors of bat distributions. The challenges of bat monitoring have constrained previous species distribution modelling efforts to temporally static presence-only approaches. Our approach extends to broader spatial and temporal scales than has been possible in the past for bats, making a substantial increase in capacity for bat conservation.

Oregon, Washington↗
USGS2015

A plan for the North American Bat Monitoring Program (NABat)

The purpose of the North American Bat Monitoring Program (NABat) is to create a continent-wide program to monitor bats at local to rangewide scales that will provide reliable data to promote effective conservation decisionmaking and the long-term viability of bat populations across the continent. This is an international, multiagency program. Four approaches will be used to gather monitoring data to assess changes in bat distributions and abundances: winter hibernaculum counts, maternity colony counts, mobile acoustic surveys along road transects, and acoustic surveys at stationary points. These monitoring approaches are described along with methods for identifying species recorded by acoustic detectors. Other chapters describe the sampling design, the database management system (Bat Population Database), and statistical approaches that can be used to analyze data collected through this program.

Report↗
USGS2014

Predicting foundation bunchgrass species abundances: Model-assisted decision-making in protected-area sagebrush steppe

Foundation species are structurally dominant members of ecological communities that can stabilize ecological processes and influence resilience to disturbance and resistance to invasion. Being common, they are often overlooked for conservation but are increasingly threatened from land use change, biological invasions, and over-exploitation. The pattern of foundation species abundances over space and time may be used to guide decision-making, particularly in protected areas for which they are iconic. We used ordinal logistic regression to identify the important environmental influences on the abundance patterns of bluebunch wheatgrass ( Pseudoroegneria spicata ), Thurber's needlegrass ( Achnatherum thurberianum ), and Sandberg bluegrass ( Poa secunda ) in protected-area sagebrush steppe. We then predicted bunchgrass abundances along gradients of topography, disturbance, and invasive annual grass abundance. We used model predictions to prioritize the landscape for implementation of a management and restoration decision-support tool. Models were fit to categorical estimates of grass cover obtained from an extensive ground-based monitoring dataset. We found that remnant stands of abundant wheatgrass and bluegrass were associated with steep north-facing slopes in higher and more remote portions of the landscape outside of recently burned areas where invasive annual grasses were less abundant. These areas represented only 25% of the landscape and were prioritized for protection efforts. Needlegrass was associated with south-facing slopes, but in low abundance and in association with invasive cheatgrass ( Bromus tectorum ). Abundances of all three species were strongly negatively correlated with occurrence of another invasive annual grass, medusahead ( Taeniatherum caput-medusae ). The rarity of priority bunchgrass stands underscored the extent of degradation and the need for prioritization. We found no evidence that insularity reduced invasibility; annual grass invasion represents a serious threat to protected-area bunchgrass communities. Our study area was entirely within the Wyoming big sagebrush ecological zone, understood to have inherently low resilience to disturbance and resistance to weed invasion. However, our study revealed important variation in abundance of the foundation species associated with resilience and resistance along the topographic-soil moisture gradient within this zone, providing an important foothold for conservation decision-making in these steppe ecosystems. We found the foundation species focus a parsimonious strategy linking monitoring to decision-making via biogeographic modeling.

Ecosphere↗
USGS2014

Accounting for false-positive acoustic detections of bats using occupancy models

1. Acoustic surveys have become a common survey method for bats and other vocal taxa. Previous work shows that bat echolocation may be misidentified, but common analytic methods, such as occupancy models, assume that misidentifications do not occur. Unless rare, such misidentifications could lead to incorrect inferences with significant management implications. 2. We fit a false-positive occupancy model to data from paired bat detector and mist-net surveys to estimate probability of presence when survey data may include false positives. We compared estimated occupancy and detection rates to those obtained from a standard occupancy model. We also derived a formula to estimate the probability that bats were present at a site given its detection history. As an example, we analysed survey data for little brown bats Myotis lucifugus from 135 sites in Washington and Oregon, USA. 3. We estimated that at an unoccupied site, acoustic surveys had a 14% chance per night of producing spurious M. lucifugus detections. Estimated detection rates were higher and occupancy rates were lower under the false-positive model, relative to a standard occupancy model. Un-modelled false positives also affected inferences about occupancy at individual sites. For example, probability of occupancy at individual sites with acoustic detections but no captures ranged from 2% to 100% under the false-positive occupancy model, but was always 100% under a standard occupancy model. 4. Synthesis and applications. Our results suggest that false positives sufficient to affect inferences may be common in acoustic surveys for bats. We demonstrate an approach that can estimate occupancy, regardless of the false-positive rate, when acoustic surveys are paired with capture surveys. Applications of this approach include monitoring the spread of White-Nose Syndrome, estimating the impact of climate change and informing conservation listing decisions. We calculate a site-specific probability of occupancy, conditional on survey results, which could inform local permitting decisions, such as for wind energy projects. More generally, the magnitude of false positives suggests that false-positive occupancy models can improve accuracy in research and monitoring of bats and provide wildlife managers with more reliable information.

Oregon, Washington↗
USGS2014

Statistical assessment on a combined analysis of GRYN-ROMN-UCBN upland vegetation vital signs

As of 2013, Rocky Mountain and Upper Columbia Basin Inventory and Monitoring Networks have multiple years of vegetation data and Greater Yellowstone Network has three years of vegetation data and monitoring is ongoing in all three networks. Our primary objective is to assess whether a combined analysis of these data aimed at exploring correlations with climate and weather data is feasible. We summarize the core survey design elements across protocols and point out the major statistical challenges for a combined analysis at present. The dissimilarity in response designs between ROMN and UCBN-GRYN network protocols presents a statistical challenge that has not been resolved yet. However, the UCBN and GRYN data are compatible as they implement a similar response design; therefore, a combined analysis is feasible and will be pursued in future. When data collected by different networks are combined, the survey design describing the merged dataset is (likely) a complex survey design. A complex survey design is the result of combining datasets from different sampling designs. A complex survey design is characterized by unequal probability sampling, varying stratification, and clustering (see Lohr 2010 Chapter 7 for general overview). Statistical analysis of complex survey data requires modifications to standard methods, one of which is to include survey design weights within a statistical model. We focus on this issue for a combined analysis of upland vegetation from these networks, leaving other topics for future research. We conduct a simulation study on the possible effects of equal versus unequal probability selection of points on parameter estimates of temporal trend using available packages within the R statistical computing package. We find that, as written, using lmer or lm for trend detection in a continuous response and clm and clmm for visually estimated cover classes with “raw” GRTS design weights specified for the weight argument leads to substantially different results and/or computational instability. However, when only fixed effects are of interest, the survey package (svyglm and svyolr) may be suitable for a model-assisted analysis for trend. We provide possible directions for future research into combined analysis for ordinal and continuous vital sign indictors.

Natural Resource Report↗
USGS2013

Woodland dynamics at the northern range periphery: A challenge for protected area management in a changing world

Managers of protected natural areas increasingly are confronted with novel ecological conditions and conflicting objectives to preserve the past while fostering resilience for an uncertain future. This dilemma may be pronounced at range peripheries where rates of change are accelerated and ongoing invasions often are perceived as threats to local ecosystems. We provide an example from City of Rocks National Reserve (CIRO) in southern Idaho, positioned at the northern range periphery of pinyon-juniper (P-J) woodland. Reserve managers are concerned about P-J woodland encroachment into adjacent sagebrush steppe, but the rates and biophysical variability of encroachment are not well documented and management options are not well understood. We quantified the rate and extent of woodland change between 1950 and 2009 based on a random sample of aerial photo interpretation plots distributed across biophysical gradients. Our study revealed that woodland cover remained at approximately 20% of the study area over the 59-year period. In the absence of disturbance, P-J woodlands exhibited the highest rate of increase among vegetation types at 0.37% yr −1 . Overall, late-successional P-J stands increased in area by over 100% through the process of densification (infilling). However, wildfires during the period resulted in a net decrease of woody evergreen vegetation, particularly among early and mid-successional P-J stands. Elevated wildfire risk associated with expanding novel annual grasslands and drought is likely to continue to be a fundamental driver of change in CIRO woodlands. Because P-J woodlands contribute to regional biodiversity and may contract at trailing edges with global warming, CIRO may become important to P-J woodland conservation in the future. Our study provides a widely applicable toolset for assessing woodland ecotone dynamics that can help managers reconcile the competing demands to maintain historical fidelity and contribute meaningfully to the U.S. protected area network in a future with novel, no-analog ecosystems.

Idaho↗
USGS2012

Assessing the status and trend of bat populations across broad geographic regions with dynamic distribution models

Bats face unprecedented threats from habitat loss, climate change, disease, and wind power development, and populations of many species are in decline. A better ability to quantify bat population status and trend is urgently needed in order to develop effective conservation strategies. We used a Bayesian autoregressive approach to develop dynamic distribution models for Myotis lucifugus , the little brown bat, across a large portion of northwestern USA, using a four-year detection history matrix obtained from a regional monitoring program. This widespread and abundant species has experienced precipitous local population declines in northeastern USA resulting from the novel disease white-nose syndrome, and is facing likely range-wide declines. Our models were temporally dynamic and accounted for imperfect detection. Drawing on species–energy theory, we included measures of net primary productivity (NPP) and forest cover in models, predicting that M. lucifugus occurrence probabilities would covary positively along those gradients. Despite its common status, M. lucifugus was only detected during ∼50% of the surveys in occupied sample units. The overall naïve estimate for the proportion of the study region occupied by the species was 0.69, but after accounting for imperfect detection, this increased to ∼0.90. Our models provide evidence of an association between NPP and forest cover and M. lucifugus distribution, with implications for the projected effects of accelerated climate change in the region, which include net aridification as snowpack and stream flows decline. Annual turnover, the probability that an occupied sample unit was a newly occupied one, was estimated to be low (∼0.04–0.14), resulting in flat trend estimated with relatively high precision (SD = 0.04). We mapped the variation in predicted occurrence probabilities and corresponding prediction uncertainty along the productivity gradient. Our results provide a much needed baseline against which future anticipated declines in M. lucifugus occurrence can be measured. The dynamic distribution modeling approach has broad applicability to regional bat monitoring efforts now underway in several countries and we suggest ways to improve and expand our grid-based monitoring program to gain robust insights into bat population status and trend across large portions of North America.

Ecological Applications↗