Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Algorithms”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32Linked to original sources
USGS2021

Optimal allocation of law enforcement patrol effort to mitigate poaching activities

Poaching is a global problem causing the decline of species worldwide. Optimizing the efficiency of ranger patrols to deter poaching activity at the lowest possible cost is crucial for protecting species with limited resources. We applied decision analysis and spatial optimization algorithms to allocate efforts of ranger patrols throughout a national park. Our objective was to mitigate poaching activity at or below management risk targets for the lowest monetary cost. We examined this trade‐off by constructing a Pareto efficiency frontier using integer linear programming. We used data from a ranger‐based monitoring program in Nyungwe National Park, Rwanda. Our measure of poaching risk is based on dynamic occupancy models that account for imperfect detection of poaching activities. We found that in order to achieve a 5% reduction in poaching risk, 622 ranger patrol events (each corresponding to patrolling 1‐km 2 sites) were needed within a year at a cost of US$49,760. In order to attain a 60% reduction in poaching risk, 15,560 patrol events were needed at a cost of US \$ 1,244,800. We evaluated the trade‐off between patrol cost and poaching risk based on our model by constructing a Pareto efficiency frontier and park managers found the solution for a 50% risk reduction to be a practical trade‐off based on funding constraints (comparable to recent years) and the diminishing returns between risk mitigation and cost. This expected reduction in risk required 8,558 patrol events per year at a cost of US \$ 684,640. Our results suggest that optimal solutions could increase efficiency compared to the actual effort allocations from 2006 to 2016 in Nyungwe National Park (e.g., risk reductions of ~30% under recent budgets compared to ~50% reduction in risk under the optimal strategy). The modeling framework in this study took into account imperfect detection of poaching risk as well as the directional and conditional nature of ranger patrol events given the spatial adjacency relationships of neighboring sites and access points. Our analyses can help to improve the efficiency of ranger patrols, and the modeling framework can be broadly applied to other spatial conservation planning problems with conditional, multilevel, site selection.

Nyungwe National Park↗
USGS2022

The wildland-urban interface in the United States based on 125 million building locations

The wildland-urban interface (WUI) is the focus of many important land management issues, such as wildfire, habitat fragmentation, invasive species, and human-wildlife conflicts. Wildfire is an especially critical issue, because housing growth in the WUI increases wildfire ignitions and the number of homes at risk. Identifying the WUI is important for assessing and mitigating impacts of development on wildlands and for protecting homes from natural hazards, but data on housing development for large areas are often coarse. We created new WUI maps for the conterminous U.S. based on 125 million individual building locations, offering higher spatial precision compared to existing maps based on U.S. census housing data. Building point locations were based on a building footprint dataset from Microsoft®. We classified WUI across the conterminous U.S. at 30-m resolution using a circular neighborhood mapping algorithm with a variable radius to determine thresholds of housing density and vegetation cover. We used our maps to (1) determine the total area of the WUI and number of buildings included, (2) assess the sensitivity of WUI area included and spatial pattern of WUI maps to choice of neighborhood size, (3) assess regional differences between building-based WUI maps and census-based WUI maps, and (4) determine how building location accuracy affected WUI map accuracy. Our building-based WUI maps identified 5.6% – 18.8% of the conterminous U.S. as being in the WUI, with larger neighborhoods increasing WUI area but excluding isolated building clusters. Building-based maps identified more WUI area relative to census-based maps for all but the smallest neighborhoods, particularly in the north-central states, and large differences were attributable to high numbers of non-housing structures in rural areas. Overall WUI classification accuracy was 98.0%. For wildfire risk mapping and for general purposes, WUI maps based on the 500-m neighborhood represent the original Federal Register definition of the WUI; these maps include clusters of buildings in and adjacent to wildlands and exclude remote, isolated buildings. Our approach for mapping the WUI offers flexibility and high spatial detail, and can be widely applied to take advantage of the growing availability of high-resolution building footprint datasets and classification methods.

Conterminous United States↗
USGS2024

Declining ecological resilience and invasion resistance under climate change in the sagebrush region, United States

In water-limited dryland ecosystems of the Western United States, climate change is intensifying the impacts of heat, drought, and wildfire. Disturbances often lead to increased abundance of invasive species, in part, because dryland restoration and rehabilitation are inhibited by limited moisture and infrequent plant recruitment events. Information on ecological resilience to disturbance (recovery potential) and resistance to invasive species can aid in addressing these challenges by informing long-term restoration and conservation planning. Here, we quantified the impacts of projected future climate on ecological resilience and invasion resistance (R&R) in the sagebrush region using novel algorithms based on ecologically relevant and climate-sensitive predictors of climate and ecological drought. We used a process-based ecohydrological model to project these predictor variables and resulting R&R indicators for two future climate scenarios and 20 climate models. Results suggested widespread future R&R decreases (24%–34% of the 1.16 million km 2 study area) that are generally consistent among climate models. Variables related to rising temperatures were most strongly linked to decreases in R&R indicators. New continuous R&R indices quantified responses to climate change; particularly useful for areas without projected change in the R&R category but where R&R still may decrease, for example, some of the areas with a historically low R&R category. Additionally, we found that areas currently characterized as having high sagebrush ecological integrity had the largest areal percentage with expected declines in R&R in the future, suggesting continuing declines in sagebrush ecosystems. One limitation of these R&R projections was relatively novel future climatic conditions in particularly hot and dry areas that were underrepresented in the training data. Including more data from these areas in future updates could further improve the reliability of the projections. Overall, these projected future declines in R&R highlight a growing challenge for natural resource managers in the region, and the resulting spatially explicit datasets provide information that can improve long-term risk assessments, prioritizations, and climate adaptation efforts.

western United States↗
USGS2024

Large-scale assessment of genetic structure to assess risk of populations of a large herbivore to disease

Chronic wasting disease (CWD) can spread among cervids by direct and indirect transmission, the former being more likely in emerging areas. Identifying subpopulations allows the delineation of focal areas to target for intervention. We aimed to assess the population structure of white-tailed deer ( Odocoileus virginianus ) in the northeastern United States at a regional scale to inform managers regarding gene flow throughout the region. We genotyped 10 microsatellites in 5701 wild deer samples from Maryland, New York, Ohio, Pennsylvania, and Virginia. We evaluated the distribution of genetic variability through spatial principal component analysis and inferred genetic structure using non-spatial and spatial Bayesian clustering algorithms (BCAs). We simulated populations representing each inferred wild cluster, wild deer in each state and each physiographic province, total wild population, and a captive population. We conducted genetic assignment tests using these potential sources, calculating the probability of samples being correctly assigned to their origin. Non-spatial BCA identified two clusters across the region, while spatial BCA suggested a maximum of nine clusters. Assignment tests correctly placed deer into captive or wild origin in most cases (94%), as previously reported, but performance varied when assigning wild deer to more specific origins. Assignments to clusters inferred via non-spatial BCA performed well, but efficiency was greatly reduced when assigning samples to clusters inferred via spatial BCA. Differences between spatial BCA clusters are not strong enough to make assignment tests a reliable method for inferring the geographic origin of deer using 10 microsatellites. However, the genetic distinction between clusters may indicate natural and anthropogenic barriers of interest for management.

Maryland, New York, Ohio, Pennsylvania, West Virgi↗
USGS2024

Genetic analysis of harvest samples reveals population structure in a highly mobile generalist carnivore

Delineating wildlife population boundaries is important for effective population monitoring and management. The bobcat ( Lynx rufus ) is a highly mobile generalist carnivore that is ecologically and economically important. We sampled 1225 bobcats harvested in South Dakota, USA (2014–2019), of which 878 were retained to assess genetic diversity and infer population genetic structure using 17 microsatellite loci. We assigned individuals to genetic clusters ( K ) using spatial and nonspatial Bayesian clustering algorithms and quantified differentiation ( F ST and GST″ ) among clusters. We found support for population genetic structure at K = 2 and K = 4, with pairwise F ST and GST″ values indicating weak to moderate differentiation, respectively, among clusters. For K = 2, eastern and western clusters aligned closely with historical bobcat management units and were consistent with a longitudinal suture zone for bobcats previously identified in the Great Plains. We did not observe patterns of population genetic structure aligning with major rivers or highways. Genetic divergence observed at K = 4 aligned roughly with ecoregion breaks and may be associated with environmental gradients, but additional sampling with more precise locational data may be necessary to validate these patterns. Our findings reveal that cryptic population structure may occur in highly mobile and broadly distributed generalist carnivores, highlighting the importance of considering population structure when establishing population monitoring programs or harvest regulations. Our study further demonstrates that for elusive furbearers, harvest can provide an efficient, broad-scale sampling approach for genetic population assessments.

South Dakota↗
USGS2024

Treed Gaussian processes for animal movement modeling

Wildlife telemetry data may be used to answer a diverse range of questions relevant to wildlife ecology and management. One challenge to modeling telemetry data is that animal movement often varies greatly in pattern over time, and current continuous-time modeling approaches to handle such nonstationarity require bespoke and often complex models that may pose barriers to practitioner implementation. We demonstrate a novel application of treed Gaussian process (TGP) modeling, a Bayesian machine learning approach that automatically captures the nonstationarity and abrupt transitions present in animal movement. The machine learning formulation of TGPs enables modeling to be nearly automated, while their Bayesian formulation allows for the derivation of movement descriptors with associated uncertainty measures. We demonstrate the use of an existing R package to implement TGPs using the familiar Markov chain Monte Carlo algorithm. We then use estimated movement trajectories to derive movement descriptors that can be compared across individuals and populations. We applied the TGP model to a case study of lesser prairie-chickens ( Tympanuchus pallidicinctus ) to demonstrate the benefits of TGP modeling and compared distance traveled and residence times across lesser prairie-chicken individuals and populations. For broad usability, we outline all steps necessary for practitioners to specify relevant movement descriptors (e.g., turn angles, speed, contact points) and apply TGP modeling and trajectory comparison to their own telemetry datasets. Combining the predictive power of machine learning and the statistical inference of Bayesian methods to model movement trajectories allows for the estimation of statistically comparable movement descriptors from telemetry studies. Our use of an accessible R package allows practitioners to model trajectories and estimate movement descriptors, facilitating the use of telemetry data to answer applied management questions.

Ecology and Evolution↗
USGS2019

Spatial sampling bias and model complexity in stream-based species distribution models: A case study of Paddlefish (Polyodon spathula) in the Arkansas River basin, USA

Leveraging existing presence records and geospatial datasets, species distribution modeling has been widely applied to informing species conservation and restoration efforts. Maxent is one of the most popular modeling algorithms, yet recent research has demonstrated Maxent models are vulnerable to prediction errors related to spatial sampling bias and model complexity. Despite elevated rates of biodiversity imperilment in stream ecosystems, the application of Maxent models to stream networks has lagged, as has the availability of tools to address potential sources of error and calculate model evaluation metrics when modeling in nonraster environments (such as stream networks). Herein, we use Maxent and customized R code to estimate the potential distribution of paddlefish ( Polyodon spathula ) at a stream-segment level within the Arkansas River basin, USA, while accounting for potential spatial sampling bias and model complexity. Filtering the presence data appeared to adequately remove an eastward, large-river sampling bias that was evident within the unfiltered presence dataset. In particular, our novel riverscape filter provided a repeatable means of obtaining a relatively even coverage of presence data among watersheds and streams of varying sizes. The greatest differences in estimated distributions were observed among models constructed with default versus AIC C -selected parameterization. Although all models had similarly high performance and evaluation metrics, the AIC C -selected models were more inclusive of westward-situated and smaller, headwater streams. Overall, our results solidified the importance of accounting for model complexity and spatial sampling bias in SDMs constructed within stream networks and provided a roadmap for future paddlefish restoration efforts in the study area.

Arkansas, Colorado, Kansas, Missouri, Nebraska, Ne↗
USGS2025

Simple bagged movement models for telemetry data

Determining which statistical methods are appropriate for data is both user and data dependent and prone to change as new methodology becomes available. This process encompasses model ideation, model selection, and determining appropriate use of statistical methods. Literature on models for animal movement emerging in the past two decades has yielded a rich collection of statistical methods garnering much deserved positive attention. Among such efforts, there is limited investigation of the broader place for simple machine learning methodology in animal movement modeling. We propose a bagged (i.e., bootstrap aggregated) animal movement model using simple, off-the-shelf machine learning algorithms. The model is intuitive, retains statistical inference about characteristics of animal movement (i.e., estimated from model-based summary statistics), and only requires knowledge of elementary statistical and machine learning analysis to understand. We show by simulation that our model can provide unbiased estimates of pertinent characteristics of animal movement (e.g., daily displacement) in the presence of large and realistic location error. We believe that increasing accessible literature on simple machine learning animal movement models provides valuable pedagogical and practical support for researchers using statistical models to study animal movement.

Ecology and Evolution↗
USGS2022

Classifying behavior from short-interval biologging data: An example with GPS tracking of birds

Recent advances in digital data collection have spurred accumulation of immense quantities of data that have potential to lead to remarkable ecological insight, but that also present analytic challenges. In the case of biologging data from birds, common analytical approaches to classifying movement behaviors are largely inappropriate for these massive data sets. We apply a framework for using K -means clustering to classify bird behavior using points from short time interval GPS tracks. K -means clustering is a well-known and computationally efficient statistical tool that has been used in animal movement studies primarily for clustering segments of consecutive points. To illustrate the utility of our approach, we apply K -means clustering to six focal variables derived from GPS data collected at 1–11 s intervals from free-flying bald eagles ( Haliaeetus leucocephalus ) throughout the state of Iowa, USA. We illustrate how these data can be used to identify behaviors and life-stage- and age-related variation in behavior. After filtering for data quality, the K -means algorithm identified four clusters in >2 million GPS telemetry data points. These four clusters corresponded to three movement states: ascending, flapping, and gliding flight; and one non-moving state: perching. Mapping these states illustrated how they corresponded tightly to expectations derived from natural history observations; for example, long periods of ascending flight were often followed by long gliding descents, birds alternated between flapping and gliding flight. The K -means clustering approach we applied is both an efficient and effective mechanism to classify and interpret short-interval biologging data to understand movement behaviors. Furthermore, because it can apply to an abundance of very short, irregular, and high-dimensional movement data, it provides insight into small-scale variation in behavior that would not be possible with many other analytical approaches.

Ecology and Evolution↗
USGS2017

Detect and exploit hidden structure in fatty acid signature data

Estimates of predator diet composition are essential to our understanding of their ecology. Although several methods of estimating diet are practiced, methods based on biomarkers have become increasingly common. Quantitative fatty acid signature analysis (QFASA) is a popular method that continues to be refined and extended. Quantitative fatty acid signature analysis is based on differences in the signatures of prey types, often species, which are recognized and designated by investigators. Similarly, predator signatures may be structured by known factors such as sex or age class, and the season or region of sample collection. The recognized structure in signature data inherently influences QFASA results in important and typically beneficial ways. However, predator and prey signatures may contain additional, hidden structure that investigators either choose not to incorporate into an analysis or of which they are unaware, being caused by unknown ecological mechanisms. Hidden structure also influences QFASA results, most often negatively. We developed a new method to explore signature data for hidden structure, called divisive magnetic clustering (DIMAC). Our DIMAC approach is based on the same distance measure used in diet estimation, closely linking methods of data exploration and parameter estimation, and it does not require data transformation or distributional assumptions, as do many multivariate ordination methods in common use. We investigated the potential benefits of the DIMAC method to detect and subsequently exploit hidden structure in signature data using two prey signature libraries with quite different characteristics. We found that the existence of hidden structure in prey signatures can increase the confusion between prey types and thereby reduce the accuracy and precision of QFASA diet estimates. Conversely, the detection and exploitation of hidden structure represent a potential opportunity to improve predator diet estimates and may lead to new insights into the ecology of either predator or prey. The DIMAC algorithm is implemented in the R diet estimation package qfasar.

Ecosphere↗
USGS2017

Potential paths for male-mediated gene flow to and from an isolated grizzly bear population

For several decades, grizzly bear populations in the Greater Yellowstone Ecosystem (GYE) and the Northern Continental Divide Ecosystem (NCDE) have increased in numbers and range extent. The GYE population remains isolated and although effective population size has increased since the early 1980s, genetic connectivity between these populations remains a long-term management goal. With only ~110 km distance separating current estimates of occupied range for these populations, the potential for gene flow is likely greater now than it has been for many decades. We sought to delineate potential paths that would provide the opportunity for male-mediated gene flow between the two populations. We first developed step-selection functions to generate conductance layers using ecological, physical, and anthropogenic landscape features associated with non-stationary GPS locations of 124 male grizzly bears (199 bear-years). We then used a randomized shortest path (RSP) algorithm to estimate the average number of net passages for all grid cells in the study region, when moving from an origin to a destination node. Given habitat characteristics that were the basis for the conductance layer, movements follow certain grid cell sequences more than others and the resulting RSP values thus provide a measure of movement potential. Repeating this process for 100 pairs of random origin and destination nodes, we identified paths for three levels of random deviation (θ) from the least-cost path. We observed broad-scale concordance between model predictions for paths originating in the NCDE and those originating in the GYE for all three levels of movement exploration. Model predictions indicated that male grizzly bear movement between the ecosystems could involve a variety of routes, and verified observations of grizzly bears outside occupied range supported this finding. Where landscape features concentrated paths into corridors (e.g., because of anthropogenic influence), they typically followed neighboring mountain ranges, of which several could serve as pivotal stepping stones. The RSP layers provide detailed, spatially explicit information for land managers and organizations working with land owners to identify and prioritize conservation measures that maintain or enhance the integrity of potential areas conducive to male grizzly bear dispersal.

Montana↗
USGS2017

Macroscale patterns of synchrony identify complex relationships among spatial and temporal ecosystem drivers

Ecology has a rich history of studying ecosystem dynamics across time and space that has been motivated by both practical management needs and the need to develop basic ideas about pattern and process in nature. In situations in which both spatial and temporal observations are available, similarities in temporal behavior among sites (i.e., synchrony) provide a means of understanding underlying processes that create patterns over space and time. We used pattern analysis algorithms and data spanning 22–25 yr from 601 lakes to ask three questions: What are the temporal patterns of lake water clarity at sub‐continental scales? What are the spatial patterns (i.e., geography) of synchrony for lake water clarity? And, what are the drivers of spatial and temporal patterns in lake water clarity? We found that the synchrony of water clarity among lakes is not spatially structured at sub‐continental scales. Our results also provide strong evidence that the drivers related to spatial patterns in water clarity are not related to the temporal patterns of water clarity. This analysis of long‐term patterns of water clarity and possible drivers contributes to understanding of broad‐scale spatial patterns in the geography of synchrony and complex relationships between spatial and temporal patterns across ecosystems.

Ecosphere↗
USGS2018

Forecasting an invasive species’ distribution with global distribution data, local data, and physiological information

Understanding invasive species distributions and potential invasions often requires broad‐scale information on the environmental tolerances of the species. Further, resource managers are often faced with knowing these broad‐scale relationships as well as nuanced environmental factors related to their landscape that influence where an invasive species occurs and potentially could occur. Using invasive buffelgrass ( Cenchrus ciliaris ), we developed global models and local models for Saguaro National Park, Arizona, USA, based on location records and literature on physiological tolerances to environmental factors to investigate whether environmental relationships of a species at a global scale are also important at local scales. In addition to correlative models with five commonly used algorithms, we also developed a model using a priori user‐defined relationships between occurrence and environmental characteristics based on a literature review. All correlative models at both scales performed well based on statistical evaluations. The user‐defined curves closely matched those produced by the correlative models, indicating that the correlative models may be capturing mechanisms driving the distribution of buffelgrass. Given climate projections for the region, both global and local models indicate that conditions at Saguaro National Park may become more suitable for buffelgrass. Combining global and local data with correlative models and physiological information provided a holistic approach to forecasting invasive species distributions.

Arizona↗
USGS2019

Designing multi-scale hierarchical monitoring frameworks for wildlife to support management: A sage-grouse case study

Population monitoring is integral to the conservation and management of wildlife; yet, analyses of population demographic data rarely consider processes occurring across spatial scales, potentially limiting the effectiveness of adaptive management. Therefore, we developed a method to identify hierarchical levels of organization (i.e., populations) to define multiple spatial scales, specifically intended to help guide appropriate conservation and management actions. This approach can support mobile species with high site fidelity where surveys occur on birthing/breeding grounds or migratory stopovers. Our approach used a graphbased clustering algorithm (Spatial K’luster Analysis by Tree Edge Removal) that explicitly included habitat selection information at multiple scales and further refined with constraint-based rules. We applied these concepts to greater sage-grouse leks (breeding grounds), a species of conservation concern, in two different ecological contexts (Nevada and Wyoming, USA). The constraint-based rules accounted for inter-lek movement distances based on literature and field studies in Nevada from 2012 to 2016, included methods to support a spatially balanced monitoring design, and identified barriers to movements among leks based on resistance surfaces. We evaluated the performance of our hierarchical clusters in Nevada using independent data from radio-marked sage-grouse, and we found the finest-scaled cluster level captured ~90% of sagegrouse movements and mid-level scales captured ~97–99% of movements. We expected comparable performance for Wyoming, where we lacked radio-marked sage-grouse for an evaluation, because genetic studies estimate similar dispersal distances to our ~15 km inter-lek movement distance in Nevada. For sage-grouse and other mobile species with high site fidelity, our approach to defining these frameworks could prove valuable for conservation and management applications, such as improving estimation of scale-dependent population trends and guiding the prescription of management actions at spatial scales that align with identified threats. Specific to sage-grouse, our analysis sets the stage for designing a monitoring framework that relies on comparison of short- and long-term population trends across our defined spatial scales and identifies and disentangles factors driving local (e.g., habitat quality) and regional (e.g., climate) population changes, thereby supporting scale-dependent management and research needs for adaptive management practices.

Nevada, Wyoming↗
USGS2021

Spatial network clustering reveals elk population structure and local variation in prevalence of chronic wasting disease

Spatial organization plays prominent roles in disease transmission, genetics, and demography of wildlife populations and is therefore an important consideration not only for wildlife management, but also for inference about populations and processes. We used hierarchical agglomerative clustering of a spatial graph network to partition Wind Cave National Park (WICA) into five regions used by 163 female elk ( Cervus elaphus ) marked with global positioning system collars during 2005–08 and 2011–13. We grouped elk based on differential use of the five regions, developed a priori models for inter-group variation in the occurrence of chronic wasting disease (CWD), and used Akaike's information criterion to compare models and stratify regions. Previous descriptions of elk population structure, which have been based on social contact or overlap of individual ranges, have distinguished spatially disjunct population subsets. Constructing hierarchical partitions of the landscape enabled us to also discern and describe overlapping and nested subsets. During 2016–18, apparent park-wide prevalence of CWD was 0.18 (90% CI = [0.146, 0.182]); however, prevalence within three spatial strata used primarily by different elk ranged from 0.03 ([0.008, 0.074]) to 0.29 ([0.211, 0.375]). In context with published estimates of recruitment, predation, and anthropogenic mortality, such differences in prevalence equate to increasing local abundance of elk in southwestern WICA, stable to declining abundance in the west/northwest, and rapidly declining abundance in the east. Despite the modest size of WICA (11,357 ha), park-wide averages conflate effects of elk distribution and disease, obscuring spatial patterns with profound implications for study and management of elk and CWD. Graph networks have been used widely in ecology to describe such phenomena as social relationships, connectivity of habitat patches, animal movements, and the spread of disease. Extension to partitioning of geographic range is straightforward but entails different considerations. We discuss allocation of sampling effort, construction of an initial partition, specification of a model for graph cohesion, selection of a clustering algorithm, and identification of useful partitions.

North Dakota↗
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

Quantifying flow and nonflow management impacts on an endangered fish by integrating data, research, and expert opinion

Managers charged with recovering endangered species in regulated river segments often have limited flexibility to alter flow regimes and want estimates of the expected population benefits associated with both flow and nonflow management actions. Disentangling impacts on different life stages from concurrently applied actions is essential for determining the effectiveness of each action, but difficult without models that integrate multiple information sources. Here, we develop and fit an integrated population model for endangered Rio Grande Silvery Minnow ( Hybognathus amarus ) in the Middle Rio Grande, New Mexico. We integrate catch per unit effort monitoring data collected during 2002–2018 with population estimates, data collected during rescue of minnow from drying pools, habitat availability estimates, laboratory results, releases of hatchery reared minnow, and expert opinion. We use expert elicitation to develop a larval carrying capacity index as an informed proxy for the complex interactions among flow, habitat, and life history in this species. We evaluate the model using out-of-sample forecasts of 2019 and 2020, develop an algorithm to identify supplemental water releases that maximize benefits to the minnow, and quantify the effectiveness of various actions. Experts generally agreed on the duration and timing of flow requirements and disagreed regarding the importance of different magnitudes. The integrated model with the larval carrying capacity index outperformed two alternative models in forecasting catch in 2019 and 2020. The model estimates that minnow abundance varied by more than three orders of magnitude between 2002 and 2018 and that in a few years recruitment was limited by spawner abundance. Evaluation of the expected benefits of flow and nonflow management actions to fall population abundance across different years suggests that efficient addition of water to the base hydrograph is the most effective action in most, but not all years. Many actions are effective only under certain hydrologic and population conditions and the effectiveness of different actions varies in different sections of the study area. Widespread water extraction and river regulation combined with periodic drought and ongoing climate change may necessitate creative management of federally listed fish species in arid systems informed by thorough analyses of management effectiveness.

Ecosphere↗
USGS2022

Models combining multiple scales of inference capture hydrologic and climatic drivers of riparian tree distributions

Predicting species geographic distributions is key to managing invasive species, conserving biodiversity, and understanding species' environmental requirements. Species distribution models (SDMs) commonly focus on climatic predictors, but other environmental factors can also be essential, particularly for species with specialized habitats defined by hydrologic, topographic, or edaphic conditions (e.g., riparian, wetland, alpine, coastal, serpentine). Here, we demonstrate a novel approach for capturing strong effects of both hydrologic and climatic predictors in SDMs for riparian plants, by merging analyses targeted at environmental drivers within riparian ecosystems and across the western USA (3.8 × 10 6 km 2 ). We developed presence-background SDMs from five algorithms for three invasive riparian trees ( Tamarix ramossisima / chinensis [saltcedar], Elaeagnus angustifolia [Russian olive], and Ulmus pumila [Siberian elm]) and three native Populus spp. (cottonwoods). We used separate background datasets to develop models with different spatial scales of inference: (1) spatially filtered random points to represent available habitat across the study area and (2) target-group points from Salix (willow) occurrences to represent available riparian habitat. Random-background models captured hydrologic drivers of riparian tree distributions relative to the largely upland western USA, whereas Salix -background models captured climatic drivers within the context of riparian ecosystems. Combining predictions from the two backgrounds identified hydrologically suitable habitats within climatically suitable regions, resulting in fewer false “absences” than either background alone, improving predictions over previous SDMs, and providing more complete information to guide management decisions. Surprisingly, the predicted habitat for U. pumila , a newly recognized riparian invader, was as or more extensive than Populus deltoides / fremontii , T. ramossisima / chinensis , and E. angustifolia , the most common riparian tree complexes in the western USA. Watersheds constituting 20% of U. pumila predicted habitat contained no occurrence records, indicating high risk of future and unrecognized invasions. Combining models from random and ecosystem-specific target-group backgrounds may improve SDMs for species from many specialized habitats, providing a method to link predicted distributions to localized geographic features while capturing broad-scale climatic requirements.

western United States↗