Modelling the Distribution and Habitat Suitability of the European Wildcat (Felis silvestris) in North-Western Spain and Its Conservation Implications.

Human activities have resulted in severe habitat degradation and fragmentation at a global scale. Despite this scenario, some carnivore species that adapted to the new conditions are expanding, leading to close coexistence with humans and the emergence of potential conflicts. In this work, we used a European wildcat (Felis silvestris) observations database of more than 350 sightings over 17 years in NW Spain to build suitability models based on environmental, topographic, climatic, and human impact variables. MaxEnt was used to analyse the availability of suitable habitats for the species at a regional scale. Our results showed that less than one third of the suitable area for the species had confirmed wildcat presence. Elevation, the percentage of forested area, and footpath density were the three main variables conditioning wildcat presence, with the first two variables having positive effects and footpath density negatively affecting wildcat presence. The selection of high areas and forest areas by the species seems to be related to food availability, while the avoidance of footpaths seems to be related to the fact that main mortality causes are linked to human disturbances. The results enhance the understanding of the European wildcat ecology and provide insight into potential management plans to ensure the conservation of one of the main populations of the species throughout its range.

Body length changes for Atlantic salmon (Salmo salar) over five decades exhibit weak spatial synchrony over a broad latitudinal gradient.

Understanding the factors that drive spatial synchrony among populations or species is important for management and recovery of populations. The range-wide declines in Atlantic salmon (Salmo salar) populations may be the result of broad-scale changes in the marine environment. Salmon undergo rapid growth in the ocean; therefore changing marine conditions may affect body size and fecundity estimates used to evaluate whether stock reference points are met. Using a dataset that spanned five decades, 172,268 individuals, and 19 rivers throughout Eastern Canada, we investigated the occurrence of spatial synchrony in changes in the body size of returning wild adult Atlantic salmon. Body size was then related to conditions in the marine environment (i.e., climate indices, thermal habitat availability, food availability, density-dependence, and fisheries exploitation rates) that may act on all populations during the ocean feeding phase of their life cycle. Body size increased during the 1980s and 1990s for salmon that returned to rivers after one (1SW) or two winters at sea (2SW); however, significant changes were only observed for 1SW and/or 2SW in some mid-latitude and northern rivers (10/13 rivers with 10 of more years of data during these decades) and not in southern rivers (0/2), suggesting weak spatial synchrony across Eastern Canada. For 1SW salmon in nine rivers, body size was longer when fisheries exploitation rates were lower. For 2SW salmon, body size was longer when suitable thermal habitat was more abundant (significant for 3/8 rivers) and the Atlantic Multidecadal Oscillation was higher (i.e., warmer sea surface temperatures; significant for 4/8 rivers). Overall, the weak spatial synchrony and variable effects of covariates on body size across rivers suggest that changes in Atlantic salmon body size may not be solely driven by shared conditions in the marine environment. Regardless, body size changes may have consequences for population management and recovery through the relationship between size and fecundity.

Vegetation structure from LiDAR explains the local richness of birds across Denmark.

Classic ecological research into the determinants of biodiversity patterns emphasised the important role of three-dimensional (3D) vegetation heterogeneity. Yet, measuring vegetation structure across large areas has historically been difficult. A growing focus on large-scale research questions has caused local vegetation heterogeneity to be overlooked compared with more readily accessible habitat metrics from, for example, land cover maps. Using newly available 3D vegetation data, we investigated the relative importance of habitat and vegetation heterogeneity for explaining patterns of bird species richness and composition across Denmark (42,394 km2 ). We used standardised, repeated point counts of birds conducted by volunteers across Denmark alongside metrics of habitat availability from land-cover maps and vegetation structure from rasterised LiDAR data (10 m resolution). We used random forest models to relate species richness to environmental features and considered trait-specific responses by grouping species by nesting behaviour, habitat preference and primary lifestyle. Finally, we evaluated the role of habitat and vegetation heterogeneity metrics in explaining local bird assemblage composition. Overall, vegetation structure was equally as important as habitat availability for explaining bird richness patterns. However, we did not find a consistent positive relationship between species richness and habitat or vegetation heterogeneity; instead, functional groups displayed individual responses to habitat features. Meanwhile, habitat availability had the strongest correlation with the patterns of bird assemblage composition. Our results show how LiDAR and land cover data complement one another to provide insights into different facets of biodiversity patterns and demonstrate the potential of combining remote sensing and structured citizen science programmes for biodiversity research. With the growing coverage of LiDAR surveys, we are witnessing a revolution of highly detailed 3D data that will allow us to integrate vegetation heterogeneity into studies at large spatial extents and advance our understanding of species' physical niches.

Individual differences in habitat selection mediate landscape level predictions of a functional response.

Predicting future space use by animals requires models that consider both habitat availability and individual differences in habitat selection. The functional response in habitat selection posits animals adjust their habitat selection to availability, but population-level responses to availability may differ from individual responses. Generalized functional response (GFR) models account for functional responses by including fixed effect interactions between habitat availability and selection. Population-level resource selection functions instead account for individual selection responses to availability with random effects. We compared predictive performance of both approaches using a functional response in elk (Cervus canadensis) selection for mixed forest in response to road proximity, and avoidance of roads in response to mixed forest availability. We also investigated how performance changed when individuals responded differently to availability from the rest of the population. Individual variation in road avoidance decreased performance of both models (random effects: ß = 0.69, 95% CI 0.47, 0.91; GFR: ß = 0.38, 95% CI 0.05, 0.71). Changes in individual road and forest availability affected performance of neither model, suggesting individual responses to availability different from the functional response mediated performance. We also found that overall, both models performed similarly for predicting mixed forest selection (F1, 58 = 0.14, p = 0.71) and road avoidance (F1, 58 = 0.28, p = 0.60). GFR estimates were slightly better, but its larger number of covariates produced greater variance than the random effects model. Given this bias-variance trade-off, we conclude that neither model performs better for future space use predictions.

Demographic changes in Pleistocene sea turtles were driven by past sea level fluctuations affecting feeding habitat availability.

Pleistocene environmental changes are generally assumed to have dramatically affected species' demography via changes in habitat availability, but this is challenging to investigate due to our limited knowledge of how Pleistocene ecosystems changed through time. Here, we tracked changes in shallow marine habitat availability resulting from Pleistocene sea level fluctuations throughout the last glacial cycle (120-14 thousand years ago; kya) and assessed correlations with past changes in genetic diversity inferred from genome-wide SNPs, obtained via ddRAD sequencing, in Caribbean hawksbill turtles, which feed in coral reefs commonly found in shallow tropical waters. We found sea level regression resulted in an average 75% reduction in shallow marine habitat availability during the last glacial cycle. Changes in shallow marine habitat availability correlated strongly with past changes in hawksbill turtle genetic diversity, which gradually declined to ~1/4th of present-day levels during the Last Glacial Maximum (LGM; 26-19 kya). Shallow marine habitat availability and genetic diversity rapidly increased after the LGM, signifying a population expansion in response to warming environmental conditions. Our results suggest a positive correlation between Pleistocene environmental changes, habitat availability and species' demography, and that demographic changes in hawksbill turtles were potentially driven by feeding habitat availability. However, we also identified challenges associated with disentangling the potential environmental drivers of past demographic changes, which highlights the need for integrative approaches. Our conclusions underline the role of habitat availability on species' demography and biodiversity, and that the consequences of ongoing habitat loss should not be underestimated.

Does stress mess with rodents' heads? Influence of habitat amount and genetic factors in mandible fluctuating asymmetry in South American water rats (Nectomys squamipes, Sigmodontinae) from Brazilian Atlantic rainforest remnants.

Loss of developmental stability can lead to deviations from bilateral symmetry (i.e. Fluctuating Asymmetry - FA), and is thought to be caused by environmental and genetic factors associated with habitat loss and stress. Therefore, levels of FA might be a valuable tool to monitor wild populations if FA serves as an indicator of exposure to stress due to impacts of habitat loss and fragmentation. In studies examining FA and habitat fragmentation, FA levels are often explained by loss of genetic variation, though few studies have addressed FA's use as indicator of environmental impact. Here, we investigated whether habitat loss, genetic variation, and/or inbreeding affect the developmental instability in Brazilian Atlantic forest populations of a Neotropical water rat (Nectomys squamipes). We sampled individuals from eight sites within Atlantic forest remnants with different amounts of available forest habitat and assessed FA levels with geometric morphometric techniques using adult mandibles. We used observed heterozygosity (Ho) and inbreeding coefficient (Fis), from seven microsatellite markers, as a proxy of genetic variation at individual and population levels. Populations were not significantly different for shape or size FA levels. Furthermore, interindividual variation in both shape and size FA levels and interpopulational differences in size FA levels were best explained by chance. However, habitat amount was negatively associated with both interpopulational variance and average shape FA levels. This association was stronger in populations living in areas with <28% of forest cover, which presented higher variance and higher average FA, suggesting that Nectomys squamipes might have a tolerance threshold to small availability of habitat. Our work is one of the first to use FA to address environmental stress caused by habitat loss in small mammal populations from a Neotropical biome. We suggest that shape FA might serve as a conservation tool to monitor human impact on natural animal populations.

Epifaunal invertebrate assemblages associated with branching Pocilloporids in Moorea, French Polynesia.

Reef-building corals can harbour high abundances of diverse invertebrate epifauna. Coral characteristics and environmental conditions are important drivers of community structure of coral-associated invertebrates; however, our current understanding of drivers of epifaunal distributions is still unclear. This study tests the relative importance of the physical environment (current flow speed) and host quality (e.g., colony height, surface area, distance between branches, penetration depth among branches, and background partial mortality) in structuring epifaunal communities living within branching Pocillopora colonies on a back reef in Moorea, French Polynesia. A total of 470 individuals belonging to four phyla, 16 families and 39 genera were extracted from 36 Pocillopora spp. colonies. Decapods were the most abundant epifaunal organisms (accounting for 84% of individuals) found living in Pocillopora spp. While coral host characteristics and flow regime are very important, these parameters were not correlated with epifaunal assemblages at the time of the study. Epifaunal assemblages associated with Pocillopora spp. were consistent and minimally affected by differences in host characteristics and flow regime. The consistency in abundance and taxon richness among colonies (regardless of habitat characteristics) highlighted the importance of total habitat availability. With escalating effects of climate change and other localized disturbances, it is critical to preserve branching corals to support epifaunal communities.

Planning forest restoration within private land holdings with conservation co-benefits at the landscape scale.

Forest loss is mainly due to the conversion of forest to agriculture, mostly in private lands. Forest restoration is a global priority, yet restoration targets are ambitious and budget-limited. Therefore, assessing the outcome of alternative decisions on land-use within private lands is paramount to perform cost-effective restoration. We present a novel framework that incorporates spatial planning for forest restoration within private lands with conservation co-benefits at the landscape scale. As a case study, we used three real landscapes of 10.000 ha with differing amounts of forest cover in the Atlantic Forest region of Brazil, and three hypothetical animal species with different dispersal abilities. We estimated the total amount of forest that landholders must restore to comply with the Native Vegetation Protection Law, which requires landholders to reforest 20% of their land within a 20-year time frame. We compared the cost-effectiveness of five restoration strategies based on the improvement in habitat availability and restoration costs. The most cost-effective strategy depends on a landscape's initial amount of forest cover and the species of concern. We revealed that spatial planning for restoration in private lands increased habitat availability up to 12 times more than random restoration, which was always the least cost-effective strategy. Cost-effective large-scale restoration in Brazil depends on public policies that assist landholders to comply with the law and on prioritizing areas for restoration within private lands. We show that by adding habitat availability as target in spatial prioritization, benefits for biodiversity can be hastened at low additional cost, even in real world scenarios with severe spatial constraints. Despite constraints, spatially planned restoration for law compliance in Brazil increased landscape permeability by creating corridors and stepping stones. Our framework should be used to plan restoration in Brazilian private lands and can be customized for other regions worldwide.

Evaluating the effects of landscape structure on the recovery of an invasive vertebrate after population control.

CONTEXT: Effective landscape control of invasive species is context-dependent due to the interplay between the landscape structure, local population dynamics, and metapopulation processes. We use a modelling approach incorporating these three elements to explore the drivers of recovery of populations of invasive species after control. OBJECTIVES: We aim to improve our understanding of the factors influencing the landscape-level control of invasive species. METHODS: We focus on the case study of invasive brushtail possum (Trichosurus vulpecula) control in New Zealand. We assess how 13 covariates describing the landscape, patch, and population features influence the time of population recovery to a management density threshold of two possums/ha. We demonstrate the effects of those covariates on population recovery under three scenarios of population growth: logistic growth, strong Allee effects, and weak Allee effects. RESULTS: Recovery times were rapid regardless of the simulated population dynamics (average recovery time < 2 years), although populations experiencing Allee effects took longer to recover than those growing logistically. Our results indicate that habitat availability and patch area play a key role in reducing times to recovery after control, and this relationship is consistent across the three simulated scenarios. CONCLUSIONS: The control of invasive possum populations in patchy landscapes would benefit from a patch-level management approach (considering each patch as an independent management unit), whereas simple landscapes would be better controlled by taking a landscape-level view (the landscape as the management unit). Future research should test the predictions of our models with empirical data to refine control operations.

River network connectivity and fish diversity.

Frequent and severe disruptions of natural river flows associated with human activities significantly alter hydrological connectivity in large river networks, with deleterious effects on fish diversity. Understanding the relationship between fish diversity and river network connectivity is fundamental to ensuring species persistence, ecosystem integrity, and human well-being. Here, we provide a review of the mechanisms by which river network connectivity (RNC) affects fish diversity. We review the relationships between forms, systems and types of RNC and fish diversity, based on more than 100 previous studies. In summary, sustaining RNC promotes fish diversity in longitudinal and lateral axes, and species sorting, dispersal dynamics, and habitat availability are the main factors driving the distribution of fish diversity, followed by nutrition and trophic dynamics. Our work highlights the effects of RNC on fish diversity, and provides a mechanistic understanding of how RNC affects fish diversity across river basins, thus providing scientific guidance for protecting fish biodiversity and improving the health of river network ecosystems.

Spatial aggregation of aquatic habitats affects oviposition patterns in Aedes mosquitoes.

Colonization, including oviposition, is an important driver of population and community dynamics both within and across habitat patches. Most research has focused on the roles of habitat availability or quality on colonization and its outcomes. However, the spatial distribution of habitats also likely affects these processes. We conducted field experiments in Georgia, USA, using clustered and dispersed arrays of equal numbers of oviposition patches to investigate how patch aggregation influenced oviposition by Aedes mosquitoes. We tested the effects of aggregation on: (1) the total number of eggs an array received, (2) the proportion of patches within an array that received eggs, and (3) the number of eggs per colonized patch. We compared results to predictions from three models (Field of Dreams, Propagule Redirection, and Excess Attraction), which vary in the degree to which arrays attract colonists and apportion those colonists among patches. Clustered arrays received 22% more eggs than dispersed arrays, with clustered patches significantly more likely to receive eggs. At the species level, A. albopictus responded more to clustering than did A. triseriatus. These results are inconsistent with Propagule Redirection, but support the Excess Attraction and Field of Dreams models. Although clustered arrays occupied a relatively small area, they attracted at least as many ovipositing mosquitoes as did dispersed arrays. However, the number of eggs per colonized patch did not differ between clustered and dispersed arrays. Therefore, density dependence among larvae, and hence the production of adult mosquitoes on a per-patch basis, should be similar in dispersed and clustered landscapes.

Functional responses in habitat selection: clarifying hypotheses and interpretations.

A fundamental challenge in habitat ecology and management is understanding the mechanisms generating animal distributions. Studies of habitat selection provide a lens into such mechanisms, but are often limited by unrealistic assumptions. For example, most studies assume that habitat selection is constant with respect to the availability of resources, such that habitat use remains proportional to availability. To the contrary, a growing body of work has shown the fallacy of this assumption, indicating that animals modify their behavior depending on the context at broader scales. This has been termed a functional response in habitat selection. Furthermore, a diversity of methods is employed to model functional responses in habitat selection, with little attention to how methodology might affect scientific and conservation conclusions. Here, we first review the conceptual and statistical foundations of methods currently used to model functional responses and clarify the ecological tests evaluated within each approach. We then use a combination of simulated and empirical data sets to evaluate the similarities and differences among approaches. Importantly, we identified multiple statistical issues with the most widely applied approaches to understand functional responses, including: (1) a complex and important role of random- or individual-level intercepts in adjusting individual-level regression coefficients as resource availability changes and (2) a sensitivity of results to poorly informed individual-level coefficients estimated for animals with low availability of a given resource. Consequently, we provide guidance on applying approaches that are insensitive to these issues with the goal of advancing our understanding of animal habitat ecology and management. Finally, we characterize the management implications of assuming similarity between the current approaches to model functional responses with two empirical data sets of federally threatened species: Canada lynx (Lynx canadensis) in the United States and woodland caribou (Rangifer tarandus caribou) in Canada. Collectively, our assessment helps clarify the similarities and differences among current approaches and, therefore, assists the integration of functional responses into the mainstream of habitat ecology and management.

Genetic diversity and gene flow of the threatened Brazilian endemic parrotfish Scarus trispinosus (Valenciennes, 1840).

The greenback parrotfish, Scarus trispinosus, is the largest herbivorous fish inhabiting Southwestern Atlantic reefs, and was recently included in the IUCN red list of threatened species as endangered due to the overexploitation of their populations. The aim of this work was to evaluate the existence of structured populations (i.e. genetic unities) along a coast of approximately 2000 km of the NE Brazilian coast. The transferability of 17 primers synthesized for Scarus rubroviolaceus was tested for S. trispinosus and five transferable loci were validated and used. Two localities within the Abrolhos Bank, off the Central Brazilian coast (Corumbau and Caravelas) and in close proximity to the MPA, which encompasses the largest remnants of the S. trispinosus population, exhibited higher levels of genetic richness. Remaining locations, Pernambuco, Porto Seguro and Rio Grande do Norte exhibited lower genetic diversity. We found no genetic differences among sampled localities however, when those samples were gathered into latitudinal groups (northern vs southern) a subtle but significant genetic substructuring was revealed. It is proposed that a combination of high local individual admixture favoured by habitat connectivity drived genetic homogeneity at regional scales while larval dispersal contributed to heterogeneities observed at large scales maintaining gene flow through oceanographic currents.

Node self-connections and metapopulation persistence: reply to Saura (2018).

Saura () claims that studies using the Probability of Connectivity metric (PC) had already demonstrated the importance of including node self-connections in network metrics. As originally defined and used, PC cannot test the importance of self-connections. However, with key terms redefined, PC could be a useful tool in future work.

Node self-connections in network metrics.

Zamborain-Mason et al. (Ecol. Lett., 20, 2017, 815-831) state that they have newly proposed network metrics that account for node self-connections. Network metrics incorporating node self-connections, also referred to as intranode (intrapatch) connectivity, were however already proposed before and have been widely used in a variety of conservation planning applications.

No place like home? A test of the natal habitat-biased dispersal hypothesis in Scandinavian wolves.

Natal dispersal is an important mechanism for the viability of populations. The influence of local conditions or experience gained in the natal habitat could improve fitness if dispersing individuals settle in an area with similar habitat characteristics. This process, defined as 'natal habitat-biased dispersal' (NHBD), has been used to explain distribution patterns in large carnivores, but actual studies evaluating it are rare. We tested whether grey wolf Canis lupus territory establishment was influenced by the habitat characteristics of the natal territory using the long-term monitoring of the Scandinavian wolf population. We paired the locations of natal and established territories, accounted for available habitats along the dispersing route, and compared their habitat characteristics for 271 wolves during 1998-2012. Wolves with the shortest dispersal distances established in natal-like habitat types more than expected by chance, whereas wolves that dispersed longer distances did not show NHBD. The pattern was consistent for male and female wolves, with females showing more NHBD than males. Chances to detect NHBD increased with the size of habitat defined as available. This highlights the importance of considering the biological characteristics of the studied species when defining habitat availability. Our methodological approach can prove useful to inform conservation and management to identify habitats to be selected by reintroduced or naturally expanding populations.

Habitat availability drives the distribution-abundance relationship in phytophagous true bugs in managed grasslands.

The nearly universal positive relationship between the distribution and abundance of species has been explained by several hypotheses but hitherto no consensus has been reached. Here, we used monitoring data of 105 phytophagous true bug species (Heteroptera) from 150 grassland sites over six years to test how (1) range position, (2) resource use, (3) resource availability, (4) density-dependent habitat selection, (5) metapopulation dynamics, and (6) habitat dispersal affect the distribution-abundance relationship. For the use in a confirmatory path analysis, we constructed causal pathways representing the hypothesized relationships and tested them separately and in a combined analysis. Our results show that the distribution-abundance relationship in phytophagous true bugs is driven by habitat-availability. An increasing local density of the host-plants increases the distribution of the species in the landscape, which in turn increases their local abundance. Thereby habitat availability facilitates dispersal success. We conclude that local abundance of herbivores facing habitat destruction could decline owing to a decrease in population dynamics between sites at the landscape scale. Finally, our results underline the potential of confirmatory path analysis for testing competing hypotheses.

Pulsed food resources, but not forest cover, determine lifetime reproductive success in a forest-dwelling rodent.

The relative contributions of habitat and food availability on fitness may provide evidence for key habitat features needed to safeguard population persistence. However, defining habitat quality for a species can be a complex task, especially if knowledge on the relationship between individual performance and habitat quality is lacking. Here, we determined the relative importance of the availability of suitable forest habitat, body mass and food from masting tree species on female lifetime reproductive success (LRS) of Siberian flying squirrels (Pteromys volans). We calculated LRS of 500 female flying squirrels based on a 22-year-long longitudinal dataset of two populations from western Finland. We assessed with generalised additive models the potential effects of availability of suitable habitat and cumulative lifetime availability of food from masting tree species on female LRS, longevity and fecundity. On a reduced dataset, we evaluated the importance of female winter body mass and conducted a piecewise path analysis to determine how variables were connected. According to generalised additive models female longevity, fecundity and LRS were mainly determined by variation in cumulative lifetime availability of food from masting alder and birch. Instead, habitat and body mass had a smaller role. The path analysis indicated that lifetime food availability had a direct effect on longevity and fecundity, and these had an equal effect on LRS at both study sites. Our results on LRS show that the occurrence of tree masting events during a female flying squirrel's lifetime has a profoundly larger effect on LRS than the cover of suitable forest habitat. Furthermore, this study emphasises the importance of both fecundity and longevity, and the indirect effects of food availability via those components, as determinants of lifetime fitness in female flying squirrels.

Shallow rocky nursery habitat for fish: Spatial variability of juvenile fishes among this poorly protected essential habitat.

Coastal nursery habitats are essential for the renewal of adult fish populations. We quantified the availability of a coastal nursery habitat (shallow heterogeneous rocky bottoms) and the spatial variability of its juvenile fish populations along 250km of the Catalan coastline (France and Spain). Nurseries were present in 27% of the coastline, but only 2% of them benefited from strict protection status. For nine taxa characteristic of this habitat, total juvenile densities varied significantly between nursery sites along the coastline, with the highest densities being found on the northern sites. Recruitment level (i.e. a proxy of nursery value) was not explained by protection level, but it was moderately and positively correlated with an anthropization index. Patterns of spatial variations were taxa-specific. Exceptional observations of four juveniles of the protected grouper Epinephelus marginatus were recorded. Our data on habitat availability and recruitment levels provides important informations which help to focus MPA management efforts.

Persistence at distributional edges: Columbia spotted frog habitat in the arid Great Basin, USA.

A common challenge in the conservation of broadly distributed, yet imperiled species is understanding which factors facilitate persistence at distributional edges, locations where populations are often vulnerable to extirpation due to changes in climate, land use, or distributions of other species. For Columbia spotted frogs (Rana luteiventris) in the Great Basin (USA), a genetically distinct population segment of conservation concern, we approached this problem by examining (1) landscape-scale habitat availability and distribution, (2) water body-scale habitat associations, and (3) resource management-identified threats to persistence. We found that areas with perennial aquatic habitat and suitable climate are extremely limited in the southern portion of the species' range. Within these suitable areas, native and non-native predators (trout and American bullfrogs [Lithobates catesbeianus]) are widespread and may further limit habitat availability in upper- and lower-elevation areas, respectively. At the water body scale, spotted frog occupancy was associated with deeper sites containing abundant emergent vegetation and nontrout fish species. Streams with American beaver (Castor canadensis) frequently had these structural characteristics and were significantly more likely to be occupied than ponds, lakes, streams without beaver, or streams with inactive beaver ponds, highlighting the importance of active manipulation of stream environments by beaver. Native and non-native trout reduced the likelihood of spotted frog occupancy, especially where emergent vegetation cover was sparse. Intensive livestock grazing, low aquatic connectivity, and ephemeral hydroperiods were also negatively associated with spotted frog occupancy. We conclude that persistence of this species at the arid end of its range has been largely facilitated by habitat stability (i.e., permanent hydroperiod), connectivity, predator-free refugia, and a commensalistic interaction with an ecosystem engineer. Beaver-induced changes to habitat quality, stability, and connectivity may increase spotted frog population resistance and resilience to seasonal drought, grazing, non-native predators, and climate change, factors which threaten local or regional persistence.