Geology Reports⌕ Search

Geology topics

B. Gardner

Publications and source records attributed to B. Gardner.

9 recordsLinked to original sources
USGS2024

Merging integrated population models and individual-based models to project population dynamics of recolonizing species

Recolonizing species exhibit unique population dynamics, namely dispersal to and colonization of new areas, that have important implications for management. A resulting challenge is how to simultaneously model demographic and movement processes so that recolonizing species can be accurately projected over time and space. We introduce a framework for spatially explicit projection modeling that harnesses the rigorous parameter estimation made possible by an integrated population model (IPM) and the flexible movement modeling made possible by an individual-based model (IBM). Our framework has two components: [1] a Bayesian IPM-driven age- and state-structured population model that governs the population state process and estimation of demographic rates, and [2] an IBM-driven spatial model that allows for the projection of dispersal and habitat colonization. We applied this model framework to estimate current and project future dynamics of gray wolves ( Canis lupus ) in Washington State, USA. We used data from 74 telemetered wolves and yearly pup and pack counts to parameterize the model, and then projected statewide dynamics over 50 years. Mean population growth was 1.29 (95 % Bayesian Credible Interval = 1.26–1.33) during initial recolonization from 2009 to 2020 and decreased to 1.02 (95 % Prediction Interval = 0.98–1.04) in the projection period (2021–2070). Our results suggest that gray wolves have an ~100 % probability of colonizing the last of Washington State's three specified recovery regions by 2030, regardless of alternative assumptions about how dispersing wolves select new territories. Our spatially explicit projection model can be used to project the dynamics of any species for which spatial spread is an important driver of population dynamics.

Washington↗
USGS2022

Integrated animal movement and spatial capture–recapture models: Simulation, implementation, and inference

Over the last decade, spatial capture–recapture (SCR) models have become widespread for estimating demographic parameters in ecological studies. However, the underlying assumptions about animal movement and space use are often not realistic. This is a missed opportunity because interesting ecological questions related to animal space use, habitat selection, and behavior cannot be addressed with most SCR models, despite the fact that the data collected in SCR studies — individual animals observed at specific locations and times — can provide a rich source of information about these processes and how they relate to demographic rates. We developed SCR models that integrated more complex movement processes that are typically inferred from telemetry data, including a simple random walk, correlated random walk (i.e., short-term directional persistence), and habitat-driven Langevin diffusion. We demonstrated how to formulate, simulate from, and fit these models with standard SCR data using data-augmented Bayesian analysis methods. We evaluated their performance through a simulation study, in which we varied the detection, movement, and resource selection parameters. We also examined different numbers of sampling occasions and assessed performance gains when including auxiliary location data collected from telemetered individuals. Across all scenarios, the integrated SCR movement models performed well in terms of abundance, detection, and movement parameter estimation. We found little difference in bias for the simple random walk model when reducing the number of sampling occasions from T = 25 to T = 15. We found some bias in movement parameter estimates under several of the correlated random walk scenarios, but incorporating auxiliary location data improved parameter estimates and significantly improved mixing during model fitting. The Langevin movement model was able to recover resource selection parameters from standard SCR data, which is particularly appealing because it explicitly links the individual-level movement process with habitat selection and population density. We focused on closed population models, but the movement models developed here can be extended to open SCR models. The movement process models could also be easily extended to accommodate additional “building blocks” of random walks, such as central tendency (e.g., territoriality) or multiple movement behavior states, thereby providing a flexible and coherent framework for linking animal movement behavior to population dynamics, density, and distribution.

Ecology↗
USGS2021

Evaluation of camera trap-based abundance estimators for unmarked populations

Estimates of species abundance are critical to understand population processes and to assess and select management actions. However, capturing and marking individuals for abundance estimation, while providing robust information, can be economically and logistically prohibitive, particularly for species with cryptic behavior. Camera traps can be used to collect data at temporal and spatial scales necessary for estimating abundance, but the use of camera traps comes with limitations when target species are not uniquely identifiable (i.e., “unmarked”). Abundance estimation is particularly useful in the management of invasive species, with herpetofauna being recognized as some of the most pervasive and detrimental invasive vertebrate species. However, the use of camera traps for these taxa presents additional challenges with relevancy across multiple taxa. It is often necessary to use lures to attract animals in order to obtain sufficient observations, yet lure attraction can influence species’ landscape use and potentially induce bias in abundance estimators. We investigated these challenges and assessed the feasibility of obtaining reliable abundance estimates using camera-trapping data on a population of invasive brown treesnakes ( Boiga irregularis ) in Guam. Data were collected using camera traps in an enclosed area where snakes were subject to high-intensity capture–recapture effort, resulting in presumed abundance of 116 snakes (density = 23/ha). We then applied spatial count, random encounter and staying time, space to event, and instantaneous sampling estimators to photo-capture data to estimate abundance and compared estimates to our presumed abundance. We found that all estimators for unmarked populations performed poorly, with inaccurate or imprecise abundance estimates that limit their usefulness for management in this system. We further investigated the sensitivity of these estimators to the use of lures (i.e., violating the assumption that animal behavior is unchanged by sampling) and camera density in a simulation study. Increasing the effective distances of a lure (i.e., lure attraction) and camera density both resulted in biased abundance estimates. Each estimator rarely recovered truth or suffered from convergence issues. Our results indicate that, when limited to unmarked estimators and the use of lures, camera traps alone are unlikely to produce abundance estimates with utility for brown treesnake management.

Ecological Applications↗
USGS2010

Estimating black bear density using DNA data from hair snares

DNA-based mark-recapture has become a methodological cornerstone of research focused on bear species. The objective of such studies is often to estimate population size; however, doing so is frequently complicated by movement of individual bears. Movement affects the probability of detection and the assumption of closure of the population required in most models. To mitigate the bias caused by movement of individuals, population size and density estimates are often adjusted using ad hoc methods, including buffering the minimum polygon of the trapping array. We used a hierarchical, spatial capturerecapture model that contains explicit components for the spatial-point process that governs the distribution of individuals and their exposure to (via movement), and detection by, traps. We modeled detection probability as a function of each individual's distance to the trap and an indicator variable for previous capture to account for possible behavioral responses. We applied our model to a 2006 hair-snare study of a black bear (Ursus americanus) population in northern New York, USA. Based on the microsatellite marker analysis of collected hair samples, 47 individuals were identified. We estimated mean density at 0.20 bears/km 2 . A positive estimate of the indicator variable suggests that bears are attracted to baited sites; therefore, including a trap-dependence covariate is important when using bait to attract individuals. Bayesian analysis of the model was implemented in WinBUGS, and we provide the model specification. The model can be applied to any spatially organized trapping array (hair snares, camera traps, mist nests, etc.) to estimate density and can also account for heterogeneity and covariate information at the trap or individual level. ?? The Wildlife Society.

Journal of Wildlife Management↗
USGS2010

Predicting species distributions from checklist data using site-occupancy models

Aim: (1) To increase awareness of the challenges induced by imperfect detection, which is a fundamental issue in species distribution modelling; (2) to emphasize the value of replicate observations for species distribution modelling; and (3) to show how 'cheap' checklist data in faunal/floral databases may be used for the rigorous modelling of distributions by site-occupancy models. Location: Switzerland. Methods: We used checklist data collected by volunteers during 1999 and 2000 to analyse the distribution of the blue hawker, Aeshna cyanea (Odonata, Aeshnidae), a common dragonfly in Switzerland. We used data from repeated visits to 1-ha pixels to derive 'detection histories' and apply site-occupancy models to estimate the 'true' species distribution, i.e. corrected for imperfect detection. We modelled blue hawker distribution as a function of elevation and year and its detection probability of elevation, year and season. Results: The best model contained cubic polynomial elevation effects for distribution and quadratic effects of elevation and season for detectability. We compared the site-occupancy model with a conventional distribution model based on a generalized linear model, which assumes perfect detectability (p = 1). The conventional distribution map looked very different from the distribution map obtained using site-occupancy models that accounted for the imperfect detection. The conventional model underestimated the species distribution by 60%, and the slope parameters of the occurrence-elevation relationship were also underestimated when assuming p = 1. Elevation was not only an important predictor of blue hawker occurrence, but also of the detection probability, with a bell-shaped relationship. Furthermore, detectability increased over the season. The average detection probability was estimated at only 0.19 per survey. Main conclusions: Conventional species distribution models do not model species distributions per se but rather the apparent distribution, i.e. an unknown proportion of species distributions. That unknown proportion is equivalent to detectability. Imperfect detection in conventional species distribution models yields underestimates of the extent of distributions and covariate effects that are biased towards zero. In addition, patterns in detectability will erroneously be ascribed to species distributions. In contrast, site-occupancy models applied to replicated detection/non-detection data offer a powerful framework for making inferences about species distributions corrected for imperfect detection. The use of 'cheap' checklist data greatly enhances the scope of applications of this useful class of models. ?? 2010 Blackwell Publishing Ltd.

Journal of Biogeography↗
USGS2009

Hierarchical models for estimating density from DNA mark-recapture studies

Genetic sampling is increasingly used as a tool by wildlife biologists and managers to estimate abundance and density of species. Typically, DNA is used to identify individuals captured in an array of traps ( e. g., baited hair snares) from which individual encounter histories are derived. Standard methods for estimating the size of a closed population can be applied to such data. However, due to the movement of individuals on and off the trapping array during sampling, the area over which individuals are exposed to trapping is unknown, and so obtaining unbiased estimates of density has proved difficult. We propose a hierarchical spatial capture-recapture model which contains explicit models for the spatial point process governing the distribution of individuals and their exposure to (via movement) and detection by traps. Detection probability is modeled as a function of each individual's distance to the trap. We applied this model to a black bear (Ursus americanus) study conducted in 2006 using a hair-snare trap array in the Adirondack region of New York, USA. We estimated the density of bears to be 0.159 bears/km2, which is lower than the estimated density (0.410 bears/km2) based on standard closed population techniques. A Bayesian analysis of the model is fully implemented in the software program WinBUGS.

Ecology↗
USGS2008

Spatial and temporal statistical analysis of bycatch data: Patterns of sea turtle bycatch in the North Atlantic

Loggerhead (Caretta caretta) and leatherback (Dermochelys coriacea) sea turtle distributions and movements in offshore waters of the western North Atlantic are not well understood despite continued efforts to monitor, survey, and observe them. Loggerhead and leatherback sea turtles are listed as endangered by the World Conservation Union, and thus anthropogenic mortality of these species, including fishing, is of elevated interest. This study quantifies spatial and temporal patterns of sea turtle bycatch distributions to identify potential processes influencing their locations. A Ripley's K function analysis was employed on the NOAA Fisheries Atlantic Pelagic Longline Observer Program data to determine spatial, temporal, and spatio-temporal patterns of sea turtle bycatch distributions within the pattern of the pelagic fishery distribution. Results indicate that loggerhead and leatherback sea turtle catch distributions change seasonally, with patterns of spatial clustering appearing from July through October. The results from the space-time analysis indicate that sea turtle catch distributions are related on a relatively fine scale (30-200 km and 1-5 days). The use of spatial and temporal point pattern analysis, particularly K function analysis, is a novel way to examine bycatch data and can be used to inform fishing practices such that fishing could still occur while minimizing sea turtle bycatch. ?? 2008 NRC.

Canadian Journal of Fisheries and Aquatic Sciences↗
USGS2008

Hierarchical modeling of bycatch rates of sea turtles in the western North Atlantic

Previous studies indicate that the locations of the endangered loggerhead Caretta caretta and critically endangered leatherback Dermochelys coriacea sea turtles are influenced by water temperatures, and that incidental catch rates in the pelagic longline fishery vary by region. We present a Bayesian hierarchical model to examine the effects of environmental variables, including water temperature, on the number of sea turtles captured in the US pelagic longline fishery in the western North Atlantic. The modeling structure is highly flexible, utilizes a Bayesian model selection technique, and is fully implemented in the software program WinBUGS. The number of sea turtles captured is modeled as a zero-inflated Poisson distribution and the model incorporates fixed effects to examine region-specific differences in the parameter estimates. Results indicate that water temperature, region, bottom depth, and target species are all significant predictors of the number of loggerhead sea turtles captured. For leatherback sea turtles, the model with only target species had the most posterior model weight, though a re-parameterization of the model indicates that temperature influences the zero-inflation parameter. The relationship between the number of sea turtles captured and the variables of interest all varied by region. This suggests that management decisions aimed at reducing sea turtle bycatch may be more effective if they are spatially explicit. ?? Inter-Research 2008.

Endangered Species Research↗