Geology Reports⌕ Search

Geology topics

J. Andrew Royle

Publications and source records attributed to J. Andrew Royle.

At least 55 records · Page 3Linked to original sources
USGS2020

Migratory behavior and winter geography drive differential range shifts of eastern birds in response to recent climate change

Over the past half century, migratory birds in North America have shown divergent population trends relative to resident species, with the former declining rapidly and the latter increasing. The role that climate change has played in these observed trends is not well understood, despite significant warming over this period. We used 43 y of monitoring data to fit dynamic species distribution models and quantify the rate of latitudinal range shifts in 32 species of birds native to eastern North America. Since the early 1970s, species that remain in North America throughout the year, including both resident and migratory species, appear to have responded to climate change through both colonization of suitable area at the northern leading edge of their breeding distributions and adaption in place at the southern trailing edges. Neotropical migrants, in contrast, have shown the opposite pattern: contraction at their southern trailing edges and no measurable shifts in their northern leading edges. As a result, the latitudinal distributions of temperate-wintering species have increased while the latitudinal distributions of neotropical migrants have decreased. These results raise important questions about the mechanisms that determine range boundaries of neotropical migrants and suggest that these species may be particularly vulnerable to future climate change. Our results highlight the potential importance of climate change during the nonbreeding season in constraining the response of migratory species to temperature changes at both the trailing and leading edges of their breeding distributions. Future research on the interactions between breeding and nonbreeding climate change is urgently needed.

Proceedings of the National Academy of Sciences of↗
USGS2020

Acoustic space occupancy: Combining ecoacoustics and lidar to model biodiversity variation and detection bias across heterogeneous landscapes

There is global interest in quantifying changing biodiversity in human-modified landscapes. Ecoacoustics may offer a promising pathway for supporting multi-taxa monitoring, but its scalability has been hampered by the sonic complexity of biodiverse ecosystems and the imperfect detectability of animal-generated sounds. The acoustic signature of a habitat, or soundscape, contains information about multiple taxa and may circumvent species identification, but robust statistical technology for characterizing community-level attributes is lacking. Here, we present the Acoustic Space Occupancy Model, a flexible hierarchical framework designed to account for detection artifacts from acoustic surveys in order to model biologically relevant variation in acoustic space use among community assemblages. We illustrate its utility in a biologically and structurally diverse Amazon frontier forest landscape, a valuable test case for modeling biodiversity variation and acoustic attenuation from vegetation density. We use complementary airborne lidar data to capture aspects of 3D forest structure hypothesized to influence community composition and acoustic signal detection. Our novel analytic framework permitted us to model both the assembly and detectability of soundscapes using lidar-derived estimates of forest structure. Our empirical predictions were consistent with physical models of frequency-dependent attenuation, and we estimated that the probability of observing animal activity in the frequency channel most vulnerable to acoustic attenuation varied by over 60%, depending on vegetation density. There were also large differences in the biotic use of acoustic space predicted for intact and degraded forest habitats, with notable differences in the soundscape channels predominantly occupied by insects . This study advances the utility of ecoacoustics by providing a robust modeling framework for addressing detection bias from remote audio surveys while preserving the rich dimensionality of soundscape data, which may be critical for inferring biological patterns pertinent to multiple taxonomic groups in the tropics. Our methodology paves the way for greater integration of remotely sensed observations with high-throughput biodiversity data to help bring routine, multi-taxa monitoring to scale in dynamic and diverse landscapes.

Mato Grosso↗
USGS2020

Integrating side-scan sonar and acoustic telemetry to estimate the annual spawning run size of Atlantic sturgeon in the Hudson River

There is considerable interest in evaluating the status and trends of sturgeon populations, yet many traditional approaches to estimating the abundance of fishes are intractable due to their biology and rarity. Side-scan sonar has recently emerged as an effective tool for censusing sturgeon in rivers, yet challenges remain for censusing open populations that may visit specific habitats over periods of time (e.g., spawning runs). We use a hierarchical model to integrate side-scan sonar with acoustic telemetry, to estimate the proportion of a spawning run fitted with acoustic tags (12%; 95% CRI = 8-16%) and extrapolate to the total run size in 2014. Our investigation represents a novel approach to generating run size estimates in a large river and provides the first estimate of Atlantic sturgeon spawning run size for the Hudson River (N ̂ = 466; 95% CRI = 310-745) since the fishery moratorium in the 1990’s. Our estimate suggests that the Hudson River holds one of the largest contemporary populations of Atlantic sturgeon, but also indicates that it remains sharply depleted relative to virgin conditions.

Hudson River↗
USGS2020

Consequences of ignoring group association in spatial capture-recapture analysis

Many models in population ecology, including spatial capture–recapture (SCR) models, assume that individuals are distributed and detected independently of one another. In reality, this is rarely the case – both antagonistic and gregarious relationships lead to non-independent spatial configurations, with territorial exclusion at one end of the spectrum and group-living at the other. Previous simulation studies suggest that grouping has limited impact on the outcome of SCR analyses. However, group associations entail not only spatial clustering of activity centers but also coordinated space use by group members, potentially impacting both ecological and observation processes underlying SCR analysis. We simulated SCR scenarios with different strengths of aggregation (clustering of individuals into groups with shared activity centers) and cohesion (synchronization of detection patterns of members of a group). We then fit SCR models to the simulated data sets and evaluated the effect of aggregation and cohesion on parameter estimates. Low to moderate aggregation and cohesion did not impact the bias and precision of estimates of density and the scale parameter of the detection function. However, non-independence between individuals led to high levels of overdispersion. Overdispersion strongly decreased the coverage of confidence intervals around parameter estimates, thereby increasing the probability of erroneous predictions. Our results indicate that SCR models are robust to moderate levels of aggregation and cohesion. Nonetheless, spatial dependence between individuals can lead to false inference. We recommend that practitioners 1) test for the presence of overdispersion in SCR data caused by aggregation and cohesion, and, if necessary, 2) correct their variance estimates using the overdispersion factor ĉ . Approaches for doing both are described in this paper. We also urge the development of SCR models that incorporate spatial associations between individuals not only to account for overdispersion but also to obtain quantitative information about social aspects of study populations.

Wildlife Biology↗
USGS2020

Modeling population dynamics with count data

In this chapter, we describe models of open populations that are subject to change over time due to additions and subtractions. Additions may be in the form of recruitment and immigration, and subtractions may be in the form of mortality, emigration, or both. Conceptually, these models are described by the Birth-Immigration-Death-Emigration (BIDE) model of population dynamics (Conroy and Carroll, 2009). In most cases, we will not formally distinguish between the two types of additions or of subtractions (birth/immigration or death/emigration), although sometimes this may be possible depending on the timescale of the study, spatial structure, and specific model assumptions (Zhao et al., 2017; see Section 2.10). In addition, distinguishing the different dynamic processes may also become possible in the presence of auxiliary data on some demographic rates, in the context of integrated population models (IPMs, Besbeas et al., 2002; see also Chapter 10). One type of open model, which allows for temporal variation in abundance but not explicit dynamics, is the simple model of temporary emigration (Kendall et al., 1997), which supposes that population size Nt changes randomly among open (primary) periods t, as a Binomial realization from some larger superpopulation. Over short timescales, this simple model may provide a sensible description of variation in population size over time.

Book chapter↗
USGS2020

Modeling false positives

Many of the models we are concerned with included explicit descriptions of false negative errors. However, false positive errors can also be commin in practice, especially in citizen science applications where observer skill is highly variable. In addition, new methods which determine detection based on statistical classification or machine learning methods are also prone to false positive errors which must be accounted for. An early treatment of the false positive detection problem by Royle & Link (2006) recognized that false positive errors can be accommodated by a mixture model for detection probability: one value of detection at occupied sites and another non-zero value at unoccupied sites. This model has been extended greatly in recent years to include more informative data about false positives including validation or confirmation data (Miller et al. 2011) and multiple detection methods, among others. A new frontier for the application of false positives models lies in the use of modern technologies such as bioacoustics for efficient automated monitoring. For these technologies to realize their promise there must be improvements in automated processing of the vast quantities of output produced. Statistical classification methods (machine learning) are fallible and necessarily produce false positive detections. Therefore models which account for this process are necessary (Chambert et al. 2017). It stands to reason that false positives will need to be accounted for in other new technologies that rely on automated digital processing, including eDNA, genetic barcoding, and automated detection in remote camera studies. We devise a new occupancy model that integrates data from bioacoustics sampling with an occupancy model. This integrated model allows occupancy probability to inform species classification of samples and vice versa bioacoustics detection data inform occupancy. We provide a proof of concept for this new model in this chapter. As the core hierarchical model for the false positives models covered in this chapter are just ordinary occupancy models, extension of the ideas to open systems poses no technical challenges. We provide a suite of illustrations of these extensions. Perhaps the most prominent mechanism that leads to false positive errors it he mis-classification of species detections, or the confusion of one species for another. Very little work has been done on developing models based on this mechanistic understanding although Chambert et al. (2018) develop this idea as a 2-species occupancy model with error. We believe one important area of future research is to extend these ideas to truly multi-species systems.

Book chapter↗
USGS2019

Occupancy Patterns of Breeding American Black Ducks

Occupancy patterns can assist with the determination of habitat limitation during breeding or wintering periods and can help guide population and habitat management efforts. American black ducks ( Anas rubripes ; black ducks) are thought to be limited by habitat and food availability during the winter, but breeding sites may also limit the size or growth potential of the population. The Canadian Wildlife Service conducts an annual breeding waterfowl survey that we used to explore the hypothesis that black duck carrying capacity is limited by wetlands available for breeding in Québec, Canada. We applied single‐visit, multi‐species occupancy models to the 1990–2015 population survey data to determine if there was evidence the black duck population was limited by breeding habitat. Using a dynamic (multi‐season) occupancy modeling approach, we estimated latent occupancy (occupancy accounting for imperfect detection) of black ducks and then used latent occupancy estimates to derive occupancy, colonization, and extirpation rates. We jointly modeled the occupancy dynamics of black ducks and other duck species in wetlands where both species were present. Throughout the duration of the survey, 44% of wetlands were never observed to be occupied by black ducks. Occupancy models showed wetland size was positively associated with occupancy at the first time step (initial occupancy) and colonization. All 2‐species models indicated initial black duck occupancy, persistence (continued occupancy), and colonization were positively associated with the presence of a second species. Colonization rate over the 26‐year period ranged from 7% to 27% across all models. Extirpation rates were similar and were constant through time within each model. Low occupancy rates, combined with approximately equal colonization and extirpation rates, suggest there are available wetlands for breeding black ducks in their core breeding area. If breeding habitats are not saturated, this suggests migration or wintering areas may be more limiting to black duck population abundance.

Ontario, Quebec↗
USGS2019

Modeling spatially and temporally complex range dynamics when detection is imperfect

Species distributions are determined by the interaction of multiple biotic and abiotic factors, which produces complex spatial and temporal patterns of occurrence. As habitats and climate change due to anthropogenic activities, there is a need to develop species distribution models that can quantify these complex range dynamics. In this paper, we develop a dynamic occupancy model that uses a spatial generalized additive model to estimate non-linear spatial variation in occupancy not accounted for by environmental covariates. The model is flexible and can accommodate data from a range of sampling designs that provide information about both occupancy and detection probability. Output from the model can be used to create distribution maps and to estimate indices of temporal range dynamics. We demonstrate the utility of this approach by modeling long-term range dynamics of 10 eastern North American birds using data from the North American Breeding Bird Survey. We anticipate this framework will be particularly useful for modeling species’ distributions over large spatial scales and for quantifying range dynamics over long temporal scales.

Scientific Reports↗
USGS2019

Incorporating citizen science data in spatially explicit integrated population models

Information about population abundance, distribution, and demographic rates is critical for understanding a species’ ecology and for effective conservation and management. To collect data over large spatial and temporal extents for such inferences, especially for species with low densities or wide distributions, citizen science can be an efficient approach. Integrated models have also emerged as an important methodology to estimate population parameters by combining multiple types of data, including citizen science data. We developed a spatially explicit integrated model that combines opportunistically collected presence–absence (PA) data, commonly collected in citizen science efforts, with systematically collected spatial capture–recapture (SCR) data, which are often limited to small spatial and temporal extents. We conducted single and multi‐season simulations with parameters informed by North American black bear ( Ursus americanus ) populations, to evaluate the influence of varying amounts of opportunistic PA data collected at larger spatial and temporal extents on the estimation of population‐level parameters. Integrating opportunistic PA data increased the precision and accuracy of posterior estimates of abundance, and survival and recruitment rates. In some cases, adding PA locations improved abundance estimates more than increasing PA detection probability. Posterior estimates were as precise and unbiased as when higher quality, but sparse, SCR data were available. We also applied the integrated model to SCR and citizen science PA data collected on black bears in New York, with results consistent with our simulations. Our findings indicate that citizen science in integrated models can be a cost‐efficient way to improve estimates of population parameters and increase the spatiotemporal extent of inference. Continued developments with integrated models and citizen science data will offer additional ways to improve our understanding of population structure and demographics.

Ecology↗
USGS2019

Reserve design to optimize functional connectivity and animal density

Ecological distance-based spatial capture–recapture models (SCR) are a promising approach for simultaneously estimating animal density and connectivity, both of which affect spatial population processes and ultimately species persistence. We explored how SCR models can be integrated into reserve-design frameworks that explicitly acknowledge both the spatial distribution of individuals and their space use resulting from landscape structure. We formulated the design of wildlife reserves as a budget-constrained optimization problem and conducted a simulation to explore 3 different SCR-informed optimization objectives that prioritized different conservation goals by maximizing the number of protected individuals, reserve connectivity, and density-weighted connectivity. We also studied the effect on our 3 objectives of enforcing that the space-use requirements of individuals be met by the reserve for individuals to be considered conserved (referred to as home-range constraints). Maximizing local population density resulted in fragmented reserves that would likely not aid long-term population persistence, and maximizing the connectivity objective yielded reserves that protected the fewest individuals. However, maximizing density-weighted connectivity or preemptively imposing home-range constraints on reserve design yielded reserves of largely spatially compact sets of parcels covering high-density areas in the landscape with high functional connectivity between them. Our results quantify the extent to which reserve design is constrained by individual home-range requirements and highlight that accounting for individual space use in the objective and constraints can help in the design of reserves that balance abundance and connectivity in a biologically relevant manner.

Conservation Biology↗
USGS2019

oSCR: A spatial capture–recapture R package for inference about spatial ecological processes

Spatial capture–recapture (SCR) methods have become widely applied in ecology. The immediate adoption of SCR is due to the fact that it resolves some major criticisms of traditional capture–recapture methods related to heterogeneity in detectabililty, and the emergence of new technologies (e.g. camera traps, non‐invasive genetics) that have vastly improved our ability to collection spatially explicit observation data on individuals. However, the utility of SCR methods reaches far beyond simply convenience and data availability. SCR presents a formal statistical framework that can be used to test explicit hypotheses about core elements of population and landscape ecology, and has profound implications for how we study animal populations. In this software note, we describe the technical basis and analytical workflow of oSCR, an R package for analyzing spatial encounter history data using a multi‐session sex‐structured likelihood. The impetus for developing oSCR was to create an accessible and transparent analysis tool that allows users to conveniently and intuitively formulate statistical models that map directly to fundamental processes of interest in spatial population ecology (e.g. space use, resource selection, density and connectivity). We have placed an emphasis on creating a transparent and accessible code base that is coupled with a logical workflow that we hope stimulates active participation in further technical developments.

Ecography↗
USGS2019

Integrated modeling reveals shifts in waterfowl population dynamics under climate change

1. Climate change has been identified as one of the most important drivers of wildlife populations. The development of appropriate conservation strategies relies on reliable predictions of population responses to climate change, which require in-depth understanding of the complex relationships between climate and population dynamics through density dependent demographic processes. Integrated population models (IPMs) are a type of modeling approach that unify the analyses of demography and abundance data, providing opportunities to understand and predict population demography and dynamics under climate change. 2. In this study we developed dynamic N-mixture models for large scale population estimates, which became an important component of the IPM we used in data analysis. We then analyzed four decades (1974-2014) of Mallard ( Anas platyrhynchos ) breeding population survey, band-recovery, and climate data covering a large spatial extent from North American prairies through boreal habitat to Alaska. Our goals were to examine the complex relationships among climate, density dependent processes, waterfowl population demography and dynamics, identify the key demographic parameters that are sensitive to climate change and are influential to population growth, and forecast population responses to climate change. 3. Our results revealed the interactive effects of temperature and density dependent processes on Mallard recruitment and to a less extent apparent survival. We also found that recruitment explained more variance of population growth than apparent survival. We then forecasted a decrease in Mallard breeding population density in the Northern Prairie Potholes and an increase in Mallard breeding population density in the northern part of our study area, indicating potential shifts in Mallard population dynamics under future climate change. 4. Synthesis and applications Different strategies need to be considered across regions to conserve waterfowl populations under climate change. Strategies that facilitate recruitment are essential for high-density populations that are relatively vulnerable to climate change. By contrast, low-density populations are relatively resilient to climate change and their habitats may serve as future climate refugia. Adaptive management is essential for evaluating management consequences. Our modelling framework approach can be easily adapted for other species and thus has wide applications in ecology and conservation.

Ecography↗
USGS2019

Genetic tagging in the Anthropocene: Scaling ecology from alleles to ecosystems

The Anthropocene is an era of marked human impact on the world. Quantifying these impacts 51 has become central to understanding the dynamics of coupled human-natural systems, resource52 dependent livelihoods, and biodiversity conservation. Ecologists are facing growing pressure to 53 quantify the size, distribution, and trajectory of wild populations in a cost-effective and socially54 acceptable manner. Genetic tagging, combined with modern computational and genetic analyses, 55 is an under-utilized tool to meet this demand, especially for wide-ranging, elusive, sensitive, and 56 low-density species. Genetic tagging studies are now revealing unprecedented insight into the 57 mechanisms that control the density, trajectory, connectivity and human-wildlife conflict for 58 populations over vast spatial scales. Here we outline the application of, and ecological inferences 59 from, new analytical techniques applied to genetically-tagged individuals, contrast this approach 60 with conventional methods, and describe how genetic tagging can be better applied to address 61 outstanding questions in ecology. We provide example analyses using a long-term genetic 62 tagging dataset of grizzly bears in the Canadian Rockies. The genetic tagging toolbox is a 63 powerful and overlooked ensemble that ecologists and conservation biologists can leverage to 64 generate evidence and meet the challenges of the Anthropocene.

Ecological Applications↗
USGS2019

Spatial capture–recapture for categorically marked populations with an application to genetic capture–recapture

Recently introduced unmarked spatial capture–recapture (SCR), spatial mark–resight (SMR), and 2‐flank spatial partial identity models (SPIMs) extend the domain of SCR to populations or observation systems that do not always allow for individual identity to be determined with certainty. For example, some species do not have natural marks that can reliably produce individual identities from photographs, and some methods of observation produce partial identity samples as is the case with remote cameras that sometimes produce single‐flank photographs. Unmarked SCR, SMR, and SPIM share the feature that they probabilistically resolve the uncertainty in individual identity using the spatial location where samples were collected. Spatial location is informative of individual identity in spatially structured populations because a sample is more likely to have been produced by an individual living near the trap where it was recorded than an individual living further away from the trap. Further, the level of information about individual identity that a spatial location contains is related to two key ecological concepts, population density and home range size, which we quantify using a proposed Identity Diversity Index (IDI). We show that latent and partial identity SCR models produce imprecise and biased density estimates in many high IDI scenarios when data are sparse. We then extend the unmarked SCR model to incorporate categorical, partially identifying covariates, which reduce the level of uncertainty in individual identity, increasing the reliability and precision of density estimates, and allowing reliable density estimation in scenarios with higher IDI values and with more sparse data. We illustrate the performance of this “categorical SPIM” via simulations and by applying it to a black bear data set using microsatellite loci as categorical covariates, where we reproduce the full data set estimates with only slightly less precision using fewer loci than necessary for confident individual identification. We then discuss how the categorical SPIM can be applied to other wildlife sampling scenarios such as remote camera surveys, where natural or researcher‐applied partial marks can be observed in photographs. Finally, we discuss how the categorical SPIM can be added to SMR, 2‐flank SPIM, or other latent identity SCR models.

Ecosphere↗
USGS2019

Dynamic N-mixture models with temporal variability in detection probability

In theory parameters of dynamic N-mixture models can be estimated with multiple years of data without the robust design under the assumption of constant detection probability. However, such an assumption can rarely be met in long-term studies, and the consequences of violating this assumption in the inferences of dynamic N-mixture models have not been assessed. In this study we used simulation studies to evaluate inferences of the original dynamic N-mixture model and two of its spatial extensions in the face of temporal variability in detection probability. We first evaluated the dynamic N-mixture models when detection probability that varied temporally was wrongly treated as a constant. We then evaluated if the robust design was necessary for dynamic N-mixture models to provide valid parameter estimates when detection probability was correctly assumed to vary temporally. Our results showed that, when detection probability that varied temporally was wrongly treated as a constant, biases were introduced in the parameter estimates of dynamic N-mixture models. When detection probability was correctly assumed to vary temporally, the models could provide valid parameter estimates with the robust design. The model could also provide valid parameter estimates when detection probability was a random effect, even without the robust design. Based on our results, we strongly recommended considering temporal variability in detection probability when using dynamic N-mixture models to analyze long-term data and adopting the robust design in long-term surveys. Our work here is not only useful for data analysis but also important for research design, and thus are relevant to a wide range of studies.

Ecological Modelling↗
USGS2018

Observer-free experimental evaluation of habitat and distance effects on the detection of anuran and bird vocalizations

Acoustic surveys of vocalizing animals are conducted to determine density, distribution, and diversity. Acoustic surveys are traditionally performed by human listeners, but automated recording devices (ARD) are becoming increasingly popular. Signal strength decays, or attenuates, with increasing distance between source and receiver and some habitat types may differentially increase attenuation beyond the effects of distance alone. These combined effects are rarely accounted for in acoustic monitoring programs. We evaluated the performance of three playback devices and three ARD models using the calls of six anurans, six birds, and four pure tones. Based on these evaluations, we determined the optimal playback and recording devices. Using these optimal devices, we broadcast and recorded vocalizations in five habitat types along 1,000 m transects. We used generalized linear models to test for effects of habitat, distance, species, environmental, and landscape variables. We predicted detection probabilities for each vocalization, in each habitat type, from 0 to 1,000 m. Among playback devices, only a remote predator caller simulated vocalizations consistently. Differences of ~10 dB were observed among ARDs. For all species, we found differences in detectability between open and closed canopy habitats. We observed large differences in predicted detection probability among species in each habitat type, as well as along 1,000 m transects. Increases in temperature, barometric pressure, and wind speed significantly decreased detection probability. However, aside from differences among species, habitat, and distance, topography impeding a line‐of‐sight between sound source and receiver had the greatest negative influence on detections. Our results suggest researchers should model the effects of habitat, distance, and frequency on detection probability when performing acoustic surveys. To optimize survey design, we recommend pilot measurements among varying habitats.

Ecology and Evolution↗