Habitats and Diversity of Subterranean Macroscopic Freshwater Invertebrates: Main Gaps and Future Trends

: Caves are the best studied aquatic subterranean habitat, but there is a wide variety of these habitats, ranging in depth below the surface and size of the spaces (pore or habitat size). Both factors are important in setting limits to species composition and richness. In addition to caves, among the most important shallow aquatic subterranean habitats are the hyporheal (underﬂow of rivers and streams), the hypotelminorheal (very superﬁcial drainages with water exiting in seeps), epikarst, and calcrete aquifers. Although it is little studied, both body size and species composition in the di ﬀ erent habitats is di ﬀ erent. Because of high levels of endemism and di ﬃ culty in access, no subterranean habitats are well sampled, even caves. However, there are enough data for robust generalizations about some geographic patterns. Individual hotspot caves are concentrated in the Dinaric region of southern Europe, and overall, tropical regions have fewer obligate aquatic cave dwellers (stygobionts). In all subterranean aquatic habitats, regional diversity is much higher than local diversity, but local diversity (especially single cave diversity) may be a useful predictor of regional species richness. In Europe there is a ridge of high aquatic subterranean species richness basically extending east from the French–Spanish border. Its cause may be either high productivity or that long-term temperature oscillations are at a minimum. With increased collecting and analysis, global and continental trends should become clearer.


Introduction
The existence of eyeless and depigmented animals in the darkness of cave streams has been known since at least 1537 [1]. Of course, caves themselves have been known since at least the Paleolithic [2]. Groundwater also has a long history of human knowledge and use, and the first known well dates back to the Neolithic [3]. On the biological side, in 1907 Racovitza [4] mentioned studies from the 1890s on eyeless, depigmented species from artesian wells in Texas and wells in the Canterbury Plain of New Zealand.
It is probably Racovitza [4] who first pointed to the potential advantages of considering all terrestrial and freshwater aphotic habitats and their inhabitants together. Earlier, North American neo-Lamarckians had a strong interest in the evolution of eyelessness due to its connection to evolution by disuse (e.g., Packard [5]). While there was a general recognition of the unity of the subterranean domain, especially by European scientists, e.g., Ginet and Decou [6], it was not until the publication of the groundbreaking book "Groundwater Ecology" [7] that at least the aquatic side of subterranean biology became a distinct discipline with its own tenets.
In this review, we consider three areas. The first is the range of aquatic subterranean habitats and what unites and divides them. There have been two overarching classifications of freshwater subterranean habitats. One, due to Botosaneanu [8] who divided aquatic subterranean habitats (he included the marine interstitial) into two branches-"milieu perméables en grande" and "milieu perméables en petit". He included more than 30 habitats, and the two most familiar are caves and the hyporheic, the underflow of streams and rivers. The other is due to Culver and Pipan [9,10] in which they attach equal, if not greater, weight to the vertical depth of the habitat compared to habitat size. We reassess these divisions and revisit the range of aquatic subterranean habitats.
The second area we review is the vexing question: how complete is our knowledge of the subterranean biota? In addition to the nearly universal shortfall of taxonomists to describe new species that have already been collected (the Linnean impediment), there are regions of the Earth's surface that have been little sampled for aquatic subterranean fauna, especially in the tropics (the Wallacean impediment), and there are subterranean habitats that are poorly sampled due to sampling difficulties (the Racovitzan impediment [11]). All of these problems are especially severe in aquatic subterranean habitats because of the high levels of endemism of the species [12,13]. What can we say about overall fauna patterns in light of this uncertainty?
The third area we review is: what is the global pattern of richness of aquatic subterranean invertebrates? After reviewing the pattern of individual hotspot sites globally, we review some of the potential "fixes" to incomplete data, especially functions such as those of Chao [14] that estimate missing species, as well as assessing the potential impact of missing data (e.g., Zagmajster et al. [15]).

Aquatic Subterranean Habitats
Botosaneanu's [8] compilation of subterranean habitats, produced as an addendum to his extensive review of the subterranean fauna is an appropriate place to begin (Table 1). We separated out from his dichotomy of "perméable en grand" and "perméable en petit" the categories of springs, which he included in "perméable en petit". Springs are both habitats themselves (ecotones between surface and subterranean waters), but also convenient collecting points for subterranean species from more inaccessible (often unknown) subterranean habitats. These habitats may have fine or coarse sediments and may be shallow or deep.
Botosaneanu suggested only species limited to subterranean habitats be considered, but in practice it makes more sense to consider all species in springs (e.g., [16,17]), especially since not only species limited to subterranean habitats are blind and eyeless [18]. Springs themselves can be classified in a number of different ways, the oldest being based on discharge rate [19]. Other classifications, which can be quite elaborate (see Springer and Stevens [20]), are based on characteristics of the hydraulic head, geomorphologic structure, and water quality and temperature [21]. All of these are more elaborate than Botosaneanu's [8] (Table 1). He listed springs as porous habitats [8] but we give them a separate category. Perhaps the reason that speleobiologists have not taken up more elaborate subdivisions of springs is that the fauna of springs often has few, if any, species showing the characteristics of subterranean life such as reduced eyes and pigment (see Botosaneanu [22]). However, this is not always the case. Dumnicka and Galas [23] show that a significant fraction of the subterranean fauna of Poland can be found in springs. Certainly, the classification of springs from the point of view of the groundwater fauna needs more attention. For example, Dumnicka et al. [17] show that substrate in the spring has a major effect on faunal composition. Slow-flowing waters and lakes; 3.
Small pools on clay or mud; 5.
Dripping or percolation water.
There can be more elaborate subdivisions. For example, Poulson [30] distinguishes several types of cave streams, such as shallow streams and moderately deep master shaft drain streams. While no doubt useful in the context of Mammoth Cave, where he worked, their generality seems very limited. Each of the cave subdivisions has its own strengths and weaknesses, and it seems counterproductive to a priori decide how detailed such a subdivision should be.
Small pore habitats, categorized by Botosaneanu [8] (Table 1), were largely divided on the basis of water flow and salinity. An interesting feature of the classification is the hypotelminorheal-a persistent wet spot, a kind of perched aquifer, fed by subsurface water in a slight depression in an area of low to moderate slope, rich in organic matter, underlain by an impermeable clay layer typically less than 50 cm below the surface-see Culver and Pipan [9] and Meštrov [31]. It appears in the classification both as a porous habitat and as a spring (called seepage springs by Keany et al. [32]). If more were known about the drainage area of all springs, this duality of classification should occur for all spring types. The fauna collected at seepage springs seems to be a mixture of species primarily found in the hypotelminorheal and species primarily found in the seepage spring itself [33].
Just as caves dominate large pore (diameter) habitats, the hyporheal dominates small pore (diameter) habitats, at least in terms of the amount of research done. Orghidan [34], (see Käser [35] for an English translation) coined the term and defined it as the zone of interstitial spaces constituted by the sediments of the stream bed. To our knowledge there have not been formally proposed subdivisions but even at the scale of a few meters there are differences, especially due to upwelling and downwelling [36]. In different streams, the connections between the hyporheic, the groundwater zone, and underlying impermeable substratum vary and have a profound effect on the composition and abundance of the fauna. Malard et al. [36] list five cases but do not formally name them:

2.
Hyporheic zone created only by advected channel water (no ground water).

3.
Hyporheic zone created by advection by both channel water and ground water.

4.
Hyporheic zone created only by infiltration of channel water beneath the stream bed (no parafluvial flow).

5.
A perched hyporheic zone created only by infiltration of channel water beneath the stream bed.
Their fifth category would appear to include seepage springs (outlet of the hypotelminorheic).
As knowledge of subterranean habitats has grown, the dichotomy between large and small pore habitats has grown increasingly problematic. Habitats such as calcrete aquifers and the hypotelminorheic are likely intermediate in pore size, at least based on the size of the organisms found in these habitats. Although epikarst and percolating water are part of karst, they are small pore habitats, based on the size of the inhabitants [10]. There is accumulating evidence [24] that the size of the habitat spaces has an impact on body size. Their model of how this works is shown in Figure 1 and the overall pattern in relationship to habitat categories is shown in Figure 2.
Given these size differences in both habitat and organism, one would expect the community composition in different habitats to be quite different. Dumnicka et al. [37], basing on 280 records for interstitial habitat, 150 for cave waters and 50 for wells, looked at this question for the Polish groundwater annelids fauna. Somewhat surprisingly, relatively little separation of habitats was found, although there was a tendency for the interstitial fauna to differ from the cave fauna along the first axis of their correspondence analysis (  Hypothetical relationship between pore size (habitat diameter) in subterranean habitats and body size, with selective forces indicated by arrows. Below a minimum (mi), there is not sufficient space for animals to occur without burrowing. Above a maximum (mj), body size is likely constrained by other factors, such as phylogenetic and structural constraints. The relationship need not be linear but is presented as such for simplicity. From Pipan and Culver [24]. Used with permission of the National Speleological Society (www.caves.org). Hypothetical relationship between pore size (habitat diameter) in subterranean habitats and body size, with selective forces indicated by arrows. Below a minimum (m i ), there is not sufficient space for animals to occur without burrowing. Above a maximum (m j ), body size is likely constrained by other factors, such as phylogenetic and structural constraints. The relationship need not be linear but is presented as such for simplicity. From Pipan and Culver [24]. Used with permission of the National Speleological Society (www.caves.org). and body size, with selective forces indicated by arrows. Below a minimum (mi), there is not sufficient space for animals to occur without burrowing. Above a maximum (mj), body size is likely constrained by other factors, such as phylogenetic and structural constraints. The relationship need not be linear but is presented as such for simplicity. From Pipan and Culver [24]. Used with permission of the National Speleological Society (www.caves.org).

Figure 2.
Histograms of body lengths of inhabitants (bottom) and log of pore size (top) for different subterranean habitats. From Pipan and Culver [24], modified. Log pore (habitat) size was used because the range of habitat sizes ranged over several magnitudes. The MSS (milieu souterrain superficiel) is the habitat of interconnected cracks and crevices of scree slopes, especially covered ones. Used with permission of the National Speleological Society (www.caves.org). An interesting question is how many species are found in the different subterranean habitats in a region. In Poland, which was largely glaciated, Dumnicka and Galas [23] demonstrate that there are more records per sample in interstitial waters and wells than in caves. In fact, there are more records per sample in springs than in caves ( Figure 4). As far as we can determine, these are the only data of this type available. It would be interesting to compare the Polish data with data from an unglaciated area, areas which generally have a much richer cave fauna.
Culver and Pipan [10,39], Pipan and Culver [18] and Blatnik et al. [40] suggested that there was a third major category of subterranean habitats, close to the surface and with intermediate-sized habitat spaces-shallow subterranean habitats (see Table 1). In their book-length treatment of the topic, Culver and Pipan [10] emphasized the vertical division-the distance from the surface and defined shallow subterranean habitats as occurring less than 10 m from the surface. They suggested several unifying features: Close surface-subsurface connections (except for calcrete aquifers); • Availability of organic carbon and nutrients; • Generally small habitat (pore) size.
Habitats include the hypotelminorheal, epikarst, hyporheal, and calcrete aquifers. Halse [25] calls calcrete aquifers deep subterranean habitats but gives no criteria for this choice. Calcretes in the Pilbara region are deeper than 10 m [13] while those in the Yilgarn are often less than 10 m [41]. Of course, all caves with natural entrances occur in part at depths of less than 10 m. Culver and Pipan [10] included only habitats that were typically less than 10 m in depth. Whether calcretes are included or not, the grouping of habitats into shallow versus deep has proven to be frequently used, although it remains to be seen how useful the distinction is.
composition in different habitats to be quite different. Dumnicka et al. [37], basing on 280 records for interstitial habitat, 150 for cave waters and 50 for wells, looked at this question for the Polish groundwater annelids fauna. Somewhat surprisingly, relatively little separation of habitats was found, although there was a tendency for the interstitial fauna to differ from the cave fauna along the first axis of their correspondence analysis ( Figure 3) (Supplementary Materials). The well fauna, collected in various kinds of wells, overlapped broadly with both. Hahn and Fuchs [38] found a similar pattern for the German fauna. It is curious that there have not been more analyses of this question, and it is worth pursuing.  Table S1. From Dumnicka et al. [37], modified.
An interesting question is how many species are found in the different subterranean habitats in a region. In Poland, which was largely glaciated, Dumnicka and Galas [23] demonstrate that there are more records per sample in interstitial waters and wells than in caves. In fact, there are more records per sample in springs than in caves ( Figure 4). As far as we can determine, these are the only  Table S1. From Dumnicka et al. [37], modified.
Water 2020, 12, x FOR PEER REVIEW 7 of 17 data of this type available. It would be interesting to compare the Polish data with data from an unglaciated area, areas which generally have a much richer cave fauna. Culver and Pipan [10,39], Pipan and Culver [18] and Blatnik et al. [40] suggested that there was a third major category of subterranean habitats, close to the surface and with intermediate-sized habitat spaces-shallow subterranean habitats (see Table 1). In their book-length treatment of the topic, Culver and Pipan [10] emphasized the vertical division-the distance from the surface and defined shallow subterranean habitats as occurring less than 10 m from the surface. They suggested

The Struggle to Measure Aquatic Subterranean Biodiversity
Aquatic subterranean biodiversity surveys share with most other invertebrate surveys both a Linnean shortfall (not enough taxonomists) and a Wallacean shortfall (not all areas sampled). However, there are a number of taxonomists, especially in Europe, actively working on major subterranean groups, including amphipods, copepods, isopods, and oligochaetes. With respect to the Wallacean shortfall, tropical caves, if not other subterranean habitats, are being increasingly well sampled (e.g., [42][43][44]).
However, there are a number of aspects of the subterranean fauna that make unbiased sampling very difficult. Nearly all the analyses and datasets are based on those species limited to subterranean habitats-stygobionts for aquatic habitats. However, in most subterranean habitats there are species that show little morphological modification for subterranean life and/or are also found in surface habitats [18]. Such stygophiles are part of the subterranean community and reproduce in subterranean habitats. Analysis of stygobionts is an analysis of the highly evolved component of subterranean communities specialized and limited to subterranean habitats, not the entire subterranean community.
A second issue is the omnipresence of undescribed species. Some of these species that have been collected are awaiting taxonomic analysis, some have been observed but not collected, and some have not even been observed but are considered likely to be present in a region by expert taxonomists. Obviously the last two are fraught with the likelihood of overestimation or exaggeration of numbers [45]. Bolded letters spell out PASCALIS and it is usually written this way. Collected, undescribed species present a special problem. If they are ignored, then some regions, like Brazil-see Trajano and Bichuette [46]-will appear to be species-poor when in fact they are species rich. One possible solution to the problem is to apply a "discount rate". For example, Culver et al. [45] report that of 19 species listed as undescribed by Holsinger et al. [47] in their enumeration of the West Virginia cave fauna, six turned out to be previously described species. This results in a discount rate of 0.68. The omnipresence of undescribed species can also lead to inappropriate comparisons between regions based on data taken from different time frames-what Culver et al. [45] term the fallacy of provincialism.
A third and related issue is the high frequency of local endemism, often single site endemism. An example of this is from the large-scale European study of subterranean biodiversity, PASCALIS (Protocols for the Assessment and Conservation of Aquatic Life in the Subsurface) [48]. In this study of all known stygobionts from six European countries, there was a strong negative relationship between number of species and number of 0.2 • × 0.2 • grid cells occupied by a species (Figure 5). Of the 930 described stygobionts, 396 were known from a single grid cell.
Water 2020, 12, 2170 8 of 16 Water 2020, 12, x FOR PEER REVIEW 8 of 17 A second issue is the omnipresence of undescribed species. Some of these species that have been collected are awaiting taxonomic analysis, some have been observed but not collected, and some have not even been observed but are considered likely to be present in a region by expert taxonomists. Obviously the last two are fraught with the likelihood of overestimation or exaggeration of numbers [45].Bolded letters spell out PASCALIS and it is usually written this way. Collected, undescribed species present a special problem. If they are ignored, then some regions, like Brazil-see Trajano and Bichuette [46]-will appear to be species-poor when in fact they are species rich. One possible solution to the problem is to apply a "discount rate". For example, Culver et al. [45] report that of 19 species listed as undescribed by Holsinger et al. [47] in their enumeration of the West Virginia cave fauna, six turned out to be previously described species. This results in a discount rate of 0.68. The omnipresence of undescribed species can also lead to inappropriate comparisons between regions based on data taken from different time frames-what Culver et al. [45] term the fallacy of provincialism.
A third and related issue is the high frequency of local endemism, often single site endemism. An example of this is from the large-scale European study of subterranean biodiversity, PASCALIS (Protocols for the Assessment and Conservation of Aquatic Life In the Subsurface) [48]. In this study of all known stygobionts from six European countries, there was a strong negative relationship between number of species and number of 0.2° × 0.2° grid cells occupied by a species (Figure 5). Of the 930 described stygobionts, 396 were known from a single grid cell. This means that unless all cells are sampled, many single cell endemics will be missed. In practice it appears that the pattern of species richness is unaffected if single cell and other narrow endemics are not included [50]. More generally, missing species do not seem to affect geographic This means that unless all cells are sampled, many single cell endemics will be missed. In practice it appears that the pattern of species richness is unaffected if single cell and other narrow endemics are not included [50]. More generally, missing species do not seem to affect geographic pattern. Deharveng et al. [49] measured these missing species using a jackknife procedure that resamples and takes into account the number of narrow endemics. Their results, by country, are shown in Table 2.  This means that unless all cells are sampled, many single cell endemics will be missed. In practice it appears that the pattern of species richness is unaffected if single cell and other narrow endemics are not included [50]. More generally, missing species do not seem to affect geographic pattern. Deharveng et al. [49] measured these missing species using a jackknife procedure that resamples and takes into account the number of narrow endemics. Their results, by country, are shown in Table 2. The final impediment is the Racovitzan impediment, the incompleteness of sampling of subterranean habitats [11]. Caves are incompletely sampled because new caves are constantly being discovered and the number of known caves is very large (e.g., more than 10,000 in Slovenia). While sampling all caves is neither possible nor necessary, the high levels of endemism makes thorough sampling important. Nevertheless, there are some very large datasets for subterranean animals. Probably the largest is that of Zagmajster et al. [51], who assembled data for 21,000 occurrences of 1570 European aquatic subterranean species. The situation for non-cave aquatic subterranean species is particularly difficult because neither epikarst, the hyporheal, or the hypotelminorheal can be sampled easily. For the hyporheal, most samples are taken by pumping water out of the habitat through a fine mesh filter, and then sorting the samples. The Bou-Rouch pump [52,53] in many ways made the sampling of hyporheal possible. Prior to the development of a continuous filtering device [54], epikarst could only be sampled very indirectly by sampling drip pools, themselves biased samples of the organisms in dripping water [55]. The situation is even more dire for the hypotelminorheal. No good sampling device exists, and the habitat must be destructively sampled. However, Niemiller et al. [56] were able to show that eDNA of seep amphipods (genus Stygobromus) could be detected, allowing for the possibility of non-destructive sampling.
There is one way that at least the extent of undersampling can be better understood and that is to report on sites where no specimens were found. There has been a reluctance to do this, perhaps based on the sense that empty samples are failed samples. However, such reporting can be extremely informative. Dumnicka et al. [57] report those quadrats in Poland where they failed to find any water mites specialized for subterranean life, and almost all of these sites were in the glaciated areas of Poland, where the specialized subterranean fauna should be rare ( Figure 6).

Geographic Patterns of Species Richness in Aquatic Subterranean Invertebrates
The easiest, but certainly not the best, measure of species richness is the species richness at a single subterranean site, α-diversity. It is an inadequate measure because α-diversity is a small fraction of the species richness of a region, γ-diversity [58]. However, the data are much easier to accumulate since the regional analysis requires data on many sites [59]. Culver and Sket [60] published the first list, which included 20 caves and wells with 20 or more obligate subterranean species, including both aquatic and terrestrial species. Since then, the number of such sites has at least doubled, and most recently Culver and Pipan [61] published a list of sites with 25 or more stygobiotic species (Table 3).
Several points emerge. First, with the exception of the Robe River in western Australia, none of the sites are in the tropics or sub-tropics. Second, there is a concentration of hotspots in the Dinaric karst, with four of ten sites in this region.
Counts of numbers of species by country (corrected for size) give a similar picture to Table 3. There is a broad band of high species richness in the north temperate zone (China is relatively unknown) as well as in Australia (Figure 7), where aquatic subterranean diversity is largely found in calcrete aquifers, accessible only by wells [25,26].

Geographic Patterns of Species Richness in Aquatic Subterranean Invertebrates
The easiest, but certainly not the best, measure of species richness is the species richness at a single subterranean site, α-diversity. It is an inadequate measure because α-diversity is a small fraction of the species richness of a region, γ-diversity [58]. However, the data are much easier to accumulate since the regional analysis requires data on many sites [59]. Culver and Sket [60] published the first list, which included 20 caves and wells with 20 or more obligate subterranean species, including both aquatic and terrestrial species. Since then, the number of such sites has at least doubled, and most recently Culver and Pipan [61] published a list of sites with 25 or more stygobiotic species (Table 3). Several points emerge. First, with the exception of the Robe River in western Australia, none of the sites are in the tropics or sub-tropics. Second, there is a concentration of hotspots in the Dinaric karst, with four of ten sites in this region.
Counts of numbers of species by country (corrected for size) give a similar picture to Table 3. There is a broad band of high species richness in the north temperate zone (China is relatively unknown) as well as in Australia (Figure 7), where aquatic subterranean diversity is largely found in calcrete aquifers, accessible only by wells [25,26].  Calcrete aquifers are carbonate deposits that form in the vicinity of the water table as a result of evaporation of groundwater-see Culver and Pipan [10]. Aquatic subterranean species richness is very high in both the Yilgarn and Pilbara aquifers of western Australia, but it remains difficult to compare with other areas because of the distinctive and different way researchers have measured species richness in Australia [45,63].
More detailed analyses at the global scale are not available, but there is more detailed information available for Europe [48]. Zagmajster et al. [51] showed that there is a ridge of high species richness in southern Europe (Figure 8), one that corresponds to a similar ridge in terrestrial species richness [59]. The explanation for this ridge may be that it corresponds to a ridge of primary productivity or that long term temperature oscillations are at a minimum. Calcrete aquifers are carbonate deposits that form in the vicinity of the water table as a result of evaporation of groundwater-see Culver and Pipan [10]. Aquatic subterranean species richness is very high in both the Yilgarn and Pilbara aquifers of western Australia, but it remains difficult to compare with other areas because of the distinctive and different way researchers have measured species richness in Australia [45,63].
More detailed analyses at the global scale are not available, but there is more detailed information available for Europe [48]. Zagmajster et al. [51] showed that there is a ridge of high species richness in southern Europe (Figure 8), one that corresponds to a similar ridge in terrestrial species richness [59]. The explanation for this ridge may be that it corresponds to a ridge of primary productivity or that long term temperature oscillations are at a minimum.
The above analyses of the patterns of aquatic subterranean species richness is an analysis of only the specialists for cave life. e.g., stygobionts. Many cave streams, especially in glaciated areas, have functioning, reproducing communities but without any specialized species. A good example of this is the aquatic fauna in Swildon's Hole, a nearly 10 km long cave in the Mendip Hills of the United Kingdom. Knight [64] found 38 taxa in the main stream, mainly Trichoptera and Diptera. There were no stygobionts, yet a functioning community. The situation in the tributaries was different where three of ten species were specialized for subterranean life.
Very little is known about the geographic pattern of species richness in non-cave aquatic subterranean habitats. The PASCALIS project [48] included extensive sampling of the hyporheal, but only eight relatively small regions were included in the study (Figure 8): Wallonia (Belgium), Jura (France), Rousillon (France), Cantabria (Spain), Padano-Alpine region (Italy), Slovenia, Rhône valley (France), and Garonne (France). Malard et al. [58] showed that species richness was highest in Slovenia, followed by the Rhone and the Garonne (Figure 9). Slovenia and the Garonne are on the ridge of aquatic species richness (Figure 8). One of the features of porous aquifers in general and hyporheal in particular is its fine scale heterogeneity, which is evident in Figure 9.
More detailed analyses at the global scale are not available, but there is more detailed information available for Europe [48]. Zagmajster et al. [51] showed that there is a ridge of high species richness in southern Europe (Figure 8), one that corresponds to a similar ridge in terrestrial species richness [59]. The explanation for this ridge may be that it corresponds to a ridge of primary productivity or that long term temperature oscillations are at a minimum. The above analyses of the patterns of aquatic subterranean species richness is an analysis of only the specialists for cave life. e.g., stygobionts. Many cave streams, especially in glaciated areas, have functioning, reproducing communities but without any specialized species. A good example of this is the aquatic fauna in Swildon's Hole, a nearly 10 km long cave in the Mendip Hills of the United Kingdom. Knight [64] found 38 taxa in the main stream, mainly Trichoptera and Diptera. There were no stygobionts, yet a functioning community. The situation in the tributaries was different where three of ten species were specialized for subterranean life.
Very little is known about the geographic pattern of species richness in non-cave aquatic subterranean habitats. The PASCALIS project [48] included extensive sampling of the hyporheal, but only eight relatively small regions were included in the study ( Figure 8): Wallonia (Belgium), Jura (France), Rousillon (France), Cantabria (Spain), Padano-Alpine region (Italy), Slovenia, Rhône valley (France), and Garonne (France). Malard et al. [58] showed that species richness was highest in Slovenia, followed by the Rhone and the Garonne (Figure 9). Slovenia and the Garonne are on the ridge of aquatic species richness (Figure 8). One of the features of porous aquifers in general and hyporheal in particular is its fine scale heterogeneity, which is evident in Figure 9. Pipan et al. [65] did a small regional study of epikarst copepods in Slovenia. They analyzed 81 drips from13 caves in three karst areas of Slovenia. As with Malard et al.'s analysis, they found small Pipan et al. [65] did a small regional study of epikarst copepods in Slovenia. They analyzed 81 drips from13 caves in three karst areas of Slovenia. As with Malard et al.'s analysis, they found small scale differences and that α diversity (within drip species richness) was small compared to β diversity. Of the 18 species, only three were accounted for by α diversity, three by differences within a cave, six by differences between caves in a region, and 18 by differences among regions ( Figure 10). When the data are viewed in another way, one that emphasizes the occurrence of "hotspot" drips, a different pattern emerges.
Water 2020, 12, x FOR PEER REVIEW 13 of 17 scale differences and that α diversity (within drip species richness) was small compared to β diversity. Of the 18 species, only three were accounted for by α diversity, three by differences within a cave, six by differences between caves in a region, and 18 by differences among regions ( Figure 10). When the data are viewed in another way, one that emphasizes the occurrence of "hotspot" drips, a different pattern emerges. Figure 10. Relative contribution of within drip (α-diversity), among drip, among cave, and among region diversity (all β-diversity) to the overall diversity of 30 epikarst copepod species, relative to random expectation. From Pipan et al. [65].
A few drips contribute a disproportionate share of species diversity. The most species-rich drip in the Dinaric karst has 10 species and the entire Županova jama has 13 species, so this drip contributes 40 percent of the species diversity known from the entire Dinaric karst! The task of assessing epikarst species diversity would be considerably simplified if we had a method of determining which drips were hotspots prior to sampling, but we do not.

Conclusions
Obligate aquatic subterranean organisms (stygobionts) occur in a variety of subterranean habitats that vary both in depth and in pore size. This fauna is relatively well known in Europe and North America, and poorly known elsewhere. Sampling this fauna, particularly the non-cave fauna remains difficult, resulting in a Racovitzan shortfall for these habitats. Species richness is highest in mid-temperate latitudes and the Dinaric karst in southern Europe is a hotspot of stygobiotic species richness. Endemism is high and consequently β diversity is much higher than α diversity.
With increased amounts of data and new analytical tools, both continental and global patterns of species richness, and its explanation, should become clearer.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Table S1. Full names of species cited in Figure 3 and the text.  Table S1. Full names of species cited in Figure 3 and in the text. Figure 10. Relative contribution of within drip (α-diversity), among drip, among cave, and among region diversity (all β-diversity) to the overall diversity of 30 epikarst copepod species, relative to random expectation. From Pipan et al. [65].
A few drips contribute a disproportionate share of species diversity. The most species-rich drip in the Dinaric karst has 10 species and the entire Županova jama has 13 species, so this drip contributes 40 percent of the species diversity known from the entire Dinaric karst! The task of assessing epikarst species diversity would be considerably simplified if we had a method of determining which drips were hotspots prior to sampling, but we do not.

Conclusions
Obligate aquatic subterranean organisms (stygobionts) occur in a variety of subterranean habitats that vary both in depth and in pore size. This fauna is relatively well known in Europe and North America, and poorly known elsewhere. Sampling this fauna, particularly the non-cave fauna remains difficult, resulting in a Racovitzan shortfall for these habitats. Species richness is highest in mid-temperate latitudes and the Dinaric karst in southern Europe is a hotspot of stygobiotic species richness. Endemism is high and consequently β diversity is much higher than α diversity.
With increased amounts of data and new analytical tools, both continental and global patterns of species richness, and its explanation, should become clearer.