Predicting the Occurrence of Cave-Inhabiting Fauna Based on Features of the Earth Surface Environment

Mountain caves of the central region of Veracruz: A vertebrate biodiversity reservoir in a Neotropical hotspot

The mountain region of central Veracruz, Mexico hosts a large system of karst and volcanic caves that are unexplored. In particular, the vertebrates that inhabit these subterranean ecosystems are unknown. This study evaluated the diversity of mammals, birds, reptiles, amphibians, and fish in three environments (euphotic, disphotic, and aphotic) of 16 caves of different geological origin (12 karst caves and 4 volcanic caves) distributed along an altitudinal gradient (300-2400 m a.s.l.). We found a richness of 242 vertebrate species (184 birds, 30 mammals, 15 reptiles, 12 amphibians, and 1 fish) and an abundance of a total of 11,323 individuals (4,969 mammals, 6,483 birds, 36 reptiles, 27 amphibians, and 5 fish). The richness of all vertebrate classes was higher in karst than in volcanic caves. Vertebrate diversity was also higher at mid-altitudes between 600-899 m a.s.l. Diversity varied between environments, where bird and reptile richness was higher in the euphotic environment, while mammal and amphibian diversity was higher in the aphotic environment. The similarity in the composition of vertebrate species does not depend on the distance between karstic and volcanic caves. Volcanic and karst caves shared on average up to 70% and 55% of vertebrate species, which indicates that only 30% and 45% of species, respectively, is different in each cave type. Given the vulnerability and fragility of these subterranean ecosystems, as well as the important diversity that they contain, we recommend including the caves of the central region of Veracruz in the conservation agenda of local governments and communities. Community-based conservation can help ensure the presence of vertebrate species in the caves of this region.

Extending beyond individual caves: a graph theory approach broadening conservation priorities in Amazon iron ore caves

The Amazon is renowned worldwide for its biological significance, but it also harbors substantial mineral reserves. Among these, the ferruginous geosystems of the region are critical for iron ore extraction, accounting for 10% of Brazil's export revenue. Additionally, this region holds a significant speleological heritage with more than 1,000 caves. However, cave conservation efforts are often in conflict with land use, necessitating mediation through environmental regulations. While conservation decisions typically consider only the caves' characteristics, such an approach fails to account for the interactions among cave communities and their surrounding landscape. This poses a challenge to reserve design for cave conservation purposes. To address this issue, we assessed the predictors that influence the similarity among cave communities, suggesting the use of this parameter as a proxy for subterranean connectivity. Applying graph theory, we proposed a tool to aid in the selection of priority caves for conservation purposes. Our study involved the sampling of invertebrates in 69 iron ore caves and analyzing 28 environmental variables related to these subterranean habitats and adjacent landscape. Our analysis revealed that landscape and habitat characteristics are more important than geographical distance in determining patterns of similarity among caves. Our graph approach highlighted densely interconnected clusters based on similarity. However, specific caves stood out for harboring exclusive fauna and/or exhibiting habitat specificity, making them unique in the study area. Thus, we recommend prioritizing cave clusters for conservation, assembling both singular caves and others that influence them. It is crucial to note that protocols for the protection of subterranean biodiversity must consider measures that encompass both the caves and the surrounding landscape. Our methodology provides insights into the connectivity among caves, identifies existing groups, highlights singular (or unique) cavities that require preservation, and recognizes those influencing these unique habitats. This methodological advancement is crucial for the development of better conservation policies for the speleological heritage in areas under constant economic pressure.

Perspectives and pitfalls in preserving subterranean biodiversity through protected areas

Subterranean ecosystems (comprising terrestrial, semi-aquatic, and aquatic components) are increasingly threatened by human activities; however, the current network of surface-protected areas is inadequate to safeguard subterranean biodiversity. Establishing protected areas for subterranean ecosystems is challenging. First, there are technical obstacles in mapping three-dimensional ecosystems with uncertain boundaries. Second, the rarity and endemism of subterranean organisms, combined with a scarcity of taxonomists, delays the accumulation of essential biodiversity knowledge. Third, establishing agreements to preserve subterranean ecosystems requires collaboration among multiple actors with often competing interests. This perspective addresses the challenges of preserving subterranean biodiversity through protected areas. Even in the face of uncertainties, we suggest it is both timely and critical to assess general criteria for subterranean biodiversity protection and implement them based on precautionary principles. To this end, we examine the current status of European protected areas and discuss solutions to improve their coverage of subterranean ecosystems.

Lithology and disturbance drive cavefish and cave crayfish occurrence in the Ozark Highlands ecoregion

Diverse communities of groundwater-dwelling organisms (i.e., stygobionts) are important for human wellbeing; however, we lack an understanding of the factors driving their distributions, making it difficult to protect many at-risk species. Therefore, our study objective was to determine the landscape factors related to the occurrence of cavefishes and cave crayfishes in the Ozark Highlands ecoregion, USA. We sampled cavefishes and cave crayfishes at 61 sampling units using both visual and environmental DNA surveys. We then modeled occurrence probability in relation to lithology and human disturbance while accounting for imperfect detection. Our results indicated that occurrence probability of cave crayfishes was negatively associated with human disturbance, whereas there was a weak positive relationship between cavefish occurrence and disturbance. Both cavefishes and cave crayfishes were more likely to occur in limestone rather than dolostone lithology. Our results indicate structuring factors are related to the distribution of these taxa, but with human disturbance as a prevalent modifier of distributions for cave crayfishes. Limiting human alteration near karst features may be warranted to promote the persistence of some stygobionts. Moreover, our results indicate current sampling efforts are inadequate to detect cryptic species; therefore, expanding sampling may be needed to develop effective conservation actions.

Identifying priority areas for bat conservation in the Western Ghats mountain range, peninsular India

Understanding patterns of species distribution and diversity plays a vital role in biodiversity conservation. Such documentation is frequently lacking for bats, which are relatively little studied and often threatened. The Western Ghats biodiversity hotspot in peninsular India is a bat hotspot with 63 species. We conducted a comprehensive bat survey across the southern Western Ghats and used maximum entropy modeling (MaxEnt) to model the potential distribution of 37 bat species for which sufficient data were available. We generated binary maps of each species using species-specific thresholds to estimate suitable habitat areas and overlaid binary maps of species to produce bat hotspots (we use the term “bat hotspot” for regions that were suitable for more than 25 bat species). We also estimated species richness across protected area networks in the southern Western Ghats to assess the level of protection. The highest levels of species richness were found mainly along the southmost Periyar–Agastyamalai landscape. The study also identified a 1,683 km2 area of potential bat hotspot and 726 km2 (43%) of the total bat hotspots are currently within the protected area network. However, more than 50% of suitable habitats for each of the 37 species remain unprotected. Therefore, conservation decisions are needed to take into account both bat hotspots and species with restricted distributions.

An Overview of Subterranean Biodiversity Hotspots

Riding a wave of interest in biodiversity patterns in surface-dwelling communities, in 2000, Culver and Sket [...]

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.

Extreme genetic diversity among springtails (Collembola) in subterranean calcretes of arid Australia

The subterranean islands hypothesis for calcretes of the Yilgarn region in Western Australia applies to many stygobitic (subterranean-aquatic) species that are "trapped" evolutionarily within isolated aquifers due to their aquatic lifestyles. In contrast, little is known about the distribution of terrestrial-subterranean invertebrates associated with the calcretes. We used subterranean Collembola from the Yilgarn calcretes to test the hypothesis that troglobitic species, those inhabiting the subterranean unsaturated (non-aquatic) zone of calcretes, are also restricted in their distribution and represent reciprocally monophyletic and endemic lineages. We used the barcoding fragment of the mtDNA cytochrome c oxidase subunit 1 (COI) gene from 183 individuals to reconstruct the phylogenetic history of the genus Pseudosinella Schäffer (Collembola, Lepidocyrtidae) from 10 calcretes in the Yilgarn. These calcretes represent less than 5% of the total possible calcretes in this region, yet we show that their diversity for subterranean Collembola comprises a minimum of 25 new species. Regionally, multiple levels of diversity exist in Pseudosinella, indicative of a complex evolutionary history for this genus in the Yilgarn. These species have probably been impacted by climatic oscillations, facilitating their dispersal across the landscape. The results represent a small proportion of the undiscovered diversity in Australia's arid zone.

Hotspots of Subterranean Biodiversity

Worldwide, caves and groundwater habitats harbor thousands of species modified and limited to subterranean habitats in karst. Data are concentrated in Europe and USA, where a number of detailed analyses have been performed. Much less is known with respect to global patterns due to a lack of data. This special issue will focus on and discuss the global patterns of individual hotspot caves and groundwater habitats.

Geographically structured genetic diversity in the cave beetle Darlingtonea kentuckensis Valentine, 1952 (Coleoptera, Carabidae, Trechini, Trechina)

Cave beetles of the eastern USA are one of many poorly studied groups of insects and nearly all previous work delimiting species is based solely on morphology. This study assesses genetic diversity in the monotypic cave carabid beetle genus DarlingtoneaValentine 1952, to test the relationship between putative geographical barriers to subterranean dispersal and the boundaries of genetically distinct groups. Approximately 400bp of the mitochondrial cytochrome oxidase I (COI) gene was sequenced from up to four individuals from each of 27 populations, sampled from caves along the escarpments of the Mississippian and Cumberland plateaus in eastern Kentucky, USA. The 81 individuals sequenced yielded 28 unique haplotypes. Hierarchical analyses of molecular variance (AMOVA) within and among geographically defined groups tested two a priori hypotheses of structure based on major and minor river drainages, as well as genetic distance clusters defined a posteriori from an unrooted analysis. High genetic differentiation (FST) between populations was found across analyses. The influence of isolation by distance could potentially account for much but not all of the variation found among geographically defined groups at both levels. High variability among the three northernmost genetic clusters (FCT), low variability among populations within clusters (FSC), and low within-cluster Mantel correlations indicate the importance of unidentified likely intra-karst barriers to gene flow separating closely grouped cave populations. Overall phylogeographic patterns are consistent with previous evidence of population isolation among cave systems in the region, revealing geographically structured cryptic diversity in Darlingtonea over its distribution. The landscape features considered a priori in this study were not predictive of the genetic breaks among the three northern clusters, which are genetically distinct despite their close geographic proximity.

Local- versus broad-scale environmental drivers of continental β-diversity patterns in subterranean spider communities across Europe

Macroecologists seek to identify drivers of community turnover (β-diversity) through broad spatial scales. However, the influence of local habitat features in driving broad-scale β-diversity patterns remains largely untested, owing to the objective challenges of associating local-scale variables to continental-framed datasets. We examined the relative contribution of local- versus broad-scale drivers of continental β-diversity patterns, using a uniquely suited dataset of cave-dwelling spider communities across Europe (35-70° latitude). Generalized dissimilarity modelling showed that geographical distance, mean annual temperature and size of the karst area in which caves occurred drove most of β-diversity, with differential contributions of each factor according to the level of subterranean specialization. Highly specialized communities were mostly influenced by geographical distance, while less specialized communities were mostly driven by mean annual temperature. Conversely, local-scale habitat features turned out to be meaningless predictors of community change, which emphasizes the idea of caves as the human accessible fraction of the extended network of fissures that more properly represents the elective habitat of the subterranean fauna. To the extent that the effect of local features turned to be inconspicuous, caves emerge as experimental model systems in which to study broad biological patterns without the confounding effect of local habitat features.

Contribution of rare and common species to subterranean species richness patterns

AIM

Common species contribute more to species richness patterns (SRPs) than rare species in most studies. Our aim was to test this hypothesis using a novel model system, species living exclusively in subterranean habitats. They consist of mainly rare species (small ranges), only a few of them being common (large ranges), and challenge whether rare species are less important for the development of SRPs in this environment. We separately analyzed aquatic and terrestrial species.

LOCATION

Western Balkans in southeastern Europe.

METHODS

We assembled two datasets comprising 431 beetle and 145 amphipod species, representing the model groups of subterranean terrestrial and aquatic diversity, respectively. We assessed the importance of rare and common species using the stepwise reconstruction of SRPs and subsequent correlation analyses, corrected also for the cumulative information content of the subsets based on species prevalence. We applied generalized linear regression models to evaluate the importance of rare and common species in forming SRPs. Additionally, we analyzed the contribution of rare and common species in species-rich cells.

RESULTS

Patterns of subterranean aquatic and terrestrial species richness overlapped only weakly, with aquatic species having larger ranges than terrestrial ones. Our analyses supported higher importance of common species for forming overall SRPs in both beetles and amphipods. However, in stepwise analysis corrected for information content, results were ambiguous. Common species presented a higher proportion of species than rare species in species-rich cells.

MAIN CONCLUSION

We have shown that even in habitats with the domination of rare species, it is still common species that drive SRPs. This may be due to an even spatial distribution of rare species or spatial mismatch in hotspots of rare and common species. SRPs of aquatic and terrestrial subterranean organisms overlap very little, so the conservation approaches need to be habitat specific.

Shedding light on the hidden world of subterranean fauna: A transdisciplinary research approach

Subterranean environments contain a diverse and unique obligate fauna: either aquatic living in the groundwater or terrestrial living in voids above the water table. In the arid region of the western part of the Australian continent, a particularly rich subterranean fauna coincides with a concentration of natural resource extraction operations. Since the inclusion of subterranean fauna in assessments of environmental impact in the mid-1990s, taxonomic research in Australia on this group of mainly invertebrates has grown exponentially. However, remaining knowledge gaps continue to frustrate both environmental regulators and development proponents due to high uncertainty in the decision-making process. In early 2017, the Western Australian Biodiversity Science Institute was tasked with leading the development of a research program to improve on the current state of knowledge of subterranean fauna. To balance the diverse environmental, economic and social needs of a range of stakeholders, transdisciplinary principles were applied to program development. A clear consensus on five broad focus areas to progress include: (1) data consolidation; (2) resilience to disturbance; (3) survey and sampling protocols; (4) abiotic and biotic habitat requirements; and (5) species delineation. In the context of these focus areas; we describe the research program development, reviewing the status of knowledge within each focus area, and the research initiatives to close the gaps in knowledge. We argue that, by adopting a transdisciplinary approach, the likelihood of success of the research program, as measured by the effective translation and adoption of research findings, will be maximized. This review is timely given the ever-increasing demand on groundwater systems for water extraction worldwide. A holistic understanding of the influence of anthropogenic activities on these ecosystems, and the functional role of organisms within them, will help to ensure that their health is not compromised.

Exploring the Interplay Between Local and Regional Drivers of Distribution of a Subterranean Organism

Caves are excellent model systems to study the effects of abiotic factors on species distributions due to their selective conditions. Different ecological factors have been shown to affect species distribution depending on the scale of analysis, whether regional or local. The interplay between local and regional factors in explaining the spatial distribution of cave-dwelling organisms is poorly understood. Using the troglophilic subterranean spider Artema nephilit (Araneae: Pholcidae) as a model organism, we investigated whether similar environmental predictors drive the species distribution at these two spatial scales. At the local scale, we monitored the abundance of the spiders and measured relevant environmental features in 33 caves along the Jordan Rift Valley. We then extended the analysis to a regional scale, investigating the drivers of the distribution using species distribution models. We found that similar ecological factors determined the distribution at both local and regional scales for A. nephilit. At a local scale, the species was found to preferentially occupy the outermost, illuminated, and warmer sectors of caves. Similarly, mean annual temperature, annual temperature range, and solar radiation were the most important drivers of its regional distribution. By investigating these two spatial scales simultaneously, we showed that it was possible to achieve an in-depth understanding of the environmental conditions that governs subterranean species distribution.

Conserving relics from ancient underground worlds: assessing the influence of cave and landscape features on obligate iron cave dwellers from the Eastern Amazon

The degradation of subterranean habitats is believed to represent a serious threat for the conservation of obligate subterranean dwellers (troglobites), many of which are short-range endemics. However, while the factors influencing cave biodiversity remain largely unknown, the influence of the surrounding landscape and patterns of subterranean connectivity of terrestrial troglobitic communities have never been systematically assessed. Using spatial statistics to analyze the most comprehensive speleological database yet available for tropical caves, we first assess the influence of iron cave characteristics and the surrounding landscape on troglobitic communities from the Eastern Amazon. We then determine the spatial pattern of troglobitic community composition, species richness, phylogenetic diversity, and the occurrence of frequent troglobitic species, and finally quantify how different landscape features influence the connectivity between caves. Our results reveal the key importance of habitat amount, guano, water, lithology, geomorphology, and elevation in shaping iron cave troglobitic communities. While mining within 250 m from the caves influenced species composition, increasing agricultural land cover within 50 m from the caves reduced species richness and phylogenetic diversity. Troglobitic species composition, species richness, phylogenetic diversity, and the occurrence of frequent troglobites showed spatial autocorrelation for up to 40 km. Finally, our results suggest that the conservation of cave clusters should be prioritized, as geographic distance was the main factor determining connectivity between troglobitic communities. Overall, our work sheds important light onto one of the most overlooked terrestrial ecosystems, and highlights the need to shift conservation efforts from individual caves to subterranean habitats as a whole.

Neutral colonisations drive high beta-diversity in cavernicole springtails (Collembola)

The theory of island biogeography predicts the effects of habitat isolation and size on species richness, community assembly, and the persistence of species. Various studies showed that habitat conditions and the ecology of species are also of relevance in explaining community assembly. Geographically isolated habitats like caves with rather constant environmental conditions provide models to test for the relevance of the above described variables. In this study we analysed springtails living in karst caves of the Romanian Carpathians and Dobrogea region. We considered phylogenetic relatedness, habitat and species characteristics to identify the relevant drivers of community assembly. Our data show that species richness of single caves is low. Neither phylogenetic relatedness nor habitat filtering and competitive interactions seem to shape species composition or to affect species richness. We found that glacial-interglacial cycles with subsequent range contractions and expansions might have led to independent and multiple colonisations of caves. Furthermore, single caves might have acted as refugia and thus might have provided the prerequisite for distinct evolution processes, leading to a high level of endemicity of these animal species.

Energy and speleogenesis: Key determinants of terrestrial species richness in caves

The aim of this study was to unravel the relative role played by speleogenesis (i.e., the process in which a cave is formed), landscape‐scale variables, and geophysical factors in the determination of species richness in caves. Biological inventories from 21 caves located in the southeastern Iberian Peninsula along with partial least square (PLS) regression analysis were used to assess the relative importance of the different explanatory variables. The caves were grouped according to the similarity in their species composition; the effect that spatial distance could have on similarity was also studied using correlation between matrices. The energy and speleogenesis of caves accounted for 44.3% of the variation in species richness. The trophic level of each cave was the most significant factor in PLS regression analysis, and epigenic caves (i.e., those formed by the action of percolating water) had significantly more species than hypogenic ones (i.e., those formed by the action of upward flows in confined aquifers). Dissimilarity among the caves was very high (multiple‐site β_sim = 0.92). Two main groups of caves were revealed through the cluster analysis, one formed by the western caves and the other by the eastern ones. The significant—but low—correlation found between faunistic dissimilarity and geographical distance (r = .16) disappeared once the caves were split into the two groups. The extreme beta‐diversity suggests a very low connection among the caves and/or a very low dispersal capacity of the species. In the region under study, two main factors are intimately related to the richness of terrestrial subterranean species in caves: the amount of organic material (trophic level) and the formation process (genesis). This is the first time that the history of a cave genesis has been quantitatively considered to assess its importance in explaining richness patterns in comparison with other factors more widely recognized.

Predicting the Occurrence of Cave-Inhabiting Fauna Based on Features of the Earth Surface Environment

One of the most challenging fauna to study in situ is the obligate cave fauna because of the difficulty of sampling. Cave-limited species display patchy and restricted distributions, but it is often unclear whether the observed distribution is a sampling artifact or a true restriction in range. Further, the drivers of the distribution could be local environmental conditions, such as cave humidity, or they could be associated with surface features that are surrogates for cave conditions. If surface features can be used to predict the distribution of important cave taxa, then conservation management is more easily obtained. We examined the hypothesis that the presence of major faunal groups of cave obligate species could be predicted based on features of the earth surface. Georeferenced records of cave obligate amphipods, crayfish, fish, isopods, beetles, millipedes, pseudoscorpions, spiders, and springtails within the area of Appalachian Landscape Conservation Cooperative in the eastern United States (Illinois to Virginia and New York to Alabama) were assigned to 20 x 20 km grid cells. Habitat suitability for these faunal groups was modeled using logistic regression with twenty predictor variables within each grid cell, such as percent karst, soil features, temperature, precipitation, and elevation. Models successfully predicted the presence of a group greater than 65% of the time (mean = 88%) for the presence of single grid cell endemics, and for all faunal groups except pseudoscorpions. The most common predictor variables were latitude, percent karst, and the standard deviation of the Topographic Position Index (TPI), a measure of landscape rugosity within each grid cell. The overall success of these models points to a number of important connections between the surface and cave environments, and some of these, especially soil features and topographic variability, suggest new research directions. These models should prove to be useful tools in predicting the presence of species in understudied areas.