A Dataset of Amphibian Species in U.S. National Parks

National parks and other protected areas are important for preserving landscapes and biodiversity worldwide. An essential component of the mission of the United States (U.S.) National Park Service (NPS) requires understanding and maintaining accurate inventories of species on protected lands. We describe a new, national-scale synthesis of amphibian species occurrence in the NPS system. Many park units have a list of amphibian species observed within their borders compiled from various sources and available publicly through the NPSpecies platform. However, many of the observations in NPSpecies remain unverified and the lists are often outdated. We updated the amphibian dataset for each park unit by collating old and new park-level records and had them verified by regional experts. The new dataset contains occurrence records for 292 of the 424 NPS units and includes updated taxonomy, international and state conservation rankings, hyperlinks to a supporting reference for each record, specific notes, and related fields which can be used to better understand and manage amphibian biodiversity within a single park or group of parks.

Broad-Scale Assessment of Methylmercury in Adult Amphibians

Mercury (Hg) is a toxic contaminant that has been mobilized and distributed worldwide and is a threat to many wildlife species. Amphibians are facing unprecedented global declines due to many threats including contaminants. While the biphasic life history of many amphibians creates a potential nexus for methylmercury (MeHg) exposure in aquatic habitats and subsequent health effects, the broad-scale distribution of MeHg exposure in amphibians remains unknown. We used nonlethal sampling to assess MeHg bioaccumulation in 3,241 juvenile and adult amphibians during 2017-2021. We sampled 26 populations (14 species) across 11 states in the United States, including several imperiled species that could not have been sampled by traditional lethal methods. We examined whether life history traits of species and whether the concentration of total mercury in sediment or dragonflies could be used as indicators of MeHg bioaccumulation in amphibians. Methylmercury contamination was widespread, with a 33-fold difference in concentrations across sites. Variation among years and clustered subsites was less than variation across sites. Life history characteristics such as size, sex, and whether the amphibian was a frog, toad, newt, or other salamander were the factors most strongly associated with bioaccumulation. Total Hg in dragonflies was a reliable indicator of bioaccumulation of MeHg in amphibians (R2 ≥ 0.67), whereas total Hg in sediment was not (R2 ≤ 0.04). Our study, the largest broad-scale assessment of MeHg bioaccumulation in amphibians, highlights methodological advances that allow for nonlethal sampling of rare species and reveals immense variation among species, life histories, and sites. Our findings can help identify sensitive populations and provide environmentally relevant concentrations for future studies to better quantify the potential threats of MeHg to amphibians.

Assay validation of saliva glucocorticoids in Columbia spotted frogs and effects of handling and marking

Non-invasive methods are important to the field of conservation physiology to reduce negative effects on organisms being studied. Glucocorticoid (GC) hormones are often used to assess health of individuals, but collection methods can be invasive. Many amphibians are imperiled worldwide, and saliva is a non- or semi-invasive matrix to measure GCs that has been partially validated for only four amphibian species. Validation ensures that assays are reliable and can detect changes in saliva corticosterone (sCORT) after exposure to stressors, but it is also necessary to ensure sCORT concentrations are correlated with plasma concentrations. To help validate the use of saliva in assessing CORT responses in amphibians, we captured uniquely marked Columbia spotted frogs (Rana luteiventris) on sequential days and collected baseline and stress-induced (after handling) samples. For a subset of individuals, we collected and quantified CORT in both saliva and blood samples, which have not been compared for amphibians. We tested several aspects of CORT responses and, by collecting across separate days, measured repeatability of CORT responses across days. We also evaluated whether methods common to amphibian conservation, such as handling alone or handling, clipping a toe and tagging elevated sCORT. Similar to previous studies, we show that sCORT is reliable concerning parallelism, recovery, precision and sensitivity. sCORT was weakly correlated with plasma CORT (R2 = 0.21), and we detected elevations in sCORT after handling, demonstrating biological validation. Toe clipping and tagging did not increase sCORT over handling alone, but repeated handling elevated sCORT for ~72 hours. However, sCORT responses were highly variable and repeatability was low within individuals and among capture sessions, contrary to previous studies with urinary and waterborne CORT. sCORT is a semi-invasive and rapid technique that could be useful to assess effects of anthropogenic change and conservation efforts, but will require careful study design and future validation.

Energy-related wastewater contamination alters microbial communities of sediment, water, and amphibian skin

To inform responsible energy development, it is important to understand the ecological effects of contamination events. Wastewaters, a common byproduct of oil and gas extraction, often contain high concentrations of sodium chloride (NaCl) and heavy metals (e.g., strontium and vanadium). These constituents can negatively affect aquatic organisms, but there is scarce information for how wastewaters influence potentially distinct microbiomes in wetland ecosystems. Additionally, few studies have concomitantly investigated effects of wastewaters on the habitat (water and sediment) and skin microbiomes of amphibians or relationships among these microbial communities. We sampled microbiomes of water, sediment, and skin of four larval amphibian species across a gradient of chloride contamination (0.04-17,500 mg/L Cl) in the Prairie Pothole Region of North America. We detected 3129 genetic phylotypes and 68 % of those phylotypes were shared among the three sample types. The most common shared phylotypes were Proteobacteria, Firmicutes, and Bacteroidetes. Salinity of wastewaters increased dissimilarity within all three microbial communities, but not the diversity or richness of water and skin microbial communities. Strontium was associated with lower diversity and richness of sediment microbial communities, but not those of water or amphibian skin, likely because metal deposition occurs in sediment when wetlands dry. Based on Bray Curtis distance matrices, sediment microbiomes were similar to those of water, but neither had substantial overlap with amphibian microbiomes. Species identity was the strongest predictor of amphibian microbiomes; frog microbiomes were similar but differed from that of the salamander, whose microbiome had the lowest richness and diversity. Understanding how effects of wastewaters on the dissimilarity, richness, and diversity of microbial communities also influence the ecosystem function of communities will be an important next step. However, our study provides novel insight into the characteristics of, and associations among, different wetland microbial communities and effects of wastewaters from energy production.

Empirical evidence for effects of invasive American Bullfrogs on occurrence of native amphibians and emerging pathogens

Invasive species and emerging infectious diseases are two of the greatest threats to biodiversity. American Bullfrogs (Rana [Lithobates] catesbeiana), which have been introduced to many parts of the world, are often linked with declines in native amphibians via predation and the spread of emerging pathogens such as amphibian chytrid fungus (Batrachochytrium dendrobatidis [Bd]) and ranaviruses. Although many studies have investigated the potential role of bullfrogs in the decline of native amphibians, analyses that account for shared habitat affinities and imperfect detection have found limited support for clear effects. Similarly, the role of bullfrogs in shaping the patch-level distribution of pathogens is unclear. We used eDNA methods to sample 233 sites in the southwestern USA and Sonora, Mexico (2016-2018) to estimate how the presence of bullfrogs affects the occurrence of four native amphibians, Bd, and ranaviruses. Based on two-species, dominant-subordinate occupancy models fitted in a Bayesian context, federally threatened Chiricahua Leopard Frogs (Rana chiricahuensis) and Western Tiger Salamanders (Ambystoma mavortium) were eight times (32% vs. 4%) and two times (36% vs. 18%), respectively, less likely to occur at sites where bullfrogs occurred. Evidence for the negative effects of bullfrogs on Lowland Leopard Frogs (Rana yavapaiensis) and Northern Leopard Frogs (Rana pipiens) was less clear, possibly because of smaller numbers of sites where these native species still occurred and because bullfrogs often occur at lower densities in streams, the primary habitat for Lowland Leopard Frogs. At the community level, Bd was most likely to occur where bullfrogs co-occurred with native amphibians, which could increase the risk to native species. Ranaviruses were estimated to occur at 33% of bullfrog-only sites, 10% of sites where bullfrogs and native amphibians co-occurred, and only 3% of sites where only native amphibians occurred. Of the 85 sites where we did not detect any of the five target amphibian species, we also did not detect Bd or ranaviruses; this suggests other hosts do not drive the distribution of these pathogens in our study area. Our results provide landscape-scale evidence that bullfrogs reduce the occurrence of native amphibians and increase the occurrence of pathogens, information that can clarify risks and aid the prioritization of conservation actions.

Testing whether adrenal steroids mediate phenotypic and physiologic effects of elevated salinity on larval tiger salamanders

Salinity (sodium chloride, NaCl) from anthropogenic sources is a persistent contaminant that negatively affects freshwater taxa. Amphibians can be susceptible to salinity, but some species are innately or adaptively tolerant. Physiological mechanisms mediating tolerance to salinity are still unclear, but changes in osmoregulatory hormones such as corticosterone (CORT) and aldosterone (ALDO) are prime candidates. We exposed larval barred tiger salamanders (Ambystoma mavortium) to environmentally relevant NaCl treatments (<32-4000 mg·L^-1 ) for 24 days to test effects on growth, survival, and waterborne CORT responses. Of those sampled, we also quantified waterborne ALDO from a subset. Using a glucocorticoid antagonist (RU486), we also experimentally suppressed CORT signaling of some larvae to determine if CORT mediates effects of salinity. There were no strong differences in survival among salinity treatments, but salinity reduced dry mass, snout-vent length, and body condition while increasing water content of larvae. High survival and sublethal effects demonstrated that salamanders were physiologically challenged but could tolerate the experimental concentrations. CORT signaling did not attenuate sublethal effects of salinity. Baseline and stress-induced (after an acute stressor, shaking) CORT were not influenced by salinity. ALDO was correlated with baseline CORT, suggesting it could be difficult to decouple the roles of CORT and ALDO. Future studies comparing ALDO and CORT responses of adaptively tolerant and previously unexposed populations could be beneficial to understand the roles of these hormones in tolerance to salinity. Nevertheless, our study enhances our understanding of the roles of corticosteroid hormones in mediating effects of a prominent anthropogenic stressor.

Using physiological conditions to assess current and future habitat use of a Subarctic frog

Species with especially close dependence on the environment to meet physiological requirements, such as ectotherms, are highly susceptible to the impacts of climate change. Climate change is occurring rapidly in the Subarctic and Arctic, but there is limited knowledge on ectotherm physiology in these landscapes. We investigated how environmental conditions and habitat characteristics influence the physiological conditions and habitat use of wood frogs (Rana sylvatica) in a Subarctic landscape near Churchill, Manitoba (Canada). We used plaster models to estimate water loss rates and surface body temperatures among different habitat types and at specific locations used by radio-tracked frogs. Water loss (R² = 0.67) and surface temperature (R² = 0.80) of plaster models was similar to that of live frogs. Model-based water loss rates were greater in tundra habitat than in boreal forest and ecotone habitat. Habitat use of wood frogs was strongly tied with available surface moisture and decreased water loss rates that were observed with plaster models. Environmental conditions, such as wind speed and ground temperature, explained 58% and 91% of the variation in water balance and temperature of plaster models. Maintaining physiological conditions may be challenging for semi-aquatic ectotherms in environments vulnerable to future climate change. The ability to predict physiological conditions based on environmental conditions, as demonstrated in our study, can help understand how wildlife will respond to climatic changes.

Diverse aging rates in ectothermic tetrapods provide insights for the evolution of aging and longevity

Comparative studies of mortality in the wild are necessary to understand the evolution of aging; yet, ectothermic tetrapods are underrepresented in this comparative landscape, despite their suitability for testing evolutionary hypotheses. We present a study of aging rates and longevity across wild tetrapod ectotherms, using data from 107 populations (77 species) of nonavian reptiles and amphibians. We test hypotheses of how thermoregulatory mode, environmental temperature, protective phenotypes, and pace of life history contribute to demographic aging. Controlling for phylogeny and body size, ectotherms display a higher diversity of aging rates compared with endotherms and include phylogenetically widespread evidence of negligible aging. Protective phenotypes and life-history strategies further explain macroevolutionary patterns of aging. Analyzing ectothermic tetrapods in a comparative context enhances our understanding of the evolution of aging.

Effects of salinity and a glucocorticoid antagonist, RU486, on waterborne aldosterone and corticosterone of northern leopard frog larvae

Increased salinity is an emerging contaminant of concern for aquatic taxa. For amphibians exposed to salinity, there is scarce information about the physiological effects and changes in osmoregulatory hormones such as corticosterone (CORT) and aldosterone (ALDO). Recent studies have quantified effects of salinity on CORT physiology of amphibians based on waterborne hormone collection methods, but much less is known about ALDO in iono- and osmoregulation of amphibians. We re-assayed waterborne hormone samples from a previous study to investigate effects of salinity (sodium chloride, NaCl) and a glucocorticoid receptor antagonist (RU486) on ALDO of northern leopard frog (Rana pipiens) larvae. We also investigated relationships between ALDO and CORT. Waterborne ALDO marginally decreased with increasing salinity and was, unexpectedly, positively correlated with baseline and stress-induced waterborne CORT. Importantly, ALDO increased when larvae were exposed to RU486, suggesting that RU486 may also suppress mineralocorticoid receptors or that negative feedback of ALDO is mediated through glucocorticoid receptors. Alternatively, CORT increases with RU486 treatment and might be a substrate for ALDO synthesis, which could account for increases in ALDO with RU486 treatment and the correlation between CORT and ALDO. ALDO was negatively correlated with percent water, such that larvae secreting more ALDO retained less water. Although sample sizes were limited and further validation and studies are warranted, our findings expand our understanding of adrenal steroid responses to salinization in amphibians and proposes new hypotheses regarding the co-regulation of ALDO and CORT.

Low occurrence of ranavirus in the Prairie Pothole Region of Montana and North Dakota (USA) contrasts with prior surveys

Ranaviruses are emerging pathogens that have caused mortality events in amphibians worldwide. Despite the negative effects of ranaviruses on amphibian populations, monitoring efforts are still lacking in many areas, including in the Prairie Pothole Region (PPR) of North America. Some PPR wetlands in Montana and North Dakota (USA) have been contaminated by energy-related saline wastewaters, and increased salinity has been linked to greater severity of ranavirus infections. In 2017, we tested tissues from larvae collected at 7 wetlands that ranged in salinity from 26 to 4103 mg Cl l-1. In 2019, we used environmental DNA (eDNA) to test for ranaviruses in 30 wetlands that ranged in salinity from 26 to 11754 mg Cl l-1. A previous study (2013-2014) found that ranavirus-infected amphibians were common across North Dakota, including in some wetlands near our study area. Overall, only 1 larva tested positive for ranavirus infection, and we did not detect ranavirus in any eDNA samples. There are several potential reasons why we found so little evidence of ranaviruses, including low larval sample sizes, mismatch between sampling and disease occurrence, larger pore size of our eDNA filters, temporal variation in outbreaks, low host abundance, or low occurrence or prevalence of ranaviruses in the wetlands we sampled. We suggest future monitoring efforts be conducted to better understand the occurrence and prevalence of ranaviruses within the PPR.

Distribution of tiger salamanders in northern Sonora, Mexico: comparison of sampling methods and possible implications for an endangered subspecies

Many aquatic species in the arid USA-Mexico borderlands region are imperiled, but limited information on distributions and threats often hinders management. To provide information on the distribution of the Western Tiger Salamander (Ambystoma mavortium), including the USA-federally endangered Sonoran Tiger Salamander (Ambystoma mavortium stebbinsi), we used traditional (seines, dip-nets) and modern (environmental DNA [eDNA]) methods to sample 91 waterbodies in northern Sonora, Mexico, during 2015-2018. The endemic Sonoran Tiger Salamander is threatened by introgressive hybridization and potential replacement by another sub-species of the Western Tiger Salamander, the non-native Barred Tiger Salamander (A. m. mavortium). Based on occupancy models that accounted for imperfect detection, eDNA sampling provided a similar detection probability (0.82 [95% CI: 0.56-0.94]) as seining (0.83 [0.46-0.96]) and much higher detection than dip-netting (0.09 [0.02-0.23]). Volume of water filtered had little effect on detection, possibly because turbid sites had greater densities of salamanders. Salamanders were estimated to occur at 51 sites in 3 river drainages in Sonora. These results indicate tiger salamanders are much more widespread in northern Sonora than previously documented, perhaps aided by changes in land and water management practices. However, because the two subspecies of salamanders cannot be reliably distinguished based on morphology or eDNA methods that are based on mitochondrial DNA, we are uncertain if we detected only native genotypes or if we documented recent invasion of the area by the non-native sub-species. Thus, there is an urgent need for methods to reliably distinguish the subspecies so managers can identify appropriate interventions.

Comparative Effects of Energy-Related Saline Wastewaters and Sodium Chloride on Hatching, Survival, and Fitness-Associated Traits of Two Amphibian Species

Salinity (sodium chloride [NaCl]) is a prevalent and persistent contaminant that negatively affects freshwater ecosystems. Although most studies focus on effects of salinity from road salts (primarily NaCl), high-salinity wastewaters from energy extraction (wastewaters) could be more harmful because they contain NaCl and other toxic components. Many amphibians are sensitive to salinity, and their eggs are thought to be the most sensitive life-history stage. However, there are few investigations with salinity that include eggs and larvae sequentially in long-term exposures. We investigated the relative effects of wastewaters from a large energy reserve, the Williston Basin (USA), and NaCl on northern leopard (Rana pipiens) and boreal chorus (Pseudacris maculata) frogs. We exposed eggs and tracked responses through larval stages (for 24 days). Wastewaters and NaCl caused similar reductions in hatching and larval survival, growth, development, and activity, while also increasing deformities. Chorus frog eggs and larvae were more sensitive to salinity than leopard frogs, suggesting species-specific responses. Contrary to previous studies, eggs of both species were less sensitive to salinity than larvae. Our ecologically relevant exposures suggest that accumulating effects can reduce survival relative to starting experiments with unexposed larvae. Alternatively, egg casings of some species may provide some protection against salinity. Notably, effects of wastewaters on amphibians were predominantly due to NaCl rather than other components. Therefore, findings from studies with other sources of increased salinity (e.g., road salts) could guide management of wastewater-contaminated ecosystems, and vice versa, to mitigate effects of salinization. Environ Toxicol Chem 2021;40:3137-3147. © 2021 SETAC.

Metal accumulation varies with life history, size, and development of larval amphibians

Amphibian larvae are commonly used as indicators of aquatic ecosystem health because they are susceptible to contaminants. However, there is limited information on how species characteristics and trophic position influence contaminant loads in larval amphibians. Importantly, there remains a need to understand whether grazers (frogs and toads [anurans]) and predators (salamanders) provide comparable information on contaminant accumulation or if they are each indicative of unique environmental processes and risks. To better understand the role of trophic position in contaminant accumulation, we analyzed composite tissues for 10 metals from larvae of multiple co-occurring anuran and salamander species from 20 wetlands across the United States. We examined how metal concentrations varied with body size (anurans and salamanders) and developmental stage (anurans) and how the digestive tract (gut) influenced observed metal concentrations. Across all wetlands, metal concentrations were greater in anurans than salamanders for all metals tested except mercury (Hg), selenium (Se), and zinc (Zn). Concentrations of individual metals in anurans decreased with increasing weight and developmental stage. In salamanders, metal concentrations were less correlated with weight, indicating diet played a role in contaminant accumulation. Based on batches of similarly sized whole-body larvae compared to larvae with their digestive tracts removed, our results indicated that tissue type strongly affected perceived concentrations, especially for anurans (gut represented an estimated 46-97% of all metals except Se and Zn). This suggests the reliability of results based on whole-body sampling could be biased by metal, larval size, and development. Overall, our data shows that metal concentrations differs between anurans and salamanders, which suggests that metal accumulation is unique to feeding behavior and potentially trophic position. To truly characterize exposure risk in wetlands, species of different life histories, sizes and developmental stages should be included in biomonitoring efforts.

Corticosterone mediates a growth-survival tradeoff for an amphibian exposed to increased salinity

Life-history tradeoffs are common across taxa, but growth-survival tradeoffs-usually enhancing survival at a cost to growth-are less frequently investigated. Increased salinity (NaCl) is a prevalent anthropogenic disturbance that may cause a growth-survival tradeoff for larval amphibians. Although physiological mechanisms mediating tradeoffs are seldom investigated, hormones are prime candidates. Corticosterone (CORT) is a steroid hormone that independently influences survival and growth and may provide mechanistic insight into growth-survival tradeoffs. We conducted a 24-day experiment to test effects of salinity (<32-4000 mg/L) on growth, development, survival, CORT responses, and tradeoffs among traits of larval Northern Leopard Frogs (Rana pipiens). We also experimentally suppressed CORT signaling to determine whether CORT signaling mediates effects of salinity and a growth-survival tradeoff. Increased salinity reduced survival, growth, and development. Suppressing CORT signaling in conjunction with salinity reduced survival further but also attenuated the negative effects of salinity on growth, development, and water content. CORT of control larvae increased or was stable with growth and development but decreased with growth and development for those exposed to salinity. Therefore, salinity dysregulated CORT physiology. Across all treatments, larvae that survived had higher CORT than larvae that died. By manipulating CORT signaling, we provide strong evidence that CORT physiology mediates the outcome of a growth-survival tradeoff and enhances survival. To our knowledge, this is the first study to concomitantly measure tradeoffs between growth and survival and experimentally link these changes to CORT physiology. Identifying mechanistic links between stressors and fitness-related outcomes is critical to enhance our understanding of tradeoffs.

Evaluating corticosterone as a biomarker for amphibians exposed to increased salinity and ambient corticosterone

Physiological biomarkers are commonly used to assess the health of taxa exposed to natural and anthropogenic stressors. Glucocorticoid (GC) hormones are often used as indicators of physiological stress in wildlife because they affect growth, reproduction and survival. Increased salinity from human activities negatively influences amphibians and their corticosterone (CORT; the main amphibian GC) physiology; therefore, CORT could be a useful biomarker. We evaluated whether waterborne CORT could serve as a biomarker of salt stress for three free-living amphibian species that vary in their sensitivity to salinity: boreal chorus frogs (Pseudacris maculata), northern leopard frogs (Rana pipiens) and barred tiger salamanders (Ambystoma mavortium). Across a gradient of contamination from energy-related saline wastewaters, we tested the effects of salinity on baseline and stress-induced waterborne CORT of larvae. Stress-induced, but not baseline, CORT of leopard frogs increased with increasing salinity. Salinity was not associated with baseline or stress-induced CORT of chorus frogs or tiger salamanders. Associations between CORT and salinity were also not related to species-specific sensitivities to salinity. However, we detected background environmental CORT (ambient CORT) in all wetlands and spatial variation was high within and among wetlands. Higher ambient CORT was associated with lower waterborne CORT of larvae in wetlands. Therefore, ambient CORT likely confounded associations between waterborne CORT and salinity in our analysis and possibly influenced physiology of larvae. We hypothesize that larvae may passively take up CORT from their environment and downregulate endogenous CORT. Although effects of some hormones (e.g. oestrogen) and endocrine disruptors on aquatic organisms are well described, studies investigating the occurrence and effects of ambient CORT are limited. We provide suggestions to improve collection methods, reduce variability and avoid confounding effects of ambient CORT. By making changes to methodology, waterborne CORT could still be a promising, non-invasive conservation tool to evaluate effects of salinity on amphibians.

Genomic signatures of thermal adaptation are associated with clinal shifts of life history in a broadly distributed frog

Temperature is a critical driver of ectotherm life‐history strategies, whereby a warmer environment is associated with increased growth, reduced longevity and accelerated senescence. Increasing evidence indicates that thermal adaptation may underlie such life‐history shifts in wild populations. Single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) can help uncover the molecular mechanisms of temperature‐driven variation in growth, longevity and senescence. However, our understanding of these mechanisms is still limited, which reduces our ability to predict the response of non‐model ectotherms to global temperature change.

In this study, we examined the potential role of thermal adaptation in clinal shifts of life‐history traits (i.e. life span, senescence rate and recruitment) in the Columbia spotted frog Rana luteiventris along a broad temperature gradient in the western United States.

We took advantage of extensive capture–recapture datasets of 20,033 marked individuals from eight populations surveyed annually for 14–18 years to examine how mean annual temperature and precipitation influenced demographic parameters (i.e. adult survival, life span, senescence rate, recruitment and population growth). After showing that temperature was the main climatic predictor influencing demography, we used RAD‐seq data (50,829 SNPs and 6,599 putative CNVs) generated for 352 individuals from 31 breeding sites to identify the genomic signatures of thermal adaptation.

Our results showed that temperature was negatively associated with annual adult survival and reproductive life span and positively associated with senescence rate. By contrast, recruitment increased with temperature, promoting the long‐term viability of most populations. These temperature‐dependent demographic changes were associated with strong genomic signatures of thermal adaptation. We identified 148 SNP candidates associated with temperature including three SNPs located within protein‐coding genes regulating resistance to cold and hypoxia, immunity and reproduction in ranids. We also identified 39 CNV candidates (including within 38 transposable elements) for which normalized read depth was associated with temperature.

Our study indicates that both SNPs and structural variants are associated with temperature and could eventually be found to play a functional role in clinal shifts in senescence rate and life‐history strategies in R. luteiventris. These results highlight the potential role of different sources of molecular variation in the response of ectotherms to environmental temperature variation in the context of global warming.

Accommodating the role of site memory in dynamic species distribution models

First-order dynamic occupancy models (FODOMs) are a class of state-space model in which the true state (occurrence) is observed imperfectly. An important assumption of FODOMs is that site dynamics only depend on the current state and that variations in dynamic processes are adequately captured with covariates or random effects. However, it is often difficult to understand and/or measure the covariates that generate ecological data, which are typically spatiotemporally correlated. Consequently, the non-independent error structure of correlated data causes underestimation of parameter uncertainty and poor ecological inference. Here, we extend the FODOM framework with a second-order Markov process to accommodate site memory when covariates are not available. Our modeling framework can be used to make reliable inference about site occupancy, colonization, extinction, turnover, and detection probabilities. We present a series of simulations to illustrate the data requirements and model performance. We then applied our modeling framework to 13 yr of data from an amphibian community in southern Arizona, USA. In this analysis, we found residual temporal autocorrelation of population processes for most species, even after accounting for long-term drought dynamics. Our approach represents a valuable advance in obtaining inference on population dynamics, especially as they relate to metapopulations.

Relative Toxicity and Sublethal Effects of NaCl and Energy-Related Saline Wastewaters on Prairie Amphibians

Increasing salinity in freshwater environments is a growing problem due both to the negative influences of salts on ecosystems and their accumulation and persistence in environments. Two major sources of increased salinity from sodium chloride salts (NaCl) are saline wastewaters co-produced during energy production (herein, wastewaters) and road salts. Effects of road salts have received more attention, but legacy contamination from wastewaters is widespread in some regions and spills still occur. Amphibians are sensitive to contaminants, including NaCl, because of their porous skin and osmoregulatory adaptations to freshwater. However, similarities and differences between effects of wastewaters and road salts have not been investigated. Therefore, we investigated the relative influence of wastewaters and NaCl at equivalent concentrations of chloride on three larval amphibian species that occur in areas with increased salinity. We determined acute toxicity and growth effects on Boreal Chorus Frogs (Pseudacris maculata), Northern Leopard Frogs (Rana pipiens), and Barred Tiger Salamanders (Ambystoma mavortium). We posited that wastewaters would have additive effects on amphibians compared to NaCl because wastewaters often have additional toxic heavy metals and other contaminants. For NaCl, toxicity was higher for frogs than the salamander. Toxicity of wastewaters was also similar between chorus and leopard frogs. Only chorus frog survival was lower when exposed to wastewater compared to NaCl. Mass and length of leopard and chorus frog larvae decreased with increasing salinity after only 96 hours of exposure but did not for tiger salamanders. Size of leopard frogs was lower when exposed to NaCl compared to wastewater. However, growth effects were similar between wastewater and NaCl for chorus frogs. Taken together, our results suggest that previous studies on effects of road salt could inform future studies and management of wastewater-contaminated ecosystems, and vice versa. Nevertheless, effects of road salts and wastewaters may be context-, species-, and trait-specific and require further investigations. The negative influence of salts on imperiled amphibians underscores the need to restore landscapes with increased salinity and reduce future salinization of freshwater ecosystems.

Survival estimates for the invasive American bullfrog

American bullfrogs (Lithobates catesbeianus) are significant invaders in many places and can negatively impact native species. Despite their impact and wide distribution, little is known about their demography. We used five years of capture mark-recapture data to estimate annual apparent survival of post-metamorphic bullfrogs in a population on the Buenos Aires National Wildlife Refuge in their invaded range in Arizona, U.S.A. This population is a potential source of colonists into breeding ponds used by the federally threatened Chiricahua leopard frog (L. chiricahuensis). Results from robust-design Cormack-Jolly-Seber models suggested that survival of bullfrogs was influenced by sex and precipitation but not body condition. Survival was higher for females (mean = 0.37; 95% , 0.72) than males (mean = 0.17; 95% , 0.49), and declined with reduced annual precipitation (mean = −0.36, 95%  = −2.09, 0.84). These survival estimates can be incorporated into models of population dynamics and to help predict spread of bullfrogs.

A statistical forecasting approach to metapopulation viability analysis

Conservation of at-risk species is aided by reliable forecasts of the consequences of environmental change and management actions on population viability. Forecasts from conventional population viability analysis (PVA) are made using a two-step procedure in which parameters are estimated, or elicited from expert opinion, and then plugged into a stochastic population model without accounting for parameter uncertainty. Recently developed statistical PVAs differ because forecasts are made conditional on models fitted to empirical data. The statistical forecasting approach allows for uncertainty about parameters, but it has rarely been applied in metapopulation contexts where spatially explicit inference is needed about colonization and extinction dynamics and other forms of stochasticity that influence metapopulation viability. We conducted a statistical metapopulation viability analysis (MPVA) using 11 yr of data on the federally threatened Chiricahua leopard frog (Lithobates chiricahuensis) to forecast responses to landscape heterogeneity, drought, environmental stochasticity, and management. We evaluated several future environmental scenarios and pond restoration options designed to reduce extinction risk. Forecasts over a 50-yr time horizon indicated that metapopulation extinction risk was <4% for all scenarios, but uncertainty was high. Without pond restoration, extinction risk is forecasted to be 3.9% (95% CI 0-37%) by year 2066. Restoring six ponds by increasing their hydroperiod reduced extinction risk to <1% and greatly reduced uncertainty (95% CI 0-2%). Our results suggest that managers can mitigate the impacts of drought and environmental stochasticity on metapopulation viability by maintaining ponds that hold water throughout the year and keeping them free of invasive predators. Our study illustrates the utility of the spatially explicit statistical forecasting approach to MPVA in conservation planning efforts.

Amphibian chytrid prevalence on boreal toads in SE Alaska and NW British Columbia: tests of habitat, life stages, and temporal trends

Tracking and understanding variation in pathogens such as Batrachochytrium dendrobatidis (Bd), the agent of amphibian chytridiomycosis which has caused population declines globally, is a priority for many land managers. However, relatively little sampling of amphibian communities has occurred at high latitudes. We used skin swabs collected during 2005-2017 from boreal toads Anaxyrus boreas (n = 248), in southeast Alaska (USA; primarily in and near Klondike Gold Rush National Historical Park [KLGO]) and northwest British Columbia (Canada) to determine how Bd prevalence varied across life stages, habitat characteristics, local species richness, and time. Across all years, Bd prevalence peaked in June and was >3 times greater for adult toads (37.5%) vs. juveniles and metamorphs (11.2%). Bd prevalence for toads in the KLGO area, where other amphibian species are rare or absent, was highest from river habitats (55.0%), followed by human-modified upland wetlands (32.3%) and natural upland wetlands (12.7%)-the same rank-order these habitats are used for toad breeding. None of the 12 Columbia spotted frogs Rana luteiventris or 2 wood frogs R. sylvatica from the study area tested Bd-positive, although all were from an area of low host density where Bd has not been detected. Prevalence of Bd on toads in the KLGO area decreased during 2005-2015. This trend from a largely single-species system may be encouraging or concerning, depending on how Bd is affecting vital rates, and emphasizes the need to understand effects of pathogens before translating disease prevalence into management actions.

Predicting attenuation of salinized surface- and groundwater-resources from legacy energy development in the Prairie Pothole Region

Oil and gas (energy) development in the Williston Basin, which partly underlies the Prairie Pothole Region in central North America, has helped meet U.S. energy demand for decades. Historical handling and disposal practices of saline wastewater co-produced during energy development resulted in salinization of surface and groundwater at numerous legacy energy sites. Thirty years of monitoring (1988-2018) at Goose Lake, which has been producing since the 1960s, documents long-term spatial and temporal changes in water quality from legacy energy development. Surface water quality was highly variable and decoupled from changes in groundwater quality, likely due to annual and regional climatic fluctuations. Therefore, changes in surface water-quality were not considered a reliable indicator of subsurface chloride migration. However, chloride concentrations in monitoring wells near wastewater sources exhibited systematic temporal reductions allowing for estimates of the time required for natural attenuation of groundwater to U.S. Environmental Protection Agency acute and chronic chloride toxicity benchmarks and a local background level. Point attenuation rates differed based on sediment type (outwash vs till) and yielded a range of predicted years when water-quality targets will be reached: acute - 2045 to 2113; chronic - 2069 to 2160; background - 2126 to 2275. Bulk attenuation rates from four separate years of data were used to calculate the distances chloride could migrate downgradient from the largest wastewater source. Potential distances of downgradient migration before dilution to water-quality targets decreased from 1989 to 2018: acute - 949 to 673 m; chronic - 1220 to 922 m; background - 1878 to 1525 m. Several downgradient wetlands are within these distances and will continue to receive saline contaminated groundwater for years. While these results demonstrate chloride attenuation at a legacy energy site, they also highlight the persistence of saline wastewater contamination and the need to mitigate future spills to prevent long-term salinization from energy development.

Spatial capture-recapture reveals age- and sex-specific survival and movement in stream amphibians

Life-history information sets the foundation for our understanding of ecology and conservation requirements. For many species, this information is lacking even for basic demographic rates such as survival and movement. When survival and movement estimates are available, they are often derived from mixed demographic groups and do not consider differences among life stages or sexes, which is critical, because life stages and sexes often contribute differentially to population dynamics. We used hierarchical models informed with spatial capture-mark-recapture data of Ascaphus montanus (Rocky Mountain tailed frog) in five streams and A. truei (coastal tailed frog) in one stream to estimate variation in survival and movement by sex and age, represented by size. By incorporating survival and movement into a single model, we were able to estimate both parameters with limited bias. Annual survival was similar between sexes of A. montanus [females = 0.885 (95% CI 0.614-1), males = 0.901 (0.657-1)], but was slightly higher for female A. truei [0.836 (0.560-0.993)] than for males [0.664 (0.354-0.962)]. Survival of A. montanus peaked at mid-age, suggesting that lower survival of young and actuarial senescence may influence population demographics. Our models suggest that younger A. montanus moved farther than older individuals, and that females moved farther than males in both species. Our results provide uncommon insight into age- and sex-specific rates of survival and movement that are crucial elements of life-history strategies and are important for modeling population growth and prescribing conservation actions.

Associations between environmental pollutants and larval amphibians in wetlands contaminated by energy-related brines are potentially mediated by feeding traits

Energy production in the Williston Basin, located in the Prairie Pothole Region of central North America, has increased rapidly over the last several decades. Advances in recycling and disposal practices of saline wastewaters (brines) co-produced during energy production have reduced ecological risks, but spills still occur often and legacy practices of releasing brines into the environment caused persistent salinization in many areas. Aside from sodium and chloride, these brines contain elevated concentrations of metals and metalloids (lead, selenium, strontium, antimony and vanadium), ammonium, volatile organic compounds, hydrocarbons, and radionuclides. Amphibians are especially sensitive to chloride and some metals, increasing potential effects in wetlands contaminated by brines. We collected bed sediment and larval amphibians (Ambystoma mavortium, Lithobates pipiens and Pseudacris maculata) from wetlands in Montana and North Dakota representing a range of brine contamination history and severity to determine if contamination was associated with metal concentrations in sediments and if metal accumulation in tissues varied by species. In wetland sediments, brine contamination was positively associated with the concentrations of sodium and strontium, both known to occur in oil and gas wastewater, but negatively correlated with mercury. In amphibian tissues, selenium and vanadium were associated with brine contamination. Metal tissue concentrations were higher in tadpoles that graze compared to predatory salamanders; this suggests frequent contact with the sediments could lead to greater ingestion of metal-laden materials. Although many of these metals may not be directly linked with energy development, the potential additive or synergistic effects of exposure along with elevated chloride from brines could have important consequences for aquatic organisms. To effectively manage amphibian populations in wetlands contaminated by saline wastewaters we need a better understanding of how life history traits, species-specific susceptibilities and the physical-chemical properties of metals co-occurring in wetland sediments interact with other stressors like chloride and wetland drying.

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.

Increasing connectivity between metapopulation ecology and landscape ecology

Metapopulation ecology and landscape ecology aim to understand how spatial structure influences ecological processes, yet these disciplines address the problem using fundamentally different modeling approaches. Metapopulation models describe how the spatial distribution of patches affects colonization and extinction, but often do not account for the heterogeneity in the landscape between patches. Models in landscape ecology use detailed descriptions of landscape structure, but often without considering colonization and extinction dynamics. We present a novel spatially explicit modeling framework for narrowing the divide between these disciplines to advance understanding of the effects of landscape structure on metapopulation dynamics. Unlike previous efforts, this framework allows for statistical inference on landscape resistance to colonization using empirical data. We demonstrate the approach using 11 yr of data on a threatened amphibian in a desert ecosystem. Occupancy data for Lithobates chiricahuensis (Chiricahua leopard frog) were collected on the Buenos Aires National Wildlife Refuge (BANWR), Arizona, USA from 2007 to 2017 following a reintroduction in 2003. Results indicated that colonization dynamics were influenced by both patch characteristics and landscape structure. Landscape resistance increased with increasing elevation and distance to the nearest streambed. Colonization rate was also influenced by patch quality, with semi-permanent and permanent ponds contributing substantially more to the colonization of neighboring ponds relative to intermittent ponds. Ponds that only hold water intermittently also had the highest extinction rate. Our modeling framework can be widely applied to understand metapopulation dynamics in complex landscapes, particularly in systems in which the environment between habitat patches influences the colonization process.

Heterogeneous responses of temperate-zone amphibian populations to climate change complicates conservation planning

The pervasive and unabated nature of global amphibian declines suggests common demographic responses to a given driver, and quantification of major drivers and responses could inform broad-scale conservation actions. We explored the influence of climate on demographic parameters (i.e., changes in the probabilities of survival and recruitment) using 31 datasets from temperate zone amphibian populations (North America and Europe) with more than a decade of observations each. There was evidence for an influence of climate on population demographic rates, but the direction and magnitude of responses to climate drivers was highly variable among taxa and among populations within taxa. These results reveal that climate drivers interact with variation in life-history traits and population-specific attributes resulting in a diversity of responses. This heterogeneity complicates the identification of conservation 'rules of thumb' for these taxa, and supports the notion of local focus as the most effective approach to overcome global-scale conservation challenges.

Widespread legacy brine contamination from oil production reduces survival of chorus frog larvae

Advances in drilling techniques have facilitated a rapid increase in hydrocarbon extraction from energy shales, including the Williston Basin in central North America. This area overlaps with the Prairie Pothole Region, a region densely populated with wetlands that provide numerous ecosystem services. Historical (legacy) disposal practices often released saline co-produced waters (brines) with high chloride concentrations, affecting wetland water quality directly or persisting in sediments. Despite the potential threat of brine contamination to aquatic habitats, there has been little research into its ecological effects. We capitalized on a gradient of legacy brine-contaminated wetlands in northeast Montana to conduct laboratory experiments to assess variation in survival of larval Boreal Chorus Frogs (Pseudacris maculata) reared on sediments from 3 local wetlands and a control source. To help provide environmental context for the experiment, we also measured chloride concentrations in 6 brine-contaminated wetlands in our study area, including the 2 contaminated sites used for sediment exposures. Survival of frog larvae during 46- and 55-day experiments differed by up to 88% among sediment sources (Site Model) and was negatively correlated with potential chloride exposure (Chloride Model). Five of the 6 contaminated wetlands exceeded the U.S. EPA acute benchmark for chloride in freshwater (860 mg/L) and all exceeded the chronic benchmark (230 mg/L). However, the Wetland Site model explained more variation in survival than the Chloride Model, suggesting that chloride concentration alone does not fully reflect the threat of contamination to aquatic species. Because the profiles of brine-contaminated sediments are complex, further surveys and experiments are needed across a broad range of conditions, especially where restoration or remediation actions have reduced brine-contamination. Information provided by this study can help quantify potential ecological threats and help land managers prioritize conservation strategies as part of responsible and sustainable energy development.

Additive impacts of experimental climate change increase risk to an ectotherm at the Arctic's edge

Globally, Arctic and Subarctic regions have experienced the greatest temperature increases during the last 30 years. These extreme changes have amplified threats to the freshwater ecosystems that dominate the landscape in many areas by altering water budgets. Several studies in temperate environments have examined the adaptive capacity of organisms to enhance our understanding of the potential repercussions of warming and associated accelerated drying for freshwater ecosystems. However, few experiments have examined these impacts in Arctic or Subarctic freshwater ecosystems, where the climate is changing most rapidly. To evaluate the capacity of a widespread ectotherm to anticipated environmental changes, we conducted a mesocosm experiment with wood frogs (Rana sylvatica) in the Canadian Subarctic. Three warming treatments were fully crossed with three drying treatments to simulate a range of predicted changes in wetland environments. We predicted wetland warming and drying would act synergistically, with water temperature partially compensating for some of the negative effects of accelerated drying. Across all drying regimes, a 1 °C increase in water temperature increased the odds of survival by 1.79, and tadpoles in 52-day and 64-day hydroperiod mesocosms were 4.1-4.3 times more likely to survive to metamorphosis than tadpoles in 45-day mesocosms. For individuals who survived to metamorphosis, there was only a weak negative effect of temperature on size. As expected, increased temperatures accelerated tadpole growth through day 30 of the experiment. Our results reveal that one of the dominant herbivores in Subarctic wetlands, wood frog tadpoles, are capable of increasing their developmental rates in response to increased temperature and accelerated drying, but only in an additive manner. The strong negative effects of drying on survival, combined with lack of compensation between these two environmental drivers, suggest changes in the aquatic environment that are expected in this ecosystem will reduce mean fitness of populations across the landscape.

Roles of patch characteristics, drought frequency, and restoration in long-term trends of a widespread amphibian

Despite the high profile of amphibian declines and the increasing threat of drought and fragmentation to aquatic ecosystems, few studies have examined long-term rates of change for a single species across a large geographic area. We analyzed growth in annual egg-mass counts of the Columbia spotted frog (Rana luteiventris) across the northwestern United States, an area encompassing 3 genetic clades. On the basis of data collected by multiple partners from 98 water bodies between 1991 and 2011, we used state-space and linear-regression models to measure effects of patch characteristics, frequency of summer drought, and wetland restoration on population growth. Abundance increased in the 2 clades with greatest decline history, but declined where populations are considered most secure. Population growth was negatively associated with temporary hydroperiods and landscape modification (measured by the human footprint index), but was similar in modified and natural water bodies. The effect of drought was mediated by the size of the water body: populations in large water bodies maintained positive growth despite drought, whereas drought magnified declines in small water bodies. Rapid growth in restored wetlands in areas of historical population declines provided strong evidence of successful management. Our results highlight the importance of maintaining large areas of habitat and underscore the greater vulnerability of small areas of habitat to environmental stochasticity. Similar long-term growth rates in modified and natural water bodies and rapid, positive responses to restoration suggest pond construction and other forms of management can effectively increase population growth. These tools are likely to become increasingly important to mitigate effects of increased drought expected from global climate change. Papeles de las Características del Fragmento, Frecuencia de Sequía y Restauración en las Tendencias a Largo Plazo de un Anfibio Ampliamente Distribuido.

Interactive effects of wildfire, forest management, and isolation on amphibian and parasite abundance

Projected increases in wildfire and other climate-driven disturbances will affect populations and communities worldwide, including host-parasite relationships. Research in temperate forests has shown that wildfire can negatively affect amphibians, but this research has occurred primarily outside of managed landscapes where interactions with human disturbances could result in additive or synergistic effects. Furthermore, parasites represent a large component of biodiversity and can affect host fitness and population dynamics, yet they are rarely included in studies of how vertebrate hosts respond to disturbance. To determine how wildfire affects amphibians and their parasites, and whether effects differ between protected and managed landscapes, we compared abundance of two amphibians and two nematodes relative to wildfire extent and severity around wetlands in neighboring protected and managed forests (Montana, USA). Population sizes of adult, male long-toed salamanders (Ambystoma macrodactylum) decreased with increased burn severity, with stronger negative effects on isolated populations and in managed forests. In contrast, breeding population sizes of Columbia spotted frogs (Rana luteiventris) increased with burn extent in both protected and managed protected forests. Path analysis showed that the effects of wildfire on the two species of nematodes were consistent with differences in their life history and transmission strategies and the responses of their hosts. Burn severity indirectly reduced abundance of soil-transmitted Cosmocercoides variabilis through reductions in salamander abundance. Burn severity also directly reduced C. variabilis abundance, possibly though changes in soil conditions. For the aquatically transmitted nematode Gyrinicola batrachiensis, the positive effect of burn extent on density of Columbia spotted frog larvae indirectly increased parasite abundance. Our results show that effects of wildfire on amphibians depend upon burn extent and severity, isolation, and prior land use. Through subsequent effects on the parasites, our results also reveal how changes in disturbance regimes can affect communities across trophic levels.

Rapid increases and time-lagged declines in amphibian occupancy after wildfire

Climate change is expected to increase the frequency and severity of drought and wildfire. Aquatic and moisture-sensitive species, such as amphibians, may be particularly vulnerable to these modified disturbance regimes because large wildfires often occur during extended droughts and thus may compound environmental threats. However, understanding of the effects of wildfires on amphibians in forests with long fire-return intervals is limited. Numerous stand-replacing wildfires have occurred since 1988 in Glacier National Park (Montana, U.S.A.), where we have conducted long-term monitoring of amphibians. We measured responses of 3 amphibian species to fires of different sizes, severity, and age in a small geographic area with uniform management. We used data from wetlands associated with 6 wildfires that burned between 1988 and 2003 to evaluate whether burn extent and severity and interactions between wildfire and wetland isolation affected the distribution of breeding populations. We measured responses with models that accounted for imperfect detection to estimate occupancy during prefire (0-4 years) and different postfire recovery periods. For the long-toed salamander (Ambystoma macrodactylum) and Columbia spotted frog (Rana luteiventris), occupancy was not affected for 6 years after wildfire. But 7-21 years after wildfire, occupancy for both species decreased ≥ 25% in areas where >50% of the forest within 500 m of wetlands burned. In contrast, occupancy of the boreal toad (Anaxyrus boreas) tripled in the 3 years after low-elevation forests burned. This increase in occupancy was followed by a gradual decline. Our results show that accounting for magnitude of change and time lags is critical to understanding population dynamics of amphibians after large disturbances. Our results also inform understanding of the potential threat of increases in wildfire frequency or severity to amphibians in the region.

Effects of amphibian chytrid fungus on individual survival probability in wild boreal toads

Chytridiomycosis is linked to the worldwide decline of amphibians, yet little is known about the demographic effects of the disease. We collected capture-recapture data on three populations of boreal toads (Bufo boreas [Bufo = Anaxyrus]) in the Rocky Mountains (U.S.A.). Two of the populations were infected with chytridiomycosis and one was not. We examined the effect of the presence of amphibian chytrid fungus (Batrachochytrium dendrobatidis [Bd]; the agent of chytridiomycosis) on survival probability and population growth rate. Toads that were infected with Bd had lower average annual survival probability than uninfected individuals at sites where Bd was detected, which suggests chytridiomycosis may reduce survival by 31-42% in wild boreal toads. Toads that were negative for Bd at infected sites had survival probabilities comparable to toads at the uninfected site. Evidence that environmental covariates (particularly cold temperatures during the breeding season) influenced toad survival was weak. The number of individuals in diseased populations declined by 5-7%/year over the 6 years of the study, whereas the uninfected population had comparatively stable population growth. Our data suggest that the presence of Bd in these toad populations is not causing rapid population declines. Rather, chytridiomycosis appears to be functioning as a low-level, chronic disease whereby some infected individuals survive but the overall population effects are still negative. Our results show that some amphibian populations may be coexisting with Bd and highlight the importance of quantitative assessments of survival in diseased animal populations.

Responses of pond-breeding amphibians to wildfire: short-term patterns in occupancy and colonization

Wildland fires are expected to become more frequent and severe in many ecosystems, potentially posing a threat to many sensitive species. We evaluated the effects of a large, stand-replacement wildfire on three species of pond-breeding amphibians by estimating changes in occupancy of breeding sites during the three years before and after the fire burned 42 of 83 previously surveyed wetlands. Annual occupancy and colonization for each species was estimated using recently developed models that incorporate detection probabilities to provide unbiased parameter estimates. We did not find negative effects of the fire on the occupancy or colonization rates of the long-toed salamander (Ambystoma macrodactylum). Instead, its occupancy was higher across the study area after the fire, possibly in response to a large snowpack that may have facilitated colonization of unoccupied wetlands. Naive data (uncorrected for detection probability) for the Columbia spotted frog (Rana luteiventris) initially led to the conclusion of increased occupancy and colonization in wetlands that burned. After accounting for temporal and spatial variation in detection probabilities, however, it was evident that these parameters were relatively stable in both areas before and after the fire. We found a similar discrepancy between naive and estimated occupancy of A. macrodactylum that resulted from different detection probabilities in burned and control wetlands. The boreal toad (Bufo boreas) was not found breeding in the area prior to the fire but colonized several wetlands the year after they burned. Occupancy by B. boreas then declined during years 2 and 3 following the fire. Our study suggests that the amphibian populations we studied are resistant to wildfire and that B. boreas may experience short-term benefits from wildfire. Our data also illustrate how naive presence-non-detection data can provide misleading results.