Search ARMI Database

Comparative studies of mortality in the wild are necessary to understand the evolution of aging, yet ectothermic tetrapods are under-represented in this comparative landscape despite their suitability for testing evolutionary hypotheses. We provide the first comprehensive study of aging rates and longevity across tetrapod ectotherms in the wild, utilizing data from 107 populations (77 species) of reptiles and amphibians. We test hypotheses of how thermoregulatory mode, temperature, protective phenotypes, and pace of life contribute to aging. Controlling for phylogeny and body size, ectotherms displayed a higher diversity of aging rates than endotherms, and included many groups with negligible aging. Furthermore, protective phenotypes and life-history tactics further explained macroevolutionary patterns of aging. By including ectothermic tetrapods, our comparative analyses enhance our understanding of aging evolution.

Increased salinity (sodium chloride; NaCl) is a prevalent and persistent
contaminant that negatively affects freshwater ecosystems. Although
most studies focus on effects of salinity from roads 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 to salinity and tracked
responses through larval stages (for 24 days). Wastewaters and NaCl
reduced 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.

Wetland creation is a common practice to mitigate for the loss of natural wetlands. However, there is still uncertainty about how effectively created wetlands replace habitat provided by natural wetlands. This uncertainty is due in part because post-construction monitoring of biological communities, and vertebrates especially, is rare and typically short-term (< 5 years). We estimated occupancy of 4 amphibian species in 8 created mitigation wetlands, 7 impacted wetlands, and 7 reference wetlands in the Greater Yellowstone Ecosystem in Wyoming, USA. Mitigation wetlands were created to replace wetland habitat that was lost during road construction and ranged in age from 1 to 10 years when sampled. Impacted wetlands were natural wetlands partially filled by road construction and were adjacent to a highway. We sampled for amphibian larvae during 6 summers from 2013 to 2020 and used multi-species occupancy models that estimated detection and occupancy of each of 4 amphibian species to determine how amphibian responses changed over time, especially in mitigation wetlands. Occupancy did not differ between impacted and reference wetlands for any of the 4 amphibian species. Western Toads (Anaxyrus boreas) were most common (although briefly) in created wetlands, and occupancy of Columbia Spotted Frogs (Rana luteiventris), Western Tiger Salamanders (Ambystoma mavortium), and Boreal Chorus Frogs (Pseudacris maculata) was lower in created wetlands than in impacted or reference wetlands. Wetland area was positively associated with occupancy for all 4 species and wetland vegetation cover was positively associated with Boreal Chorus Frog and Columbia Spotted Frog occupancy; these results emphasize the importance of design characteristics when planning mitigation wetlands. The link between wetland age and occupancy was complex and included threshold and quadratic relationships for three of the four species, but only Boreal Chorus Frog occupancy was still increasing slowly at the end of our study. Our results indicate created wetlands did not attain the suitability of impacted and natural wetlands for local amphibians, even several years after construction. The complicated relationships between wetland age and species-specific occupancy illustrate the importance of long-term monitoring in describing population responses to the construction of wetlands as mitigation for wetland loss.

Life-history tradeoffs are common across taxa, but growth-survival tradeoffs—which usually enhance 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 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 used 24-d trials to test effects of salinity (0 – 4000 mg/L Cl-) on growth, development, survival, CORT responses, and tradeoffs among traits of larval Northern Leopard Frogs (Rana pipiens). We also tested experimentally suppressed CORT signaling to determine whether CORT signaling mediates effects of salinity and life history tradeoffs . Increased salinity reduced survival, growth, and development. Suppressing CORT signaling in conjunction with salinity reduced survival further, but also attenuated negative effects of salinity on growth and development. 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. 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.

Salinization is harmful to amphibians and waterborne corticosterone could be a useful biomarker. Salinity was only associated with waterborne corticosterone for one of three amphibian species. Ambient corticosterone likely confounded associations and possibly influenced amphibian physiology. We provide suggestions to improve reliability of waterborne corticosterone as a biomarker of salt stress.

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 (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, selenium (Se), and zinc (Zn). Concentrations of individual metals in anurans decreased with increasing weight and developmental stage. In salamanders, which are predatory, 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 50–90% of all metals except Se and Zn). This suggest the reliability of results based on whole body sampling could be biased by metal, larval size, and development. Overall, our data suggest that metal concentrations differs among orders (anuran 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.

For species with variable phenology, it is often challenging to produce reliable estimates of population dynamics or changes in occupancy. The Arizona Toad (Anaxyrus microscaphus) is a southwestern USA endemic that has been petitioned for legal protection, but status assessments are limited by a lack of information on population trends. Also, timing and consistency of Arizona Toad breeding varies greatly, making it difficult to predict optimal survey times or effort required for detection. To help fill these information gaps, we conducted breeding season call surveys during 2013–2016 and 2019 at 86 historically occupied sites and 59 control sites across the species’ range in New Mexico. We estimated variation in mean dates of arrival and departure from breeding sites, changes in occupancy, and site-level extinction since 1959 with recently developed multi-season staggered-entry models, which relax the within-season closure assumption common to most occupancy models. Optimal timing of surveys in our study areas was approximately March 5 - March 30. Averaged across years, estimated probability of occupancy was https://0.58 (SE = 0.09) for historical sites and https://0.19 (SE = 0.08) for control sites. Occupancy increased from 2013 through 2019. Notably, even though observer error was trivial, annual detection probabilities varied from https://0.23 to https://0.75 and declined during the study; this means naïve occupancy values would have been misleading, indicating apparent declines in toad occupancy. Occupancy was lowest during the first year of the study, possibly due to changes in stream flows and conditions in many waterbodies following extended drought and recent wildfires. Although within-season closure was violated by variable calling phenology, simple multi-season models provided nearly identical estimates as staggered-entry models. Surprisingly, extinction probability was unrelated to the number of years since the first or last record at historically occupied sites. Collectively, our results suggest a lack of large, recent declines in occupancy by Arizona Toads in New Mexico, but we still lack population information from most of the species’ range.

Variation in temperature is known to influence mortality patterns in ectotherms. Even though a few experimental studies on model organisms have reported a positive relationship between temperature and actuarial senescence (i.e., the increase in mortality risk with age), how variation in climate influences senescence rate across the range of a species is still poorly understood in free-ranging animals. We filled this knowledge gap by investigating the relationships linking senescence rate, adult lifespan, and climatic conditions using long-term, capture-recapture data from multiple amphibian populations. We considered two pairs of related anuran species from the Ranidae (Rana luteiventris, Rana temporaria) and Bufonidae (Anaxyrus boreas, Bufo bufo) families, which diverged more than 100 mya and are broadly distributed in North America and Europe. Senescence rates were positively associated with mean annual temperature in all species. In addition, lifespan was negatively correlated with mean annual temperature in all species except A. boreas. In both R. luteiventris and A. boreas, mean annual precipitation and human environmental footprint both had negligible effects on senescence rates or lifespans. Overall, our findings demonstrate the critical influence of thermal conditions on mortality patterns across anuran species from temperate regions. In the current context of further global temperature increase predicted by IPCC scenarios, a widespread acceleration of aging in amphibians is expected to occur in the decades to come, which might threaten even more seriously the viability of populations and exacerbate global decline.

During December 8–10, 2020, the NPS-Greater Yellowstone Network hosted a virtual three-day amphibian and wetland science symposium entitled Linking Monitoring and Research to Amphibian Conservation in the Greater Yellowstone Ecosystem. This symposium highlighted work related to the NPS program that has monitored amphibian and wetland dynamics in Grand Teton and Yellowstone national parks since 2006. USGS-ARMI has been linked with monitoring program since the inception of the NPS monitoring program and assisted with organizing symposium presentations.

The three-day symposium served two important purposes. The first two days were intended to bring together practitioners of amphibian monitoring and park natural resource personnel with university, non-governmental scientists, and USGS scientists. Participants shared examples from their work and the current state of knowledge about research and monitoring of amphibians and their habitats. Reflecting the long history of ARMI research in the Rocky Mountains and on NPS lands, there were five presentations by ARMI scientists, and seven presentations highlighted on-going or recently completed ARMI research.

The third day was structured to provide a formal peer review of the network’s amphibian and wetland monitoring protocol. The eight reviewers, which included ARMI scientist Brian Halstead and former ARMI PhD student Paige Howell, offered recommendations to enrich and strengthen the Greater Yellowstone Network ’s ongoing monitoring efforts. In total, more than 60 people participated in the symposium and program review.

ARMI and the Greater Yellowstone Network will continue their collaboration by co-organizing a special issue in the journal Ecological Indicators. The issue will include several papers that highlight collaborative USGS-NPS research on amphibians, wetlands, and climate, as well as emerging tools and priorities for scientists and resource managers.

North America’s protected lands harbor significant biodiversity and provide habitats where species threatened by a variety of stressors in other environments can thrive. Yet disease, climate change, and other threats are not limited by land management boundaries and can interact with conditions within protected landscapes to affect sensitive populations. We examined the population dynamics of a boreal toad (Anaxyrus boreas) metapopulation at a wildlife refuge in northwestern Montana over a 16-year period (2003-2018). We used robust design capture-recapture models to estimate male population size, recruitment, and apparent survival over time and in relation to the amphibian chytrid fungus, Batrachochytrium dendrobatidis. We estimated female population size in years with sufficient captures. Finally, we examined trends in male and female toad body size and condition. We found no evidence of an effect of disease or time on male toad survival but detected a strong negative trend in recruitment of new males to the population. Estimates of male and female abundance decreased dramatically over time. Body size of males and females was inversely related to estimated population size, consistent with reduced recruitment to replace adults, but body condition of adult males was only weakly associated with abundance. Together, these results describe the demography of a near-extinction event, and point to dramatic decreases in the recruitment of new individuals to the breeding population as the cause of this decline. We surmise that processes related to the restoration of historical hydrology within the refuge adversely affected amphibian breeding habitat, and that these changes interacted with disease, life history, and other factors to restrict the recruitment of new individuals to the breeding population over time. Our results point to challenges in understanding and predicting drivers of population change and highlight that current metrics for assessing population status can have limited predictive ability.

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 measure the covariates that generate ecological data, which are often spatio-temporally 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 years of data from an amphibian community in southern Arizona, USA and find that site memory helps describe dynamic processes for most species. Our approach represents a valuable advance in obtaining inference on population dynamics, especially as they relate to metapopulations.

Data release and code for Ecological Applications paper: A statistical forecasting approach to metapopulation viability analysis

The Buenos Aires National Wildlife Refuge (BANWR) has developed a new pond designed to replicate ciénega conditions that will provide habitat for the federally threatened Chiricahua leopard frog. Reintroduced to the Refuge after extirpation in the early 2000s, the Chiricahua leopard frog has colonized areas beyond the initial reintroduction sites, now inhabiting several human-made stock tanks on the Refuge. The species’ persistence on the Refuge depends on the ongoing eradication of invasive bullfrogs and the availability of water. In desert habitats such as BANWR, the ciénega (a Spanish word for wetland) is a marsh-like habitat that is critical for numerous desert-dwelling creatures. Ciénegas on BANWR are characterized by a spring or seep that saturates the soil and allows water to pool in small areas (Hendrickson and Minckley, 1985). Ciénega vegetation includes rushes, watercress, and smartweed; they are often bordered by cottonwood and willow trees which the Refuge will plant at the new ciénega once the water has settled into the excavated depression. The new ciénega on BANWR is plumbed, meaning that water levels can be modified and maintained in response to environmental conditions that range from multi-year droughts to heavy monsoon conditions (Bezy et. al. 2007).

ARMI has collected data at BANWR since 2000, gathering information about native species (e.g., Chandler et al. 2015, Jarchow et al. 2016, , Howell et al. 2018, 2020a), invasive species such as American bullfrogs and sunfish (Suhr 2010, Howell et al. 2020b), and disease (Sigafus et al. 2014). Models assessing occupancy and movement of the Chiricahua leopard frog indicate that water availability and permanency are critical components to the its persistence at BANWR. The Refuge has used this information to make decisions about management actions such as building plumbed ponds — a non-trivial action in terms of cost and logistics. The new ciénega will not only support the Chiricahua leopard frog but will provide water for many other species that call BANWR home including the federally endangered masked bobwhite quail, great blue herons, yellow-billed cuckoos, and pronghorn.

Bezy, J., C. F. Hutchinson, and C. J. Bahre. 2007. Buenos Aires National Wildlife Refuge, 373 Arizona. Desert Plants 23:3–44.

Chandler, R., E. Muths, B. H. Sigafus, C.R. Schwalbe, C. Jarchow, and B.R. Hossack. 2015. Realizing the potential of spatially explicit metapopulation theory for predicting extinction risk. Journal of Applied Ecology. DOI: 10.1111/1365-2664.12481.

Howell, P., E. Muths, B.R. Hossack, B.H. Sigafus, and R.B. Chandler. 2018. Increasing connectivity between metapopulation ecology and landscape ecology. Ecology 99: 1119-1128.

Howell, P.E., B.R. Hossack, E. Muths, B. Sigafus, A. Chenevert- Steffler, and R. Chandler. 2020a. A statistical forecasting approach to metapopulation viability analysis. Ecological Applications 30(2), e02038

Howell, P.E. E. Muths, B.H. Sigafus, and B.R. Hossack. 2020b. Survival estimates for the invasive American bullfrog. Amphibia-Reptilia.

Hendrickson, D. A. and W. L. Minckley. 1985. Ciénegas vanishing climax communities of the American southwest. Desert Plants 6 (3): 131-175.

Sigafus, B. H., C.R. Schwalbe, B.R. Hossack, and E. Muths. 2014. Prevalence of the amphibian chytrid fungus at Buenos Aires National Wildlife Refuge, Arizona. Herpetological Review 45: 41-42.

Suhre, D. O. 2010. Dispersal and demography of the American Bullfrog (Rana catesbeiana) in a 455 semi-arid grassland. M.S. Thesis, University of Arizona, USA.

The disproportionate effects of warming for high-latitude, freshwater ecosystems has been well documented, but in some areas, changes have been further impacted by human-subsidized increases of waterfowl. To gain insight into how predicted changes in temperature and nutrient inputs might affect ecosystem function, we conducted a mesocosm experiment in the Canadian Subarctic with three levels of simulated goose enrichment and warming to measure changes in size and survival of larval wood frogs and boreal chorus frogs and primary productivity (phytoplankton and periphyton biomass). Our results highlight that the consequences of these rapid changes are non-linear and even non-intuitive, with species-specific consumer and ecosystem responses that depend on the magnitude of temperature and nutrient changes as well as community composition.

Data release for data collected from Prairie Pothole wetlands.

Data release for manuscript summarizing responses of Rocky Mountain Tailed Frogs to the 2003 Robert Fire in Glacier National Park.

Data release for acoustic monitoring of amphibians at 3 mine sites.

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.

Small waterbodies are numerically dominant in many landscapes and provide several important ecosystem services, but automated measurement of waterbodies smaller than a standard Landsat pixel (https://0.09 ha) remains challenging. To further evaluate sub-Landsat pixel techniques for estimating inundation extent of small waterbodies (basin area: https://0.06 – https://1.79 ha), we used a partial spectral unmixing method with matched filtering applied to September 1985–2018 Landsat 5 and 8 imagery from southern Arizona, USA. We estimated trends in modeled surface water area each September and evaluated the ability of several common drought indices to explain variation in mean water area. Our methods accurately classified waterbodies as dry or inundated (Landsat 5: 91.3%; Landsat 8: 98.9%) and modeled and digitized surface water areas were strongly correlated (R2 = https://0.70 – https://0.92; bias = -https://0.024 – -https://0.015 ha). Estimated surface water area was best explained by the 3-month seasonal standardized precipitation index (SPI03; July?September) and. We found a wide range of estimated relationships between drought indices (e.g., SPI vs. Palmer Drought Severity Index) and estimated water area, even for different durations of the same drought index (e. g., SPI01 vs SPI12). Mean surface area of waterbodies decreased by ~14% from September 1985 to September 2018, which matched declines in annual precipitation in the area and is consistent with broader trends of reduced inundation extent based on larger waterbodies. Estimated of surface water area and trends over time were also consistent when we limited analyses to waterbodies ? https://0.04 ha or those that varied most in size (based on CV). These results emphasize the importance of understanding local systems when relying on drought indices to infer variation in past or future surface water dynamics. Several challenges remain before widespread application of sub-pixel methods is feasible, but our results provide further evidence that partial spectral unmixing with matched filtering provides reliable measures of inundation extent of small waterbodies.

Understanding the demographic consequences of interactions among pathogens, hosts, and weather conditions is critical in determining how amphibian populations respond to disease and in identifying site-specific conservation actions that can be developed to bolster persistence of amphibian populations. We investigated population dynamics in Boreal Toads relative to abiotic (fall temperatures and snowpack) and biotic (the abundance of another anuran host and disease) characteristics of the local environment in Wyoming, USA. We used capture-recapture data and a multi-state model where state is treated as a hidden Markov process to incorporate disease state uncertainty and assess our a priori hypotheses. Our results indicate that snowpack during the coldest week of the winter is more influential to toad survival, disease transition probabilities, and the population-level prevalence of the amphibian chytrid fungus (Batrachochytrium dendrobatidis) in the spring, than temperatures in the fall or the presence of another host. As hypothesized, apparent survival at low (i.e., <25 cm) snowpack (https://0.22 [CI: 0.15–0.31]) was lower than apparent survival at high snowpack (https://90.65 [CI: 0.50–0.78]). Our findings highlight the potential for local environmental factors, like snowpack, to influence disease and host persistence, and demonstrate the ecological complexity of disease effects on population demography in natural environments. This work further emphasizes the need for improved understanding of how climate change may influence the relationships among pathogens, hosts, and their environment for wild animal populations challenged by disease.