Application of multiple-population viability analysis to evaluate species recovery alternatives

Population viability analysis (PVA) is a powerful conservation tool, but it remains impractical for many species, particularly species with multiple, broadly distributed populations for which collecting suitable data can be challenging. A recently developed method of multiple-population viability analysis (MPVA), however, addresses many limitations of traditional PVA. We built on previous development of MPVA for Lahontan cutthroat trout (LCT) (Oncorhynchus clarkii henshawi), a species listed under the U.S. Endangered Species Act, that is distributed broadly across habitat fragments in the Great Basin (U.S.A.). We simulated potential management scenarios and assessed their effects on population sizes and extinction risks in 211 streams, where LCT exist or may be reintroduced. Conservation populations (those managed for recovery) tended to have lower extinction risks than nonconservation populations (mean = 19.8% vs. 52.7%), but not always. Active management or reprioritization may be warranted in some cases. Eliminating non-native trout had a strong positive effect on overall carrying capacities for LCT populations but often did not translate into lower extinction risks unless simulations also reduced associated stochasticity (to the mean for populations without non-native trout). Sixty fish or 5–10 fish/km was the minimum reintroduction number and density, respectively, that provided near-maximum reintroduction success. This modeling framework provided crucial insights and empirical justification for conservation planning and specific adaptive management actions for this threatened species. More broadly, MPVA is applicable to a wide range of species exhibiting geographic rarity and limited availability of abundance data and greatly extends the potential use of empirical PVA for conservation assessment and planning.     

Use of ambiguous detections to improve estimates from species distribution models

As large carnivores recover throughout Europe, their distribution needs to be studied to determine their conservation status and assess the potential for human-carnivore conflicts. However, efficient monitoring of many large carnivore species is challenging due to their rarity, elusive behavior, and large home ranges. Their monitoring can include opportunistic sightings from citizens in addition to designed surveys. Two types of detection errors may occur in such monitoring schemes: false negatives and false positives. False-negative detections can be accounted for in species distribution models (SDMs) that deal with imperfect detection. False-positive detections, due to species misidentification, have rarely been accounted for in SDMs. Generally, researchers use ad hoc data-filtering methods to discard ambiguous observations prior to analysis. These practices may discard valuable ecological information on the distribution of a species. We investigated the costs and benefits of including data types that may include false positives rather than discarding them for SDMs of large carnivores. We used a dynamic occupancy model that simultaneously accounts for false negatives and positives to jointly analyze data that included both unambiguous detections and ambiguous detections. We used simulations to compare the performances of our model with a model fitted on unambiguous data only. We tested the 2 models in 4 scenarios in which parameters that control false-positive detections and true detections varied. We applied our model to data from the monitoring of the Eurasian lynx (Lynx lynx) in the European Alps. The addition of ambiguous detections increased the precision of parameter estimates. For the Eurasian lynx, incorporating ambiguous detections produced more precise estimates of the ecological parameters and revealed additional occupied sites in areas where the species is likely expanding. Overall, we found that ambiguous data should be considered when studying the distribution of large carnivores through the use of dynamic occupancy models that account for misidentification.     

Rapid,recurring, structured survey versus bioblitz for generating biodiversity data and analysis with a multispecies abundance model

A bioblitz inexpensively and quickly generates biodiversity data, but bioblitzes are often conducted with haphazard, unreplicated sampling. Results tend to be taxonomically, geographically, or temporally biased, lack metadata, and consist of lists of observed taxa that do not enable further analyses or correction for imperfect detection. A rapid, recurring, structured survey (RRSS) uses a structured sampling design and temporal and spatial replication to survey randomly selected sites on a conservation property. We participated in a loosely structured bioblitz and a subsequent RRSS at Big Canoe Creek Nature Preserve in Springville (St. Clair County), Alabama (USA) to compare observed richness derived from the 2 survey approaches. The RRSS data structure enabled us to fit models that accounted for imperfect detection to estimate abundances, occupancy probabilities, and habitat associations. The loosely structured bioblitz data could not be used in such models. We present a new integrated multispecies abundance model that we applied to avian RRSS data. Our model extension enables estimation for the community, employs data augmentation to estimate the number of undetected species, and incorporates covariates. The RRSS generated a more comprehensive and less biased list of observed taxonomic richness than the loosely structured bioblitz (e.g., 73 vs. 45 bird species and 104 vs. 63 insect families from the RRSS vs. loosely structured bioblitz, respectively). Models fit to the RRSS data identified seasonal patterns in avian community composition and allowed for estimation of habitat–occupancy relationships for insect taxa. The RRSS protocol has potential for broad transferability as a standardized, quick, and inexpensive way to inventory biodiversity and estimate ecological parameters while providing an outreach opportunity.     

Site‐Occupancy Distribution Modeling to Correct Population‐Trend Estimates Derived from Opportunistic Observations

Abstract: Species’ assessments must frequently be derived from opportunistic observations made by volunteers (i.e., citizen scientists). Interpretation of the resulting data to estimate population trends is plagued with problems, including teasing apart genuine population trends from variations in observation effort. We devised a way to correct for annual variation in effort when estimating trends in occupancy (species distribution) from faunal or floral databases of opportunistic observations. First, for all surveyed sites, detection histories (i.e., strings of detection–nondetection records) are generated. Within‐season replicate surveys provide information on the detectability of an occupied site. Detectability directly represents observation effort; hence, estimating detectablity means correcting for observation effort. Second, site‐occupancy models are applied directly to the detection‐history data set (i.e., without aggregation by site and year) to estimate detectability and species distribution (occupancy, i.e., the true proportion of sites where a species occurs). Site‐occupancy models also provide unbiased estimators of components of distributional change (i.e., colonization and extinction rates). We illustrate our method with data from a large citizen‐science project in Switzerland in which field ornithologists record opportunistic observations. We analyzed data collected on four species: the widespread Kingfisher (Alcedo atthis) and Sparrowhawk (Accipiter nisus) and the scarce Rock Thrush (Monticola saxatilis) and Wallcreeper (Tichodroma muraria). Our method requires that all observed species are recorded. Detectability was <1 and varied over the years. Simulations suggested some robustness, but we advocate recording complete species lists (checklists), rather than recording individual records of single species. The representation of observation effort with its effect on detectability provides a solution to the problem of differences in effort encountered when extracting trend information from haphazard observations. We expect our method is widely applicable for global biodiversity monitoring and modeling of species distributions.     

Longitudinal monitoring of neutral and adaptive genomic diversity in a reintroduction

Restoration programs in the form of ex-situ breeding combined with reintroductions are becoming critical to counteract demographic declines and species losses. Such programs are increasingly using genetic management to improve conservation outcomes. However, the lack of long-term monitoring of genetic indicators following reintroduction prevents assessments of the trajectory and persistence of reintroduced populations. We carried out an extensive monitoring program in the wild for a threatened small-bodied fish (southern pygmy perch, Nannoperca australis) to assess the long-term genomic effects of its captive breeding and reintroduction. The species was rescued prior to its extirpation from the terminal lakes of Australia's Murray-Darling Basin, and then used for genetically informed captive breeding and reintroductions. Subsequent annual or biannual monitoring of abundance, fitness, and occupancy over a period of 11 years, combined with postreintroduction genetic sampling, revealed survival and recruitment of reintroduced fish. Genomic analyses based on data from the original wild rescued, captive born, and reintroduced cohorts revealed low inbreeding and strong maintenance of neutral and candidate adaptive genomic diversity across multiple generations. An increasing trend in the effective population size of the reintroduced population was consistent with field monitoring data in demonstrating successful re-establishment of the species. This provides a rare empirical example that the adaptive potential of a locally extinct population can be maintained during genetically informed ex-situ conservation breeding and reintroduction into the wild. Strategies to improve biodiversity restoration via ex-situ conservation should include genetic-based captive breeding and longitudinal monitoring of standing genomic variation in reintroduced populations.     

Assessing species reintroduction sites based on future climate suitability for food resources

The reintroduction of a species that is extinct in the wild demands caution because reintroduction locations may be associated with threats, such as hunting, poor-quality habitat, and climate change. This is the case for Cyanopsitta spixii (Spix's Macaw), which has been extinct in the wild since 2000. The few living individuals were created in captivity and will be used in a reintroduction project within the species’ original distribution area, the Caatinga domain (Brazil). Because the occurrence records for this bird are old and inaccurate, we investigated the current and future environmental suitability of the 14 plant species used by C. spixii as resource. These plants are key elements for the long-term reestablishment of the species in the wild, so the use of models helps in the assessment of the effects of climate change on the availability of these resources for the species and informs selection of the best places for reintroduction. We based our models of environmental suitability on 19 bioclimatic variables and nine physical soil and topography variables. Climate projections were created for the present and for the year 2070 with an optimistic (SSP2-4.5) and a pessimistic (SSP5-8.5) climate scenario. Both future climate scenarios lead to a reduction in area of environmental suitability that overlapped for all the plant species: 33% reduction for SSP2-4.5 and 63% reduction for SSP5-8.5. If our projections materialize, climate change could thus affect the distribution of key resources, and the maintenance of C. spixii would depend on restoration of degraded areas, especially riparian forests, and the preservation of already existing natural areas. The Caatinga domain is very threatened by habitat loss and, for the success of this reintroduction project, the parties involved must act to protect the species and the resources it uses.     

Predicting species distributions from samples collected along roadsides

Predictive models of species distributions are typically developed with data collected along roads. Roadside sampling may provide a biased (nonrandom) sample; however, it is currently unknown whether roadside sampling limits the accuracy of predictions generated by species distribution models. We tested whether roadside sampling affects the accuracy of predictions generated by species distribution models by using a prospective sampling strategy designed specifically to address this issue. We built models from roadside data and validated model predictions at paired locations on unpaved roads and 200 m away from roads (off road), spatially and temporally independent from the data used for model building. We predicted species distributions of 15 bird species on the basis of point-count data from a landbird monitoring program in Montana and Idaho (U.S.A.). We used hierarchical occupancy models to account for imperfect detection. We expected predictions of species distributions derived from roadside-sampling data would be less accurate when validated with data from off-road sampling than when it was validated with data from roadside sampling and that model accuracy would be differentially affected by whether species were generalists, associated with edges, or associated with interior forest. Model performance measures (kappa, area under the curve of a receiver operating characteristic plot, and true skill statistic) did not differ between model predictions of roadside and off-road distributions of species. Furthermore, performance measures did not differ among edge, generalist, and interior species, despite a difference in vegetation structure along roadsides and off road and that 2 of the 15 species were more likely to occur along roadsides. If the range of environmental gradients is surveyed in roadside-sampling efforts, our results suggest that surveys along unpaved roads can be a valuable, unbiased source of information for species distribution models.     

Factors associated with preemptive conservation under the U.S. Endangered Species Act

In recent decades, there has been an increasing emphasis on proactive efforts to conserve species being considered for listing under the U.S. Endangered Species Act (ESA) before they are listed (i.e., preemptive conservation). These efforts, which depend on voluntary actions by public and private land managers across the species’ range, aim to conserve species while avoiding regulatory costs associated with ESA listing. We collected data for a set of social, economic, environmental, and institutional factors that we hypothesized would influence voluntary decisions to promote or inhibit preemptive conservation of species under consideration for ESA listing. We used logistic regression to estimate the association of these factors with preemptive conservation outcomes based on data for a set of species that entered the ESA listing process and were either officially listed (n = 314) or preemptively conserved (n = 73) from 1996 to 2018. Factors significantly associated with precluded listing due to preemptive conservation included high baseline conservation status, low proportion of private land across the species’ range, small total range size, exposure to specific types of threats, and species’ range extending over several states. These results highlight strategies that can help improve conservation outcomes, such as allocating resources for imperiled species earlier in the listing process, addressing specific threats, and expanding incentives and coordination mechanisms for conservation on private lands.     

How decisions about fitting species distribution models affect conservation outcomes

Species distribution models (SDMs) are increasingly used in conservation and land-use planning as inputs to describe biodiversity patterns. These models can be built in different ways, and decisions about data preparation, selection of predictor variables, model fitting, and evaluation all alter the resulting predictions. Commonly, the true distribution of species is unknown and independent data to verify which SDM variant to choose are lacking. Such model uncertainty is of concern to planners. We analyzed how 11 routine decisions about model complexity, predictors, bias treatment, and setting thresholds for predicted values altered conservation priority patterns across 25 species. Models were created with MaxEnt and run through Zonation to determine the priority rank of sites. Although all SDM variants performed well (area under the curve >0.7), they produced spatially different predictions for species and different conservation priority solutions. Priorities were most strongly altered by decisions to not address bias or to apply binary thresholds to predicted values; on average 40% and 35%, respectively, of all grid cells received an opposite priority ranking. Forcing high model complexity altered conservation solutions less than forcing simplicity (14% and 24% of cells with opposite rank values, respectively). Use of fewer species records to build models or choosing alternative bias treatments had intermediate effects (25% and 23%, respectively). Depending on modeling choices, priority areas overlapped as little as 10–20% with the baseline solution, affecting top and bottom priorities differently. Our results demonstrate the extent of model-based uncertainty and quantify the relative impacts of SDM building decisions. When it is uncertain what the best SDM approach and conservation plan is, solving uncertainty or considering alterative options is most important for those decisions that change plans the most.     

Generalized site occupancy models allowing for false positive and false negative errors

Site occupancy models have been developed that allow for imperfect species detection or "false negative" observations. Such models have become widely adopted in surveys of many taxa. The most fundamental assumption underlying these models is that "false positive" errors are not possible. That is, one cannot detect a species where it does not occur. However, such errors are possible in many sampling situations for a number of reasons, and even low false positive error rates can induce extreme bias in estimates of site occupancy when they are not accounted for. In this paper, we develop a model for site occupancy that allows for both false negative and false positive error rates. This model can be represented as a two-component finite mixture model and can be easily fitted using freely available software. We provide an analysis of avian survey data using the proposed model and present results of a brief simulation study evaluating the performance of the maximum-likelihood estimator and the naive estimator in the presence of false positive errors.     

Sentiment analysis of the Twitter response to Netflix's Our Planet documentary

The role of nature documentaries in shaping public attitudes and behavior toward conservation and wildlife issues is unclear. We analyzed the emotional content of over 2 million tweets related to Our Planet, a major nature documentary released on Netflix, with dictionary and rule-based automatic sentiment analysis. We also compared the sentiment associated with species mentioned in Our Planet and a set of control species with similar features but not mentioned in the documentary. Tweets were largely negative in sentiment at the time of release of the series. This effect was primarily linked to the highly skewed distributions of retweets and, in particular, to a single negatively valenced and massively retweeted tweet (>150,000 retweets). Species mentioned in Our Planet were associated with more negative sentiment than the control species, and this effect coincided with a short period following the airing of the series. Our results are consistent with a general negativity bias in cultural transmission and document the difficulty of evoking positive sentiment, on social media and elsewhere, in response to environmental problems.     

Assessing ecological function in the context of species recovery

Species interactions matter to conservation. Setting an ambitious recovery target for a species requires considering the size, density, and demographic structure of its populations such that they fulfill the interactions, roles, and functions of the species in the ecosystems in which they are embedded. A recently proposed framework for an International Union for Conservation of Nature Green List of Species formalizes this requirement by defining a fully recovered species in terms of representation, viability, and functionality. Defining and quantifying ecological function from the viewpoint of species recovery is challenging in concept and application, but also an opportunity to insert ecological theory into conservation practice. We propose 2 complementary approaches to assessing a species’ ecological functions: confirmation (listing interactions of the species, identifying ecological processes and other species involved in these interactions, and quantifying the extent to which the species contributes to the identified ecological process) and elimination (inferring functionality by ruling out symptoms of reduced functionality, analogous to the red-list approach that focuses on symptoms of reduced viability). Despite the challenges, incorporation of functionality into species recovery planning is possible in most cases and it is essential to a conservation vision that goes beyond preventing extinctions and aims to restore a species to levels beyond what is required for its viability. This vision focuses on conservation and recovery at the species level and sees species as embedded in ecosystems, influencing and being influenced by the processes in those ecosystems. Thus, it connects and integrates conservation at the species and ecosystem levels.     

A Multispecies Framework for Landscape Conservation Planning

Abstract: Rapidly changing landscapes have spurred the need for quantitative methods for conservation assessment and planning that encompass large spatial extents. We devised and tested a multispecies framework for conservation planning to complement single‐species assessments and ecosystem‐level approaches. Our framework consisted of 4 elements: sampling to effectively estimate population parameters, measuring how human activity affects landscapes at multiple scales, analyzing the relation between landscape characteristics and individual species occurrences, and evaluating and comparing the responses of multiple species to landscape modification. We applied the approach to a community of terrestrial birds across 25,000 km² with a range of intensities of human development. Human modification of land cover, road density, and other elements of the landscape, measured at multiple spatial extents, had large effects on occupancy of the 67 species studied. Forest composition within 1 km of points had a strong effect on occupancy of many species and a range of negative, intermediate, and positive associations. Road density within 1 km of points, percent evergreen forest within 300 m, and distance from patch edge were also strongly associated with occupancy for many species. We used the occupancy results to group species into 11 guilds that shared patterns of association with landscape characteristics. Our multispecies approach to conservation planning allowed us to quantify the trade‐offs of different scenarios of land‐cover change in terms of species occupancy.     

Evaluating institutional fit for the conservation of threatened species

Recovery and conservation of threatened species require adequate institutional responses. We tested an approach to systematically identify and measure how an institutional framework acknowledges threats and required responses for the recovery of endangered species. We measured institutional functional fit with a drivers-pressure-state-impacts-response (DPSIR) model integrated with a quantitative text mining method and qualitative analysis of statutory instruments to examine regulatory responses that support the recovery of 2 endangered species native to Australia, the bridled nailtail wallaby (Onychogalea fraenata) and the Eastern Bristlebird (Dasyornis brachypterus). The key components of the DPSIR model were present in the institutional framework at statutory and operational levels, but some institutional gaps remained in the protection and recovery of the Eastern Bristlebird, including feral predator control, weed control, and grazing management in some locations. However, regulatory frameworks varied in their geographic scope and the application and implementation of many instruments remained optional. Quantitative text mining can be used to quickly navigate a large volume of regulatory documents, but challenges remain in selection of terms, queries of co-occurrence, and interpretation of word frequency counts. To inform policy, we recommend that quantitative assessments of institutional fit be complemented with qualitative analysis and interpreted in light of the sociopolitical and institutional context.     

A robust-design formulation of the incidence function model of metapopulation dynamics applied to two species of rails

The incidence function model (IFM) uses area and connectivity to predict metapopulation dynamics. However, false absences and missing data can lead to underestimates of the number of sites contributing to connectivity, resulting in overestimates of dispersal ability and turnovers (extinctions plus colonizations). We extend estimation methods for the IFM by using a hierarchical Bayesian model to account both for false absences due to imperfect detection and for missing data due to sites not surveyed in some years. We compare parameter estimates, measures of metapopulation dynamics, and forecasts using stochastic patch occupancy models (SPOMs) among three IFM models: (1) a Bayesian formulation assuming no false absences and omitting site-year combinations with missing data; (2) a hierarchical Bayesian formulation assuming no false absences but incorporating missing data; and (3) a hierarchical Bayesian formulation allowing for imperfect detection and incorporating missing data. We fit the models to multiyear data sets of occupancy for two bird species that differ in body size and presumed dispersal ability but inhabit the same network of sites: the small Black Rail (Laterallus jamaicensis) and the medium-sized Virginia Rail (Rallus limicola). Incorporating missing data affected colonization parameters and led to lower estimates of dispersal ability for the Black Rail. Detection rates were high for the Black Rail in most years but moderate for the Virginia Rail. Incorporating imperfect detection resulted in higher occupancy and lower turnover rates for both species, with largest effects for the Virginia Rail. Forecasts using SPOMs were sensitive to both missing data and false absences; persistence in models assuming no false absences was more optimistic than from robust models. Our results suggest that incorporating false absences and missing data into the IFM can improve (1) estimates of dispersal ability and the effect of connectivity on colonization, (2) the scaling of extinction risk with patch area, and (3) forecasts of occupancy and turnover rates.     

Human disturbance and shifts in vertebrate community composition in a biodiversity hotspot

Understanding how human modification of the landscape shapes vertebrate community composition is vital to understanding the current status and future trajectory of wildlife. Using a participatory approach, we deployed the largest camera-trap network in Mesoamerica to date to investigate how anthropogenic disturbance shapes the occupancy and co-occurrence of terrestrial vertebrate species in a tropical biodiversity hotspot: the Osa Peninsula, Costa Rica. We estimated species richness in different categories of land protection with rarefaction analysis and estimated the expected occupancy with a joint species distribution model that included covariates for anthropogenic disturbance, land protection, habitat quality, and habitat availability. Areas with the most stringent land-use protections (e.g., Corcovado National Park, 24 species [95% CI 23–25]) harbored significantly more species than unprotected areas (20 species [19.7–20.3]), mainly due to a reduced presence of large-bodied species of conservation concern in unprotected areas (e.g., jaguar Panthera onca and white-lipped peccary Tayassu pecari). Small-bodied generalist species, such as opossums (Didelphidae) and armadillos (Dasypus novemcinctus), in contrast, were more common at disturbed sites, resulting in a significant difference in vertebrate community composition between sites with low and high disturbance. Co-occurrence of species was also mainly associated with response to disturbance. Similar responses to disturbance create two groups of species, those whose site-level occupancy usually increased as anthropogenic disturbance increased and those whose estimated occupancy decreased. The absence of large-bodied species entails an important loss of ecological function in disturbed areas and can hinder forest development and maintenance. Efforts to protect and restore forested landscapes are likely having a positive effect on the abundance of some threatened species. These efforts, however, must be sustained and expanded to increase connectivity and ensure the long-term viability of the wildlife community.     

Using expert knowledge to support Endangered Species Act decision-making for data-deficient species

Many questions relevant to conservation decision-making are characterized by extreme uncertainty due to lack of empirical data and complexity of the underlying ecologic processes, leading to a rapid increase in the use of structured protocols to elicit expert knowledge. Published ecologic applications often employ a modified Delphi method, where experts provide judgments anonymously and mathematical aggregation techniques are used to combine judgments. The Sheffield elicitation framework (SHELF) differs in its behavioral approach to synthesizing individual judgments into a fully specified probability distribution for an unknown quantity. We used the SHELF protocol remotely to assess extinction risk of three subterranean aquatic species that are being considered for listing under the U.S. Endangered Species Act. We provided experts an empirical threat assessment for each known locality over a video conference and recorded judgments on the probability of population persistence over four generations with online submission forms and R-shiny apps available through the SHELF package. Despite large uncertainty for all populations, there were key differences between species’ risk of extirpation based on spatial variation in dominant threats, local land use and management practices, and species’ microhabitat. The resulting probability distributions provided decision makers with a full picture of uncertainty that was consistent with the probabilistic nature of risk assessments. Discussion among experts during SHELF's behavioral aggregation stage clearly documented dominant threats (e.g., development, timber harvest, animal agriculture, and cave visitation) and their interactions with local cave geology and species’ habitat. Our virtual implementation of the SHELF protocol demonstrated the flexibility of the approach for conservation applications operating on budgets and time lines that can limit in-person meetings of geographically dispersed experts.     

The importance of incorporating functional habitats into conservation planning for highly mobile species in dynamic systems

The distribution of mobile species in dynamic systems can vary greatly over time and space. Estimating their population size and geographic range can be problematic and affect the accuracy of conservation assessments. Scarce data on mobile species and the resources they need can also limit the type of analytical approaches available to derive such estimates. We quantified change in availability and use of key ecological resources required for breeding for a critically endangered nomadic habitat specialist, the Swift Parrot (Lathamus discolor). We compared estimates of occupied habitat derived from dynamic presence‐background (i.e., presence‐only data) climatic models with estimates derived from dynamic occupancy models that included a direct measure of food availability. We then compared estimates that incorporate fine‐resolution spatial data on the availability of key ecological resources (i.e., functional habitats) with more common approaches that focus on broader climatic suitability or vegetation cover (due to the absence of fine‐resolution data). The occupancy models produced significantly (P < 0.001) smaller (up to an order of magnitude) and more spatially discrete estimates of the total occupied area than climate‐based models. The spatial location and extent of the total area occupied with the occupancy models was highly variable between years (131 and 1498 km²). Estimates accounting for the area of functional habitats were significantly smaller (2–58% [SD 16]) than estimates based only on the total area occupied. An increase or decrease in the area of one functional habitat (foraging or nesting) did not necessarily correspond to an increase or decrease in the other. Thus, an increase in the extent of occupied area may not equate to improved habitat quality or function. We argue these patterns are typical for mobile resource specialists but often go unnoticed because of limited data over relevant spatial and temporal scales and lack of spatial data on the availability of key resources. Understanding changes in the relative availability of functional habitats is crucial to informing conservation planning and accurately assessing extinction risk for mobile resource specialists.     

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. Incrementos Rápidos y Declinaciones Desfasadas en la Ocupación de Anfibios Después de un Incendio