Large-scale variation in density of an aquatic ecosystem indicator species

Evaluating density-weighted connectivity of black bears (Ursus americanus) in Glacier National Park with spatial capture-recapture models

BACKGROUND

Improved understanding of wildlife population connectivity among protected area networks can support effective planning for the persistence of wildlife populations in the face of land use and climate change. Common approaches to estimating connectivity often rely on small samples of individuals without considering the spatial structure of populations, leading to limited understanding of how individual movement links to demography and population connectivity. Recently developed spatial capture-recapture (SCR) models provide a framework to formally connect inference about individual movement, connectivity, and population density, but few studies have applied this approach to empirical data to support connectivity planning.

METHODS

We used mark-recapture data collected from 924 genetic detections of 598 American black bears (Ursus americanus) in 2004 with SCR ecological distance models to simultaneously estimate density, landscape resistance to movement, and population connectivity in Glacier National Park northwest Montana, USA. We estimated density and movement parameters separately for males and females and used model estimates to calculate predicted density-weighted connectivity surfaces.

RESULTS

Model results indicated that landscape structure influences black bear density and space use in Glacier. The mean density estimate was 16.08 bears/100 km2 (95% CI 12.52-20.6) for females and 9.27 bears/100 km2 (95% CI 7.70-11.14) for males. Density increased with forest cover for both sexes. For male black bears, density decreased at higher grizzly bear (Ursus arctos) densities. Drainages, valley bottoms, and riparian vegetation decreased estimates of landscape resistance to movement for male and female bears. For males, forest cover also decreased estimated resistance to movement, but a transportation corridor bisecting the study area strongly increased resistance to movement presenting a barrier to connectivity.

CONCLUSIONS

Density-weighed connectivity surfaces highlighted areas important for population connectivity that were distinct from areas with high potential connectivity. For black bears in Glacier and surrounding landscapes, consideration of both vegetation and valley topography could inform the placement of underpasses along the transportation corridor in areas characterized by both high population density and potential connectivity. Our study demonstrates that the SCR ecological distance model can provide biologically realistic, spatially explicit predictions to support movement connectivity planning across large landscapes.

Comparative spatially explicit approach for testing effects of soil chemicals on terrestrial wildlife bioindicator demographics

Wildlife species are often used as bioindicators to evaluate the extent and severity of environmental contamination and the effectiveness of remediation practices. A common approach for investigating population- or community-level impacts on bioindicators compares demographic parameter estimates (e.g., population size or density) between sites that were subjected to different levels of contamination. However, the traditional analytical method used in such studies is nonspatial capture-recapture, which results in conclusions about potential relationships between demographics and contaminants being inferred indirectly. Here, we extend this comparative approach to the spatially explicit framework, allowing direct estimation of said relationships and comparisons between study areas, by applying spatial capture-recapture (SCR) models to bioindicator (deer mice [Peromyscus spp.]) detection data from two study areas that were subjected to different industrial activities and remediation practices. Bioindicator density differed by 178% between the neighboring study areas, and the area with the highest soil concentrations of polychlorinated biphenyls, chromium, and zinc had the highest bioindicator density. Under the traditional nonspatial approach, we might have concluded that soil chemical levels had negligible influences on demographics. However, by modeling density as a spatial function of select chemical concentrations using SCR models, we found strong support for a positive relationship between density and soil chromium concentrations in one study area (β = 0.82), which was not masked by or associated with habitat-related metrics. To obtain reliable inferences about potential effects of environmental contamination on bioindicator demographics, we contend that a comparative spatially explicit approach using SCR ought to become standard.

Precision and bias of spatial capture-recapture estimates: A multi-site, multi-year Utah black bear case study

Spatial capture-recapture (SCR) models are powerful analytical tools that have become the standard for estimating abundance and density of wild animal populations. When sampling populations to implement SCR, the number of unique individuals detected, total recaptures, and unique spatial relocations can be highly variable. These sample sizes influence the precision and accuracy of model parameter estimates. Testing the performance of SCR models with sparse empirical data sets typical of low-density, wide-ranging species can inform the threshold at which a more integrated modeling approach with additional data sources or additional years of monitoring may be required to achieve reliable, precise parameter estimates. Using a multi-site, multi-year Utah black bear (Ursus americanus) capture-recapture data set, we evaluated factors influencing the uncertainty of SCR structural parameter estimates, specifically density, detection, and the spatial scale parameter, sigma. We also provided some of the first SCR density estimates for Utah black bear populations, which ranged from 3.85 to 74.33 bears/100 km² . Increasing total detections decreased the uncertainty of density estimates, whereas an increasing number of total recaptures and individuals with recaptures decreased the uncertainty of detection and sigma estimates, respectively. In most cases, multiple years of data were required for precise density estimates (<0.2 coefficient of variation [CV]). Across study areas there was an average decline in CV of 0.07 with the addition of another year of data. One sampled population with very high estimated bear density had an atypically low number of spatial recaptures relative to total recaptures, apparently inflating density estimates. A complementary simulation study used to assess estimate bias suggested that when <30% of recaptured individuals were spatially recaptured, density estimates were unreliable and ranged widely, in some cases to >3 times the simulated density. Additional research could evaluate these requirements for other density scenarios. Large numbers of individuals detected, numbers of spatial recaptures, and precision alone may not be sufficient indicators of parameter estimate reliability. We provide an evaluation of simple summary statistics of capture-recapture data sets that can provide an early signal of the need to alter sampling design or collect auxiliary data before model implementation to improve estimate precision and accuracy.

Nesting success of red-winged blackbirds (Agelaius phoeniceus) in marshes in an anthropogenic landscape

Recent analyses show significant population declines in many abundant avian species, especially marsh-nesting species including the red-winged blackbird (RWBL). Hypothesized causes include reduced nesting success resulting from changing land-use patterns and exposure to contaminants. Our goal was to test the hypothesis that landscape and nest characteristics as well as exposure to polychlorinated biphenyls (PCBs) correlate with nesting success. From 2008 to 2014, we measured clutch size, egg and nestling mass, hatching and fledging success and daily survival of 1293 RWBL nests from 32 marshes in the Hudson River valley of New York. Using generalized linear effect and survival models, we found that: (i) Julian date was negatively related to hatching success and clutch size but positively related to egg mass; (ii) nest height was negatively related to hatching success; (iii) nestling mass decreased with increased nest density and distance to edges; (iv) fledging success was significantly lower in nests closer to the ground that were far from water; and (v) clutch size and daily survival were higher in nests farther from water. Results showed that nesting success was correlated with variables associated with flooding, population density and predation and provided no support for the predicted negative effects of PCB exposure.

Explaining detection heterogeneity with finite mixture and non-Euclidean movement in spatially explicit capture-recapture models

Landscape structure affects animal movement. Differences between landscapes may induce heterogeneity in home range size and movement rates among individuals within a population. These types of heterogeneity can cause bias when estimating population size or density and are seldom considered during analyses. Individual heterogeneity, attributable to unknown or unobserved covariates, is often modelled using latent mixture distributions, but these are demanding of data, and abundance estimates are sensitive to the parameters of the mixture distribution. A recent extension of spatially explicit capture-recapture models allows landscape structure to be modelled explicitly by incorporating landscape connectivity using non-Euclidean least-cost paths, improving inference, especially in highly structured (riparian & mountainous) landscapes. Our objective was to investigate whether these novel models could improve inference about black bear (Ursus americanus) density. We fit spatially explicit capture-recapture models with standard and complex structures to black bear data from 51 separate study areas. We found that non-Euclidean models were supported in over half of our study areas. Associated density estimates were higher and less precise than those from simple models and only slightly more precise than those from finite mixture models. Estimates were sensitive to the scale (pixel resolution) at which least-cost paths were calculated, but there was no consistent pattern across covariates or resolutions. Our results indicate that negative bias associated with ignoring heterogeneity is potentially severe. However, the most popular method for dealing with this heterogeneity (finite mixtures) yielded potentially unreliable point estimates of abundance that may not be comparable across surveys, even in data sets with 136-350 total detections, 3-5 detections per individual, 97-283 recaptures, and 80-254 spatial recaptures. In these same study areas with high sample sizes, we expected that landscape features would not severely constrain animal movements and modelling non-Euclidian distance would not consistently improve inference. Our results suggest caution in applying non-Euclidean SCR models when there is no clear landscape covariate that is known to strongly influence the movement of the focal species, and in applying finite mixture models except when abundant data are available.

Evaluating otter reintroduction outcomes using genetic spatial capture-recapture modified for dendritic networks

Monitoring the demographics and genetics of reintroduced populations is critical to evaluating reintroduction success, but species ecology and the landscapes that they inhabit often present challenges for accurate assessments. If suitable habitats are restricted to hierarchical dendritic networks, such as river systems, animal movements are typically constrained and may violate assumptions of methods commonly used to estimate demographic parameters. Using genetic detection data collected via fecal sampling at latrines, we demonstrate applicability of the spatial capture-recapture (SCR) network distance function for estimating the size and density of a recently reintroduced North American river otter (Lontra canadensis) population in the Upper Rio Grande River dendritic network in the southwestern United States, and we also evaluated the genetic outcomes of using a small founder group (n = 33 otters) for reintroduction. Estimated population density was 0.23-0.28 otter/km, or 1 otter/3.57-4.35 km, with weak evidence of density increasing with northerly latitude (β = 0.33). Estimated population size was 83-104 total otters in 359 km of riverine dendritic network, which corresponded to average annual exponential population growth of 1.12-1.15/year since reintroduction. Growth was ≥40% lower than most reintroduced river otter populations and strong evidence of a founder effect existed 8-10 years post-reintroduction, including 13-21% genetic diversity loss, 84%-87% genetic effective population size decline, and rapid divergence from the source population (F ST accumulation = 0.06/generation). Consequently, genetic restoration via translocation of additional otters from other populations may be necessary to mitigate deleterious genetic effects in this small, isolated population. Combined with non-invasive genetic sampling, the SCR network distance approach is likely widely applicable to demogenetic assessments of both reintroduced and established populations of multiple mustelid species that inhabit aquatic dendritic networks, many of which are regionally or globally imperiled and may warrant reintroduction or augmentation efforts.

Optimal sampling design for spatial capture-recapture

Spatial capture-recapture (SCR) has emerged as the industry standard for estimating population density by leveraging information from spatial locations of repeat encounters of individuals. The precision of density estimates depends fundamentally on the number and spatial configuration of traps. Despite this knowledge, existing sampling design recommendations are heuristic and their performance remains untested for most practical applications. To address this issue, we propose a genetic algorithm that minimizes any sensible, criteria-based objective function to produce near-optimal sampling designs. To motivate the idea of optimality, we compare the performance of designs optimized using three model-based criteria related to the probability of capture. We use simulation to show that these designs outperform those based on existing recommendations in terms of bias, precision, and accuracy in the estimation of population size. Our approach, available as a function in the R package oSCR, allows conservation practitioners and researchers to generate customized and improved sampling designs for wildlife monitoring.

A noninvasive environmental monitoring tool for brominated flame-retardants (BFRs) assisted by conservation detection dogs

Fecal matter is a useful noninvasive/nondestructive media for evaluating contaminants in wildlife, as residues therein have been observed to correlate with body burdens. Conservation detection dog-handler teams can be used to optimize the acquisition of fecal samples. To build on previous work, sentinel-species' (i.e. mink (Mustela vison) and otter (Lontra canadensis)) fecal matter was opportunistically located by a detection dog team along the tri-river system of Missoula, Montana, USA. Sediments were also collected. Samples were used to develop an analytical method from fecal matter to determine habitat exposure to the brominated flame-retardants (BFRs): polybrominated diphenyl ethers (PBDEs), hexabromocyclododecane (HBCDD), 2-ethylhexyl 2, 3, 4, 5-tetrabromobenzoate (EH-TBB), di (2-ethylhexyl)-2, 3, 4, 5-tetrabromophthalate (BEH-TEBP) and decabromodiphenyl ethane (DBDPE). Sediments contained PBDEs (BDE-99 and BDE-209) and EH-TBB at detection rates of 67%, 33% and 67%, respectively. BDE-99, -209 and EH-TBB were also detected in mink and otter feces, at rates of 81%, 25% and 81%, respectively; plus BEH-TEBP at 13%. BFR levels correlated positively with human population density except along the lower Bitterroot River, where BDE-209 sediment and feces levels exceeded other sites by several orders of magnitude. Fecal matter body burden estimates indicated marginal PBDE exposure. However, exposure to their replacements, EH-TBB and BEH-TEBP, were at levels that may adversely affect healthy Mustelidae populations. Proof-of-concept was achieved; validation results were within established standards for the development of analytical methods. The established application of conservation dog-handler teams to facilitate the collection of fecal matter for BFR analysis represents a valuable, but currently underutilized environmental monitoring tool.

Investigating effects of soil chemicals on density of small mammal bioindicators using spatial capture-recapture models

Monitoring the ecological impacts of environmental pollution and the effectiveness of remediation efforts requires identifying relationships between contaminants and the disruption of biological processes in populations, communities, or ecosystems. Wildlife are useful bioindicators, but traditional comparative experimental approaches rely on a staunch and typically unverifiable assumption that, in the absence of contaminants, reference and contaminated sites would support the same densities of bioindicators, thereby inferring direct causation from indirect data. We demonstrate the utility of spatial capture-recapture (SCR) models for overcoming these issues, testing if community density of common small mammal bioindicators was directly influenced by soil chemical concentrations. By modeling density as an inhomogeneous Poisson point process, we found evidence for an inverse spatial relationship between Peromyscus density and soil mercury concentrations, but not other chemicals, such as polychlorinated biphenyls, at a site formerly occupied by a nuclear reactor. Although the coefficient point estimate supported Peromyscus density being lower where mercury concentrations were higher (β = –0.44), the 95% confidence interval overlapped zero, suggesting no effect was also compatible with our data. Estimated density from the most parsimonious model (2.88 mice/ha; 95% CI = 1.63–5.08), which did not support a density-chemical relationship, was within the range of reported densities for Peromyscus that did not inhabit contaminated sites elsewhere. Environmental pollution remains a global threat to biodiversity and ecosystem and human health, and our study provides an illustrative example of the utility of SCR models for investigating the effects that chemicals may have on wildlife bioindicator populations and communities.

Movement-assisted localization from acoustic telemetry data

BACKGROUND

Acoustic telemetry technologies are being increasingly deployed to study a variety of aquatic taxa including fishes, reptiles, and marine mammals. Large cooperative telemetry networks produce vast quantities of data useful in the study of movement, resource selection and species distribution. Efficient use of acoustic telemetry data requires estimation of acoustic source locations from detections at receivers (i.e., "localization"). Multiple processes provide information for localization estimation including detection/non-detection data at receivers, information on signal rate, and an underlying movement model describing how individuals move and utilize space. Frequently, however, localization methods only integrate a subset of these processes and do not utilize the full spatial encounter history information available from receiver arrays.

METHODS

In this paper we draw analogies between the challenges of acoustic telemetry localization and newly developed methods of spatial capture-recapture (SCR). We develop a framework for localization that integrates explicit sub-models for movement, signal (or cue) rate, and detection probability, based on acoustic telemetry spatial encounter history data. This method, which we call movement-assisted localization, makes efficient use of the full encounter history data available from acoustic receiver arrays, provides localizations with fewer than three detections, and even allows for predictions to be made of the position of an individual when it was not detected at all. We demonstrate these concepts by developing generalizable Bayesian formulations of the SCR movement-assisted localization model to address study-specific challenges common in acoustic telemetry studies.

RESULTS

Simulation studies show that movement-assisted localization models improve point-wise RMSE of localization estimates by >50% and greatly increased the precision of estimated trajectories compared to localization using only the detection history of a given signal. Additionally, integrating a signal rate sub-model reduced biases in the estimation of movement, signal rate, and detection parameters observed in independent localization models.

CONCLUSIONS

Movement-assisted localization provides a flexible framework to maximize the use of acoustic telemetry data. Conceptualizing localization within an SCR framework allows extensions to a variety of data collection protocols, improves the efficiency of studies interested in movement, resource selection, and space-use, and provides a unifying framework for modeling acoustic data.

Metal contamination of river otters in North Carolina

Aquatic apex predators are vulnerable to environmental contaminants due to biomagnification. North American river otter (Lontra canadensis) populations should be closely monitored across their range due to point and nonpoint pollution sources. Nonetheless, no information exists on environmental contaminants in the North Carolina otter population. Metals and metalloids occur naturally across the landscape, are essential for cellular function, and become toxic when concentrated unnaturally. We conducted our study across the three Furbearer Management Units (FMU) and 14 river basins of North Carolina. We determined the concentrations of arsenic, cadmium, calcium, cobalt, copper, iron, lead, magnesium, manganese, mercury, molybdenum, selenium, thallium, and zinc in liver and kidney samples from 317 otters harvested from 2009 to 2016. Arsenic, lead, and thallium samples were tested at levels below the limit of detection. With the exception of cadmium, we detected all other elements at higher levels in the liver compared with the kidney. Specifically, cadmium, cobalt, copper, iron, magnesium, manganese, mercury, molybdenum, and zinc levels differed by tissue type analyzed. Most element concentrations remained stable or increased with otter age. We detected higher levels of mercury and selenium in the Lower Pee Dee and Cape Fear river basins. River basins within the Mountain FMU were higher in cadmium, copper, iron, lead, and zinc, whereas the Coastal Plain FMU was lower in cobalt and manganese. None of the elements occurred at toxic levels. Our research establishes baseline concentration levels for North Carolina, which will benefit future monitoring efforts and provide insight into future changes in the otter population.

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 has become central to understanding the dynamics of coupled human-natural systems, resource-dependent livelihoods, and biodiversity conservation. Ecologists are facing growing pressure to quantify the size, distribution, and trajectory of wild populations in a cost-effective and socially acceptable manner. Genetic tagging, combined with modern computational and genetic analyses, is an under-utilized tool to meet this demand, especially for wide-ranging, elusive, sensitive, and low-density species. Genetic tagging studies are now revealing unprecedented insight into the mechanisms that control the density, trajectory, connectivity, and patterns of human-wildlife interaction for populations over vast spatial extents. Here, we outline the application of, and ecological inferences from, new analytical techniques applied to genetically tagged individuals, contrast this approach with conventional methods, and describe how genetic tagging can be better applied to address outstanding questions in ecology. We provide example analyses using a long-term genetic tagging dataset of grizzly bears in the Canadian Rockies. The genetic tagging toolbox is a powerful and overlooked ensemble that ecologists and conservation biologists can leverage to generate evidence and meet the challenges of the Anthropocene.