Search USGS⌕ Search

Geology topics

Robert Dorazio

Publications and source records attributed to Robert Dorazio.

At least 19 recordsLinked to original sources
USGS2019

State-space models to infer movements and behavior of fish detected in a spatial array of acoustic receivers

We developed state-space models for inferring movements and behaviors of fish implanted with acoustic transmitters and detected within a spatial array of stationary acoustic receivers. In these models fish movements and behavior switching are specified using a hidden Markov model of the changes in an individual's latent activity center. The observed number of detections at each acoustic receiver is modeled as a function of the distance between the receiver and a fish's activity center. To illustrate these models, we analyzed the detections of individual Gulf sturgeon that moved within a spatial array of acoustic receivers placed in the Suwannee River, Florida. Our models of these detections provided estimates of the locations of individual sturgeon and the periods when individuals switched between spawning and resting behaviors. We anticipate that these models will be used to analyze acoustic surveys of other species and to help design new surveys.

Canadian Journal of Fisheries and Aquatic Sciences↗
USGS2018

Time series sightability modeling of animal populations

Logistic regression models—or “sightability models”—fit to detection/non-detection data from marked individuals are often used to adjust for visibility bias in later detection-only surveys, with population abundance estimated using a modified Horvitz-Thompson (mHT) estimator. More recently, a model-based alternative for analyzing combined detection/non-detection and detection-only data was developed. This approach seemed promising, since it resulted in similar estimates as the mHT when applied to data from moose ( Alces alces ) surveys in Minnesota. More importantly, it provided a framework for developing flexible models for analyzing multiyear detection-only survey data in combination with detection/non-detection data. During initial attempts to extend the model-based approach to multiple years of detection-only data, we found that estimates of detection probabilities and population abundance were sensitive to the amount of detection-only data included in the combined (detection/non-detection and detection-only) analysis. Subsequently, we developed a robust hierarchical modeling approach where sightability model parameters are informed only by the detection/non-detection data, and we used this approach to fit a fixed-effects model (FE model) with year-specific parameters and a temporally-smoothed model (TS model) that shares information across years via random effects and a temporal spline. The abundance estimates from the TS model were more precise, with decreased interannual variability relative to the FE model and mHT abundance estimates, illustrating the potential benefits from model-based approaches that allow information to be shared across years.

PLoS ONE↗
USGS2018

Surveys of environmental DNA (eDNA): a new approach to estimate occurrence in Vulnerable manatee populations

Environmental DNA (eDNA) detection is a technique used to non-invasively detect cryptic, low density, or logistically difficult-to-study species, such as imperiled manatees. For eDNA measurement, genetic material shed into the environment is concentrated from water samples and analyzed for the presence of target species. Cytochrome b quantitative PCR and droplet digital PCR eDNA assays were developed for the 3 Vulnerable manatee species: African, Amazonian, and both subspecies of the West Indian (Florida and Antillean) manatee. Environmental DNA assays can help to delineate manatee habitat ranges, high use areas, and seasonal population changes. To validate the assay, water was analyzed from Florida’s east coast containing a high-density manatee population and produced 31564 DNA molecules l -1 on average and high occurrence (ψ) and detection (p) estimates (ψ = 0.84 [0.40-0.99]; p = 0.99 [0.95-1.00]; limit of detection 3 copies µl -1 ). Similar occupancy estimates were produced in the Florida Panhandle (ψ = 0.79 [0.54-0.97]) and Cuba (ψ = 0.89 [0.54-1.00]), while occupancy estimates in Cameroon were lower (ψ = 0.49 [0.09-0.95]). The eDNA-derived detection estimates were higher than those generated using aerial survey data on the west coast of Florida and may be effective for population monitoring. Subsequent eDNA studies could be particularly useful in locations where manatees are (1) difficult to identify visually (e.g. the Amazon River and Africa), (2) are present in patchy distributions or are on the verge of extinction (e.g. Jamaica, Haiti), and (3) where repatriation efforts are proposed (e.g. Brazil, Guadeloupe). Extension of these eDNA techniques could be applied to other imperiled marine mammal populations such as African and Asian dugongs.

Endangered Species Research↗
USGS2017

eDNAoccupancy: An R package for multi-scale occupancy modeling of environmental DNA data

In this article we describe eDNAoccupancy, an R package for fitting Bayesian, multi-scale occupancy models. These models are appropriate for occupancy surveys that include three, nested levels of sampling: primary sample units within a study area, secondary sample units collected from each primary unit, and replicates of each secondary sample unit. This design is commonly used in occupancy surveys of environmental DNA (eDNA). eDNAoccupancy allows users to specify and fit multi-scale occupancy models with or without covariates, to estimate posterior summaries of occurrence and detection probabilities, and to compare different models using Bayesian model-selection criteria. We illustrate these features by analyzing two published data sets: eDNA surveys of a fungal pathogen of amphibians and eDNA surveys of an endangered fish species.

Molecular Ecology Resources↗
USGS2017

Concepts: Integrating population survey data from different spatial scales, sampling methods, and species

Conservationists and managers are continually under pressure from the public, the media, and political policy makers to provide “tiger numbers,” not just for protected reserves, but also for large spatial scales, including landscapes, regions, states, nations, and even globally. Estimating the abundance of tigers within relatively small areas (e.g., protected reserves) is becoming increasingly tractable (see Chaps. 9 and 10), but doing so for larger spatial scales still presents a formidable challenge. Those who seek “tiger numbers” are often not satisfied by estimates of tiger occupancy alone, regardless of the reliability of the estimates (see Chaps. 4 and 5). As a result, wherever tiger conservation efforts are underway, either substantially or nominally, scientists and managers are frequently asked to provide putative large-scale tiger numbers based either on a total count or on an extrapolation of some sort (see Chaps. 1 and 2).

Book chapter↗
USGS2017

Concepts and practices: Estimating abundance of prey species using hierarchical model-based approaches

Tigers predominantly prey on large ungulate species, such as sambar (Cervus unicolor), red deer (Cervus elaphus), gaur (Bos gaurus), banteng (Bos javanicus), chital (Axis axis), muntjac (Muntiacus muntjak), wild pig (Sus scrofa), and bearded pig (Sus barbatus). The density of a tiger population is strongly correlated with the density of such prey species (Karanth et al. 2004). In the absence of direct hunting of tigers, abundance of prey in an area is the key determinant of the “carrying capacity” of that area for tigers (Chap. 2). Accurate estimates of prey abundance are often needed to assess the potential number of tigers a conservation area can support.

Book chapter↗
USGS2017

A hierarchical model for estimating the spatial distribution and abundance of animals detected by continuous-time recorders

Motivation Several spatial capture-recapture (SCR) models have been developed to estimate animal abundance by analyzing the detections of individuals in a spatial array of traps. Most of these models do not use the actual dates and times of detection, even though this information is readily available when using continuous-time recorders, such as microphones or motion-activated cameras. Instead most SCR models either partition the period of trap operation into a set of subjectively chosen discrete intervals and ignore multiple detections of the same individual within each interval, or they simply use the frequency of detections during the period of trap operation and ignore the observed times of detection. Both practices make inefficient use of potentially important information in the data. Model and data analysis We developed a hierarchical SCR model to estimate the spatial distribution and abundance of animals detected with continuous-time recorders. Our model includes two kinds of point processes: a spatial process to specify the distribution of latent activity centers of individuals within the region of sampling and a temporal process to specify temporal patterns in the detections of individuals. We illustrated this SCR model by analyzing spatial and temporal patterns evident in the camera-trap detections of tigers living in and around the Nagarahole Tiger Reserve in India. We also conducted a simulation study to examine the performance of our model when analyzing data sets of greater complexity than the tiger data. Benefits Our approach provides three important benefits: First, it exploits all of the information in SCR data obtained using continuous-time recorders. Second, it is sufficiently versatile to allow the effects of both space use and behavior of animals to be specified as functions of covariates that vary over space and time. Third, it allows both the spatial distribution and abundance of individuals to be estimated, effectively providing a species distribution model, even in cases where spatial covariates of abundance are unknown or unavailable. We illustrated these benefits in the analysis of our data, which allowed us to quantify differences between nocturnal and diurnal activities of tigers and to estimate their spatial distribution and abundance across the study area. Our continuous-time SCR model allows an analyst to specify many of the ecological processes thought to be involved in the distribution, movement, and behavior of animals detected in a spatial trapping array of continuous-time recorders. We plan to extend this model to estimate the population dynamics of animals detected during multiple years of SCR surveys.

Western Ghats, Nagarahole Tiger Reserve↗
USGS2017

Integrated species distribution models: combining presence-background data and site-occupancy data with imperfect detection

Two main sources of data for species distribution models (SDMs) are site-occupancy (SO) data from planned surveys, and presence-background (PB) data from opportunistic surveys and other sources. SO surveys give high quality data about presences and absences of the species in a particular area. However, due to their high cost, they often cover a smaller area relative to PB data, and are usually not representative of the geographic range of a species. In contrast, PB data is plentiful, covers a larger area, but is less reliable due to the lack of information on species absences, and is usually characterised by biased sampling. Here we present a new approach for species distribution modelling that integrates these two data types. We have used an inhomogeneous Poisson point process as the basis for constructing an integrated SDM that fits both PB and SO data simultaneously. It is the first implementation of an Integrated SO–PB Model which uses repeated survey occupancy data and also incorporates detection probability. The Integrated Model's performance was evaluated, using simulated data and compared to approaches using PB or SO data alone. It was found to be superior, improving the predictions of species spatial distributions, even when SO data is sparse and collected in a limited area. The Integrated Model was also found effective when environmental covariates were significantly correlated. Our method was demonstrated with real SO and PB data for the Yellow-bellied glider ( Petaurus australis ) in south-eastern Australia, with the predictive performance of the Integrated Model again found to be superior. PB models are known to produce biased estimates of species occupancy or abundance. The small sample size of SO datasets often results in poor out-of-sample predictions. Integrated models combine data from these two sources, providing superior predictions of species abundance compared to using either data source alone. Unlike conventional SDMs which have restrictive scale-dependence in their predictions, our Integrated Model is based on a point process model and has no such scale-dependency. It may be used for predictions of abundance at any spatial-scale while still maintaining the underlying relationship between abundance and area.

Methods in Ecology and Evolution↗
USGS2016

State-dependent resource harvesting with lagged information about system states

Markov decision processes (MDPs), which involve a temporal sequence of actions conditioned on the state of the managed system, are increasingly being applied in natural resource management. This study focuses on the modification of a traditional MDP to account for those cases in which an action must be chosen after a significant time lag in observing system state, but just prior to a new observation. In order to calculate an optimal decision policy under these conditions, possible actions must be conditioned on the previous observed system state and action taken. We show how to solve these problems when the state transition structure is known and when it is uncertain. Our focus is on the latter case, and we show how actions must be conditioned not only on the previous system state and action, but on the probabilities associated with alternative models of system dynamics. To demonstrate this framework, we calculated and simulated optimal, adaptive policies for MDPs with lagged states for the problem of deciding annual harvest regulations for mallards (Anas platyrhynchos) in the United States. In this particular example, changes in harvest policy induced by the use of lagged information about system state were sufficient to maintain expected management performance (e.g. population size, harvest) even in the face of an uncertain system state at the time of a decision.

PLoS ONE↗
USGS2016

Incorporating imperfect detection into joint models of communites: A response to Warton et al.

Warton et al. [1] advance community ecology by describing a statistical framework that can jointly model abundances (or distributions) across many taxa to quantify how community properties respond to environmental variables. This framework specifies the effects of both measured and unmeasured (latent) variables on the abundance (or occurrence) of each species. Latent variables are random effects that capture the effects of both missing environmental predictors and correlations in parameter values among different species. As presented in Warton et al. , however, the joint modeling framework fails to account for the common problem of detection or measurement errors that always accompany field sampling of abundance or occupancy, and are well known to obscure species- and community-level inferences.

Trends in Ecology and Evolution↗
USGS2015

Occupancy estimation for rare species using a spatially-adaptive sampling design

Summary 1. Spatially clustered populations create unique challenges for conservation monitoring programmes. Advances in methodology typically are focused on either the design or the modelling stage of the study but do not involve integration of both. 2. We integrate adaptive cluster sampling and spatial occupancy modelling by developing two models to handle the dependence induced by cluster sampling. We compare these models to scenarios using simple random sampling and traditional occupancy models via simulation and data collected on a rare plant species, Tamarix ramosissima , found in China. 3. Our simulations show a marked improvement in confidence interval coverage for the new models combined with cluster sampling compared to simple random sampling and traditional occupancy models, with greatest improvement in the presence of low detection probability and spatial correlation in occupancy. 4. Accounting for the design using the simple cluster random-effects model reduces bias considerably, and full spatial modelling reduces bias further, especially for large n when the spatial covariance parameters can be estimated reliably. Both new models build on the strength of occupancy modelling and adaptive sampling and perform at least as well, and often better, than occupancy modelling alone. 5. We believe our approach is unique and potentially useful for a variety of studies directed at patchily distributed, clustered or rare species exhibiting spatial variation.

Methods in Ecology and Evolution↗
USGS2015

Statistical models for the analysis and design of digital polymerase chain (dPCR) experiments

Statistical methods for the analysis and design of experiments using digital PCR (dPCR) have received only limited attention and have been misused in many instances. To address this issue and to provide a more general approach to the analysis of dPCR data, we describe a class of statistical models for the analysis and design of experiments that require quantification of nucleic acids. These models are mathematically equivalent to generalized linear models of binomial responses that include a complementary, log–log link function and an offset that is dependent on the dPCR partition volume. These models are both versatile and easy to fit using conventional statistical software. Covariates can be used to specify different sources of variation in nucleic acid concentration, and a model’s parameters can be used to quantify the effects of these covariates. For purposes of illustration, we analyzed dPCR data from different types of experiments, including serial dilution, evaluation of copy number variation, and quantification of gene expression. We also showed how these models can be used to help design dPCR experiments, as in selection of sample sizes needed to achieve desired levels of precision in estimates of nucleic acid concentration or to detect differences in concentration among treatments with prescribed levels of statistical power.

Analytical Chemistry↗
USGS2015

Bayesian data analysis in population ecology: motivations, methods, and benefits

During the 20th century ecologists largely relied on the frequentist system of inference for the analysis of their data. However, in the past few decades ecologists have become increasingly interested in the use of Bayesian methods of data analysis. In this article I provide guidance to ecologists who would like to decide whether Bayesian methods can be used to improve their conclusions and predictions. I begin by providing a concise summary of Bayesian methods of analysis, including a comparison of differences between Bayesian and frequentist approaches to inference when using hierarchical models. Next I provide a list of problems where Bayesian methods of analysis may arguably be preferred over frequentist methods. These problems are usually encountered in analyses based on hierarchical models of data. I describe the essentials required for applying modern methods of Bayesian computation, and I use real-world examples to illustrate these methods. I conclude by summarizing what I perceive to be the main strengths and weaknesses of using Bayesian methods to solve ecological inference problems.

Population Ecology↗
USGS2011

Occupancy and abundance of wintering birds in a dynamic agricultural landscape

Assessing wildlife management action requires monitoring populations, and abundance often is the parameter monitored. Recent methodological advances have enabled estimation of mean abundance within a habitat using presence–absence or count data obtained via repeated visits to a sample of sites. These methods assume populations are closed and intuitively assume habitats within sites change little during a field season. However, many habitats are highly variable over short periods. We developed a variation of existing occupancy and abundance models that allows for extreme spatio‐temporal differences in habitat, and resulting changes in wildlife abundance, among sites and among visits to a site within a field season. We conducted our study in sugarcane habitat within the Everglades Agricultural Area southeast of Lake Okeechobee in south Florida. We counted wintering birds, primarily passerines, within 245 sites usually 5 times at each site during December 2006–March 2007. We estimated occupancy and mean abundance of birds in 6 vegetation states during the sugarcane harvest and allowed these parameters to vary temporally or spatially within a vegetation state. Occupancy and mean abundance of the common yellowthroat ( Geothlypis trichas ) was affected by structure of sugarcane and uncultivated edge vegetation (occupancy = 1.00 [ = 0.96–1.00] and mean abundance = 7.9 [ = 3.2–19.5] in tall sugarcane with tall edge vegetation versus 0.20 [ = 0.04–0.71] and 0.22 [ = 0.04–1.2], respectively, in short sugarcane with short edge vegetation in one half of the study area). Occupancy and mean abundance of palm warblers ( Dendroica palmarum ) were constant (occupancy = 1.00, = 0.69–1.00; mean abundance = 18, = 1–270). Our model may enable wildlife managers to assess rigorously effects of future edge habitat management on avian distribution and abundance within agricultural landscapes during winter or the breeding season. The model may also help wildlife managers make similar management decisions involving other dynamic habitats such as wetlands, prairies, and even forested areas if forest management or fires occur during the field season.

Florida↗
USGS2005

Estimating size and composition of biological communities by modeling the occurrence of species

We develop a model that uses repeated observations of a biological community to estimate the number and composition of species in the community. Estimators of community-level attributes are constructed from model-based estimators of occurrence of individual species that incorporate imperfect detection of individuals. Data from the North American Breeding Bird Survey are analyzed to illustrate the variety of ecologically important quantities that are easily constructed and estimated using our model-based estimators of species occurrence. In particular, we compute site-specific estimates of species richness that honor classical notions of species-area relationships. We suggest extensions of our model to estimate maps of occurrence of individual species and to compute inferences related to the temporal and spatial dynamics of biological communities.

Journal of the American Statistical Association↗