1✉ “Emil Racoviţă” Institute of Speleology of Romanian Academy, 13 Septembrie Road, No. 13, 050711, Bucharest, Romania. University of Bucharest, Faculty of Biology, Splaiul Independentei 91-95, Bucharest, R-050095, Romania.
2“Emil Racoviţă” Institute of Speleology of Romanian Academy, 13 Septembrie Road, No. 13, 050711, Bucharest, Romania. 3University Ovidius Constanţa, Faculty of Natural Sciences, Al. Universităţii, corp B, Constanţa, Romania.
2017 - Volume: 57 Issue: 4 pages: 855-868https://doi.org/10.24349/acarologia/20174202
Exploring how community composition and diversity change across spatial and temporal scales is important for biodiversity conservation as it explains the sources of diversity and the processes that create or maintain diversity (Veech et al. 2002; Gering et al. 2003; Noda 2004). Work in this field of research using oribatid mite species has mostly focused on species-area relationships or arboreal communities (Maraun and Scheu 2000; Kaneko et al. 2005; Fischer et al. 2010). Subsoil compartments, e.g. Mesovoid Shallow Substratum (MSS) communities are less known. Described in 1980 as ”Millieu Souterrain Superficiel”, MSS is an intermediate habitat between the base of the soil and the bedrock. It is composed by a network of small cracks and voids and is commonly situated in the lower levels of scree slopes (Juberthie et al. 1980; Juberthie 1983).
The MSS’s particular environmental conditions, i.e., absence of light and photoperiod, low temperature fluctuations (Culver and Pipan 2014), very high humidity (Giachino and Vailati 2010) and bi-directional flux of organic material (Culver and Pipan 2014) make it a unique habitat that shelters and serves as a microrefuge for epigean, endogean and hypogean invertebrate species (Nitzu et al. 2014). Therefore, understanding the MSS community structure and diversity as related to spatio-temporal variation of environmental factors is important.
Recent studies revealed that the MSS is inhabited predominantly by oribatid mites (Skubała et al. 2013). Oribatid mites are distributed worldwide, including the alpine regions, and are one of the richest and most abundant of the Acari taxa in soils with high content of decaying organic matter (Krantz and Walter 2009). Oribatid abundance is influenced by environmental variables such as temperature and humidity (Culver and Pipan 2014, Mumladze et al. 2015, Pipan et al. 2011). There is limited information on oribatid mites from MSS, and the patterns and drivers of their diversity and community structure are poorly understood (Jiménez-Valverde et al. 2015). Although the interest in evaluating the importance of the MSS as a hot spot for invertebrate diversity increased in the last years (Nae and Ilie 2004; Nitzu et al. 2006; Nitzu et al. 2010; Nitzu et al. 2014; Pipan et al. 2011), there is little information on oribatid community structure and spatio-temporal dynamics in these environments (Skubała et al. 2013).
In this study we aim to investigate: i) the diversity and community structure of oribatid mites from edaphic environment, superficial MSS (i.e., the upper layer of MSS) and deep MSS (i.e., lower layer of MSS) in a sub-alpine region in the Piatra Craiului National Park, Southern Carpathians, Romania; ii) the relationships between the environmental variables, such as temperature and humidity, and diversity and community structure of oribatid mites. We discuss the results in the context of increased interest in studying MSS habitats from an ecological perspective, and we provide here a first comprehensive study on the oribatid mite communities in MSS.
The study area is located in the Piatra Craiului National Park, Southern Carpathians, one of the most important karst areas in Romania (Nitzu et al. 2014). Piatra Craiului Massif is a 20 km2 limestone ridge with more than 500 caves and diverse types of talus and scree slopes, both covered and open (see Culver and Pipan 2014).
Three sampling sites were selected: (1) Cerdacul Stanciului – a mobile limestone scree situated near Stanciului Cave; (2) Marele Grohotiș – the largest mobile nude limestone scree accumulation from Piatra Craiului Massif; and (3) Valea Seacă – a stabilized type of MSS, covered by forest (spruce and beech).
Cerdacul Stanciului and Marele Grohotiș are sub-alpine habitats, classified as "calcareous and calcashist screes of the montane to alpine levels – Thlaspietea rotundifolii" (Doniță et al. 2005) and are listed in the 8210 habitat types following Natura 2000 habitats classification. Valea Seacă is a R6111 type of habitat – Carpathian South-East communities of fixed screes with Geranium macrrorhizum, Sedum fabaria and Geranium lucidum (Doniță et al. 2005).
All Oribatida material used in this study was collected as a part of a broader study concerning the diversity of epigeal invertebrates, commonly beetles and wandering spiders (Nitzu et al. 2014). Pitfall traps were used to collect invertebrates from edaphic habitat (EDAF), and drillings for MSS (López and Oromi 2010). Pitfall traps might not be the most effective sampling technique for oribatid mites, as the traps sample surface-active invertebrates, estimating the abundance of each species as a function of its activity during the sampling period and population density in the habitat (Brown and Matthews 2016). Thus this study provides data on oribatid species activity-density, i.e., the abundance of a species in pitfall traps is an unknown function of that species’ surface activity and density in the surrounding habitat, but for simplicity we refer to oribatid mite “abundance” throughout the paper.
The MSS was sampled at two depths: 0.5 m – the superficial MSS (SMSS) and 0.75 m – the deep MSS (DMSS). The MSS at each sampling site was sampled at two altitudes: Cerdacul Stanciului at 1637 m and 1672 m, Marele Grohotiș at 1579 m and 1580 m and Valea Seacă at 1087 m and 1200 m, respectively (Table 1). However, the drilling at 1200 m was lost (represented by missing values (-) in Table 1).
To collect mites with drillings, inside of each drilling we placed a trap half filled with 70% ethanol. The upper part of the drilling was covered with a plastic lid to prevent the debris and rocks from falling inside. Each trap was emptied once a month from April to November (8 months), in 2008 and 2009. This made a total of 80 MSS samples: Cerdacul Stanciului 8 months x 2 years x 2 depths; Marele Grohotiş 8 months x 2 years x 2 depths; Valea Seacă 8 months x 1 year x 2 depths.
To collect mites from EDAF, we selected four plots, one in Cerdacul Stanciului, two at Marele Grohotiş and one plot in Valea Seacă. In each plot five pitfall traps, with 70% ethanol were set at a depth of 9 cm, covering a perimeter of 25 m2. Each of the five traps per plot was considered as one independent sample. The traps were emptied at 5-day intervals, and the material collected in each trap pooled for each month. The pitfall traps were installed during May and June in Cerdacul Stanciului and Marele Grohotiş. In Valea Seacă the pitfall traps were installed during March, May and July. This made a total of 35 EDAF samples: Cerdacul Stanciului 5 traps x 1 plot x 2 months; Marele Grohotiş 5 traps x 2 plots x 2 months; Valea Seacă 5 traps x 1 plot x 3 months (Table 1). One pitfall trap did not collect any oribatid mites and 22 pitfall traps collected at list one specimen.
The collected mites were sorted and identified to genus and species level. We used the identification keys published by van der Hammen (1952), Bernini (1978), Pérez-Iṅigo (1993, 1997), Weigmann (2006). The systematic ranking of the species was done after Subías (2004, updated in 2017). After identification the material was preserved in 70% ethanol and stored in the collection of “Emil Racoviţă” Institute of Speleology, Bucharest, Romania. Original taxonomic descriptions are not included in the References.
The temperature (T) (oC) and relative humidity (Rh) (%) were measured using a humidity LogR thermo-hygrometer “Digi-Sense” Cole – Palmer. T and Rh were measured at the ground level (hereafter soil temperature (Ts) and soil relative humidity (Rhs)) and at 0.5 m and 0.75 m depth (hereafter subsoil temperature (Tss) and subsoil relative humidity (Rhss)). Figure 1 summarizes these measurements.
The structure of mite communities in the three habitat types (EDAF, SMSS and DMSS) was accessed based on the dominance and constancy of the species and the number of the unique species. The dominance of an individual species was calculated as DO = Ni / N * 100 (%), where N = the total number of individuals in each habitat, Ni = the total number of individuals of the ith species. The species with DO ≥ 5 were considered dominant and species with DO ≥ 10 – eudominant. Constancy is the percentage of the samples in which the species occurred: C = Ls / L * 100 (%), where L = the total number of samples in each habitat, Ls = the number of samples in which the species was found. Species with C ≥ 50 were considered constant, and species with C > 75 – euconstant. The unique species are species that occur at only one site.
Because the number of samples in each habitat type was not equal and the total abundance (N) per habitat type was different (EDAF: N = 412; SMSS: N = 807; DMSS: N = 349) the data was normalized before analysis. To normalize the data, (i) for abundance analysis, we used average abundances of species, i.e., number of individuals per sample / number of samples for each habitat; and (ii) for diversity analysis, we rarefied the abundance matrix, using 300 individuals per sample. The following diversity indices were chosen to investigate the community composition and how it differs among habitat types: species richness (S), Shannon-Wiener index (H′), Simpson’s index (D) and Pielou’s (J′) index (Maguran 2004).
General Linear Mixed Models (GLMM) were applied to test whether the main community features (abundance and species richness) were related to the habitat type (EDAF, SMSS and DMSS) and to the environmental variables (Ts, Rhs, Tss and Rhss). The habitat type, Ts, Rhs, Tss and Rhss were introduced in the analysis as fixed factors, the sampling site (CS, MG and VS) was considered a random factor and month was included as a repeated measure. We built a set of four candidate models and assess the relative performance of these models using a selection technique based on Akaike’s information criterion corrected for sample size (AICc: Burnham and Anderson 2002; Johnson and Omland 2004). Models were ranked, and the one with the lowest AICc was used as the reference for calculating the AIC difference (∆i) and the likelihood of a model given the data and model weights (wi). Models within two AIC units of the AICmin were considered competitive and more plausible than others (Burnham and Anderson 2002).
The mite species assemblage relationships were modelled by applying a series of Constrained Correspondence Analyses (CCA) and the so-called “partial CCA”. They can be used to model the multivariate response of a species assemblage to a matrix of explanatory variables (Borcard et al. 1992; Legendre and Legendre 1998). The Correspondence Analysis approach is appropriate because it preserves the chi-square distance of the sample and thus correctly handles species frequency.
In particular, the following two models were considered on abundance data pooled over each month:
CCA I: A f [habitat]
CCA II: A f [environmental variables]
Where A is the species abundance matrix, the function f is the linear combination of independent variables and the operator [ ] describes the operation for partitioning out the component of variation described by the linear function within the square brackets. The independent variables included in the CCA analysis were the habitats type (EDAF, SMSS and DMSS) and the environmental variables (Ts, Rhs, Tss and Rhss). The abundance matrix was ln (x + 1) transformed to maintain normal distribution and to avoid the "arch effect" in CCA (Ter Braak 1986). The permutation procedure (based on 9999 cycles) was used to test the significance of explanatory variables in CCA for all eigenvalues (Oksanen et al. 2006).
All analyses were performed using R version 3.2.1 (R Development Core Team 2016). The abundance matrix rarefaction and the CCA were performed using vegan package (Oksanen et al. 2006). The diversity analysis was done using BiodiversityR package (Kindt 2014) and the GLMM using the lme4 package (Bates and Maechler 2010).
A total of 1568 oribatid mites belonging to 94 species, 57 genera and 28 families was collected (Table 2). The number of mites collected per sample ranged from 1 to 100, and the species richness from 1 to 21. Overall, in the three habitats, five species were eudominant, five species were dominant, three species were constant and none was euconstant (Table 2). Three species were eudominant and constant in one of the three habitats: Phthiracarus sp. in EDAF, Oribatella longispina (Berlese, 1914) in SMSS and Ceratoppia bipilis (Hermann, 1804) in DMSS. The latter species was also dominant and constant in EDAF. The most abundant species, Ceratoppia bipilis and Oribatella longispina, were present in all three habitats but had the highest abundance in SMSS. Four species were new records for Romania: Achipteria elegans (Schweizer, 1956), Eupelops plicatus (C.L. Koch, 1835), Ommatocepheus ocellatus (Michael, 1882) and Oribatella longispina Berlese, 1915 (see Nae and Ivan 2015).
Twenty-four species were common to all three habitats. On the other hand, 47 species were unique in one of the three habitats (EDAF: 27 species; SMSS: 17 species; DMSS: 3 species).The species richness was higher in edaphic habitat (EDAF) than in MSS; the diversity indices (Shannon-Wiener, Pielou’s and Simpson’s) were lower in the superficial MSS than in the deep MSS or in EDAF (Table 3).
The model selection using AIC indicated that for both abundance and species richness, only the models including the Tss and Rhss were supported (Table 4).
The permutation tests for all two CCA models applied to oribatid mite abundance matrix indicated that explanatory variables accounted for a significant portion of species distribution variation (999 permutations, P < 0.05). The first two canonical axes clarify the main community patterns (Figs 2–3: species assemblages show a gradient in species composition that is collinear with the spatial variation (Figure 2) and environmental factors (Figure 3).
The CCA of the association between oribatid mite species abundance and the habitat shows that species in the upper right quadrate – Oribatella longispina, Chamobates birulai (Kulczynski, 1902), and Pilogalumna tenuiclava (Berlese, 1908) – were associated with the SMSS (Figure 2). Four species were associated with the EDAF habitat: Achipteria coleoptrata (Linnaeus, 1758), Ceratoppia quadridentata (Haller, 1882), Cepheus dentatus (Michael, 1888) and Phthiracarus sp. (Figure 2). Oribatella sp., Hemileius initialis (Berlese, 1908) and Scheloribates laevigatus (Koch, 1835) were associated with the DMSS (Figure 2).
The temperature and relative humidity in the MSS present decreased variation in amplitude in comparison with the external temperature and relative humidity that show high fluctuations (Figs 1, 3). The CCA results for association of the oribatid mite species abundance and environmental variables showed that Ts and Rhs are the strongest determinants of Oribatida community composition (Figures 1, 3). Cepheus dentatus and Camisia biverrucata (Koch, 1839) were the species correlated with Ts, while the presence of Oribatula tibialis (Nicolet, 1855), Neoribates aurantiacus (Oudemans, 1914) and H. initialis was influenced by Rhs (Figure 3).
This study highlights the ecological importance of MSS habitats for oribatid mite communities. We found rich and diverse oribatid mite communities in MSS habitats and we report rare (A. elegans) and new records (O. ocellatus) in edaphic habitat for the study area. We also report two species new for the Romanian fauna from MSS (E. plicatus and O. longispina). A. elegans was first described by Schweizer (1956) from individuals collected in Switzerland (from spruce and larch forests with high humidity, acidic soils, lichen and moss beds and rocks) and was not reported again with new observations until now. Ommatocepheus ocellatus is frequently found on tree bark and lichens (Weigmann 2006) and is a Palearctic species (Subías 2004). Eupelops plicatus, which has Holarctic distribution (Subías 2004) and prefers forest soils or tree bark (Weigmann 2006), was collected from all three environments. Oribatella longispina was the second eudominant species in our sites with the highest abundance in MSS.
The composition of oribatid mite communities differed among the three habitats. This most likely indicates habitat specialization and explains the observed high levels of species diversity. However, the choice of the pitfall traps as sampling technique of the edaphic habitat may have influenced the type of mites that were sampled (Moreira et al. 2008) and thus the results must be carefully considered when compared to other findings.
The MSS showed the strongest evidence of habitat association (29 species were present only in the MSS), indicating dependence of oribatid mite communities on the particular environmental characteristics of the MSS. A number of authors have brought arguments in favour of the MSS as a habitat for endemic and rare species – species and genera new for science (Honciuc and Stănescu 2003; Arillo et al. 1994), and the importance of MSS as a habitat “in its own right, one with a set of unique species” (Pipan et al. 2011; Culver and Pipan 2014). However, two of the most abundant species in our study, C. bipilis and O. longispina, were found both in EDAF and MSS, although neither of these species was previously reported from the MSS. The first species is widely distributed in Holarctic (Subías 2004) and the second is well represented in East Europe (Subías 2004). Bernini (1978) reported O. longispina as a mountain species.
The observed species richness of the oribatid mites decreased from edaphic habitat through SMSS, to the DMSS (but see the Material section). The diversity indexes, however, are similar for the DMSS and edaphic habitats, all higher than for the MSS. This suggests that the MSS, especially deep MSS, plays an important role as a refuge for edaphic species, which is in line with other studies (Pipan et al. 2011; Nitzu et al. 2014) and supports the hypothesis that the MSS is a gateway in the colonization of the subterranean realm (Pipan et al. 2011).
Furthermore, environmental variables, i.e., subsoil temperature and relative humidity, were important determinants of the species richness of oribatid mites. Temperature and humidity are key factors in determining the presence of invertebrate species in edaphic and MSS habitats (Růžička et al. 1995; Růžička and Zacharda 2010). In summer, when the outside temperature increases and soil relative humidity decreases, the species enter the MSS and screes, using it as a refuge (Nitzu et al. 2014).
Environmental variables were also the major parameter influencing the structure of the oribatid mite community. Soil temperature was the most important factor determining differences in community structure followed by soil relative humidity.
The CCA results also indicated that several oribatid mite species were affected by soil temperature and relative humidity. This is in line with previous studies which indicated that the temperature and humidity shape the species community structure in these scree habitats (Blair et al. 2000, Blesić and Mitrovski 2003, Nitzu et al. 2014).
Temperature can directly and indirectly influence oribatid mite communities through its effect on moisture conditions and subsoil temperature, and it may be a key factor determining the presence of certain species (Jiménez-Valverde et al. 2015). Thus, as pointed out by other studies, the two variables are interconnected and the higher the temperature, the greater the choice of organisms for higher humidity (Madge 1964). The low amplitude of temperature variation in MSS compared with high temperature variation at soil level (Figure 1) makes the MSS a micro–refuge for epigean species when conditions in the surface habitats become too harsh (e.g., too hot) (Nitzu et al. 2014).
The main findings of our study are: (i) the edaphic and MSS habitats have rich oribatid communities that are different among the habitats; (ii) the diversity and community structure of oribatid mites are sensitive to fluctuations of environmental factors, such as temperature and relative humidity. Our results highlight the role of MSS as a reservoir of edaphic oribatid mite species. Although we could not distinguish facultative or obligate subterranean species, the fact that some species were found exclusively in MSS indicates that the MSS may offer habitat to possible trogloxene or troglophilic oribatid species (Mammola et al. 2016). Further studies are needed to confirm the role played by that the MSS in sheltering cave-dwelling mite species. However, we suggest that the MSS should be prioritized for conservation since it may be a key component in maintaining biodiversity.
We thank to Dr. Augustin Nae, from “Emil Racovitza” Institute of Speleology for collecting the material and to Dr. Otilia Ivan from the Institute of Biology, Iaşi for checking and identifying some of the oribatid species. Special thanks to Dr. Ioana Meleg from “Emil Racovitza” Institute of Speleology, to the anonymous reviewers. The study was partially supported by a grant of the Romanian National Authority for Scientific Research, CNCS – UEFISCDI, projects number PN-II-RU-TE-2014-4-1536 to RB.
Arillo A., Gil Martin J., Subías, S.L. 1994 — Oribatidos del "M.S.S" de las Islas Canarias. Poroscheloribatidae subf. n. (Acari, Oribatida) — Mém. Biospéol., 21: 1-6.
Bachmann S., Baker W.J., Brummitt N., Dransfield J., Moat J. 2004 — Elevational gradients, area and tropical island diversity: an example from the palms of New Guinea — Ecography, 27: 299–310.
Bates D., Maechler M. 2010 — Package 'lme4.' hhttp://lme4. r-forge.r-project.org/
Bernini F. 1978 — Notulae Oribatologicae XX. Il genere Oribatella in Italia (Acarida, Oribatida) — Redia, 61: 503-538.
Blair J.M., Todd T.C., Callaham Jr. M.A. 2000 — Responses of grassland soil invertebrates to natural and anthropogenic disturbances — In: Coleman D.C. and Hendrix P.F., Invertebrates as Webmasters in Ecosystems. CABI Publishing, New York, USA, pp. 43-71. doi:10.1079/9780851993942.0043
Blesić B., Mitrovski A. 2003 — Investigation of soil arthropods in the Kragujevac basin, Serbia — Acta entomologica serbica, 8: 57-62.
Borcard D., Legendre P., Drapeau P. 1992 — Partialling out the spatial component of ecological variation — Ecology, 73: 1045-1055.
Brown G.R., Matthews I.M. 2016 — A review of extensive variation in the design of pitfall traps and a proposal for a standard pitfall trap design for monitoring ground-active arthropod biodiversity — Ecology and Evolution, 6(12): 3953-3964. doi:10.1002/ece3.2176
Burnham K.P., Anderson D.R. 2002 — Model selection and multimodel interference. A practical Information – Theoretic approach, second edition — New York. U.S.A.: Springer – Verlag. pp. 515.
Culver D.C., Pipan T. 2014 — Shallow Subterranean Habitats, Ecology, Evolution and Conservation — New York, U.S.A.: Oxford University Press. pp. 258. doi:10.1093/acprof:oso/9780199646173.001.0001
Doniță N., Popescu A., Paucă – Comănescu M., Mihăilescu S., Biriș I. A. 2005 — Habitatele din România — București. Romania: Editura Tehnică Silvică. pp. 496.
Fischer B.M., Schatz H., Maraun M. 2010 — Community structure, trophic position and reproductive mode of soil and bark-living oribatid mites in an alpine grassland ecosystem — Exp. Appl. Acarol. 52(3), 221-237. doi:10.1007/s10493-010-9366-8
Gering J.C., Crist T.O., Veech J.A. 2003 — Additive partitioning of species diversity across multiple spatial scales: implications for regional conservation of biodiversity — Conserv. Biol., 17: 488-499.
Giachino P.M., Vailati D. 2010 — The Subterranean Environment. Hypogean Life, Concepts and Collecting Techniques — Verona, Italy: WBA Handbook. pp. 132.
Hammen L. van der 1952 — The Oribatei (Acari) of the Netherlands — Leiden: Brill. pp. 139.
Honciuc V., Stãnescu M. 2003 — Acarian fauna (Mesostigmata, Prostigmata, Oribatida) of Piatra Craiului National Park — Research in Piatra Craiului National Park, 1: 159-169.
Jiménez-Valverde A., Gilgado J.D., Sendra A., Pérez-Suárez G., Herrero-Borgoñón J.J., Ortuño V.M. 2015 — Exceptional invertebrate diversity in a scree slope in Eastern Spain — J. Insect Conserv., 19: 713-728. doi:10.1007/s10841-015-9794-1
Johnson J.B., Omland K.S. 2004 — Model selection in ecology and evolution — Trends Ecol. Evol., 19: 101-108. doi:10.1016/j.tree.2003.10.013
Juberthie C., Delay D., Bouillon M. 1980 — Extension du milieu souterrain en zone non calcaire: description d'un nouveau milieu et de son peuplement par les Coléoptères troglobies — Mém. Biospéol., 7: 19-52.
Juberthie C. 1983 — Le milieu souterrain: Étude et composition — Mém. Biospéol., 10: 17-65.
Kaneko N., Sugawara Y., Miyamoto T., Hasegawa M., Hiura T. 2005 — Oribatid mite community structure and tree species diversity: a link? — Pedobiologia, 49: 521-528. doi:10.1016/j.pedobi.2005.06.004
Kindt R. 2014 — BiodiversityR: GUI for biodiversity, suitability and community ecology analysis — [Internet]. Available from: ( https://cran.r-project.org/web/packages/BiodiversityR/BiodiversityR.pdf ). pp. 812.
Krantz G. W., Walter D. E. 2009 — A Manual of Acarology. Third Edition — Texas Tech University Press. pp. 807.
Legendre P., Legendre L. 1998 — Numerical Ecology — Elsevier Science. pp. 1006.
López H., Oromi P. 2010 — A pitfall trap for sampling the mesovoid shallow substratum (MSS) fauna — Speleobiology Notes, 2: 7-11.
Madge D.S. 1964 — The humidity reactions of oribatid mites — Acarologia, 6(3): 566-591.
Maguran A. E. 2004 — Measuring biological diversity — Malden: Oxford and Victoria. Blackwell Publishing. pp. 256.
Mammola S., Giachino P. M., Piano E., Jones A., Barberis M., Badino G., Isaia M. 2016 — Ecology and sampling techniques of an understudied subterranean habitat: the Milieu Souterrain Superficiel (MSS) — Sci. Nat., 103:88. doi:10.1007/s00114-016-1413-9
Maraun M., Scheu S. 2000 — The structure of oribatid mite communities (Acari, Oribatida): patterns, mechanisms and implications for future research — Ecography, 23:374-382
Moreira F.M.S, Huising E.J., Bignell D.E. 2008 — Tropical Soil Biology. Sampling and Characterization of Below-ground Biodiversity — London: Earthscan, pp. 218.
Mumladze L., Murvanidze M., Maraun M., Salakaia M. 2015 — Oribatid mite communities along an elevational gradient in Sairme gorge (Caucasus) — Exp. Appl. Acarol., 65(3): 41-51. doi:10.1007/s10493-015-9893-4
Nae A., Ilie V. 2004 — Data concerning the spider diversity (Arachnida: Aranenae) from the Closani karstic area (Oltenia, Romania), with special reference to the superficial subterranean environment — Trav. Inst. Spéol. «Émile Racovitza», 47: 31-41.
Nae I., Ivan O. 2015 — First record of three oribatid mites species (Acari: Oribatida) in the Romanian fauna — Trav. Inst. Spéol. «Émile Racovitza», 54: 33-40.
Nitzu E., Nae A., Popa I. 2006 — Eco-faunistic study on the invertebrate fauna (Araneae, Collembola and Coleoptera) from The Vârghis Gorge Natural Reserve (Eastern Carpathians, Romania), with special note on the micro-refugial role of the subterranean habitats — Trav. Inst. Spéol. «Émile Racovitza», 45-46: 31-50.
Nitzu E., Nae A., Giurginca A., Popa I. 2010 — Invertebrate communities from the mesovoid shallow substratum of the Carpatho-euxinic area: eco-faunistic and zoogeographic analysis — Trav. Inst. Spéol. «Émile Racovitza», 49: 41-79.
Nitzu E., Nae A., Bãncilã R.I., Popa I., Giurginca A., Plãiasu R. 2014 —Scree habitats: ecological function, species conservation and spatial-temporal variation in the arthropod community — Syst. Biodivers., 12: 65-75. doi:10.1080/14772000.2013.878766
Noda T. 2004 — Spatial hierarchical approach in community ecology: a way beyond high context-dependency and low predictability in local phenomena — Popul. Ecol., 46: 105-117.
Oksanen J., Blanchet F.G., Kindt R., Legendre L., Minchin P.R., O'Hara R.B., Simpson G.L., Solymos P., Stevens M.H.H., Wagner H. 2006 — Community Ecology Package. R package version 2.0. — [Internet] Available from: ( https://cran.r-project.org/web/packages/vegan/index.html ).
Pérez-Iṅigo C. 1993 — Acari, Oribatei, Poronota — In: Ramos, M.A. (Ed.), Fauna Iberica, 3, Museo National de Ciencias Naturales. Madrid, pp. 320.
Pérez-Iṅigo C. 1997 — Acari: Oribatei, Gymnonota — In: Ramos, M.A. (Ed.), Fauna Iberica, 9, Museo Nacional de Ciencias Naturales, Madrid, pp. 374.
Pipan T., López H., Oromí P., Polak S., Culver D. C. 2011 — Temperature variation and the presence of troglobionts in terrestrial shallow subterranean habitats — Journal of natural history, 45: 253-273. doi:10.1080/00222933.2010.523797
R Core Team 2014 — R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Internet update 2016, available from ( http://www.R-project.org/ .
Růžička V., Hajer J., Zacharda M. 1995 — Arachnid population patterns in underground cavities of a stony debris field (Araneae, Opiliones, Pseudoscorpionidea, Acari: Prostigmata, Rhagidiidae) — Pedobiologia, 39: 42-51.
Růžička V., Zacharda M. 2010 — Variation and diversity of spiders assemblages along a thermal gradient in scree slopes and adjacent cliffs — Pol. J. Ecol., 58: 361-369.
Schweizer J. 1956 — Die Landmilben des Schweizerischen Nationalparkes 3. Teil: Sarcoptiformes Reuter 1909 — In: Ergebnisse der wissenschaftlichen Untersuchungen des schweizerischen Nationalparks (N.F.). Bd. 5. Druck Lüdin AG. Liestal. pp. 215-377.
Skubała P., Dethier M., Madej G., Solarz K., Makoł J., Kaźmierski A. 2013 — How many mite species dwell in subterranean habitats? A survey of Acari in Belgium — Zool. Anz., 252: 307-318. doi:10.1016/j.jcz.2012.09.001
Subías L.S. 2004 — Listado sistemático, sinonímico y biogeográfico de los ácaros oribátidos (Acariformes: Oribatida) del mundo (excepto fósiles) — Graellsia, 60: 3-305. Internet update 2017, pp. 593, available from: ( http://escalera.bio.ucm.es/usuarios/bba/cont/docs/RO_1.pdf ).
Ter Braak J.F.C. 1986 — Canonical Correspondence Analysis: a new eigenvector technique for Multivariate Direct Gradient Analysis — Ecology, 67: 1167-1179.
Veech J.A., Summerville K.S., Crist T.O., Gering J.C. 2002 —The additive partitioning of species diversity: recent revival of an old idea —Oikos, 99: 3-9.
Weigmann G. 2006 — Hornmilben (Oribatida) — In: Dahl, F. series founder. Die Tierwelt Deutschlands, part 76. Goecke andandand Evers, Keltern, pp. 520.