Joint estimation of habitat dynamics and species interactions: disturbance reduces co-occurrence of non-native predators with an endangered toad

Native amphibian toxin reduces invasive crayfish feeding with potential benefits to stream biodiversity

BACKGROUND

Biodiversity is generally reduced when non-native species invade an ecosystem. Invasive crayfish, Procambarus clarkii, populate California freshwater streams, and in the Santa Monica Mountains (Los Angeles, USA), their introduction has led to trophic cascades due to omnivorous feeding behavior and a rapid rate of population growth. The native California newt, Taricha torosa, possesses a neurotoxin, tetrodotoxin (TTX), that affects freshwater animal behavior. Given P. clarkii has a limited evolutionary history with TTX, we hypothesized that TTX may affect crayfish feeding behaviors. To determine if TTX affects P. clarkii behavior, we measured cumulative movement and various feeding behaviors of P. clarkii exposed to (i) waterborne, ecologically realistic concentrations of TTX (~ 3.0 × 10- 8 moles/L), (ii) an anuran chemical cue to account for intraguild cues, or (iii) a T. torosa chemical cue with quantitated TTX in it (~ 6.2 × 10- 8 moles/L).

RESULTS

We found that the presence of TTX in any form significantly reduced crayfish movement and decreased the amount of food consumed over time. Crayfish responses to the anuran treatment did not significantly differ from controls.

CONCLUSION

Our laboratory results show that naturally occurring neurotoxin from native California newts limits invasive crayfish foraging and feeding rates, which may play a role in preserving local stream ecosystems by limiting invasive crayfish behaviors that are detrimental to biodiversity.

Research on Space Occupancy, Activity Rhythm and Sexual Segregation of White-Lipped Deer (Cervus albirostris) in Forest Habitats of Jiacha Gorge on Yarlung Zangbo River Basin Based on Infrared Camera Technology

The white-lipped deer (Cervus albirostris) is a rare and endangered species found in the Qinghai-Tibet Plateau in China. To understand the space occupancy, activity rhythm, and sexual segregation of the white-lipped deer, 24,096 effective photos and 827 effective videos were captured using infrared cameras from February 2020 to January 2022. The ecology and behavior of the white-lipped deer in Jiacha Gorge were studied in more detail using site occupancy models, relative abundance index, and other technologies and methods. The results show that The occupancy predicted by the model exceeds or approaches 0.5. The occupancy increases with greater altitude and with larger EVI values, while the detection rate increases with altitude only during spring and decreases with EVI values only in summer. The daily activity peaks for white-lipped deer were observed from 7:00 to 11:00 and 17:00 to 22:00, with annual activity peaks occurring from April to June and from September to November. From July to the following January, white-lipped deer mostly move in mixed-sex groups, while during the remainder of the year, they predominantly associate with individuals of the same sex. Climate, vegetation coverage, food resources, and human disturbance collectively influenced the behavior and habitat utilization of white-lipped deer. The foundational research conducted on white-lipped deer over the past two years is expected to enhance the basic understanding of white-lipped deer in the Qinghai-Tibet Plateau and contribute to future protection and management decisions.

Disentangling direct and indirect effects of extreme events on coastal wetland communities

One of the primary ways in which climate change will impact coastal freshwater wetlands is through changes in the frequency, intensity, timing and distribution of extreme weather events. Disentangling the direct and indirect mechanisms of population- and community-level responses to extreme events is vital to predicting how species composition of coastal wetlands will change under future conditions. We extended static structural equation modelling approaches to incorporate system dynamics in a multi-year multispecies occupancy model to quantify the effects of extreme weather events on a coastal freshwater wetland system. We used data from an 8-year study (2009-2016) on St. Marks National Wildlife Refuge in Florida, USA, to quantify species-specific and community-level changes in amphibian and fish occupancy associated with two flooding events in 2012 and 2013. We examine how physical changes to the landscape, including potential changes in salinity and increased wetland connectivity, may have contributed to or exacerbated the effects of these extreme weather events on the biota of isolated coastal wetlands. We provide evidence that the primary effects of flooding on the amphibian community were through indirect mechanisms via changes in the composition of the sympatric fish community that may have had lethal (i.e. through direct predation) or non-lethal (i.e. through direct or indirect competitive interactions) effects. In addition, we have shown that amphibian species differed in their sensitivity to direct flooding effects and indirect changes in the fish community and wetland-specific conductance, which led to variable responses across the community. These effects led to the overall decline in amphibian species richness from 2009 to 2016, suggesting that wetland-breeding amphibian communities on St. Marks National Wildlife Refuge may not be resilient to predicted changes in coastal disturbance regimes because of climate change. Understanding both direct and indirect effects, as well as species interactions, is important for predicting the effects of a changing climate on individual species, communities and ecosystems.

Range‐wide persistence of the endangered arroyo toad ( Anaxyrus californicus ) for 20+ years following a prolonged drought

Prolonged drought due to climate change has negatively impacted amphibians in southern California, U.S.A. Due to the severity and length of the current drought, agencies and researchers had growing concern for the persistence of the arroyo toad (Anaxyrus californicus), an endangered endemic amphibian in this region. Range‐wide surveys for this species had not been conducted for at least 20 years. In 2017–2020, we conducted collaborative surveys for arroyo toads at historical locations. We surveyed 88 of the 115 total sites having historical records and confirmed that the arroyo toad is currently extant in at least 61 of 88 sites and 20 of 25 historically occupied watersheds. We did not detect toads at almost a third of the surveyed sites but did detect toads at 18 of 19 specific sites delineated in the 1999 Recovery Plan to meet one of four downlisting criteria. Arroyo toads are estimated to live 7–8 years, making populations susceptible to prolonged drought. Drought is estimated to increase in frequency and duration with climate change. Mitigation strategies for drought impacts, invasive aquatic species, altered flow regimes, and other anthropogenic effects could be the most beneficial strategies for toad conservation and may also provide simultaneous benefits to several other native species that share the same habitat.

Impacts of a Non-indigenous Ecosystem Engineer, the American Beaver (Castor canadensis), in a Biodiversity Hotspot

Non-native species having high per capita impacts in invaded communities are those that modulate resource availability and alter disturbance regimes in ways that are biologically incompatible with the native biota. In areas where it has been introduced by humans, American beaver (Castor canadensis) is an iconic example of such species due to its capacity to alter trophic dynamics of entire ecosystems and create new invasional pathways for other non-native species. The species is problematic in several watersheds within the Southern California-Northern Baja California Coast Ecoregion, a recognized hotspot of biodiversity, due to its ability to modify habitat in ways that favor invasive predators and competitors over the region's native species and habitat. Beaver was deliberately introduced across California in the mid-1900s and generally accepted as non-native to the region up to the early 2000s; however, articles promoting the idea that beaver may be a natural resident have gained traction in recent years, due in large part to the species' charismatic nature rather than by presentation of sound evidence. Here, we discuss the problems associated with beaver disturbance and its effects on conserving the region's native fauna and flora. We refute arguments underlying the claim that beaver is native to the region, and review paleontological, zooarchaeological, and historical survey data from renowned field biologists and naturalists over the past ~160 years to show that no evidence exists that beaver arrived by any means other than deliberate human introduction. Managing this ecosystem engineer has potential to reduce the richness and abundance of other non-native species because the novel, engineered habitat now supporting these species would diminish in beaver-occupied watersheds. At the same time, hydrologic functionality would shift toward more natural, ephemeral conditions that favor the regions' native species while suppressing the dominance of the most insidious invaders.

What is the effect of poaching activity on wildlife species?

Poaching is a pervasive threat to wildlife, yet quantifying the direct effect of poaching on wildlife is rarely possible because both wildlife and threat data are infrequently collected concurrently. In this study, we used poaching data collected through the Management Information System (MIST) and wildlife camera trap data collected by the Tropical Ecology Assessment and Monitoring (TEAM) network from 2014 to 2017 in Volcanoes National Park, Rwanda. We implemented co-occurrence multi-season occupancy models that accounted for imperfect detection to investigate the effect of poaching on initial occupancy, colonization, and extinction of five mammal species. Specifically, we focused on two species of conservation concern (mountain gorilla [Gorilla beringei beringei] and golden monkey [Cercopithecus mitis kandti]), and three species targeted by poachers (black-fronted duiker [Cephalophus nigrifrons], bushbuck [Tragelaphus scriptus], and African buffalo [Syncerus caffer]). We found that the probability of local extinction was highest in sites with poaching activity for golden monkey and bushbuck. In addition, the probability of initial occupancy for golden monkey was highest in sites without poaching activity. We only found weak evidence of effects of poaching on parameters governing the occupancy dynamics of the other species. All species showed evidence of poaching presence affecting the probability of detection of the wildlife species. This is the first study to our knowledge to combine direct threat observations from ranger-based monitoring data with camera trap wildlife observations to quantify the effect of poaching on wildlife. Given the widespread collection of ranger-based monitoring and camera trap data, our approach is broadly applicable to numerous protected areas and has the potential to significantly improve conservation management. Specifically, the relationship between poaching activity and wildlife population dynamics can be combined with information on the relationship between ranger patrols and poaching activity to develop models useful for making wise decisions about ranger patrol deployment.

Estimation of breeding probability can make monitoring data more revealing: a case study of amphibians

Monitoring programs serve to detect trends in the distribution and abundance of species. To do so, monitoring programs often use static state variables. Dynamic state variables that describe population dynamics might be more valuable because they allow for a mechanistic understanding of the processes that lead to population trends. We fit multistate occupancy models to data from a country-wide multispecies amphibian occupancy monitoring program and estimated occupancy and breeding probabilities. If breeding probabilities are determinants of occupancy dynamics, then they may serve in monitoring programs as state variables that describe dynamic processes. The results showed that breeding probabilities were low and that a large proportion of the populations had to be considered to be non-breeding populations (i.e., populations where adults are present but no breeding occurs). For some species, the majority of populations were non-breeding populations. We found that non-breeding populations have lower persistence probabilities than populations where breeding occurs. Breeding probabilities may thus explain trends in occupancy but they might also explain other ecological phenomena, such as the success of invasive species, which had high breeding probabilities. Signs of breeding, i.e., the presence of eggs and larvae, were often hard to detect. Importantly, non-breeding populations also had low detection probabilities, perhaps because they had lower abundances. We suggest that monitoring programs should invest more in the detection of life history stages indicative of breeding, and also into the detection of non-breeding populations. We conclude that breeding probability should be used as a state variable in monitoring programs because it can lead to deeper insights into the processes driving occupancy dynamics.

A note on investigating co-occurrence patterns and dynamics for many species, with imperfect detection and a log-linear modeling parameterization

Patterns in, and the underlying dynamics of, species co-occurrence is of interest in many ecological applications. Unaccounted for, imperfect detection of the species can lead to misleading inferences about the nature and magnitude of any interaction. A range of different parameterizations have been published that could be used with the same fundamental modeling framework that accounts for imperfect detection, although each parameterization has different advantages and disadvantages.We propose a parameterization based on log-linear modeling that does not require a species hierarchy to be defined (in terms of dominance) and enables a numerically robust approach for estimating covariate effects.Conceptually, the parameterization is equivalent to using the presence of species in the current, or a previous, time period as predictor variables for the current occurrence of other species. This leads to natural, "symmetric," interpretations of parameter estimates.The parameterization can be applied to many species, in either a maximum likelihood or Bayesian estimation framework. We illustrate the method using camera-trapping data collected on three mesocarnivore species in South Texas.

The impact of climate change induced alterations of streamflow and stream temperature on the distribution of riparian species

Distributions of riparian species will likely shift due to climate change induced alterations in temperature and rainfall patterns, which alter stream habitat. Spatial forecasting of suitable habitat in projected climatic conditions will inform management interventions that support wildlife. Challenges in developing forecasts include the need to consider the large number of riparian species that might respond differently to changing conditions and the need to evaluate the many different characteristics of streamflow and stream temperature that drive species-specific habitat suitability. In particular, in dynamic environments like streams, the short-term temporal resolution of species occurrence and streamflow need to be considered to identify the types of conditions that support various species. To address these challenges, we cluster species based on habitat characteristics to select habitat representatives and we evaluate regional changes in habitat suitability using short-term, temporally explicit metrics that describe the streamflow and stream temperature regime. We use stream-specific environmental predictors rather than climatic variables. Unlike other studies, the stream-specific environmental predictors are generated from the time that species were observed in a particular reach, in addition to long term trends, to evaluate habitat preferences. With species occurrence data from local monitoring surveys and streamflow and stream temperature modeled from downscaled Coupled Model Intercomparison Project - Phase 5 (CMIP5) climate projections, we predict change in habitat suitability at the end-of-century. The relative importance of hydrology and stream temperature varied by cluster. High altitudinal, cold water species' distributions contracted, while lower elevation, warm water species distributions expanded. Modeling with short-term temporally explicit environmental metrics did produce different end-of-century projections than using long-term averages for some of the representative species. These findings can help wildlife managers prioritize conservation efforts, manage streamflow, initiate monitoring of species in vulnerable clusters, and address stressors, such as passage barriers, in areas projected to be suitable in future climate conditions.

Historical museum collections and contemporary population studies implicate roads and introduced predatory bullfrogs in the decline of western pond turtles

The western pond turtle (WPT), recently separated into two paripatrically distributed species (Emys pallida and Emys marmorata), is experiencing significant reductions in its range and population size. In addition to habitat loss, two potential causes of decline are female-biased road mortality and high juvenile mortality from non-native predatory bullfrogs (Rana catesbeiana). However, quantitative analyses of these threats have never been conducted for either species of WPT. We used a combination of historical museum samples and published and unpublished field studies shared with us through personal communications with WPT field researchers (B. Shaffer, P. Scott, R. Fisher, C. Brown, R. Dagit, L. Patterson, T. Engstrom, 2019, personal communications) to quantify the effect of roads and bullfrogs on WPT populations along the west coast of the United States. Both species of WPT shift toward increasingly male biased museum collections over the last century, a trend consistent with increasing, female-biased road mortality. Recent WPT population studies revealed that road density and proximity were significantly associated with increasingly male-biased sex ratios, further suggesting female-biased road mortality. The mean body size of museum collections of E. marmorata, but not E. pallida, has increased over the last 100 years, consistent with reduced recruitment and aging populations that could be driven by invasive predators. Contemporary WPT population sites that co-occur with bullfrogs had significantly greater average body sizes than population sites without bullfrogs, suggesting strong bullfrog predation on small WPT hatchlings and juveniles. Overall, our findings indicate that both species of WPT face demographic challenges which would have been difficult to document without the use of both historical data from natural history collections and contemporary demographic field data. Although correlational, our analyses suggest that female-biased road mortality and predation on small turtles by non-native bullfrogs are occurring, and that conservation strategies reducing both may be important for WPT recovery.

Estimation of metapopulation colonization rates from disturbance history and occurrence-pattern data

Occurrence patterns of many sessile species in dynamic landscapes are not in equilibrium due to their slow rates of metapopulation colonization and extinction. Colonization-extinction data enable the estimation of colonization rates for such species, but collecting the necessary data may require long waiting times between sampling years. Methods for estimating colonization rates of nonequilibrium metapopulations from single occurrence-pattern data have so far relied on additional data on patch ages and on past patch connectivities. We present an approach where metapopulation colonization rates are estimated from occurrence-pattern data and from disturbance history data that inform of past patch dynamics and that can be collected together with occurrence-pattern data. We estimated parameter values regulating patch and metapopulation dynamics by simulating patch network and metapopulation histories that result in present-like patch network configurations and metapopulation occurrence patterns. We tested our approach using occurrence-pattern data of the epiphytic lichen Lobaria pulmonaria in Fennoscandian forests, and fire-scar data that inform of the 400-yr history of fires and host tree dynamics in the same landscapes. The estimated model parameters were similar to estimates obtained using colonization-extinction data. The projected L. pulmonaria occupancy into the future also agreed with the respective projections that were made using the model estimated from colonization-extinction data. Our approach accelerates the estimation of metapopulation colonization rates for sessile species that are not in metapopulation equilibrium with the current landscape structure.

Drought-mediated extinction of an arid-land amphibian: insights from a spatially explicit dynamic occupancy model

Understanding how natural and anthropogenic processes affect population dynamics of species with patchy distributions is critical to predicting their responses to environmental changes. Despite considerable evidence that demographic rates and dispersal patterns vary temporally in response to an array of biotic and abiotic processes, few applications of metapopulation theory have sought to explore factors that explain spatiotemporal variation in extinction or colonization rates. To facilitate exploring these factors, we extended a spatially explicit model of metapopulation dynamics to create a framework that requires only binary presence-absence data, makes few assumptions about the dispersal process, and accounts for imperfect detection. We apply this framework to 22 yr of biannual survey data for lowland leopard frogs, Lithobates yavapaiensis, an amphibian that inhabits arid stream systems in the southwestern United States and northern Mexico. Our results highlight the importance of accounting for factors that govern temporal variation in transition probabilities, as both extinction and colonization rates varied with hydrologic conditions. Specifically, local extinctions were more frequent during drought periods, particularly at sites without reliable surface water. Colonization rates increased when larval and dispersal periods were wetter than normal, which increased the probability that potential emigrants metamorphosed and reached neighboring sites. Extirpation of frogs from all sites in one watershed during a period of severe drought demonstrated the influence of site-level features, as frogs persisted only in areas where most sites held water consistently and where the amount of sediment deposited from high-elevation wildfires was low. Application of our model provided novel insights into how climate-related processes affected the distribution and population dynamics of an arid-land amphibian. The approach we describe has application to a wide array of species that inhabit patchy environments, can improve our understanding of factors that govern metapopulation dynamics, and can inform strategies for conservation of imperiled species.

Quantifying climate sensitivity and climate-driven change in North American amphibian communities

Changing climate will impact species' ranges only when environmental variability directly impacts the demography of local populations. However, measurement of demographic responses to climate change has largely been limited to single species and locations. Here we show that amphibian communities are responsive to climatic variability, using >500,000 time-series observations for 81 species across 86 North American study areas. The effect of climate on local colonization and persistence probabilities varies among eco-regions and depends on local climate, species life-histories, and taxonomic classification. We found that local species richness is most sensitive to changes in water availability during breeding and changes in winter conditions. Based on the relationships we measure, recent changes in climate cannot explain why local species richness of North American amphibians has rapidly declined. However, changing climate does explain why some populations are declining faster than others. Our results provide important insights into how amphibians respond to climate and a general framework for measuring climate impacts on species richness.

Ecological correlates of the spatial co-occurrence of sympatric mammalian carnivores worldwide

The composition of local mammalian carnivore communities has far-reaching effects on terrestrial ecosystems worldwide. To better understand how carnivore communities are structured, we analysed camera trap data for 108 087 trap days across 12 countries spanning five continents. We estimate local probabilities of co-occurrence among 768 species pairs from the order Carnivora and evaluate how shared ecological traits correlate with probabilities of co-occurrence. Within individual study areas, species pairs co-occurred more frequently than expected at random. Co-occurrence probabilities were greatest for species pairs that shared ecological traits including similar body size, temporal activity pattern and diet. However, co-occurrence decreased as compared to other species pairs when the pair included a large-bodied carnivore. Our results suggest that a combination of shared traits and top-down regulation by large carnivores shape local carnivore communities globally.

Host-pathogen metapopulation dynamics suggest high elevation refugia for boreal toads

Emerging infectious diseases are an increasingly common threat to wildlife. Chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), is an emerging infectious disease that has been linked to amphibian declines around the world. Few studies exist that explore amphibian-Bd dynamics at the landscape scale, limiting our ability to identify which factors are associated with variation in population susceptibility and to develop effective in situ disease management. Declines of boreal toads (Anaxyrus boreas boreas) in the southern Rocky Mountains are largely attributed to chytridiomycosis but variation exists in local extinction of boreal toads across this metapopulation. Using a large-scale historic data set, we explored several potential factors influencing disease dynamics in the boreal toad-Bd system: geographic isolation of populations, amphibian community richness, elevational differences, and habitat permanence. We found evidence that boreal toad extinction risk was lowest at high elevations where temperatures may be suboptimal for Bd growth and where small boreal toad populations may be below the threshold needed for efficient pathogen transmission. In addition, boreal toads were more likely to recolonize high elevation sites after local extinction, again suggesting that high elevations may provide refuge from disease for boreal toads. We illustrate a modeling framework that will be useful to natural resource managers striving to make decisions in amphibian-Bd systems. Our data suggest that in the southern Rocky Mountains high elevation sites should be prioritized for conservation initiatives like reintroductions.

Longevity and population age structure of the arroyo southwestern toad (Anaxyrus californicus) with drought implications

The arroyo southwestern toad is a specialized and federally endangered amphibian endemic to the coastal plains and mountains of central and southern California and northwestern Baja California. It is largely unknown how long these toads live in natural systems, how their population demographics vary across occupied drainages, and how hydrology affects age structure. We used skeletochronology to estimate the ages of adult arroyo toads in seven occupied drainages with varying surface water hydrology in southern California. We processed 179 adult toads with age estimates between 1 and 6 years. Comparisons between skeletochronological ages and known ages of PIT tagged toads showed that skeletochronology likely underestimated toad age by up to 2 years, indicating they may live to 7 or 8 years, but nonetheless major patterns were evident. Arroyo toads showed sexual size dimorphism with adult females reaching a maximum size of 12 mm greater than males. Population age structure varied among the sites. Age structure at sites with seasonally predictable surface water was biased toward younger individuals, which indicated stable recruitment for these populations. Age structures at the ephemeral sites were biased toward older individuals with cohorts roughly corresponding to higher rainfall years. These populations are driven by surface water availability, a stochastic process, and thus more unstable. Based on our estimates of toad ages, climate predictions of extreme and prolonged drought events could mean that the number of consecutive dry years could surpass the maximum life span of toads making them vulnerable to extirpation, especially in ephemeral freshwater systems. Understanding the relationship between population demographics and hydrology is essential for predicting species resilience to projected changes in weather and rainfall patterns. The arroyo toad serves as a model for understanding potential responses to climatic and hydrologic changes in Mediterranean stream systems. We recommend development of adaptive management strategies to address these threats.

Spatio-temporal interactions facilitate large carnivore sympatry across a resource gradient

Species within a guild vary their use of time, space and resources, thereby enabling sympatry. As intra-guild competition intensifies, such behavioural adaptations may become prominent. We assessed mechanisms of facilitating sympatry among dhole (Cuon alpinus), leopard (Panthera pardus) and tiger (Panthera tigris) in tropical forests of India using camera-trap surveys. We examined population-level temporal, spatial and spatio-temporal segregation among them across four reserves representing a gradient of carnivore and prey densities. Temporal and spatial overlaps were higher at lower prey densities. Combined spatio-temporal overlap was minimal, possibly due to chance. We found fine-scale avoidance behaviours at one high-density reserve. Our results suggest that: (i) patterns of spatial, temporal and spatio-temporal segregation in sympatric carnivores do not necessarily mirror each other; (ii) carnivores are likely to adopt temporal, spatial, and spatio-temporal segregation as alternative mechanisms to facilitate sympatry; and (iii) carnivores show adaptability across a gradient of resource availability, a driver of inter-species competition. We discuss behavioural mechanisms that permit carnivores to co-occupy rather than dominate functional niches, and adaptations to varying intensities of competition that are likely to shape structure and dynamics of carnivore guilds.

When habitat matters: Habitat preferences can modulate co-occurrence patterns of similar sympatric species

Disentangling the role of competition in regulating the distribution of sympatric species can be difficult because species can have different habitat preferences or time use that introduce non-random patterns that are not related to interspecific interactions. We adopted a multi-step approach to systematically incorporate habitat preferences while investigating the co-occurrence of two presumed competitors, morphologically similar, and closely related ground-dwelling birds: the brown tinamou (Crypturellus obsoletus) and the tataupa tinamou (C. tataupa). First, we used single-species occupancy models to identify the main landscape characteristics affecting site occupancy, while accounting for detection probability. We then used these factors to control for the effect of habitat while investigating species co-occurrence. In addition, we investigated species present-time partitioning by measuring the degree of overlap in their activity time. Both species were strictly diurnal and their activity time highly overlapped (i.e., the species are not present-time partitioning). The distribution of the two species varied across the landscape, and they seemed to occupy opposite portions of the study area, but co-occurrence models and species interaction factors suggested that the tinamous have independent occupancy and detection. In addition, co-occurrence models that accounted for habitat performed better than models without habitat covariates. The observed co-occurrence pattern is more likely related to habitat preferences, wherein species segregated by elevation. These results provide evidence that habitat characteristics can play a bigger role than interspecific interactions in regulating co-existence of some species. Therefore, exploring habitat preferences while analyzing co-occurrence patterns is essential, in addition to being a feasible approach to achieve more accurate estimation of parameters reflecting species interactions. Occupancy models can be a valuable tool in such modeling.

Jaguar interactions with pumas and prey at the northern edge of jaguars' range

We present the first study that evaluates jaguar-puma interactions in the arid lands of northern Mexico, where jaguars have their northernmost breeding population and both predators are persecuted for livestock depredation. We tested whether jaguars are the dominant species in this unique ecosystem, where: (1) pumas outnumber jaguars, (2) pumas are better adapted to arid environments, and (3) jaguars and pumas are of similar size. We analyzed four years of data with two approaches; a two species conditional occupancy model and an activity patterns analysis. We used camera location and prey presence as covariates for jaguar and puma detection and presence probabilities. We also explored overlap in activities of predators and prey. Where both species were detected, peccary presence was positively correlated with both jaguar and puma presence, whereas in areas where jaguars were detected but pumas were not, deer presence explained the probability of jaguar presence. We found that both predators were more likely to co-occur together than to be found independently, and so we rejected the hypothesis that jaguars were the dominant species in our study area. Predators were mainly nocturnal and their activity patterns overlapped by 60%. Jaguar, as compared with puma, overlapped more with deer and calves; puma overlapped with calves more than with other prey, suggesting a preference. We believe exploring predator relationships at different scales may help elucidate mechanisms that regulate their coexistence.

Life history plasticity does not confer resilience to environmental change in the mole salamander (Ambystoma talpoideum)

Plasticity in life history strategies can be advantageous for species that occupy spatially or temporally variable environments. We examined how phenotypic plasticity influences responses of the mole salamander, Ambystoma talpoideum, to disturbance events at the St. Marks National Wildlife Refuge (SMNWR), FL, USA from 2009 to 2014. We observed periods of extensive drought early in the study, in contrast to high rainfall and expansive flooding events in later years. Flooding facilitated colonization of predatory fishes to isolated wetlands across the refuge. We employed multistate occupancy models to determine how this natural experiment influenced the occurrence of aquatic larvae and paedomorphic adults and what implications this may have for the population. We found that, in terms of occurrence, responses to environmental variation differed between larvae and paedomorphs, but plasticity (i.e. the ability to metamorphose rather than remain in aquatic environment) was not sufficient to buffer populations from declining as a result of environmental perturbations. Drought and fish presence negatively influenced occurrence dynamics of larval and paedomorphic mole salamanders and, consequently, contributed to observed short-term declines of this species. Overall occurrence of larval salamanders decreased from 0.611 in 2009 to 0.075 in 2014 and paedomorph occurrence decreased from 0.311 in 2009 to 0.121 in 2014. Although variation in selection pressures has likely maintained this polyphenism previously, our results suggest that continued changes in environmental variability and the persistence of fish in isolated wetlands could lead to a loss of paedomorphosis in the SMNWR population and, ultimately, impact regional persistence in the future.

Species interactions and the effects of climate variability on a wetland amphibian metacommunity

Disentangling the role that multiple interacting factors have on species responses to shifting climate poses a significant challenge. However, our ability to do so is of utmost importance to predict the effects of climate change on species distributions. We examined how populations of three species of wetland-breeding amphibians, which varied in life history requirements, responded to a six-year period of extremely variable precipitation. This interval was punctuated by both extensive drought and heavy precipitation and flooding, providing a natural experiment to measure community responses to environmental perturbations. We estimated occurrence dynamics using a discrete hidden Markov modeling approach that incorporated information regarding habitat state and predator-prey interactions. This approach allowed us to measure how metapopulation dynamics of each amphibian species was affected by interactions among weather, wetland hydroperiod, and co-occurrence with fish predators. The pig frog, a generalist, proved most resistant to perturbations, with both colonization and persistence being unaffected by seasonal variation in precipitation or co-occurrence with fishes. The ornate chorus frog, an ephemeral wetland specialist, responded positively to periods of drought owing to increased persistence and colonization rates during periods of low-rainfall. Low probabilities of occurrence of the ornate chorus frog in long-duration wetlands were driven by interactions with predators due to low colonization rates when fishes were present. The mole salamander was most sensitive to shifts in water availability. In our study area, this species never occurred in short-duration wetlands and persistence probabilities decreased during periods of drought. At the same time, negative effects occurred with extreme precipitation because flooding facilitated colonization of fishes to isolated wetlands and mole salamanders did not colonize wetlands once fishes were present. We demonstrate that the effects of changes in water availability depend on interactions with predators and wetland type and are influenced by the life history of each of our species. The dynamic species occurrence modeling approach we used offers promise for other systems when the goal is to disentangle the complex interactions that determine species responses to environmental variability.

Quantifying population declines based on presence-only records for red-list assessments

Accurate trend estimates are necessary for understanding which species are declining and which are most in need of conservation action. Imperfect species detection may result in unreliable trend estimates because this may lead to the overestimation of declines. Because many management decisions are based on population trend estimates, such biases could have severe consequences for conservation policy. We used an occupancy-modeling framework to estimate detectability and calculate nationwide population trends for 14 Swiss amphibian species both accounting for and ignoring imperfect detection. Through the application of International Union for Conservation of Nature Red List criteria to the different trend estimates, we assessed whether ignoring imperfect detection could affect conservation policy. Imperfect detection occurred for all species and detection varied substantially among species, which led to the overestimation of population declines when detectability was ignored. Consequently, accounting for imperfect detection lowered the red-list risk category for 5 of the 14 species assessed. We demonstrate that failing to consider species detectability can have serious consequences for species management and that occupancy modeling provides a flexible framework to account for observation bias and improve assessments of conservation status. A problem inherent to most historical records is that they contain presence-only data from which only relative declines can be estimated. A move toward the routine recording of nonobservation and absence data is essential if conservation practitioners are to move beyond this toward accurate population trend estimation.

Estimating Effects of Species Interactions on Populations of Endangered Species

Global change causes community composition to change considerably through time, with ever-new combinations of interacting species. To study the consequences of newly established species interactions, one available source of data could be observational surveys from biodiversity monitoring. However, approaches using observational data would need to account for niche differences between species and for imperfect detection of individuals. To estimate population sizes of interacting species, we extended N-mixture models that were developed to estimate true population sizes in single species. Simulations revealed that our model is able to disentangle direct effects of dominant on subordinate species from indirect effects of dominant species on detection probability of subordinate species. For illustration, we applied our model to data from a Swiss amphibian monitoring program and showed that sizes of expanding water frog populations were negatively related to population sizes of endangered yellow-bellied toads and common midwife toads and partly of natterjack toads. Unlike other studies that analyzed presence and absence of species, our model suggests that the spread of water frogs in Central Europe is one of the reasons for the decline of endangered toad species. Thus, studying population impacts of dominant species on population sizes of endangered species using data from biodiversity monitoring programs should help to inform conservation policy and to decide whether competing species should be subject to population management.

Movement, demographics, and occupancy dynamics of a federally-threatened salamander: evaluating the adequacy of critical habitat

Critical habitat for many species is often limited to occupied localities. For rare and cryptic species, or those lacking sufficient data, occupied habitats may go unrecognized, potentially hindering species recovery. Proposed critical habitat for the aquatic Jollyville Plateau salamander (Eurycea tonkawae) and two sister species were delineated based on the assumption that surface habitat is restricted to springs and excludes intervening stream reaches. To test this assumption, we performed two studies to understand aspects of individual, population, and metapopulation ecology of E. tonkawae. First, we examined movement and population demographics using capture-recapture along a spring-influenced stream reach. We then extended our investigation of stream habitat use with a study of occupancy and habitat dynamics in multiple headwater streams. Indications of extensive stream channel use based on capture-recapture results included frequent movements of >15 m, and high juvenile abundance downstream of the spring. Initial occupancy of E. tonkawae was associated with shallow depths, maidenhair fern presence and low temperature variation (indicative of groundwater influence), although many occupied sites were far from known springs. Additionally, previously dry sites were three times more likely to be colonized than wet sites. Our results indicate extensive use of stream habitats, including intermittent ones, by E. tonkawae. These areas may be important for maintaining population connectivity or even as primary habitat patches. Restricting critical habitat to occupied sites will result in a mismatch with actual habitat use, particularly when assumptions of habitat use are untested, thus limiting the potential for recovery.

Integrating Multiple Distribution Models to Guide Conservation Efforts of an Endangered Toad

Species distribution models are used for numerous purposes such as predicting changes in species' ranges and identifying biodiversity hotspots. Although implications of distribution models for conservation are often implicit, few studies use these tools explicitly to inform conservation efforts. Herein, we illustrate how multiple distribution models developed using distinct sets of environmental variables can be integrated to aid in identification sites for use in conservation. We focus on the endangered arroyo toad (Anaxyrus californicus), which relies on open, sandy streams and surrounding floodplains in southern California, USA, and northern Baja California, Mexico. Declines of the species are largely attributed to habitat degradation associated with vegetation encroachment, invasive predators, and altered hydrologic regimes. We had three main goals: 1) develop a model of potential habitat for arroyo toads, based on long-term environmental variables and all available locality data; 2) develop a model of the species' current habitat by incorporating recent remotely-sensed variables and only using recent locality data; and 3) integrate results of both models to identify sites that may be employed in conservation efforts. We used a machine learning technique, Random Forests, to develop the models, focused on riparian zones in southern California. We identified 14.37% and 10.50% of our study area as potential and current habitat for the arroyo toad, respectively. Generally, inclusion of remotely-sensed variables reduced modeled suitability of sites, thus many areas modeled as potential habitat were not modeled as current habitat. We propose such sites could be made suitable for arroyo toads through active management, increasing current habitat by up to 67.02%. Our general approach can be employed to guide conservation efforts of virtually any species with sufficient data necessary to develop appropriate distribution models.

The effects of hydropattern and predator communities on amphibian occupancy

Complex, interactive ecological constraints regulate species distributions, and understanding these factors is crucial for predicting species persistence. We used occupancy analysis, which corrects for imperfect detection, to test the importance of abiotic and biotic habitat and landscape factors on probability of occupancy by Boreal Chorus Frog (Pseudacris maculata (Agassiz, 1850)) tadpoles. We hypothesized that hydropattern and predators are primarily important because they affect desiccation and predation risk and can interact in ways difficult to predict. We surveyed 62 wetland sites across an elevational gradient in Colorado, USA, and modeled patterns in P. maculata occupancy. Tadpoles were most frequently present in intermediate-length hydropattern systems with lower desiccation risk and no predatory fish because of occasional drying. Pseudacris maculata occupancy had a strong negative relationship with fish presence, while tadpoles, odonate larvae, and Barred Tiger Salamanders (Ambystoma mavortium mavortium Baird, 1850) frequently co-occurred. Dry seasonal conditions will likely result in fewer intermediate-length hydropattern ponds available for amphibian breeding. We hypothesize that this will force P. maculata to breed in habitats with fish. As habitats shrink, predators that co-occur with P. maculata are expected to concentrate in the remaining habitat and increase predation risk for developing tadpoles (assuming predators are similarly constricted in their habitat use as amphibians are).