Next Article in Journal
Plant Diversity Patterns and Conservation Implications under Climate-Change Scenarios in the Mediterranean: The Case of Crete (Aegean, Greece)
Next Article in Special Issue
Are the Neglected Tipuloidea Crane Flies (Diptera) an Important Component for Subterranean Environments?
Previous Article in Journal
Spatial Segregation between Red Foxes (Vulpes vulpes), European Wildcats (Felis silvestris) and Domestic Cats (Felis catus) in Pastures in a Livestock Area of Northern Spain
Previous Article in Special Issue
Cave Communities: From the Surface Border to the Deep Darkness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Abiotic Community Constraints in Extreme Environments: Epikarst Copepods as a Model System

1
ZRC SAZU Karst Research Institute, Novi trg 2, SI-1000 Ljubljana, Slovenia
2
Departments of Biology and Statistics, University of Florida, MCC Statistical Consulting LLC, Gainesville, FL 32611, USA
3
Department of Environmental Science, American University, 4400 Massachusetts Ave. NW, Washington, DC 20016, USA
*
Author to whom correspondence should be addressed.
Diversity 2020, 12(7), 269; https://doi.org/10.3390/d12070269
Submission received: 9 June 2020 / Revised: 30 June 2020 / Accepted: 2 July 2020 / Published: 7 July 2020
(This article belongs to the Special Issue Cave Communities: From the Surface Border to the Deep Darkness)

Abstract

:
The general hypothesis that the overall presence or absence of one or more species in an extreme habitat is determined by physico-chemical factors was investigated using epikarst copepod communities as a model system, an example of an extreme environment with specialized, often rare species. The relationship between the presence or absence of epikarst copepods from drips in six Slovenian caves and 12 physico-chemical factors (temperature, conductivity, pH, Ca2+, Na+, K+, Mg2+, NH4+, and Cl, NO2, NO3, and SO42−) was explored. Statistical analyses included principal components analysis, logistic mixed models, stepwise logistic multivariate regression, classification trees, and random forests. Parametric statistical analyses demonstrated the overall importance of two variables—temperature and conductivity. The more flexible statistical approaches, namely categorical trees and random forests, indicate that temperature and concentrations of Ca2+ and Mg2+ were important. This may be because they are essential nutrients or, at least in the case of Ca2+, its importance in molting. The correlation of Cl and NO3 with copepod abundance may be due to unmeasured variables that vary at the scale of individual cave, but in any case, the values have an anthropogenic component. This contrasts with factors important in individual species’ niche separation, which overlap with the community parameters only for NO3.

Graphical Abstract

1. Introduction

1.1. Background

Extreme environments, such as polar ice caps, hypersaline lakes, acidic peatlands, hyper-arid deserts, groundwater aquifers, and caves, often harbor a specialized fauna that is both numerically rare and patchily distributed. Rarity is typical in extreme environments for which resource scarcity is at least part of what makes them extreme. This resource scarcity is, with the exception of chemoautotrophic sites [1] and sites very near the surface [2], a universal characteristic of aphotic habitats, which lacks primary productivity. In addition, sites with extremes of temperature, such as polar ice caps, or aridity, such as hyper-arid deserts, are also typically very low primary productivity habitats. Sampling rare species and populations of such extreme environments presents unique statistical challenges itself [3], and analysis of ecological relationships among species and between species and their environment is especially difficult given the low numbers of individuals of any one species. Whilst classic ecological techniques like Canonical Correspondence Analysis [4] and Outlying Mean Index [5] address the relationship between species occurrences and abiotic parameters, the spatial distribution of taxa (i.e., occupancy) is still an unexplored question. This question is interesting because it addresses the question of the overall constraints to adaptation to extreme environments.
In extreme environments, there is often a suite of species that show specialized physiological, morphological, and behavioral adaptations to the environment. There are numerous examples in other extreme environments, such as stomatal modification of desert dwelling plants [6], and heat-resistant enzymes in species living in thermal vents on the ocean floor [7]. Convergent adaptation in extreme environments is perhaps best developed, or at least best studied, in the numerous species found in caves and other subterranean habitats. Species specialized for and isolated in subterranean environments typically have elaborated extra-optic sensory structures as well as reduced or absent eyes and pigment [8]. This convergent morphology is an expression of the constraints imposed by aphotic environments. Such communities of specialists, whether in thermal vents, caves, or other extreme environments, are often cited as examples of convergence and parallelism [9,10,11,12,13]. The study of the whole community of morphologically convergent species, rather than individual species, has yielded insights into not only adaptation but also community function (e.g., [14,15]).

1.2. Epikarst as a Model Extreme Environment

Epikarst is the shallow part of karst areas, where stress release, climate, tree roots, and geological processes (especially dissolution of rock) fracture and enlarge rock joints and cracks, creating a more porous zone than the underlying limestone [16]. It is directly beneath the soil, and retains water from infiltration, typically in a 3–10 m thick zone (Figure 1). It is subject to periodic drying and to periodic flushing of water that dislodge the plankton sized organisms that inhabit it [17]. Epikarst is aphotic, low in organic carbon [18] with habitat spaces that limit size and shape of animals occupying the space [19].
Copepods are the predominant group found in epikarst throughout the world, but other crustacean taxa occur regularly, including Amphipoda, Ostracoda, and Syncarida [2]. Sampling of the epikarst fauna is necessarily indirect, and the least biased samples are those that capture the invertebrates found in epikarst water dripping into caves [20,21]. Due to this indirect sampling, most environmental information is about the physical and chemical characteristics of the dripping water. A universal feature of such epikarst samples is that a significant number of samples have no animals (e.g., [20,22,23,24]).

1.3. Aims and Goals

The question we wish to examine is whether these rare, specialized communities are constrained in their spatial distribution by physical and chemical factors. It is common that not all extreme environments of a particular kind, e.g., caves, harbor a specialized fauna. In any region, a number of caves have no specialized fauna, although caves without species are not always reported (but see [25]). The absence of a specialized fauna from a particular site may be because of barriers to dispersal to the site, historical accident, or because the physico-chemical conditions at the site are unsuitable for the fauna in general.
Our hypothesis is that differences in chemical and physical factors among extreme habitats provide an explanation for the presence or absence of a specialized fauna from a particular site. To test this, we analyzed an extensive data set on the occurrence of copepods specialized for life in a shallow subterranean habitat—epikarst, sampled in central Slovenia [20,26]. Pipan [26] took > 250 monthly chemical and biological samples from dripping epikarst water in six caves in central Slovenia. Slightly more than half of these samples lacked copepods. The question then is how the physico-chemical parameters differed between the two kinds of samples (those with and without fauna), and whether there was a biological explanation for the difference. These epikarst copepod communities are especially appropriate for analysis because they are diverse and with variable densities among and within subterranean patches [27]. Thus, we can not only elucidate the important community-wide parameters, but we can also compare these parameters with the individual species parameters that differentiate their niches.

2. Methods and Materials

2.1. Study Area and Sampling Procedures

The study area is in the western Dinaric Mountains of Slovenia. Six caves were studied—Črna jama, Dimnice, Pivka jama, Postojnska jama, Škocjanske jame, and Županova jama (see [20]).
Direct sampling of epikarst habitats is not possible due to the small size of the cavities and their inaccessibility. Thus, the epikarst water and epikarst fauna must be sampled indirectly by taking samples of the percolation water that drips directly from the ceiling. To do this, water from trickles was directed through a funnel into a collecting container. On two sides, the plastic containers have 5 × 5 cm openings covered with a 60 μm diameter mesh to retain animals in the container. At monthly intervals, the water and the organisms in the containers were collected and preserved in 4% formalin. This sample is likely biased in favor or smaller organisms since they are more easily dislodged by water currents [1,2], but it is less biased than pool samples because there is differential survival and reproduction in pools [2]. For the purpose of this study, only the presence or absence of adult copepods was used for analysis. Juveniles were not included because adults should be a better representation of a permanent population [1,2]. In addition, a series of chemical and physical measurements were made at the time of the biological sampling.
Temperature and conductivity were measured in situ by a conductivity meter (WTWTM Model LF 91); pH was measured by a pH meter (WTWTM Model 323). Water samples for ionic analysis were collected and stored in plastic containers and kept at 4 °C prior to measurement. An ion chromatograph (MetrohmTM Model 761 Compact IC) was used to analyze the concentrations of major ions. Chloride, nitrite, nitrate, phosphate, and sulphate were determined using an anion separation column. Sodium, ammonium, potassium, calcium, and magnesium were analyzed using a cation separation column. Concentrations of PO43− were always below the level of detectability of 0.05 mg/L, and therefore not included in the analysis. Detection limits of ions ranged from 50 to 100 μg/L.
Sampling (Table 1) began in Postojnska jama, but when five months of sampling produced relatively few copepods, sampling was switched to other parts of the Postojna-Planina Cave System (Črna jama and Pivka jama). Sampling in three other caves in the region (Dimnice, Škocjanske jame, and Županova jama) was performed for the whole period. Five drips were sampled in each cave, except in Postojnska jama, where 10 drips were sampled. The copepod fauna in the studied drips shows no temporal correlation for each drip [21], and samples within and among drips are treated as independent.

2.2. Statistical Methods

Summary statistics and exploratory analyses.—Stratified means and standard errors were calculated for each physico-chemical variable where strata were identified as each distinct combination of cave and month. In order to explore the relationships among variables, principal components for each cave were obtained for the same set of variables, in this case ignoring the sampling site effect. Thus, we assumed that the relationships among variables are the same regardless of the cave within which the samples were collected. A mixed model with the interaction of month and cave as a fixed effect and sampling site (drip) within cave as a random effect was run for each physico-chemical variable to assess temporal pattern within each cave. The residuals from each model of the physico-chemical variables were checked to verify that the assumptions of the mixed model were met. In some cases, transformations were required to obtain normality and homogeneous variance while in some of the other variables, the model was modified to allow unequal residual variances among months. Where relevant, pairwise differences of means were tested using the Tukey–Kramer method for controlling the experiment-wise error rate.
Modeling the probability of copepod presence in a sample. We ran an all possible subsets logistic regression analysis in which every possible subset of the physico-chemical variables was used as the predictor set. The models were ordered by number of variables and value of the χ2 statistic for the overall test of significance. From these the variables identified in the most parsimonious model with the highest χ2 value, they were chosen for a logistic mixed model for predicting probability of copepod presence with month fixed and sampling site (drip) within cave as a random effect was run.
One concern is that the actual relationship between the explanatory variables and the probability of copepod presence may be different than that assumed for a logistic regression model. In the logistic regression model all variables are assumed to be additive in their effect on the probability of presence, and that the relationship can be described by a logistic curve. In reality, a variable may have a threshold effect such that there is a discontinuous jump in probability at some level of the explanatory variable, or several variables could interact in their effect on the probability. To address these issues, the non-parametric random forest and classification tree approach was used [28,29,30]. The recursive partitioning approach known as random forests was used to provide individual measures of variable importance [29]. This utilizes a bootstrap method that allows for the detection of the relative importance of co-varying physico-chemical measurements. The basic approach is that sets of five randomly selected variables (out of our total of 12) are used to determine a classification tree. This is repeated many times (in our case 2000). For each tree, the important variables are recorded. This allows variables that may be masked by correlation with other predictor variables to be identified as important predictors. The best predictors from random forests were then used to construct classification trees. The predictor space is recursively partitioned into a set of rectangular areas such that samples with (or without) copepods are grouped. A probability of presence of copepods is predicted within each partition. The process is repeated until no improvement in prediction is possible. Confusion matrices showing the accuracy of prediction of the models were calculated. All computations were done in either R (cran.r-project.org), SAS© v.9.2 (SAS Institute, Cary, NC, USA), or JMP® v9.0.0 (SAS Institute, Cary, NC, USA).

3. Results

3.1. What Are the Physical and Chemical Conditions Present in Epikarst Water?

Table 2 summarizes the basic statistical properties of the stratified samples of physico-chemical variables. Overall, the physico-chemical values are typical for carbonate and calcium rich and slightly alkaline waters found in karst regions (Table 2). As is typical of carbonate waters, calcium dominated the cations, and pH was slightly above neutrality. Conductivity was high, a reflection of the relatively high concentrations of Ca2+. The median temperature (8.9 °C) was very similar to the mean annual temperature for the region, approximately 9.0 °C.
Table 3 shows the correlations among the 12 variables. There are four correlations greater than 0.5: K+ and NO3; Na+ and NO3; Mg2+ and SO42−; and Na+ and Cl. The major cations, Ca2+, K+, Mg2+, and Na+, are positively correlated, most of them significantly so (Table 3). Overall, 19 of 66 correlations had absolute values greater than 0.25 and 32 of 66 were statistically significant.
To investigate the differences among caves and dates for environmental variables, PCA analysis, by cave, is shown in Figure 2. Three sites (Pivka jama, Postojnska jama, and Škocjanske jame) have no distinct clusters of points but with outliers (Pivka jama with one and Škocjanske jame with two). The other three caves (Črna jama, Dimnice and Županova jama) have two or three clusters of points, but except for Črna jama, there was no clustering of PCA scores by date. In Črna jama values for October and November formed a distinct cluster, along the first principal component. Loadings on the first component of Črna jama were highest for Mg2+, Ca2+, K+, NH4+, Cl, and SO42−. There was variation in factor loadings among caves. In general, cation and anion loadings were positive for component 1.
Loadings on the second component were more variable. As Postojnska jama was sampled at different months (April through October) than Črna jama and Pivka jama (October through February, Table 1), it is interesting to compare these three PCA plots. In Postojnska jama, temperature has a very different loading than in Črna jama and Pivka jama, but this pattern is not repeated for other variables. In Postojnska jama, the first two PCA axes accounted lowest percent of the variance (43.6%) and in Pivka jama, the first two axes account for the highest percent of the variance (65.2%).
Of the 252 samples taken at various times in the six caves, 136 (54.0%) had adult copepods. For none of the six caves were there significant monthly differences in the presence or absence of copepods (Table 4). Therefore, the connection between physico-chemical factors and the presence or absence of copepods can be considered without the confounding effects of cave and time interactions.
In contrast, for at least one cave, average monthly values for every physico-chemical parameter showed at least one significant monthly difference (Table 4, Figure 3). There were differences among caves, and the most prominent of these anomalous patterns was Županova jama. It was unique in showing no temporal differences in either Mg2+ or Na+, and shared with Postojnska jama the lack of temporal differences in temperature and Cl (Table 4). Postojnska jama was unique in showing temporal differences in NH4+. With the exception of temperature (Figure 3, upper panel), none of the variables showed a yearly cyclical pattern. From April to November, temperature remained near the yearly high for each cave, and then, in all except Županova jama, fell several degrees. In the case of Postojnska jama, there were no data during the winter months. More typical was the pattern shown by conductivity (lower panel of Figure 3), with occasional erratic changes of highs and lows, except for Črna jama which had no significant monthly differences. There was no overall seasonality with respect to the appearance of highs and lows, nor was there any consistent spacing of highs and lows.

3.2. Which Variables Could Be Important in Limiting Copepod Distribution?

The results of the all possible subsets regression analysis are shown in Table 5. No model with more than six variables gave a significant improvement in fit over the six variable model which included temperature, conductivity, K+, NO2, Na+, and Cl (Table 5). When a seventh variable is added (Ca2+), the corresponding increase in χ2 is less than one. When the random effect of sampling site within a cave was included in the model with these explanatory variables, only conductivity and temperature show a significant effect, both being significant at p = 0.022.
A more general and realistic approach, which assumes neither independence or monotonicity of predictor variables is a random forest, showing the factors affecting the presence or absence of epikarst copepods is a random forest (Figure 4). As expected, temperature and conductivity are the most important, as was the case for the logistic regression. However, four other variables seem likely to play a role in determining copepod distribution, listed in decreasing order of importance: Cl, NO3, Ca2+, and Mg2+. How these factors influence presence or absence of epikarst copepods can be determined by a classification tree using these six factors. Three of the six variables are shared with all possible subsets regression (Table 5)—temperature, conductivity, and Cl.
The best model from the classification tree (Table 6) is shown in Figure 5 and summarized in Table 6. Overall, the tree classified 192 of 252 samples (76%) correctly with regard to the presence or absence of copepods, and had a generalized R2 of 0.45. The presence of copepods was correctly predicted 69% (94 of 136) of the time and the absence of copepods was correctly predicted 84% (98 of 116) of the time. The first branch was for temperature and 57 of 136 occurrences of copepods could be predicted solely on the basis of temperature being less than 8.2 °C (Table 6). There was no further branching on the left but for temperatures greater than 8.2 °C, Ca2+ concentration was the next partition with high concentrations also being associated with the presence of copepods. The remainder of the partitions (Figure 5) involved Cl (twice), NO3, Mg2+, and Ca2+ (again). Of the eight terminal leafs, three were associated with absence, and the most important of these (81 of 116 samples without copepods) were samples with high temperature, intermediate Ca2+, low Cl, low Mg2+, and high NO3 (Table 6).

4. Discussion

4.1. Patterns of Variation of Physico-Chemical Factors

The overall physico-chemical characteristics of epikarst drip water are that of water in contact with carbonate rock, with resulting high levels of conductivity, a predominance of Ca2+ cations, slightly basic pH, and temperatures close to the mean annual air temperature (Table 2) [31].
The differences among samples were caused by a combination of location (Figure 2), seasonality (Table 4), and water retention time. We did not measure residence time of water in this study, but previous studies in Postojnska jama [32], found that residence time of water varied from 2.5 months to over a year. Similarly Kluge et al. [33] found retention times of one to three years in three German caves. Three caves—Dimnice, Postojnska jama, and Škocjanske jame—have relatively thick ceilings [20,26], so that, all things being equal, residence time of water may be greater in these caves. If residence time is correlated with presence (and abundance) of copepods, then differences in residence time will result in differences among the caves. The underlying geology of the caves is nearly identical, with all caves formed in Jurassic limestones, although there may be some dolomite in Županova jama [34].
There are factors that vary on the scale of the projected area of a cave onto the surface that influence the physico-chemical variables that we measured. For example, sites in Pivka jama are high in NO3 [20]. We suspect, but cannot demonstrate, that this is the result of the presence of a campground and associated structures [35]. Other potentially anthropogenically causes include high Cl in Dimnice [20], which is also correlated with NO3 (Table 3), and may indicate water from wells with high concentrations of both [36]. There is no commercial activity and minimal agricultural activity in the area, but elevated levels of Cl nonetheless suggest anthropogenic inputs because of the relative rarity of naturally occurring Cl. High Cl concentrations may be from animal waste, human wastewater, or other sources such as salt licks.
Given variation in residence time of water underground, it is not surprising that there is almost no seasonality, but high temporal variability with respect to physico-chemical variables. One manifestation of the variability of residence time is the response to the flow rate of drips to rainfall events, according to Kogovšek [32]. After a dry period or drought, the flow rate of drips does not increase after a rainfall event, but once the epikarst layer is saturated with water, the increase of flow can be in a matter of hours. Of course, the water exiting the drip is not the precipitation water itself, but some of the water previously stored. That is, rainfall has a piston effect on the water in epikarst [37,38].
Temperature showed a rather different pattern than the chemical variables (Figure 3). The differences in temperature among caves are likely the result of details of the differences in sampling times (Table 1), climate differences (especially with respect to Škocjanske jame, which is in region of more Mediterranean climate than the others), and vertical distance from the surface. Except for Postojnska jama, for which there were no winter samples (Table 1), all caves showed a seasonal pattern of reduced temperatures from January to April (Figure 3, upper panel). In some cases, this difference was small. In Županova jama, the monthly least squares means (see Figure 3), varied less than 1 °C, but the coldest months were between January and March. While these differences were not statistically significant, the seasonal pattern in Črna jama, Dimnice, Pivka jama, and Škocjanske jame was statistically significant, and the temperature differences were greater. Kogovšek [32] also found a seasonal pattern in temperature for two drips in Postojnska jama, although the range was less than 1 °C. The relationship between temperature and residence time is a complex one (see [39,40,41]) for a more detailed discussion, but the temperature of resident water eventually approaches the mean annual temperature [39]. Temperature also reflects season, as seen in this study. The lower temperatures in winter are indicative of other events, such as changes in evapo-transpiration and precipitation.

4.2. The Epikarst Copepod Physico-Chemical Niche

A useful beginning point of the analysis is to see which individual physico-chemical factors can account for presence or absence of epikarst copepods in a sample. If the effect of cave is included as a covariate, only temperature and conductivity were significant predictors of the presence or absence of copepods in a sample. However, when a non-parametric multivariate approach, using variables shown to be important in a random forest analysis (Figure 4) in a classification tree (Figure 5, Table 6), was employed, the most important correlates of copepod presence were:
  • Temperature,
  • Ca2+,
  • Mg2+,
  • Cl,
  • NO3.
Conductivity itself was unimportant in the classification tree, presumably because some its major components, especially Ca2+ and Mg2+, were exposed by the random forest analysis, which teased apart correlated variables (see Table 3). Correlation analysis indicated that conductivity was not pairwise correlated with any of the cations (Table 3). The results of the classification tree argue that the relationship of conductivity with the other variables is more complicated than simple pairwise relationships. In fact, a regression tree (not shown) using conductivity as the response variable and the other five variables as predictors indicated that conductivity is explained in order by Ca2+, Cl, NO3, and Mg2+ which is the same order that these variables enter the tree for predicting presence of copepods.
Copepods tend to be in samples with lower temperatures (Table 6), and so the connection between temperature and copepod presence is likely a seasonal one. It is very likely that the relationship between occurrence probability and temperature is driven by copepods getting washed out of epikarst in greater numbers in winter because flow rates are greater. We did not measure flow rates but Kogovšek [32] continuously monitored discharge from two epikarst drips in Postojnska jama for a period of more than two years, and demonstrated that discharge rates were highest in the winter. Rouch [42] observed a similar pattern of copepod drift from a karst spring.
Of the other variables shown to be important in the classification tree, Ca2+ has a strong connection with the biology of copepods. It is critical in the molting process, and some subterranean crustacean species, like the amphipod Gammarus minus, are limited to carbonate springs [43]. Additionally, of interest is that copepods tend not to be found in water with temperature greater than 8.2 °C and Ca2+ concentrations greater than 57.5 mg/L. Water in epikarst can be supersaturated with respect to Ca2+ (part of the mechanism of deposition of CaCO3 in caves (e.g., stalactites)) and this may cause physiological problems for animals in this water. While carbonate geochemists have long focused on the Ca2+ -HCO3 system, we suggest it deserves more attention from biologists working in the same systems. Mg2+ is also a critical nutrient [44], and it is possible that it is limiting in some contexts.
While Mg2+ concentration may or may not be biologically significant, it seems likely that the correlation of epikarst copepod abundance and Cl is due to some other unmeasured variable, one that varies at the scale of cave. Cl concentrations in Dimnice are twice as high (5.54 ± 1.12 mg/L) as in any other cave [26]. We suspect that it is not Cl that is important but some other unmeasured factor.
The correlation with NO3 is perhaps also the result of some other unmeasured variable, but nitrate is also biologically important. It is a frequent contaminant of karst aquifers, resulting from agricultural runoff, septic tanks, and perhaps atmospheric deposition [45]. There are few studies of the nitrogen cycle in caves or epikarst, but available evidence suggests that it is not a limiting nutrient [15,46]. However, there are still a number of puzzling aspects of the nitrogen cycle (see, [45]), such as whether nitrogen fixation occurs in caves. If not, Barton [47] points out it is likely to be limiting in some circumstances.
If we take classification trees as the most general approach to the understanding of the connections of copepod occurrence to physico-chemical parameters, then we have the following factors, listed in order of importance:
  • Temperature, which is likely a reflection of flow velocities rather than community structure.
  • Calcium and perhaps magnesium ions, which are important, both as essential nutrients and in molting.
  • Anthropogenically augmented ions—Cl and perhaps NO3—may indicate contamination from upgradient well water, or they may be surrogates for particular epikarst sites, where some unmeasured variable is important.
The multi-faceted statistical approach, combined with an emphasis on the overall community rather than individual species, has made some sense of the complex patterns of variation of physico-chemical variables.

4.3. The Relationship between Community Physico-Chemical Niche and Individual Physico-Chemical Niches

Pipan [20,26] and Pipan et al. [27] analyzed the same data from a different perspective, one that emphasized niche separation among the epikarst copepod species. Using the same variables with the addition of ceiling thickness, they found that the following parameters were significant factors in distinguishing individual species in a Canonical Correspondence Analysis (CCA): ceiling thickness; NO3; K+; and Na+. Of these, only NO3 was important on a community-wide basis. Thus, the parameters by which the species are separated are, for the most part, distinct from the parameters that predict the presence or absence of one or more species (see Figure 5).
An example of individual niche analysis for NO3 is shown in Figure 6. This is reflected in the presence of some species, especially Brycocamptus dacicus, Bryocamptus n.sp., Moraria varica, and Maraenobiotus brucei, only in high NO3 concentrations and only in Pivka jama (Figure 6). Whatever the source, it points to nitrate as an important factor in organizing communities (see [44]). What is also apparent in Figure 6 is the difficulty in separating or even characterizing the physico-chemical niches of the different species. Of the 27 species found in drips, there were data on more than 100 individuals for only three species, and only an additional four species had more than 10 individuals for which nitrate data were available.
The data we used were not originally collected for the purpose of elucidating physico-chemical niches, and this is the case for many ecological and bioinventory studies. In many of these studies data collected on basic water chemistry remains unconnected and often unanalyzed with respect to the organisms being studied. The results of our study suggest that a careful field study exploring the impact of variation in the physical and chemical characteristics of water on the likelihood of copepods being present may yield additional insights into the forces that control aquatic community structure and dynamics. These impacts could be implicit such as when variables act as surrogates for other factors (such as temperature) or explicit, i.e., with direct effects (such as is likely the case with Ca2+). There are other constraints on the epikarst copepod community, especially the lack of light and low levels of organic carbon, that are a formidable barrier for not just copepods, but any species to survive in epikarst. Thus, the physico-chemical constraints suggested by the classification tree (Figure 5), are not absolute, but constraints in the context of no light and little organic carbon, among the extreme conditions that characterize epikarst in general.

5. Conclusions

A focus on general community occupancy of extreme habitats may ultimately yield insights into environmental constraints on even groups that seem ubiquitous and unconstrained by environmental factors, like copepods. There were physico-chemical differences between samples with fauna and without fauna. The physico-chemical factors that constrained epikarst copepod communities in this study were, with the exception of NO3, different than those that separated the physico-chemical conditions in which different species were found. Some of the community constraints were likely surrogates for age of the epikarst water; and some, like NO3 and Ca2+, may be important biological requirements themselves. In general, the study of such community constraints is made difficult by co-variation of explanatory parameters, and the technique of random forests and classification trees is a useful way out of this dilemma. The study of community constraints in extreme habitats is a fruitful field for further inquiry.

Author Contributions

Conceptualization, T.P. and D.C.C.; methodology, T.P. and D.C.C.; software and formal analysis, M.C.C.; writing—original draft preparation, review and editing, T.P., D.C.C. and M.C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a Mellon Fund grant from the College of Arts and Sciences of American University for D.C.C. and T.P. was funded by EU H2020 project eLTER, RI-SI-LifeWatch and Slovenian Research Agency, P6–0119.

Acknowledgments

Work on this project by DCC was supported by a Mellon Fund grant from the College of Arts and Sciences of American University. Work on this project by TP was supported by the Karst Research Programme P6–0119, RI-SI-LifeWatch, and EU H2020 project eLTER.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Engel, A.S. Chemolithoautotrophy. In Encyclopedia of Caves, 3rd ed.; White, W.B., Culver, D.C., Pipan, T., Eds.; Academic Press: London, UK, 2019; pp. 267–276. [Google Scholar]
  2. Culver, D.C.; Pipan, T. Shallow Subterranean Habitats. Ecology, Evolution, and Conservation; Oxford University Press: Oxford, UK, 2014. [Google Scholar]
  3. Christman, M.C. Sequential sampling for rare or geographically clustered populations. In Sampling Rare or Elusive Species; Thompson, W.L., Ed.; Island Press: Washington, DC, USA, 2004; pp. 134–148. [Google Scholar]
  4. Ter Braak, C.J.F.; Verdonschot, P.J.M. Canonical correspondence analysis and related multivariate methods in aquatic ecology. Aquat. Sci. 1995, 57, 1015–1621. [Google Scholar] [CrossRef]
  5. Doledec, S.; Chessel, D.; Gimaret-Carpentier, C. Niche separation in community analysis: A new method. Ecology 2000, 81, 2914–2927. [Google Scholar] [CrossRef]
  6. Ward, D. The Biology of Deserts; Oxford University Press: Oxford, UK, 2009. [Google Scholar]
  7. Van Dover, C. The Ecology of Deep-Sea Hydrothermal Vents; Princeton University Press: Princeton, NJ, USA, 2000. [Google Scholar]
  8. Culver, D.C.; Pipan, T. The Biology of Caves and Other Subterranean Habitats, 2nd ed.; Oxford University Press: Oxford, UK, 2019. [Google Scholar]
  9. Christiansen, K.A. Convergence and parallelism in Cave Entomobryinae. Evolution 1961, 15, 288–301. [Google Scholar] [CrossRef]
  10. Wiens, J.J.; Chippendale, P.T.; Hillis, D.M. When are phylogenetic analyses misled by convergence? A case study in Texas Cave Salamanders. Syst. Biol. 2003, 52, 501–514. [Google Scholar] [CrossRef]
  11. Pipan, T.; Culver, D.C. Convergence and divergence in the subterranean realm: A reassessment. Biol. J. Linn. Soc. 2012, 107, 1–14. [Google Scholar] [CrossRef]
  12. Protas, M.E.; Jeffery, W.R. Evolution and development in cave animals: From fish to crustaceans. Wiley Interdiscip. Rev. Dev. Biol. 2013, 1, 823–845. [Google Scholar] [CrossRef] [Green Version]
  13. Culver, D.C.; Pipan, T. Shifting paradigms of the evolution of cave life. Acta Carsologica 2015, 44, 415–425. [Google Scholar] [CrossRef]
  14. Venarsky, M.P.; Huntsman, B.M.; Huryn, A.D.; Benstead, J.P.; Kuhajda, B.R. Quantitative food web analysis supports the energy limitation hypothesis in cave stream ecosystems. Oecologia 2014, 176, 859–869. [Google Scholar] [CrossRef]
  15. Venarsky, M.P.; Huntsman, B.M. Food webs in caves. In Cave Ecology; Moldovan, O.T., Kováč, L., Halse, S., Eds.; Springer: Gland, Switzerland, 2018; pp. 309–330. [Google Scholar]
  16. Bakalowicz, M. Epikarst. In Encyclopedia of Caves, 3rd ed.; White, W.B., Culver, D.C., Pipan, T., Eds.; Academic Press: London, UK, 2019; pp. 394–398. [Google Scholar]
  17. Pipan, T.; Culver, D.C. Copepod distribution as an indicator of epikarst system connectivity. Hydrogeol. J. 2007, 15, 817–822. [Google Scholar] [CrossRef]
  18. Simon, K.S.; Pipan, T.; Culver, D.C. A conceptual model of the flow and distribution of organic carbon in caves. J. Caves Karst Stud. 2007, 69, 279–284. [Google Scholar]
  19. Pipan, T.; Culver, D.C. The unity and diversity of the subterranean realm with respect to invertebrate body size. J. Caves Karst Stud. 2017, 79, 1–9. [Google Scholar] [CrossRef]
  20. Pipan, T. Epikarst—A Promising Habitat; Založba ZRC: Ljubljana, Slovenia, 2005. [Google Scholar]
  21. Pipan, T.; Holt, N.; Culver, D.C. How to protect a diverse, poorly known, inaccessible fauna: Identification of source and sink habitats in the epikarst. Aquat. Conserv. Mar. Freshw. Ecosyst. 2010, 20, 748–755. [Google Scholar] [CrossRef]
  22. Meleg, I.N.; Fiers, F.; Moldovan, O.T. Assessing copepod (Crustacea: Copepoda) species richness at different spatial scales in Northwestern Romanian caves. Subterr. Biol. 2011, 9, 103–112. [Google Scholar] [CrossRef]
  23. Meleg, I.N.; Moldovan, O.T.; Iepure, S.; Fiers, F.; Brad, T. Diversity patterns of fauna in dripping water of caves from Transylvania. Ann. Limnol. Int. J. Lim. 2011, 47, 185–197. [Google Scholar] [CrossRef]
  24. Papi, F.; Pipan, T. Ecological studies of an epikarst community in Snežna jama na planini Arto—An ice cave in North Central Slovenia. Acta Carsol. 2011, 40, 505–513. [Google Scholar] [CrossRef]
  25. Fong, D.W.; Culver, D.C.; Hobbs, H.H., III; Pipan, T. The Invertebrate Cave Fauna of West Virginia, 2nd ed.; West Virginia Speleological Survey: Barrackville, WV, USA, 2007. [Google Scholar]
  26. Pipan, T. Ekologija Ceponožnih Rakov (Crustacea: Copepoda) v Prenikajoči Vodi Izbranih Kraških jam. Ph.D. Thesis, University of Ljubljana, Ljubljana, Slovenia, May 2003. [Google Scholar]
  27. Pipan, T.; Blejec, A.; Brancelj, A. Multivariate analysis of copepod assemblages in epikarstic waters of some Slovenian caves. Hydrobiologia 2006, 559, 213–223. [Google Scholar] [CrossRef]
  28. Strobl, C.; Boulestein, A.L.; Kneib, T.; Augustin, T.; Zeileis, A. Conditional variable importance for random forests. BMC Bioinform. 2008, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Strobl, C.; Malley, J.; Tutz, G. An introduction to recursive partitioning: Rationale, application, and characteristics of classification and regression trees, bagging, and random forests. Psychol. Methods 2009, 14, 323–348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Nicodemus, K.K.; Malley, J.D.; Strobl, C.; Ziegler, A. The behavior of random forest permutation-based variable importance under predictor correlation. BMC Bioinfom. 2010, 11. [Google Scholar] [CrossRef]
  31. Herman, J.S. Water chemistry in caves. In Encyclopedia of Caves, 3rd ed.; White, W.B., Culver, D.C., Pipan, T., Eds.; Academic Press: London, UK, 2019; pp. 1136–1143. [Google Scholar]
  32. Kogovšek, J. Characteristics of Percolation Through the Karst Vadose Zone; ZRC Publishing: Ljubljana, Slovenia, 2010. [Google Scholar]
  33. Kluge, T.; Reichelman, D.F.C.; Wieser, M.; Spötl, C.; Sültenfuß, J.; Schröder Ritzrau, A.; Niggemann, S.; Aeschbach-Hertig, W. Dating cave drip water by tritium. J. Hydrol. 2010, 394, 396–406. [Google Scholar] [CrossRef]
  34. Gospodarič, R. Speleogeological data about Taborska jama and its vicinity. Acta Carsol. 1987, 16, 19–34. [Google Scholar]
  35. Kogovšek, J. Natural purification of sanitary sewage during the vertical percolation in Pivka jama. Acta Carsol. 1987, 16, 123–139. [Google Scholar]
  36. Mercado, A. Nitrate and chloride pollution of aquifers; a regional study with the aid of a single cell model. Water Resour. Res. 1976, 12, 731–747. [Google Scholar] [CrossRef]
  37. Williams, P.W. The role of the epikarst in karst and cave hydrogeology: A review. Int. J. Speleol. 2008, 37, 1–10. [Google Scholar] [CrossRef] [Green Version]
  38. Berthelin, R.; Hartmann, A. The shallow subsurface of karst systems: Review and directions. In Eurokarst 2018, Besançon: Advances in the Hydrogeology and Karst and Carbonate Reservoirs; Bertrand, S., Denimal, M., Steinmann, R.P., Eds.; Springer Nature: Gland, Switzerland, 2020; pp. 61–68. [Google Scholar]
  39. Badino, G. Underground meteorology—“What’s the weather underground? ” Acta Carsol. 2010, 39, 427–448. [Google Scholar] [CrossRef]
  40. Covington, M.D.; Luhmann, A.J.; Gabrovšek, F.; Wicks, C.M. Mechanisms of heat exchange between water and rock in karst conduits. Water Resour. Res. 2011, 47. [Google Scholar] [CrossRef] [Green Version]
  41. Covington, M.D.; Perne, M. Consider a cylindrical cave: A physicists’ view of cave and karst science. Acta Carsol. 2015, 44, 363–380. [Google Scholar] [CrossRef] [Green Version]
  42. Rouch, R. Recherches sur lex eaux souterraines—12—Le système karstique du Baget. 1. La phénomène d’hémorragie’ au niveau de l’exutoire principal. Ann. Spéléol. 1970, 25, 665–709. [Google Scholar]
  43. Glazier, D.S.; Horne, M.T.; Lehman, M.E. Abundance, body composition and reproductive output of Gammarus minus (Crustacea: Amphipoda) in ten cold springs differing in pH and ionic content. Freshw. Biol. 1992, 28, 149–163. [Google Scholar] [CrossRef]
  44. Katz, B.G. Nitrate contamination of karst groundwater. In Encyclopedia of Caves, 3rd ed.; White, W.B., Culver, D.C., Pipan, T., Eds.; Academic Press: London, UK, 2019; pp. 756–760. [Google Scholar]
  45. Simon, K.S.; Benfield, E.F. Ammomium retention and whole-stream metabolism in cave streams. Hydrobiologia 2002, 482, 31–39. [Google Scholar] [CrossRef]
  46. Pasternak, K.; Kocot, J.; Horecka, A. Biochemistry of magnesium. J. Elem. 2010, 15, 601–616. [Google Scholar] [CrossRef]
  47. Barton, H.A. Starving artists: Bacterial oligotrophic heterotrophy in caves. In Microbial Life of Cave Systems; Engel, A.S., Ed.; De Gruyter: Berlin, Germany, 2015; pp. 79–104. [Google Scholar]
Figure 1. Conceptual model of epikarst. Gray arrows indicate the direction of slow water flow and black arrows are faster flow paths. From [18]. Used with permission of ZRC SAZU, Založba ZRC.
Figure 1. Conceptual model of epikarst. Gray arrows indicate the direction of slow water flow and black arrows are faster flow paths. From [18]. Used with permission of ZRC SAZU, Založba ZRC.
Diversity 12 00269 g001
Figure 2. PCA analysis by cave. On the left hand panels, individual samples are shown. For Črna jama, the ellipse encloses all samples taken during October and November. It was the only such cluster of sampling dates for the six caves. See text for details. On the right hand panels the eigenvectors for the physical and chemical variables are shown.
Figure 2. PCA analysis by cave. On the left hand panels, individual samples are shown. For Črna jama, the ellipse encloses all samples taken during October and November. It was the only such cluster of sampling dates for the six caves. See text for details. On the right hand panels the eigenvectors for the physical and chemical variables are shown.
Diversity 12 00269 g002
Figure 3. Temporal pattern of temperature (upper panel) and conductivity (lower panel). Vertical bars are 95% confidence intervals. For temperature, Črna jama, Dimnice, Pivka jama, and Škocjanske jame had significant temporal variation. For conductivity, Dimnice, Pivka jama, Postojnska jama, Škocjanske jame, and Županova jama had significant temporal variation.
Figure 3. Temporal pattern of temperature (upper panel) and conductivity (lower panel). Vertical bars are 95% confidence intervals. For temperature, Črna jama, Dimnice, Pivka jama, and Škocjanske jame had significant temporal variation. For conductivity, Dimnice, Pivka jama, Postojnska jama, Škocjanske jame, and Županova jama had significant temporal variation.
Diversity 12 00269 g003
Figure 4. Bar graph of the relative importance of different variables in generating classification trees, when five variables are repeatedly drawn at random.
Figure 4. Bar graph of the relative importance of different variables in generating classification trees, when five variables are repeatedly drawn at random.
Diversity 12 00269 g004
Figure 5. Classification tree for presence or absence of epikarst copepods. The size of the black rectangle is the relative number of empty samples and the size of the gray rectangle is the relative number of samples with copepods. See Table 6 for details.
Figure 5. Classification tree for presence or absence of epikarst copepods. The size of the black rectangle is the relative number of empty samples and the size of the gray rectangle is the relative number of samples with copepods. See Table 6 for details.
Diversity 12 00269 g005
Figure 6. Box and whiskers plot of the log NO3 concentrations (mg/L) for different copepod species, which are arranged by their abundance (indicated by the vertical lines). The rectangles enclose the middle 50% of the data, the line across each rectangle is the group median, the ‘‘whiskers’’ are the minimum and maximum values. Species names are Par-ty2: Parastenocaris sp.2; Spe-inf: Speocyclopes infernus; Par-nol: Parastenocaris nolli alpine; Ela-cve: Elaphoidella cvetkae; Mar-bru: Maraenobiotus brucei; Bry-bal: Brycamptus balcanicus; Bry-dac: Brycamptus dacicus; Sty-nsp: Stygepactophanes n.sp.; Ela-sta: Elaphoidella stammeri; Mor-var: Moraria varica; Bry-pyr: Bryocamptus pyrenaicus; Par-ty3: Parastenocaris sp. 3; Mar-pop: Moraria poppei; Par-ty1: Parastenocaris sp. 1; Nit-sp: Nitocrella n.sp.; Ela-kie: Elaphoidella kieferi; Mor-spA: Moraria sp. B; Ela-nsp: Elaphoidella n.sp.; Mor-dum: Morarioipsis dumonti; Dia-lan: Diacyclops languidoides; Ela-mil: Elaphoidella millennii; Mor-sco: Morariopsis scontenophila; Bry-nsp: Bryocamptus n.sp.; Mor-sta: Moraria stankovitchi; Mor-spB: Moraria sp. B; Par-and: Parastenocaris andreji; Phy-vig: Phyllognathopus viguieri.
Figure 6. Box and whiskers plot of the log NO3 concentrations (mg/L) for different copepod species, which are arranged by their abundance (indicated by the vertical lines). The rectangles enclose the middle 50% of the data, the line across each rectangle is the group median, the ‘‘whiskers’’ are the minimum and maximum values. Species names are Par-ty2: Parastenocaris sp.2; Spe-inf: Speocyclopes infernus; Par-nol: Parastenocaris nolli alpine; Ela-cve: Elaphoidella cvetkae; Mar-bru: Maraenobiotus brucei; Bry-bal: Brycamptus balcanicus; Bry-dac: Brycamptus dacicus; Sty-nsp: Stygepactophanes n.sp.; Ela-sta: Elaphoidella stammeri; Mor-var: Moraria varica; Bry-pyr: Bryocamptus pyrenaicus; Par-ty3: Parastenocaris sp. 3; Mar-pop: Moraria poppei; Par-ty1: Parastenocaris sp. 1; Nit-sp: Nitocrella n.sp.; Ela-kie: Elaphoidella kieferi; Mor-spA: Moraria sp. B; Ela-nsp: Elaphoidella n.sp.; Mor-dum: Morarioipsis dumonti; Dia-lan: Diacyclops languidoides; Ela-mil: Elaphoidella millennii; Mor-sco: Morariopsis scontenophila; Bry-nsp: Bryocamptus n.sp.; Mor-sta: Moraria stankovitchi; Mor-spB: Moraria sp. B; Par-and: Parastenocaris andreji; Phy-vig: Phyllognathopus viguieri.
Diversity 12 00269 g006
Table 1. Summary of sampling coverage of the six study caves, by month. An X indicates that all of the drips in the cave were sampled for physico-chemical parameters, at least once. Data from [26].
Table 1. Summary of sampling coverage of the six study caves, by month. An X indicates that all of the drips in the cave were sampled for physico-chemical parameters, at least once. Data from [26].
MonthPostojnska JamaPivka JamaČrna JamaŠkocjanske JameDimniceŽupanova Jama
Jan XXXXX
Feb XXX
Mar XX
AprX
MayX XXX
JunX XXX
JulX XXX
AugX
SepX XXX
OctXXXXXX
Nov XXXXX
Dec XXXXX
Table 2. Basic statistics for measured variables, using stratified estimates (caves and months define the strata). There were 35 strata and 252 observations. All concentrations are in mg/L. Conductivity is in μS/cm. Temperature is in °C.
Table 2. Basic statistics for measured variables, using stratified estimates (caves and months define the strata). There were 35 strata and 252 observations. All concentrations are in mg/L. Conductivity is in μS/cm. Temperature is in °C.
StatisticTemperatureConductivitypHCa2+K+Mg2+Na+NH4+ClNO2NO3SO42−
Mean8.16356.647.7937.710.470.931.360.0962.340.00302.846.18
S.E. of Mean0.1064.6410.0120.7470.0130.0250.0860.0050.2190.00040.3600.135
Median8.91365.007.7340.600.410.910.950.0841.440.0000.785.95
S.E. of Median0.0434.910.0121.470.0130.0270.0430.0020.050.0010.0830.15
Table 3. Pairwise correlations among the variables, using data from all caves and all dates. Values in bold are greater than 0.5 and values marked with an asterisk are statistically significant (p < 0.05).
Table 3. Pairwise correlations among the variables, using data from all caves and all dates. Values in bold are greater than 0.5 and values marked with an asterisk are statistically significant (p < 0.05).
NO3NO2NH4+SO42−K+Ca2+Na+Mg+Cl°CpH
NO31
NO20.202 *1
NH4+0.140 *0.0651
SO42−0.0750.0250.0611
K+0.500 *0.316 *0.327 *0.261 *1
Ca2+0.194 *−0.0340.0340.492 *0.282 *1
Na+0.622 *0.198 *0.264 *0.183 *0.470 *0.210 *1
Mg2+0.247 *0.0230.261 *0.588 *0.420 *0.461 *0.406 *1
Cl0.480 *0.0300.162 *0.0230.0880.157 *0.687 *0.216 *1
Temperature−0.1130.094−0.0290.072−0.062−0.271 *−0.208 *0.022−0.258 *1
pH0.121−0.009−0.035−0.0740.156 *−0.0230.058−0.1050.021−0.0391
Conductivity0.294 *0.055−0.023−0.228 *−0.0390.1210.005−0.0530.085−0.047−0.123
Table 4. General linear mixed model type III tests of fixed effect of month/cave combination on presence/absence of copepods and on the chemical and physical variables. NO2 was not included because it was not possible to normalize the variable. The list of caves with monthly differences includes all caves with at least one significant difference among months. Some of the results are graphically displayed in Figure 3.
Table 4. General linear mixed model type III tests of fixed effect of month/cave combination on presence/absence of copepods and on the chemical and physical variables. NO2 was not included because it was not possible to normalize the variable. The list of caves with monthly differences includes all caves with at least one significant difference among months. Some of the results are graphically displayed in Figure 3.
Covariance Parameter Estimates (SE) for the Random Effects
VariableTransformationSampling Site within CaveResidualNumerator dfDenominator dfF-ValuepCaves with Monthly Differences
Presence of copepodsBinary2.050 (1.193)na371790.740.86None
TemperatureNone1.603 (0.462)1.070 (0.113)43178.111.64<0.0001All caves except Postojnska jama and Županova jama
ConductivityNone3595.6 (1026.4)1935.9 (204.8)43177.85.22<0.0001All caves except Črna jama
pHNone0.00566 (0.00239)0.0263 (0.0028)43179.43.31<0.0001Pivka jama, Postojnska jama, Županova jama
Ca2+None17.90 (7.37)*4396.516.58<0.0001All caves
K+None0.00923 (0.00408)**4393.53.00<0.0001Pivka jama, Postojnska jama, Županova jama
Mg2+None0.0838 (0.0257)0.0772 (0.0082)43176.25.26<0.0001All caves except Županova jama
Na+Log0.232 (0.075)0.320 (0.035)43176.44.71<0.0001All caves except Županova jama
NH4+None0***43105.82.090.0012Postojnska jama
ClSquare root0.258 (0.071)0.0841 (0.0089)43177.34.66<0.0001All caves except Postojnska jama and Županova jama
NO3Log0.888 (0.242)0.232 (0.025)43177.17.77<0.0001All caves
SO42−None1.449 (0.501)3.126 (0.301)43177.29.25<0.0001All caves
* Unequal variances among months, with residuals ranging from 8.95 to 338.01; ** Unequal variances among months, with residuals ranging from 0.0064 to 0.0776; *** Unequal variances among months, with residuals ranging from 0.000245 to 0.00370.
Table 5. Best subset of physico-chemical data for predicting presence of copepods for each of k = 1,…7 explanatory variables using logistic regression. The six variable model was chosen for further analysis because addition of an additional variable added less than one to the χ2 value.
Table 5. Best subset of physico-chemical data for predicting presence of copepods for each of k = 1,…7 explanatory variables using logistic regression. The six variable model was chosen for further analysis because addition of an additional variable added less than one to the χ2 value.
Number of Variablesχ2Variables
112.978Temperature
218.265Temperature, Conductivity
322.426Temperature, Conductivity, K+
424.350Temperature, Conductivity, K+, NO2
526.117Temperature, Conductivity, K+, NO2, Na+
627.782Temperature, Conductivity, K+, NO2, Na+, Cl
728.441Temperature, Conductivity, K+, NO2, Na+, Cl, Ca2+
Table 6. Terminal “leafs” of classification tree (see Figure 5). Conductivity differences were not in the optimal classification tree. Blanks indicate that no dichotomy for that variable exists for that leaf. Concentrations are in mg/L and temperature is in °C.
Table 6. Terminal “leafs” of classification tree (see Figure 5). Conductivity differences were not in the optimal classification tree. Blanks indicate that no dichotomy for that variable exists for that leaf. Concentrations are in mg/L and temperature is in °C.
TemperatureCa2+ClNO3Mg2+Frequency with CopepodsFrequency without Copepodsn
<8.2 0.6680.33285
≥8.2≥57.5 0.7580.24222
≥8.2<15.3≥2.11 0.8430.1579
≥8.2<15.3<2.11 0.4100.59022
≥8.2<57.5 and ≥15.3 <0.21 0.5870.41315
≥8.2<57.5 and ≥15.3 ≥0.21≥1.250.4470.55311
≥8.2<57.5 and ≥15.3≥2.85≥0.21<1.250.5560.4445
≥8.2<57.5 and ≥15.3<2.85≥0.21<1.250.0520.94883

Share and Cite

MDPI and ACS Style

Pipan, T.; Christman, M.C.; Culver, D.C. Abiotic Community Constraints in Extreme Environments: Epikarst Copepods as a Model System. Diversity 2020, 12, 269. https://doi.org/10.3390/d12070269

AMA Style

Pipan T, Christman MC, Culver DC. Abiotic Community Constraints in Extreme Environments: Epikarst Copepods as a Model System. Diversity. 2020; 12(7):269. https://doi.org/10.3390/d12070269

Chicago/Turabian Style

Pipan, Tanja, Mary C. Christman, and David C. Culver. 2020. "Abiotic Community Constraints in Extreme Environments: Epikarst Copepods as a Model System" Diversity 12, no. 7: 269. https://doi.org/10.3390/d12070269

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop