Geology topics

J. Andrew Royle

Publications and source records attributed to J. Andrew Royle.

At least 109 records · Page 6Linked to original sources

USGS2015

Book review: Analysis of capture–recapture data

Analysis of Capture–Recapture Data by McCrea and Morgan is an excellent, easy to read monograph about capture–recapture models. In this book, the authors provide a concise overview of traditional closed population capture–recapture models (Models M0, Mb, Mh, etc.), individual covariate models, and open population models such as the Cormack–Jolly–Seber, Jolly–Seber models, multi-state models, and more recent developments such as occupancy models, state-space models, and integrated population models. The authors write “In this book we aim to cover the many modern developments in the area of capture–recapture and related models, and to set them in historical context of relevant research over the past 100 years.” The book does a good job of achieving this objective. And it is a very easy to read because it is well organized and the writing is clear and concise. I would recommend this book as a reference for the quantitative ecologist or statistician interested in knowing what’s out there. Review info: Analysis of Capture–Recapture Data. By Rachel S. McCrea and Byron J. T. Morgan, 2015. ISBN: 978-14398-6590, xx+292 pp.

Journal of Agricultural, Biological, and Environme↗

USGS2015

Likelihood analysis of spatial capture-recapture models for stratified or class structured populations

We develop a likelihood analysis framework for fitting spatial capture-recapture (SCR) models to data collected on class structured or stratified populations. Our interest is motivated by the necessity of accommodating the problem of missing observations of individual class membership. This is particularly problematic in SCR data arising from DNA analysis of scat, hair or other material, which frequently yields individual identity but fails to identify the sex. Moreover, this can represent a large fraction of the data and, given the typically small sample sizes of many capture-recapture studies based on DNA information, utilization of the data with missing sex information is necessary. We develop the class structured likelihood for the case of missing covariate values, and then we address the scaling of the likelihood so that models with and without class structured parameters can be formally compared regardless of missing values. We apply our class structured model to black bear data collected in New York in which sex could be determined for only 62 of 169 uniquely identified individuals. The models containing sex-specificity of both the intercept of the SCR encounter probability model and the distance coefficient, and including a behavioral response are strongly favored by log-likelihood. Estimated population sex ratio is strongly influenced by sex structure in model parameters illustrating the importance of rigorous modeling of sex differences in capture-recapture models.

Ecosphere↗

USGS2015

Comparing spatial capture–recapture modeling and nest count methods to estimate orangutan densities in the Wehea Forest, East Kalimantan, Indonesia

Accurate information on the density and abundance of animal populations is essential for understanding species' ecology and for conservation planning, but is difficult to obtain. The endangered orangutan ( Pongo spp.) is an example; due to its elusive behavior and low densities, researchers have relied on methods that convert nest counts to orangutan densities and require substantial effort for reliable results. Camera trapping and spatial capture–recapture (SCR) models could provide an alternative but have not been used for primates. We compared density estimates calculated using the two methods for orangutans in the Wehea Forest, East Kalimantan, Indonesia. Camera trapping/SCR modeling produced a density estimate of 0.16 ± 0.09–0.29 indiv/km 2 , and nest counts produced a density estimate of 1.05 ± 0.18–6.01 indiv/km 2 . The large confidence interval of the nest count estimate is probably due to high variance in nest encounter rates, indicating the need for larger sample size and the substantial effort required to produce reliable results using this method. The SCR estimate produced a very low density estimate and had a narrower but still fairly wide confidence interval. This was likely due to unmodeled heterogeneity and small sample size, specifically a low number of individual captures and recaptures. We propose methodological fixes that could address these issues and improve precision. A comparison of the overall costs and benefits of the two methods suggests that camera trapping/SCR modeling can potentially be a useful tool for assessing the densities of orangutans and other elusive primates, and warrant further investigation to determine broad applicability and methodological adjustments needed.

East Kalimantan↗

USGS2015

An open-population hierarchical distance sampling model

Modeling population dynamics while accounting for imperfect detection is essential to monitoring programs. Distance sampling allows estimating population size while accounting for imperfect detection, but existing methods do not allow for direct estimation of demographic parameters. We develop a model that uses temporal correlation in abundance arising from underlying population dynamics to estimate demographic parameters from repeated distance sampling surveys. Using a simulation study motivated by designing a monitoring program for island scrub-jays ( Aphelocoma insularis ), we investigated the power of this model to detect population trends. We generated temporally autocorrelated abundance and distance sampling data over six surveys, using population rates of change of 0.95 and 0.90. We fit the data generating Markovian model and a mis-specified model with a log-linear time effect on abundance, and derived post hoc trend estimates from a model estimating abundance for each survey separately. We performed these analyses for varying number of survey points. Power to detect population changes was consistently greater under the Markov model than under the alternatives, particularly for reduced numbers of survey points. The model can readily be extended to more complex demographic processes than considered in our simulations. This novel framework can be widely adopted for wildlife population monitoring.

Ecology↗

USGS2014

Modelling non-Euclidean movement and landscape connectivity in highly structured ecological networks

Movement is influenced by landscape structure, configuration and geometry, but measuring distance as perceived by animals poses technical and logistical challenges. Instead, movement is typically measured using Euclidean distance, irrespective of location or landscape structure, or is based on arbitrary cost surfaces. A recently proposed extension of spatial capture-recapture (SCR) models resolves this issue using spatial encounter histories of individuals to calculate least-cost paths (ecological distance: Ecology , 94 , 2013, 287) thereby relaxing the Euclidean assumption. We evaluate the consequences of not accounting for movement heterogeneity when estimating abundance in highly structured landscapes, and demonstrate the value of this approach for estimating biologically realistic space-use patterns and landscape connectivity. We simulated SCR data in a riparian habitat network, using the ecological distance model under a range of scenarios where space-use in and around the landscape was increasingly associated with water (i.e. increasingly less Euclidean). To assess the influence of miscalculating distance on estimates of population size, we compared the results from the ecological and Euclidean distance based models. We then demonstrate that the ecological distance model can be used to estimate home range geometry when space use is not symmetrical. Finally, we provide a method for calculating landscape connectivity based on modelled species-landscape interactions generated from capture-recapture data. Using ecological distance always produced unbiased estimates of abundance. Explicitly modelling the strength of the species-landscape interaction provided a direct measure of landscape connectivity and better characterised true home range geometry. Abundance under the Euclidean distance model was increasingly (negatively) biased as space use was more strongly associated with water and, because home ranges are assumed to be symmetrical, produced poor characterisations of home range geometry and no information about landscape connectivity. The ecological distance SCR model uses spatially indexed capture-recapture data to estimate how activity patterns are influenced by landscape structure. As well as reducing bias in estimates of abundance, this approach provides biologically realistic representations of home range geometry, and direct information about species-landscape interactions. The incorporation of both structural (landscape) and functional (movement) components of connectivity provides a direct measure of species-specific landscape connectivity.

Methods in Ecology and Evolution↗

USGS2014

Small mammal use of native warm-season and non-native cool-season grass forage fields

Recent emphasis has been put on establishing native warm-season grasses for forage production because it is thought native warm-season grasses provide higher quality wildlife habitat than do non-native cool-season grasses. However, it is not clear whether native warm-season grass fields provide better resources for small mammals than currently are available in non-native cool-season grass forage production fields. We developed a hierarchical spatially explicit capture-recapture model to compare abundance of hispid cotton rats ( Sigmodon hispidus ), white-footed mice ( Peromyscus leucopus ), and house mice ( Mus musculus ) among 4 hayed non-native cool-season grass fields, 4 hayed native warm-season grass fields, and 4 native warm-season grass-forb ("wildlife") fields managed for wildlife during 2 summer trapping periods in 2009 and 2010 of the western piedmont of North Carolina, USA. Cotton rat abundance estimates were greater in wildlife fields than in native warm-season grass and non-native cool-season grass fields and greater in native warm-season grass fields than in non-native cool-season grass fields. Abundances of white-footed mouse and house mouse populations were lower in wildlife fields than in native warm-season grass and non-native cool-season grass fields, but the abundances were not different between the native warm-season grass and non-native cool-season grass fields. Lack of cover following haying in non-native cool-season grass and native warm-season grass fields likely was the key factor limiting small mammal abundance, especially cotton rats, in forage fields. Retention of vegetation structure in managed forage production systems, either by alternately resting cool-season and warm-season grass forage fields or by leaving unharvested field borders, should provide refugia for small mammals during haying events.

Wildlife Society Bulletin↗

USGS2014

Estimating landscape resistance to dispersal

Dispersal is an inherently spatial process that can be affected by habitat conditions in sites encountered by dispersers. Understanding landscape resistance to dispersal is important in connectivity studies and reserve design, but most existing methods use resistance functions with cost parameters that are subjectively chosen by the investigator. We develop an analytic approach allowing for direct estimation of resistance parameters that folds least cost path methods typically used in simulation approaches into a formal statistical model of dispersal distributions. The core of our model is a frequency distribution of dispersal distances expressed as least cost distance rather than Euclidean distance, and which includes terms for feature-specific costs to dispersal and sex (or other traits) of the disperser. The model requires only origin and settlement locations for multiple individuals, such as might be obtained from mark–recapture studies or parentage analyses, and maps of the relevant habitat features. To evaluate whether the model can estimate parameters correctly, we fit our model to data from simulated dispersers in three kinds of landscapes (in which resistance of environmental variables was categorical, continuous with a patchy configuration, or continuous in a trend pattern). We found maximum likelihood estimators of resistance and individual trait parameters to be approximately unbiased with moderate sample sizes. We applied the model to a small grizzly bear dataset to demonstrate how this approach could be used when the primary interest is in the prediction of costs and found that estimates were consistent with expectations based on bear ecology. Our method has important practical applications for testing hypotheses about dispersal ecology and can be used to inform connectivity planning efforts, via the resistance estimates and confidence intervals, which can be used to create a data-driven resistance surface.

Landscape Ecology↗

USGS2014

Estimating migratory connectivity of birds when re-encounter probabilities are heterogeneous

Understanding the biology and conducting effective conservation of migratory species requires an understanding of migratory connectivity – the geographic linkages of populations between stages of the annual cycle. Unfortunately, for most species, we are lacking such information. The North American Bird Banding Laboratory (BBL) houses an extensive database of marking, recaptures and recoveries, and such data could provide migratory connectivity information for many species. To date, however, few species have been analyzed for migratory connectivity largely because heterogeneous re-encounter probabilities make interpretation problematic. We accounted for regional variation in re-encounter probabilities by borrowing information across species and by using effort covariates on recapture and recovery probabilities in a multistate capture–recapture and recovery model. The effort covariates were derived from recaptures and recoveries of species within the same regions. We estimated the migratory connectivity for three tern species breeding in North America and over-wintering in the tropics, common (Sterna hirundo), roseate (Sterna dougallii), and Caspian terns (Hydroprogne caspia). For western breeding terns, model-derived estimates of migratory connectivity differed considerably from those derived directly from the proportions of re-encounters. Conversely, for eastern breeding terns, estimates were merely refined by the inclusion of re-encounter probabilities. In general, eastern breeding terns were strongly connected to eastern South America, and western breeding terns were strongly linked to the more western parts of the nonbreeding range under both models. Through simulation, we found this approach is likely useful for many species in the BBL database, although precision improved with higher re-encounter probabilities and stronger migratory connectivity. We describe an approach to deal with the inherent biases in BBL banding and re-encounter data to demonstrate that this large dataset is a valuable source of information about the migratory connectivity of the birds of North America.

Ecology and Evolution↗

USGS2014

Trap configuration and spacing influences parameter estimates in spatial capture-recapture models

An increasing number of studies employ spatial capture-recapture models to estimate population size, but there has been limited research on how different spatial sampling designs and trap configurations influence parameter estimators. Spatial capture-recapture models provide an advantage over non-spatial models by explicitly accounting for heterogeneous detection probabilities among individuals that arise due to the spatial organization of individuals relative to sampling devices. We simulated black bear ( Ursus americanus) populations and spatial capture-recapture data to evaluate the influence of trap configuration and trap spacing on estimates of population size and a spatial scale parameter, sigma, that relates to home range size. We varied detection probability and home range size, and considered three trap configurations common to large-mammal mark-recapture studies: regular spacing, clustered, and a temporal sequence of different cluster configurations (i.e., trap relocation). We explored trap spacing and number of traps per cluster by varying the number of traps. The clustered arrangement performed well when detection rates were low, and provides for easier field implementation than the sequential trap arrangement. However, performance differences between trap configurations diminished as home range size increased. Our simulations suggest it is important to consider trap spacing relative to home range sizes, with traps ideally spaced no more than twice the spatial scale parameter. While spatial capture-recapture models can accommodate different sampling designs and still estimate parameters with accuracy and precision, our simulations demonstrate that aspects of sampling design, namely trap configuration and spacing, must consider study area size, ranges of individual movement, and home range sizes in the study population.

PLoS ONE↗

USGS2014

Inferences about population dynamics from count data using multi-state models: A comparison to capture-recapture approaches

Wildlife populations consist of individuals that contribute disproportionately to growth and viability. Understanding a population's spatial and temporal dynamics requires estimates of abundance and demographic rates that account for this heterogeneity. Estimating these quantities can be difficult, requiring years of intensive data collection. Often, this is accomplished through the capture and recapture of individual animals, which is generally only feasible at a limited number of locations. In contrast, N-mixture models allow for the estimation of abundance, and spatial variation in abundance, from count data alone. We extend recently developed multistate, open population N-mixture models, which can additionally estimate demographic rates based on an organism's life history characteristics. In our extension, we develop an approach to account for the case where not all individuals can be assigned to a state during sampling. Using only state-specific count data, we show how our model can be used to estimate local population abundance, as well as density-dependent recruitment rates and state-specific survival. We apply our model to a population of black-throated blue warblers ( Setophaga caerulescens ) that have been surveyed for 25 years on their breeding grounds at the Hubbard Brook Experimental Forest in New Hampshire, USA. The intensive data collection efforts allow us to compare our estimates to estimates derived from capture–recapture data. Our model performed well in estimating population abundance and density-dependent rates of annual recruitment/immigration. Estimates of local carrying capacity and per capita recruitment of yearlings were consistent with those published in other studies. However, our model moderately underestimated annual survival probability of yearling and adult females and severely underestimates survival probabilities for both of these male stages. The most accurate and precise estimates will necessarily require some amount of intensive data collection efforts (such as capture–recapture). Integrated population models that combine data from both intensive and extensive sources are likely to be the most efficient approach for estimating demographic rates at large spatial and temporal scales.

New Hampshire↗

USGS2014

A hierarchical model combining distance sampling and time removal to estimate detection probability during avian point counts

Imperfect detection during animal surveys biases estimates of abundance and can lead to improper conclusions regarding distribution and population trends. Farnsworth et al. (2005) developed a combined distance-sampling and time-removal model for point-transect surveys that addresses both availability (the probability that an animal is available for detection; e.g., that a bird sings) and perceptibility (the probability that an observer detects an animal, given that it is available for detection). We developed a hierarchical extension of the combined model that provides an integrated analysis framework for a collection of survey points at which both distance from the observer and time of initial detection are recorded. Implemented in a Bayesian framework, this extension facilitates evaluating covariates on abundance and detection probability, incorporating excess zero counts (i.e. zero-inflation), accounting for spatial autocorrelation, and estimating population density. Species-specific characteristics, such as behavioral displays and territorial dispersion, may lead to different patterns of availability and perceptibility, which may, in turn, influence the performance of such hierarchical models. Therefore, we first test our proposed model using simulated data under different scenarios of availability and perceptibility. We then illustrate its performance with empirical point-transect data for a songbird that consistently produces loud, frequent, primarily auditory signals, the Golden-crowned Sparrow ( Zonotrichia atricapilla ); and for 2 ptarmigan species ( Lagopus spp.) that produce more intermittent, subtle, and primarily visual cues. Data were collected by multiple observers along point transects across a broad landscape in southwest Alaska, so we evaluated point-level covariates on perceptibility (observer and habitat), availability (date within season and time of day), and abundance (habitat, elevation, and slope), and included a nested point-within-transect and park-level effect. Our results suggest that this model can provide insight into the detection process during avian surveys and reduce bias in estimates of relative abundance but is best applied to surveys of species with greater availability (e.g., breeding songbirds).

The Auk↗

USGS2014

Modeling structured population dynamics using data from unmarked individuals

The study of population dynamics requires unbiased, precise estimates of abundance and vital rates that account for the demographic structure inherent in all wildlife and plant populations. Traditionally, these estimates have only been available through approaches that rely on intensive mark–recapture data. We extended recently developed N -mixture models to demonstrate how demographic parameters and abundance can be estimated for structured populations using only stage-structured count data. Our modeling framework can be used to make reliable inferences on abundance as well as recruitment, immigration, stage-specific survival, and detection rates during sampling. We present a range of simulations to illustrate the data requirements, including the number of years and locations necessary for accurate and precise parameter estimates. We apply our modeling framework to a population of northern dusky salamanders ( Desmognathus fuscus ) in the mid-Atlantic region (USA) and find that the population is unexpectedly declining. Our approach represents a valuable advance in the estimation of population dynamics using multistate data from unmarked individuals and should additionally be useful in the development of integrated models that combine data from intensive (e.g., mark–recapture) and extensive (e.g., counts) data sources.

Ecology↗

USGS2014

Reply to Efford on ‘Integrating resource selection information with spatial capture-recapture’

1. We proposed ( Methods in Ecology and Evolution , 2013, 4) a model for combining telemetry data with spatial capture–recapture (SCR) data that was vigorously criticized by Efford ( Methods in Ecology and Evolution , 2014, 000, 000). Efford's main claim was that our encounter probability model was incorrect, and therefore our R code and simulation results were wrong. 2. In fact, our encounter probability model is correct under the Poisson point process model that we used as a basis for our integrated model. On the other hand, the basis for Efford's claims clearly rest on the assumption of an alternative model which, while possibly useful, is distinct from that analysed in Royle et al. (Methods in Ecology and Evolution, 2013, 4). 3. A key point of Royle et al. ( Methods in Ecology and Evolution , 2013, 4) was that active resource selection induces heterogeneity in encounter probability which, if unaccounted for, should bias estimates of population size or density. The models of Royle et al. ( Methods in Ecology and Evolution , 2013, 4) and Efford ( Methods in Ecology and Evolution , 2014, 000, 000) merely amount to alternative models of resource selection, and hence varying amounts of heterogeneity in encounter probability.

Methods in Ecology and Evolution↗

USGS2014

Northeast regional and state trends in anuran occupancy from calling survey data (2001-2011) from the North American Amphibian Monitoring Program

We present the first regional trends in anuran occupancy from North American Amphibian Monitoring Program (NAAMP) data from 11 northeastern states using an 11 years of data. NAAMP is a long-term monitoring program where observers collect data at assigned random roadside routes using a calling survey technique. We assessed occupancy trends for 17 species. Eight species had statistically significant regional trends, of these seven were negative (Anaxyrus fowleri, Acris crepitans, Pseudacris brachyphona, Pseudacris feriarum-kalmi complex, Lithobates palustris, Lithobates pipiens, and Lithobates sphenocephalus) and one was positive (Hyla versicolor-chrysoscelis complex). We also assessed state level trends for 101 species/state combinations, of these 29 showed a significant decline and nine showed a significant increase in occupancy.

Herpetological Conservation and Biology↗

USGS2013

Hierarchical spatial capture-recapture models: Modeling population density from stratified populations

Capture–recapture studies are often conducted on populations that are stratified by space, time or other factors. In this paper, we develop a Bayesian spatial capture–recapture (SCR) modelling framework for stratified populations – when sampling occurs within multiple distinct spatial and temporal strata. We describe a hierarchical model that integrates distinct models for both the spatial encounter history data from capture–recapture sampling, and also for modelling variation in density among strata. We use an implementation of data augmentation to parameterize the model in terms of a latent categorical stratum or group membership variable, which provides a convenient implementation in popular BUGS software packages. We provide an example application to an experimental study involving small-mammal sampling on multiple trapping grids over multiple years, where the main interest is in modelling a treatment effect on population density among the trapping grids. Many capture–recapture studies involve some aspect of spatial or temporal replication that requires some attention to modelling variation among groups or strata. We propose a hierarchical model that allows explicit modelling of group or strata effects. Because the model is formulated for individual encounter histories and is easily implemented in the BUGS language and other free software, it also provides a general framework for modelling individual effects, such as are present in SCR models.

Methods in Ecology and Evolution↗

USGS2013

Estimating true instead of apparent survival using spatial Cormack-Jolly-Seber models

Survival is often estimated from capture–recapture data using Cormack–Jolly–Seber (CJS) models, where mortality and emigration cannot be distinguished, and the estimated apparent survival probability is the product of the probabilities of true survival and of study area fidelity. Consequently, apparent survival is lower than true survival unless study area fidelity equals one. Underestimation of true survival from capture–recapture data is a main limitation of the method. We develop a spatial version of the CJS model that allows estimation of true survival. Besides the information about whether a specific individual was encountered at a given occasion, it is often recorded where the encounter occurred. Thus, information is available about the fraction of dispersal that occurs within the study area, and we use it to model dispersal and estimate true survival. Our model is formulated hierarchically and consists of survival, dispersal and observation submodels, assuming that encounters are possible anywhere within a study area. In a simulation study, our new spatial CJS model produced accurate estimates of true survival and dispersal behaviour for various sizes and shapes of the study area, even if emigration is substantial. However, when the information about dispersal is scarce due to low survival, low recapture probabilities and high emigration, the estimators are positively biased. Moreover, survival estimates are sensitive to the assumed dispersal kernel. We applied the spatial CJS model to a data set of adult red-backed shrikes ( Lanius collurio ). Apparent survival of males ( c . 0·5) estimated with the CJS model was larger than in females ( c . 0·4), but the application of the spatial CJS model revealed that both sexes had similar survival probabilities ( c . 0·6). The mean breeding dispersal distance in females was c . 700 m, while males dispersed only c . 250 m between years. Spatial CJS models enable study of dispersal and survival independent of study design constraints such as imperfect detection and size of the study area provided that some of the dispersing individuals remain in the study area. We discuss possible extensions of our model: alternative dispersal models and the inclusion of covariates and of a habitat suitability map.

Methods in Ecology and Evolution↗

USGS2013

A hierarchical nest survival model integrating incomplete temporally varying covariates

Nest success is a critical determinant of the dynamics of avian populations, and nest survival modeling has played a key role in advancing avian ecology and management. Beginning with the development of daily nest survival models, and proceeding through subsequent extensions, the capacity for modeling the effects of hypothesized factors on nest survival has expanded greatly. We extend nest survival models further by introducing an approach to deal with incompletely observed, temporally varying covariates using a hierarchical model. Hierarchical modeling offers a way to separate process and observational components of demographic models to obtain estimates of the parameters of primary interest, and to evaluate structural effects of ecological and management interest. We built a hierarchical model for daily nest survival to analyze nest data from reintroduced whooping cranes (Grus americana) in the Eastern Migratory Population. This reintroduction effort has been beset by poor reproduction, apparently due primarily to nest abandonment by breeding birds. We used the model to assess support for the hypothesis that nest abandonment is caused by harassment from biting insects. We obtained indices of blood-feeding insect populations based on the spatially interpolated counts of insects captured in carbon dioxide traps. However, insect trapping was not conducted daily, and so we had incomplete information on a temporally variable covariate of interest. We therefore supplemented our nest survival model with a parallel model for estimating the values of the missing insect covariates. We used Bayesian model selection to identify the best predictors of daily nest survival. Our results suggest that the black fly Simulium annulus may be negatively affecting nest survival of reintroduced whooping cranes, with decreasing nest survival as abundance of S. annulus increases. The modeling framework we have developed will be applied in the future to a larger data set to evaluate the biting-insect hypothesis and other hypotheses for nesting failure in this reintroduced population; resulting inferences will support ongoing efforts to manage this population via an adaptive management approach. Wider application of our approach offers promise for modeling the effects of other temporally varying, but imperfectly observed covariates on nest survival, including the possibility of modeling temporally varying covariates collected from incubating adults.

Wisconsin↗

USGS2013

Inference for finite-sample trajectories in dynamic multi-state site-occupancy models using hidden Markov model smoothing

Ecologists and wildlife biologists increasingly use latent variable models to study patterns of species occurrence when detection is imperfect. These models have recently been generalized to accommodate both a more expansive description of state than simple presence or absence, and Markovian dynamics in the latent state over successive sampling seasons. In this paper, we write these multi-season, multi-state models as hidden Markov models to find both maximum likelihood estimates of model parameters and finite-sample estimators of the trajectory of the latent state over time. These estimators are especially useful for characterizing population trends in species of conservation concern. We also develop parametric bootstrap procedures that allow formal inference about latent trend. We examine model behavior through simulation, and we apply the model to data from the North American Amphibian Monitoring Program.

Environmental and Ecological Statistics↗