applied sciences

Article The Chemical Properties Influencing the Oribatid (; ) Abundance and Diversity in Coal Ash Basin Vicinage

Zuzana Feketeová 1,* , Barbara Mangová 2,* and Malvína Ciernikovˇ á 1

1 Department of Soil Science, Faculty of Natural Sciences, Comenius University in Bratislava, Mlynská Dolina, Ilkoviˇcova6, 84215 Bratislava, Slovakia; [email protected] 2 Institute of Zoology, Slovak Academy of Sciences, Dúbravská cesta 9, 84104 Bratislava, Slovakia * Correspondence: [email protected] (Z.F.); [email protected] (B.M.)

Abstract: The samples of two technogenic sediments (MOS, coal ash sediment; MOD, a mixture of ash and natural soil) and two natural (MOM, meadow; MOF—forest; both Cambisol Dystric) in the vicinage of the coal ash basin were studied. We evaluated risk element concentrations and select-ed chemical and microbiological parameters to determine their influence on the structure of Oribatida. High concentrations of various toxic elements, alkaline pH, and low hu-midity negatively affected the abundance of oribatid . The microbial indicator values showed that the soil microbial community formed in technogenic sediments could effectively use organic carbon. However, considering the wide C/N ratio of the substrates, the process of soil organic matter (SOM) decomposition was slowed down, and thus nutrients were less available for the mites, which could be among the reasons why we did not find any individual of oribatid mite in MOS. In MOD, however, we found representatives of

 three species, but only Tectocepheus velatus sarekensis established an abundant community at the highly  contaminated site. Anthropogenic pressure resulted in the selection of r-strategists, which became

Citation: Feketeová, Z.; Mangová, B.; dominant and reduced the whole community’s species’ diversity. Therefore, we consider it a suitable Cierniková,ˇ M. The Soil Chemical indicator of improper human intervention in the ecosystem. Properties Influencing the Oribatid Mite (Acari; Oribatida) Abundance Keywords: oribatid mites; metal pollution; Tectocepheus velatus sarekensis; bioindicators; diversity and Diversity in Coal Ash Basin Vicinage. Appl. Sci. 2021, 11, 3537. https://doi.org/10.3390/ app11083537 1. Introduction Global changes in climate and land use are affecting the of soil Academic Editor: Sung-Deuk Choi in many ecosystems [1]. Oribatid mites represent species-rich and highly abundant soil mites inhabiting the soil’s upper layers [2]. There is increasing interest in their reaction to Received: 28 March 2021 environmental conditions such as metal pollution because of their essential role in detrital Accepted: 13 April 2021 Published: 15 April 2021 food webs. They are actively involved in the soil decomposition process, increasing soil fertility by recycling nutrients out of dead organic matter. The low reproduction and long

Publisher’s Note: MDPI stays neutral life-cycle rates suggest that oribatid mites may be sensitive to environmental change [3]. with regard to jurisdictional claims in They are not passive inhabitants of ecosystems; instead, they are strong interactors. Despite published maps and institutional affil- being generally sensitive to disturbance, oribatid mites show a considerable variation in iations. sensitivity among taxonomic groups. The various scientific papers declare that the number of species is declining toward the more significant disturbance, and, at the same time, the composition of their entire community is changing [4,5]. The large number of factors determine the distribution and abundance of soil mites. Among the abiotic factors, soil pH and soil organic matter (SOM) represent the most significant drivers concerning soil Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. pH variation on soil arthropods’ presence [6,7] and arthropods’ role in SOM degradation. This article is an open access article Highly increased levels of metals, wide C/N ratio, and pH values were recognized as the distributed under the terms and structuring forces that influence the distribution, diversity, and richness of oribatid species, conditions of the Creative Commons whereas the soil reaction (pH) influences these differences fundamentally [8]. Changes Attribution (CC BY) license (https:// in the value of the soil reaction mostly come from pollution from industrial production creativecommons.org/licenses/by/ and mining. Its value affects the soil’s chemical, physical, and biological properties and 4.0/). the solubility of organic carbon and increases the bioavailability of various toxic elements,

Appl. Sci. 2021, 11, 3537. https://doi.org/10.3390/app11083537 https://www.mdpi.com/journal/applsci Appl. Sci. 2021, 11, 3537 2 of 14

mainly metals. Soils around heating power are often enriched in metals due to ash blowout of refuse heaps, and toxic metals tend to accumulate in high concentrations in the topsoil near point sources [9]. Oribatid mites proved to have a high capacity to colonize postindustrial dumps [10]. With an increasing metal load of the soil, sensitive species may be replaced by resistant ones without changing the total number of mites. One of the critical reasons for reducing the abundance and species’ richness of Oribatida at highly contaminated sites may be explained by toxic effects through the uptake of toxic metal-loaded microorganisms, espe- cially fungi [11]. In general, oribatids and their communities are relatively tolerant of toxic metal exposure [12,13]. Many studies have measured their responses to environmental pollution at the community level [14,15]. The studies proved that heavy pollution exposure decreases both abundance and species’ diversity, but species vary in their sensitivity to pollution, and some species may even benefit from moderate pollution levels [14,16]. Orib- atida includes some species with resistance to metals to maintain their populations in even highly contaminated soils [15–17]. The findings mentioned above suggest that the oribatid mites have a great potential to be used as bioindicators of changing soil environment. Therefore, it is crucial to understand how the mite community responds to ongoing soil ecosystems’ changes. Soil microarthropods play a vital role in the nutrient cycle due to the litter’s fragmentation, increasing the microbial action area and stimulating microorganisms’ growth [18]. They accelerate the circulation of elements within ecosystems, either directly or indirectly, by stimulating soil microbes and contributing to soil fertility maintenance [19]. Their essential role in litter decay as detritivores suggests that any change in their community structure in response to global environmental change is likely to have further consequences for ecosystem functioning. The hypotheses were as follows: (1) Anthropogenic pressure results in the selection of r-strategists, which will become dominant and reduce the whole community’s species’ diversity. (2) Oribatids can be considered appropriate indicators of improper human intervention in the soil ecosystem, reflecting extremely low abundance compared to their abundance in undisturbed . This study’s main goal was to characterize the oribatid mite community, which inhabits the coal ash basin’s extreme environment and its close surrounding. Knowing the biological parameters and the properties of soil and contaminants acting in situ is crucial for proper ecological risk assessment. Concerning coal ash disposal sites, representative data about toxicity on oribatid mites are lacking. Therefore, the aim was to detect the influence of environmental load caused by coal processing to the Oribatida community as a change in species’ occurrence. Based on the hypotheses we made, we targeted identifying species with negative responses under stress and those prospering in deteriorating living conditions.

2. Materials and Methods 2.1. Study Area This investigation aimed to study the species’ composition and diversity of Oribatid mites and microbiological activity in a power ash pond, the dam, and surrounding meadow and forest soil located in Zvolen—Môt’ová, Slovakia (Central Europe), and the relationship between biological parameters of such soils and their chemical properties (Figure1). The ash pond’s dam is 10 m high, and the whole area reaches approximately 7 hectares. The sedimentation of ash sediment with permanent storage in the basin was carried out from 1991 to 2005. This storage represents more than 320,000 m3 of stored ash. At present, it is operated in a cyclical manner—the period of flotation and ash export alternates [20]. Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 14

(Figure 1). The ash pond’s dam is 10 m high, and the whole area reaches approximately 7 hectares. The sedimentation of ash sediment with permanent storage in the basin was car- Appl. Sci. 2021, 11, 3537 ried out from 1991 to 2005. This storage represents more than 320,000 m³ of stored ash. At 3 of 14 present, it is operated in a cyclical manner—the period of flotation and ash export alter- nates [20].

FigureFigure 1. 1.Study Study area area inin the map of of Slovakia Slovakia..

IndividualIndividual sites sites at at which which samples were were taken taken represent represent four four biotopes biotopes with withgradually gradually decreasingdecreasing exposure exposure to to the the contaminant.contaminant. They They include include (1) (1)ash ash pond pond sediment—ash sediment—ash (MOS), (MOS), (2) forming Technosol (MOD)—a heterogeneous mixture of ash and natural soil, (3) (2) forming Technosol (MOD)—a heterogeneous mixture of ash and natural soil, (3) nearby nearby meadow (MOM) and surrounding forest soil (MOF), both representing soil type meadow (MOM) and surrounding forest soil (MOF), both representing soil type Dystric Cam- Dystric Cambisol according to the World Reference Base for Soil Resources WRB Classifi- bisolcation according System to [21]. the The World sites Reference were not Baseseparated for Soil from Resources each other WRB by any Classification physical barrier System [21]. The(Figure sites were 2). The not ash separated sediment fromwas sampled each other in two by lines any physicalof a dry surface barrier of (Figure the impound-2). The ash sedimentment: the was ash sampled deposit’s in central two lines part and of a the dry dam. surface Sites of MOM the impoundment: and MOF represent the the ash clos- deposit’s Appl. Sci. 2021, 11, x FOR PEER REVIEWcentral est natural part and biotopes the dam. where Sites the MOM concentrations and MOF of representany risk elements the closest do not natural exceed biotopes the stat- where4 of 14

theutory concentrations limit. of any risk elements do not exceed the statutory limit.

FigureFigure 2. 2. TheThe sampling sampling area area around around the the ash ash pond pond ( (leftleft above above):): MOS MOS (ash (ash sediment) sediment) and and MOD MOD (forming (forming Technosol), Technosol), ( (rightright above), the eggs of Charadrius dubius in the technogenic sediment of the pond’s dam, (left below) MOM (meadow), (right above), the eggs of Charadrius dubius in the technogenic sediment of the pond’s dam, (left below) MOM (meadow), (right below) MOF (forest). below) MOF (forest). 2.2. Soil Sampling and Analyses We took two sets of samples (one set for microbiological analyses, one set for soil fauna extraction purpose) from each study site from the 0–10-cm depth interval. Soil sam- ples were obtained using soil cores: five soil cores per field, from which eight samples (with a total volume of 0.20 dm³ each) were combined to give one composite soil sample. Afterwards, the samples were placed in sterile plastic bags and put into the cooling box. The fresh soil samples were taken to the laboratory immediately. Once in the laboratory, one part of the samples was ground and passed through a 2-mm sieve; the other part was put in the Tullgren apparatus. Before the necessary chemical and physical analysis proce- dures, the samples were air-dried. The soil samples intended for microbiological analyses were kept field moist and stored until the analyses were performed. We measured the soil moisture content by the gravimetric method. A subsample of a fresh, sieved composite sample was weighed, oven-dried (at 105 °C) until there was no further mass loss, then reweighed. The moisture content in this paper was expressed as a mass of water per mass of dry soil. Soil pH was determined with a glass electrode (Sentix 980) in suspension with distilled water and 1 M KCl, respectively, using 1:2.5 weight ratio (20 g of sample, 50 mL of H₂O or KCl). Electrical conductivity of the ash sediment was measured in the water- saturated paste with TetraCon 925 conductivity probe. The total amount of organic carbon and nitrogen was determined by a method based on the Pregl–Dumas principle using the NCHSO “FLASH 2000” Elemental Element Analyzer from Thermo Scientific. Samples were analyzed for the content of potentially toxic elements by a combination of induc- tively coupled plasma (ICP) emission spectrometry and ICP mass spectrometry. Samples were pretreated by hot acid digestion. Content of 40 elements was determined (Mo, Cu, Pb, Zn, Ag, Ni, Co, Mn, Fe, As, U, Th, Sr, Cd, Sb, Bi, V, Ca, P, La, Cr, Mg, Ba, Ti, Al, Na, K, W, Zr, Ce, Sn, Y, Nb, Ta, Be, Sc, Li, S, Rb, Hf). The analysis of potentially toxic elements in studied samples was carried out in ACME Analytical Laboratories Ltd. (Vancouver, BC, Canada). Appl. Sci. 2021, 11, 3537 4 of 14

2.2. Soil Sampling and Analyses We took two sets of samples (one set for microbiological analyses, one set for soil fauna extraction purpose) from each study site from the 0–10-cm depth interval. Soil samples were obtained using soil cores: five soil cores per field, from which eight samples (with a total volume of 0.20 dm3 each) were combined to give one composite soil sample. Afterwards, the samples were placed in sterile plastic bags and put into the cooling box. The fresh soil samples were taken to the laboratory immediately. Once in the laboratory, one part of the samples was ground and passed through a 2-mm sieve; the other part was put in the Tullgren apparatus. Before the necessary chemical and physical analysis procedures, the samples were air-dried. The soil samples intended for microbiological analyses were kept field moist and stored until the analyses were performed. We measured the soil moisture content by the gravimetric method. A subsample of a fresh, sieved composite sample was weighed, oven-dried (at 105 ◦C) until there was no further mass loss, then reweighed. The moisture content in this paper was expressed as a mass of water per mass of dry soil. Soil pH was determined with a glass electrode (Sentix 980) in suspension with distilled water and 1 M KCl, respectively, using 1:2.5 weight ratio (20 g of sample, 50 mL of H2O or KCl). Electrical conductivity of the ash sediment was measured in the water-saturated paste with TetraCon 925 conductivity probe. The total amount of organic carbon and nitrogen was determined by a method based on the Pregl–Dumas principle using the NCHSO “FLASH 2000” Elemental Element Analyzer from Thermo Scientific. Samples were analyzed for the content of potentially toxic elements by a combination of inductively coupled plasma (ICP) emission spectrometry and ICP mass spectrometry. Samples were pretreated by hot acid digestion. Content of 40 elements was determined (Mo, Cu, Pb, Zn, Ag, Ni, Co, Mn, Fe, As, U, Th, Sr, Cd, Sb, Bi, V, Ca, P, La, Cr, Mg, Ba, Ti, Al, Na, K, W, Zr, Ce, Sn, Y, Nb, Ta, Be, Sc, Li, S, Rb, Hf). The analysis of potentially toxic elements in studied samples was carried out in ACME Analytical Laboratories Ltd. (Vancouver, BC, Canada).

2.3. Mite Extraction and Species’ Identification The Berlese–Tullgren method [22] was used as a method of mite extraction. In the laboratory, ethanol-preserved specimens were cleared in lactic acid (60%) and placed on glass slides with Swann’s solution. Individuals were determined using light microscopy (300×–675×) to the species’ level utilizing keys and original species’ descriptions [23–40].

2.4. Microbiological Analyses Microbial biomass carbon was measured using the fumigation extraction procedure described by Vance et al. [41]. Soil respiration rate was measured as carbon dioxide evolution over 24 h, according to Schinner [42]. We calculated the microbial indicators, microbial metabolic quotient (qCO2), as the amount of basal carbon dioxide produced per unit microbial biomass carbon [43] and substrate availability index (SAI) as the potential for basal production of carbon dioxide ratio according to Cheng et al. [44]. Microbial quotient (qMic) was determined from the microbial biomass carbon to total organic carbon ratio MBC/Corg [45]. This paper’s microbial measurements are the mean of triplicate determinations and expressed on a soil-oven, dry-weight basis.

2.5. Statistical Analyses The statistical estimation was carried out with PAST version 4.03 [46]. Abundance (A) is expressed by the number of individuals in dm2 of the soil surface (ind./dm2). Species’ number (SN) in a sample was defined as the number of species to 1.6 dm2 of soil surface. The species’ dominance (SD) was expressed as the percentage of a particular species’ individuals to the total number of oribatid mites. According to Tischler [47], a modified categorization into five dominance classes was used to identify dominance. Diversity Appl. Sci. 2021, 11, 3537 5 of 14

(Shannon–Weaver’s (H’)) (S—total number of species; pi—the proportion of species i relative to the total number of species) [48],

S H = − ∑ pilnpi (1) j=1

Margalef Richness Index in Biodiversity [49](S—total number of species; n—total number of individuals in the sample),

M = (S − 1)/Log(n) (2)

and equitability (J) (H—Shannon diversity; N—species number) [50] indices were used to evaluate the species’ richness and evenness of oribatocoenoses.

eH S = (3) N All relationships between soil chemical and biological variables were evaluated using Pearson’s correlation coefficients. Differences at p < 0.05 level were deemed significant. Hierarchical clustering (UPGMA: unweighted pair group method with arithmetic mean) was used to visualize the relationship between individual species and coenoses. A matrix of Bray–Curtis distances between coenoses was calculated from the original data matrix based on individual species’ dominance values.

3. Results The main polluting factor on the site was the exceeding concentration of various risk elements, among which eight elements (As, Ba, Cd, Co, Cr, Cu, Ni, Zn) exceeded the permitted limits in the ash sediment MOS and its dam MOD, while arsenic and copper concentrations reached the indication value for remediation (Table1). The MOM and MOF areas represented natural soils (both Cambisol Dystric), in which the concentrations of risk elements did not exceed the reference value.

Table 1. Concentrations of risk elements’ metals in technogenic substrates MOS and MOD; rv, reference value (µg.g−1) represents the maximum permissible concentration of a hazardous substance in the soil set by Act No. 220/2004 Protection and Utilization of Agricultural Soil by Ministry of , Environment and Regional Development of the Slovak Republic.

As Cu Ni Co Ba Cr Zn Cd MOS 185 99.8 105.9 39 763 176 113 <0.5 MOD 79 529.5 120 50 1897 219 217 1.9 rv 29 36 35 20 500 130 140 0.8

The technogenic sediment of the basin (MOS) and the dam (MOD) showed an alkaline soil reaction. In contrast, the meadow (MOM) and the nearby forest’s pH (MOF) showed weakly acidic values. The results of electrical conductivity (1008 µS.cm−1) measurement in the ash sediment suggested that the risk of mobilization of detected toxic elements from the tailings and the adverse effects on biota was not significant. The percentage of organic carbon was high in all monitored samples; the lowest value was found in the MOD sample. The carbon-to-nitrogen ratio, an essential indicator of SOM decomposition rate, was extensive (43–45) in the MOS and MOD samples; in the other samples, this ratio was meaningfully lower (Table2). We also measured significantly lower humidity in the MOS (11.76 %) and MOD (13.64 %) samples compared to all monitored samples. Microbial biomass carbon (MBC) values were very high in all examined samples and ranged from 2.37 to 2.68 mg C/g dry matter. Both basal and potential respiration values were low in MOS and MOD samples. We recorded slightly higher basal respiration in the MOM and MOF, but the potential respiration values were many times higher than in the MOS and Appl. Sci. 2021, 11, 3537 6 of 14

MOD samples (Table2). The microbial quotient (qMic) refers to the microbial activity potential to mineralize organic matter in the soil. This quotient’s lower values are usually associated with the reduced availability of organic material for the microorganisms. The microbial quotient was low in the samples (Table2). On the contrary, the metabolic quotient (qCO2) values were surprisingly lower in MOS and MOD than in the meadow and forest samples (Table2). The metabolic quotient expresses the efficiency of carbon utilization by the microbial community. Its values increase when the microbial community in the soil is exposed to stress. On the other hand, low qCO2 values signal optimal carbon consumption by the soil microorganisms. We found a bit higher values of this parameter in the MOM and MOF. The substrate availability index (SAI), indicating better microorganisms’ ability to utilize the readily available organic substances provided by the technogenic substrate, was lower in MOS and MOD compared to natural soils MOM and MOF, which was entirely unexpected (Table2).

Table 2. Chemical and microbiological properties of soils (MOS, ash sediment; MOD, heterogeneous mixture of ash and

natural soil; MOM, meadow; MOF, forest); pH (H2O/KCL), negative log of the hydrogen ion concentration; w, water content in the sample; Corg, total organic carbon; Ntot, total organic nitrogen; MBC, microbial biomass carbon (dry weight); CO2B, basal respiration (per 24 h); CO2P, potential respiration (per 24 h); SAI, substrate availability index (CO2P/CO2B); qCO2, metabolic quotient 100× (CO2B/MBC); qMic, microbial quotient (MBC/Corg); C/N, total organic carbon-and-nitrogen ratio.

pH pH w Corg Ntot MBC CO2B CO2P SAI qCO2 qMIC C/N (H2O) (KCl) mg mg mg mg %%% −1 −1 −1 −1 C.g CO2.g CO2.g CO2.g MOS 9.17 9.04 11.76 2.63 0.06 2.68 0.05 0.06 1.29 1.74 1.02 43.83 0.32 0.44 1.41 0.17 <0.01 0.06 0.02 0.03 0.66 0.81 3.59 7.30 MOD 8.92 8.71 13.64 1.37 0.03 2.58 0.04 0.12 2.84 1.67 1.88 45.67 0.41 0.25 1.48 0.61 <0.01 0.15 0.03 0.06 0.47 0.46 0.82 8.40 MOM 5.90 5.29 26.68 2.48 0.14 2.37 0.26 1.57 6.03 11.00 0.96 17.71 0.28 0.30 3.22 1.08 <0.01 0.14 0.04 0.21 1.04 0.85 0.13 5.52 MOF 5.56 5.04 28.72 3.82 0.28 2.66 0.18 1.52 8.32 6.85 0.70 13.64 0.32 0.29 4.66 1.11 0.09 0.15 0.04 0.22 0.14 2.53 0.08 3.52

A total of 32 species of oribatid mites were recorded (MOF 62.5%, MOM 74.9%, MOD 9.4%, MOS 0.0%). In the sample MOS, we did not find any representative of the oribatid mites. In the sample MOD, we recorded only three species. Opiella (Moritzoppia) translamellata (Willmann, 1923) and Poecilochthonius italicus (Berlese, 1910) were found in the sample MOD in very low abundance, while the eurytopic species Tectocepheus velatus sarekensis (Trägårdh, 1910) was numerous. Oribatocoenoses MOF and MOM were formed mainly by species of moist meadow and forest habitats (Table3).

Table 3. The list of Oribatida species found at each of three sampling sites near the coal ash basin (listed in alphabetical ). No individual (*) was found in the sample MOS. SN, species number; A, abundance; D, diversity.

MOD MOM MOF SPCIES Abr. SN A D SN A D SN A D Ceratoppia quadridentata (Haller, CerQua C * * * 5 3.13 1.27 * * * 1882) Ceratozetes mediocris Berlese, CerMed A * * * 56 35.00 14.21 51 31.86 15.45 1908 Chamobates cuspidatus (Michael, ChaCus B * * * * * * 1 0.63 0.30 1884) Damaeus sp. DamSp B * * * * * * 1 0.63 0.30 obvia (Berlese, 1913) GalObv B * * * 4 2.50 1.02 7 4.38 2.12 Gustavia microcephala (Nicolet, GusMic A * * * 16 10.00 4.06 17 10.63 5.15 1855) Appl. Sci. 2021, 11, 3537 7 of 14

Table 3. Cont.

MOD MOM MOF SPCIES Abr. SN A D SN A D SN A D Hermaniella dolosa Grandjean, HerDol B * * * * * * 4 2.50 1.21 1931 rufulus C.L. HypRuf A * * * 58 36.25 14.73 15 9.38 4.55 Koch, 1840 Liebstadia pannonica (Willmann, LiePan A * * * * * * 82 51.25 24.85 1951) Metabelba papillipes (Nicolet, MetPap C * * * 6 3.75 1.52 1 0.63 0.30 1855) Nothrus borrusicus Sellnick, NotBor C * * * 1 0.63 0.25 * * * 1928 Nothrus palustris C.L. Koch, NotPal C * * * 8 5.00 2.03 3 1.88 0.91 1839 Opiella (Rhinoppia) nasuta OppNas C * * * 9 5.63 2.28 * * * (Moritz, 1965) Opiella (Rhinoppia) similifallax OppSim A * * * 23 14.38 5.84 * * * (Subias et Rodriguez, 1986) Opiella (Moritzoppia) OppTra B 2 1.25 1.15 * * * 12 7.50 3.64 translamelata (Willmann, 1923) Opiella uliginosa (Willmann, OppUli B * * * * * * 12 7.50 3.64 1919) Phthiracarus ferrugineus (C.L. PhtFer B * * * * * * 3 1.88 0.91 Koch, 1841) Phthiracarus laevigatus (C.L. PhtLae C * * * 1 0.63 0.25 * * * Koch, 1844) Pilogalumna crassiclava (Berlese, PilCra C * * * 1 0.63 0.25 * * * 1914) Pilogalumna tenuiclava (Berlese, PilTen B * * * * * * 8 5.00 2.42 1908) Platynothrus peltifer (C.L. Koch, PlaPel A * * * 1 0.63 0.25 36 22.50 10.91 1839) Podoribates longipes Berlese, PodLon C * * * 6 3.75 1.52 * * * 1908 Poecilochthonius italicus PoeIta− 1 0.63 0.57 * * * * * * (Berlese, 1910) Punctoribates punctum (C.L. PunPun C * * * 4 2.50 1.02 4 2.50 1.21 Koch, 1839) Ramusella (Insculptoppia) furcata RamFur A * * * 51 31.87 12.94 67 41.88 20.31 (Willmann, 1928) ardua (C.L. Koch, RhyArd C * * * 3 1.88 0.76 * * * 1841) Scheloribates laevigatus (C.L. SchLae A * * * 68 42.50 17.26 3 1.88 0.91 Koch, 1836) Scheloribates pallidulus (C.L. SchPal A * * * 54 33.75 13.71 * * * Koch, 1841) Steganacarus (Atropacarus) SteStr C * * * 2 1.25 0.51 1 0.63 0.30 striculus (C.L. Koch, 1835) Suctobelbella sarekensis SucSar C * * * 8 5.00 2.03 * * * (Forsslund, 1941) Tectocepheus velatus sarekensis TecSar− 171 106.88 98.28 5 3.13 1.27 * * * Trägårdh, 1910 Xenillus tegeocranus (Hermann, XenTeg C * * * 4 2.50 1.02 2 1.25 0.61 1804) 174 108.75 100 394 246.25 100 330 206.25 100

Overall, 37.5% of the species showed a tendency to come freely into the neighboring habi- tat and create numerically stable and comparable communities except Platynothrus peltifer (C.L. Koch, 1839), which showed a strong preference for MOF, and species (C.L. Koch, 1840) and Scheloribates laevigatus (C.L. Koch, 1836), preferring MOM. We surprisingly found O. (M.) translamelata, a moisture-loving species, in the MOD sample, but as it belongs to the species that prefer forest ecosystems, its long-term remaining in the meadow is questionable. Its co-occurrence with Poecilochthonius italicus (Berlese, 1910), which prefers the habitat’s meso-xerophilic nature, indicates one (or both) of these two species’ accidental occurrence in the MOD mite community. The MOM and MOF oribatid Appl. Sci. 2021, 11, 3537 8 of 14

communities showed comparable diversity (H’ = 2.44: 2.24; M = 4.87: 4.22) and equitability (J = 0.78: 0.75) representing all dominant categories. The low values of diversity (H’ = 0.10; M = 0.44) and equitability (J = 0.09) of the MOD community together with the uneven distri- bution of dominances in the community and the fact that no individual was found in the MOS substrate indicate extremely disturbed soils close to the ash pond. Despite the higher values of diversity in the natural soils at MOM and MOF sites than in the ash pond and its dam, we cannot consider the mite community’s state optimal. Mangová et al. [51] reported much higher diversity (H’ = 4.60: 4.19; M = 11.75: 9.74), even though there was lower equitability (J = 0.33: 0.27) in the various meadow and forest soils across Northern Slovakia. Similar values of H’ diversity were recorded for ruderal habitats and less managed urban parks (2.72: 2.27), where lower values of equitability (J = 0.48: 0.41) and higher values of diversity were also recorded (M = 12.90: 10.73) [52]. The oribatocoenoses of MOM and MOF can, thus, be considered as balanced, species-poorer communities. The correlation analysis outlined that the abundance and diversity of oribatid mites increase with decreasing the substrate’s pH. There was a positive correlation between mite abundance and diversity and soil moisture, and the same was observed with soil microbial respiration. Interestingly, the highest was the soil carbon-to-nitrogen rate and the lowest was mite diversity and abundance. We also found a negative correlation with mite abundance and increased arsenic, nickel, cobalt, and chromium concentrations. The risk elements demonstrably negatively affected mite diversity, primarily cobalt and nickel (Table4).

Table 4. The correlation between chemical and microbiological attributes and oribatid mites’ diversity

and abundance. A, mite abundance; D, mite diversity; pH (H2O/KCl), negative log of the hydrogen ion concentration; w, water content in the sample; MBC, microbial biomass carbon (dry weight); SMA, soil microbial activity (respiration); C/N, total organic carbon-and-nitrogen ratio. Significance labels: * p < 0.01; ** p < 0.05; ns, non-significant.

AD pH −0.92 ** −0.98 * w 0.84 * 0.91 ** MBC −0.66 ** −0.58 ** SMA 0.88 ** 0.96 * C/N −0.86 ** −0.96 * As −0.96 * −0.86 ** Cu ns −0.59 ** Ni −0.84 ** −0.96 * Co −0.73 ** −0.91 ** Ba ns −0.63 * Cr −0.65 ** −0.85 **

The interrelationships between species and coenoses, shown by cluster analysis (Figure3), divided the oribatid mite community into three basic clusters. The Cluster A included species reaching an eudominant and dominant position in the MOM and MOF sites, predominantly freely passing between sites, except for P.peltifer, O. (R.) similifallax, S. laevigatus, and Scheloribates pallidulus (C.L. Koch, 1841), which achieved a high dominance in only one of these sites. The Cluster B included species bound to the MOF site, preferring a forest-habitat type, of which only O. (M.) Translamellata occurred in a recent position at the MOD site. Grouping C represented a recent to subrecent species bound mainly to the meadow character of the MOM locality. Species Punctoribates punctum (C.L. Koch, 1839) and Xenillus tegeocranus (Hermann, 1804) had a similar position MOF locality. P.italicus and T.v. sarekensis, occurring at the MOD site, did not belong to any of the mentioned clusters; moreover, their joint occurrence appeared to be accidental (Figure3). Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 14

The interrelationships between species and coenoses, shown by cluster analysis (Fig- ure 3), divided the oribatid mite community into three basic clusters. The Cluster A in- cluded species reaching an eudominant and dominant position in the MOM and MOF sites, predominantly freely passing between sites, except for P. peltifer, O. (R.) similifallax, S. laevigatus, and Scheloribates pallidulus (C.L. Koch, 1841), which achieved a high domi- nance in only one of these sites. The Cluster B included species bound to the MOF site, preferring a forest-habitat type, of which only O. (M.) Translamellata occurred in a recent position at the MOD site. Grouping C represented a recent to subrecent species bound mainly to the meadow character of the MOM locality. Species Punctoribates punctum (C.L. Koch, 1839) and Xenillus tegeocranus (Hermann, 1804) had a similar position MOF locality. Appl. Sci. 2021, 11, 3537 P. italicus and T.v. sarekensis, occurring at the MOD site, did not belong to any of the men- 9 of 14 tioned clusters; moreover, their joint occurrence appeared to be accidental (Figure 3).

Figure 3. Dendrogram resulting from a clustering procedure using the Bray–Curtis method and

showing the similarities of the oribatid mite community at the study’s four sites. Figure 3. Dendrogram resulting from a clustering procedure using the Bray–Curtis method and showing4. the Discussion similarities of the oribatid mite community at the study’s four sites. Reduction of the abundance and species’ richness of oribatid mites in soils close to 4. Discussionindustrial factories is a well-known phenomenon. The effects of metal pollution from different sources on the community structure of soil arthropods have been demonstrated by many researchers [53–56]. The authors underlined that oribatids generally seemed to be more susceptible to toxic metals than most of the other arthropods [12,57–60]. A generally observed pattern was that communities were impoverished with increasing metal concentration since most species decrease in abundance, although some species may survive even the most polluted sites [61,62]. We can confirm this claim as we detected very low mite abundance at contaminated sites with one dominating species, T. v. sarekensis. The reduced number and low abundance of some species in the oribatid community at the MOD might indicate that this community is resistant, with reduced self-regulation possibilities. There are more reasons why oribatids’ abundance and diversity were strongly affected at the polluted sites. First, many pollutants quickly accumulate in the soil, which is the living space for many oribatid species. Second, fungal hyphae, known to accumulate metals effectively, are the primary food source for many species [63]. In MOS and MOD Appl. Sci. 2021, 11, 3537 10 of 14

samples, we repeatedly tried to cultivate microscopic fungi, which are considered a required component of juveniles’ diet of Oribatida, but without success. The pond’s and dam’s sediment did not favor the development of microscopic fungi, which also reflected in the abundance of oribatid mites. The MOS and MOD samples showed meaningfully lower diversity and equitability than meadow and forest samples. On the other hand, even though the diversity values were much higher in the natural soil than in the technogenic sediments, we can consider the oribatocoenosis at this site relatively poor. The abundance of mites and species’ number of Oribatida were also highly reduced in the MOS and MOD, gradually increasing toward the natural, unpolluted sites, MOM and MOF. To a certain degree, the soil protects mites against pollutants and, therefore, whole absence of Oribatida was observed only in highly polluted ash sediments. Our study revealed that the recorded values of pollutants (especially the concentrations of As, Cu, Ni, Co) influenced the mite communities’ structure in the ash pond and its dam (Table4). T. v. sarekensis was found to resist a wide range of toxic metal pollution. It was the only species that thrived in high abundance in heavily polluted samples of forming Technosol MOD. According to Siepel [63], it was the only oribatid species observed as comparatively abundant at highly polluted sites. T. v. sarekensis is considered an opportunistic herbo- fungivorous species. It is also well known for belonging to a group of parthenogenetic r-selected organisms, which can recolonize areas where sexual reproduction may be hin- dered by environmental conditions [64]. Its life strategy and high tolerance allow it to occupy free niches after the displacement of other species and acquire a highly dominant position in the community [4,65–68]. Most oribatid species are, although defined as K- selected organisms with low dispersion ability [69]. They cannot adapt to new conditions, escape easily from the stress of disturbance, or quickly recolonize disturbed habitats. Only a few eurypotent species from oribatocoenoses, usually Oppiidaea and , can survive in permanently stressed biotopes. The technogenic substrates often form an unfavorable microhabitat for mites because of the low amount of SOM, low humidity, and wide C/N range. At optimal humidity, the oribatids can survive even in poor substrates with a minimum SOM content; in substrates with a wide C/N range, they are limited mainly by the lack of SOM immobilized by soil microorganisms. We found a low (none) oribatid abundance and diversity in samples of MOS and MOD, where the C/N range was extensive even though the amounts of MBC and respiration suggested a good functioning microbial community. Moreover, microbial indicators’ values showed that despite the reduced organic matter available to the microorganisms, the microbial community could efficiently utilize the organic carbon, even more than in both unpolluted sites. The wide range of C/N, on the other hand, suggests that in the nitrogen-poor substrates, MOS and MOD decomposition processes were slowed down due to microbial immobilization of organic matter. This could be the other reason for the total absence or reduced numbers of mites in these kinds of substrates. Manu [17] demonstrated that, even though the heavy metal pollution in the study site was high, the total number of species and numerical abundance reached high values due to the high humidity and acidity of the substrate, which provided a more favorable habitat. We also confirmed a positive correlation between humidity and mite abundance (r = 0.84, p < 0.01). The abundance was meaningfully lower in MOD samples, where we measured significantly lower substrate moisture than in the natural soils. Therefore, it seems that the availability of soil water highly influences the performance of soil fauna. Depending on the species, different strategies can be used to cope with low soil-water potential. Soil microarthropods may be resistant to desiccation by reducing the water permeability of their integument. This strategy has been adopted by many epedaphic (surface-dwelling) species, whereas euedaphic (soil-dwelling) species have little or no control of evaporation across the integument [70]. In general, soil moisture’s positive effect on soil communities’ abundance has been emphasized [71,72]. Holmstrup et al. [72] presented a study that showed that the abundance of microarthropods in copper-contaminated soil was stimulated by enforced drought in a field experiment. It was expected that the Appl. Sci. 2021, 11, 3537 11 of 14

combination of copper and drought would harm soil microarthropods. We confirmed a negative correlation between mite diversity and copper (Table4), which, connected with the low humidity of the technogenic substrates MOS and MOD, could have had a considerable effect on the mite community. Soil acidity unlikely affects oribatid mite communities directly since most oribatid mite species tolerate even very acidic conditions [62,73]. Furthermore, soil microarthropods’ community structure has been relatively resistant to substantial soil pH changes [74]. The soil pH affects the colonization of forests by soil macrodecomposers, which affects oribatid mites via changes in the thickness of the litter layer since soil acidity may affect calcium availability. A negative correlation was observed between pH and mite abundance and diversity (Table4). Although we cannot ignore the impact of drought and various toxic metals on the mite community, our observations showed that, with increasing pH, the diversity of oribatid decreases. However, in the case of MOS and MOD samples, we are inclined to believe that the very low diversity values were caused by a combination of critical factors, e.g., higher pH, low humidity, and high concentrations of some risk elements.

5. Conclusions The results obtained in this study provide additional knowledge of the community structure of oribatid mites in the disturbed soil environment around a coal ash basin. The results confirm that the impoundment’s technogenic sediment and its dam contain high concentrations of various toxic elements, which, in connection with alkaline pH, low substrate humidity, and wide C/N range, has a negative effect on the abundance of oribatid mites. Even though the substrate’s microbial community seems to utilize the accessible organic carbon efficiently and survive without energy loss, the organic compound stays immobilized for other soil organisms. This suggestion is supported by changes in oribatid mite community composition, which was negatively reflected in their abundance and percentage of dominances in acarocoenoses. Despite such a substrate’s unsuitable properties, some oribatid species can withstand critical metal concentrations and establish an abundant community. Compared to the community’s composition in natural soils, there was an observed reduction in species’ diversity, exchanging long-lived species for short-lived ones. T. v. sarekensis was found to survive in a disturbed habitat in a small population and ultimately increase its abundance. In conclusion, we confirmed that high values of diversity and equitability occurring in natural soil habitats decrease significantly in localities with regular anthropogenic intervention. When protecting soil ecosystems, there must be the biodiversity of soil , emphasizing oribatid mites. The differential response of T. v. sarekensis to adverse conditions makes it a suitable indicator for assessing the degree of ecosystem vulnerability.

Author Contributions: Conceptualization, Z.F. and B.M.; methodology, Z.F., B.M. and M.C.;ˇ software, Z.F. and B.M.; validation, Z.F. and B.M.; formal analysis, Z.F. and B.M. investigation, Z.F. and B.M.; resources, Z.F. and B.M.; data curation, Z.F. and B.M.; writing—original draft preparation, Z.F. writing—review and editing, Z.F., B.M. and M.C.;ˇ visualization, Z.F. and B.M.; supervision, Z.F. and B.M.; project administration, B.M. and M.C.;ˇ funding acquisition, B.M. and M.C.ˇ All authors have read and agreed to the published version of the manuscript. Funding: This research was funded by Slovak Scientific Grant Agency VEGA (Vedecká grantová agentúra) projects 2/0111/18, 1/0658/19, 2/0096/19, and 1/0712/20. Institutional Review Board Statement: The authors declare that the accepted principles of ethical conduct have been followed strictly. Since, in the case of soil mite study, no ethical issues are required so far, our knowledge is clear. Informed Consent Statement: Not applicable. Data Availability Statement: Not applicable. Conflicts of Interest: The authors declare no conflict of interest. Appl. Sci. 2021, 11, 3537 12 of 14

References 1. Prasannakumar, N.R.; Kumar, K.P. Impact of Climate Change on Arthropod Diversity. In Arthropod Diversity and Conservation in the Tropics and Sub-Tropics; Chakravathy, A.K., Sridhara, S., Eds.; Springer Science + Bussines Media: Singapore, 2016; pp. 1–18. [CrossRef] 2. Siepel, H.; de Ruiter-Dijkman, E.M. Feeding guilds of oribatid mites based on their carbohydrase activities. Soil. Biol. Biochem. 1993, 25, 1491–1497. [CrossRef] 3. Behan-Pelletier, V.M. Oribatid mite biodiversity in agroecosystems: Role for bioindication. Agric. Ecosyst. Environ. 1999, 74, 411–423. [CrossRef] 4. Caruso, T.; Migliorini, M. A new formulation of the geometric series with applications to oribatid (Acari, Oribatida) species assemblages from human-disturbed Mediterranean areas. Ecol. Model. 2006, 195, 402–406. [CrossRef] 5. Caruso, T.; Migliorini, M. Micro-arthropod communities under human disturbance: Is taxonomic aggregation a valuable tool for detecting multivariate change? Evidence from Mediterranean soil oribatid coenoses. Acta Oecologica 2006, 30, 46–53. [CrossRef] 6. van Straalen, N.M. Evaluation of bioindicator systems derived from soil arthropod communities. Appl. Soil. Ecol. 1998, 9, 429–437. [CrossRef] 7. van Straalen, N.M.; Verhoef, H.A. The Development of a Bioindicator System for Soil Acidity Based on Arthropod pH Preferences. J. Appl. Ecol. 1997, 34, 217–232. [CrossRef] 8. Shevchenko, O.S.; Kolodochka, L.A. Species complexes of the oribatid mites (, Oribatei) in soils of urban street lawns with different pollution rates. Vestn. Zool. 2013, 47, 49–52. [CrossRef] 9. Frankovská, J.; Kordík, J.; Slaninka, I. Atlas of Remediation Methods for Contaminated Sites; State Geological Institute of Dionýz Štúr: Bratislava, Slovakia, 2010; p. 360. ISBN 978-80-89343-39-3. (In Slovak) 10. Murvanidze, M.; Mumladze, L.; Arabuli, T.; Kvavadze, E. Oribatid mite colonization of sand and manganese tailing sites. Acarologia 2013, 53, 203–215. [CrossRef] 11. Khan, A.G.; Kuek, C.; Chaudhry, T.M.; Khoo, C.S.; Hayes, W.J. Role of plants, mycorrhizae and phytochelators in heavy metal contaminated land remediation. Chemosphere 2000, 41, 197–207. [CrossRef] 12. Van Straalen, N.M.; Butovsky, R.O.; Pokarzhevskii, A.D.; Zaitsev, A.S.; Verhoef, S.C. Metal concentrations in soil and invertebrates in the vicinity of a metallurgical factory near Tula (Russia). Pedobiologia 2001, 45, 451–466. [CrossRef] 13. Köhler, H.R.; Alberti, G.; Seniczak, S.; Seniczak, A. Lead-induced hsp70 and hsp60 pattern transformation and leg malformation during postembryonic development in the oribatid mite, Archegozetes longisetosus Aoki. Comp. Biochem. Physiol. C-Toxicol. Pharmacol. 2005, 141, 398–405. [CrossRef][PubMed] 14. Seniczak, S.; Dabrowski, J.; Dlugosz, J. Effect of copper smelting air pollution on the mites (Acari) associated with young Scots pine forests polluted by a copper smelting works at Giogow, Poland.1. Arboreal mites. Water Air Soil. Pollut. 1997, 94, 71–84. [CrossRef] 15. Skubala, P.; Kafel, A. Oribatid mite communities and metal bioaccumulation in oribatid species (Acari, Oribatida) along the heavy metal gradient in forest ecosystems. Environ. Pollut. 2004, 132, 51–60. [CrossRef] 16. Visser, S. Management of microbial processes in surface mined land reclamation in Western Canada. In Soil Reclamation Processes, Microbial Analysis and Applications; Tate, R.L., Klein, D.A., Eds.; Marcel-Dekker: New York, NY, USA, 1985; pp. 203–241. 17. Manu, M.; Honciuc, V.; Neagoe, A.; Băncilă, R.I.; Iordache, V.; Onete, M. Soil mite communities (Acari: , Oribatida) as bioindicators for environmental conditions from polluted soils. Sci. Rep. 2019, 9, 20250. [CrossRef][PubMed] 18. Madej, G.; Kozub, M. Possibilities of using soil microarthropods, with emphasis on mites (Arachnida, Acari, Mesostigmata) in assessment of successional stages in a reclaimed coal mine dump (Pszów, S Poland). Biol. Lett. 2014, 51, 19–36. [CrossRef] 19. Berg, M.P.; De Ruiter, P.C.; Didden, W.; Janssen, M.P.M.; Schouten, A.J.; Verhoef, H.A. Community food web, decomposition and nitrogen mineralisation in a stratified Scots pine forest soil. Oikos 2001, 94, 130–142. [CrossRef] 20. Masaroviˇcová, M.; Slávik, I.; Koval’ková, J. Komplexný Monitoring Odkalíska SR (ˇcast’ 6); STU: Bratislava, Slovkia, 2008; p. 57. 21. IUSS Working Group WRB. World Reference Base for Soil Resources 2014, Update 2015. International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; World Soil Resources Reports; FAO: Rome, Italy, 2015; p. 106. 22. Tullgren, A. En enkel apparat för automatikst vittjande av sällgods. Tijdschr Entomol. 1917, 45, 97–100. 23. Ermilov, S.G. Laboratory Cultivation of Oribatid Mites of the Superfamily Crotonioidea (Acari, Oribatida) with the Purpose of Studying Their Development; Vector TiS. Press: Nizhny Novgorod, Russia, 2008; p. 54. 24. Ermilov, S.G. Morphology of juvenile stages of Gustavia microcephala (Acari, Oribatida, Gustaviidae). Acarina 2010, 18, 73–78. 25. Ermilov, S.G.; Ryabinin, N.A.; Anichkin, A.E. Morphology of juvenile instars of two oribatid mite species of the family Hermanni- idae (Acari, Oribatida). Zool Zh. 2012, 91, 657–668. [CrossRef] 26. Ermilov, S.G.; Weigmann, G.; Tolstikov, A.V. Morphology of adult and juvenile instars of Galumna obvia (Acari, Oribatida, ), with discussion of its taxonomic status. ZooKeys 2013, 357, 11–28. [CrossRef] 27. Ghiljarov, M.S.; Krivolutsky, D.A. Opredelitel Obitajushchikh v Pochve Kleshchei, Sarcoptiformes; Nauka press: Moscow, Russia, 1975; p. 492. 28. Krivolutsky, D.A.; Lebrun, P.; Kunst, M.; Akimov, I.A.; Bayartogtokh, B.; Vasiliu, N.; Golosova, L.D.; Grishina, L.G.; Karppinen, E.; Kramnoy, V.J.; et al. Oribatid Mites, Morphology, Development, Phylogeny, , Methods of Study and Characteristics of the Model Species Nothrus Palustris. C.L. Koch; Nauka Press: Moscow, Russia, 1995; p. 240. Appl. Sci. 2021, 11, 3537 13 of 14

29. Kunst, M. Nadkohorta Pancíˇrníci—Oribatei. In KlíˇcZvíˇreny CSSRˇ ; Daniel, M., Cernˇ ý, V., Eds.; Academia: Praha, Czechoslovakia, 1971; pp. 531–580. 30. Olszanowski, Z. A monograph of the and Camisidae of Poland (Acari: Oribatida: Crotonoidea). Int. J. Invertebr. Taxon. 1996, 5, 1–201. 31. Pavlitshenko, P.G. A guide to the Ceratozetoid mites (Oribatei, Ceratozetoidea) of Ukraine; Natsional'noi Akademii nauk Ukrainy: Kiew, Ukraine, 1994; p. 142. 32. Pfingstl, T.; Krisper, G. No difference in the juveniles of two Tectocepheus species (Acari: Oribatida, ). Acarologia 2011, 51, 199–218. [CrossRef] 33. Seniczak, S. Morphology of juvenile stages of some (Acarina, Oribatei) I. Pedobiologia 1975, 15, 249–261. 34. Seniczak, S. The morphology of juvenile stages of moss mites of the family (Acari, Oribatei) I. Acta Zool. Cracov. 1980, 24, 487–500. 35. Seniczak, S. The morphology of the juvenile stages of moss mites of the family Scheloribatidae (Acari, Oribatida) II. Ann. Zool. 1990, 43, 299–308. 36. Seniczak, S.; Norton, R.A.; Seniczak, A. Morphology of Eniochthonius minutissimus (Berlese, 1904) and Hypochthonius rufulus C.L. Koch, 1835 (Acari: Oribatida: Hypochthonioidea). Ann. Zool. 2009, 59, 373–386. [CrossRef] 37. Seniczak, S.; Seniczak, A.; Kaczmarek, S.; Zelazna, E. Systematic status of Oribatula Berlese, 1895 (Acari: Oribatida: ) in the light of the ontogeny of three species. Int. J. Acarol. 2012, 38, 664–680. [CrossRef] 38. Seniczak, A.; Seniczak, S. Morphology of juvenile stages of Pilogalumna crassiclava (Berlese, 1914) and P. ornatula Grandjean, 1956 (Acari: Oribatida: Galumnidae). Ann. Zool. 2007, 57, 841–850. [CrossRef] 39. Seniczak, S.; Seniczak, A. Morphology of three species of Crotonioidea Thorell, 1876 (Acari: Oribatida), and relations between some genera. Zool. Anz. 2009, 248, 195–211. [CrossRef] 40. Subías, L.S.; Balogh, P. Identification keys to the genera of Oppiidae Grandjean, 1951 (Acari: Oribatei). Acta. Zool. Hung. 1989, 35, 355–412. 41. Vance, E.; Brookes, P.; Jenkinson, D. An extraction method for measuring soil microbial biomass C. Soil Biol. Biochem. 1987, 19, 703–707. [CrossRef] 42. Schinner, F.; Öhlinger, R.; Kandeler, E. Bodenbiologische Arbeitsmethoden; Metzler, J.B., Ed.; Springer: Berlin, Germany, 1991; p. 385. 43. Dilly, O.; Munch, J.-C. Ratios between estimates of microbial biomass content and microbial activity in soils. Biol. Fertil. Soils 1998, 27, 374–379. [CrossRef] 44. Cheng, W.; Zhang, Q.; Coleman, D.C.; Carroll, C.R.; Hoffman, C.A. Is available carbon limiting microbial respiration in the rhizosphere? Soil Biol. Biochem. 1996, 28, 1283–1288. [CrossRef] 45. Anderson, T.-H. Microbial eco-physiological indicators to asses soil quality. Agric. Ecosyst. Environ. 2003, 98, 285–293. [CrossRef] 46. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST–Paleontological statistics software package for education and data analysis. Paleontol. Electron. 2001, 4, 9. [CrossRef] 47. Tischler, W. Synökologie der Landtiere; Gustav Fischer Verlag: Stuttgart, Germany, 1955; p. 414. 48. Shannon, C.E.; Weaver, W. The Mathematical Theory of Communication; The University of Illinois Press: Urbana, IL, USA, 1964; p. 146. 49. Margalef, R. Information theory in ecology. Gen. Syst. 1958, 3, 36–71. 50. Sheldon, A.L. Equitability indices: Dependence on the species count. Ecology 1969, 50, 466–467. [CrossRef] 51. Mangová, B.; Krumpál, M.; L’uptáˇcik,P. Pancierniky (Acari: Oribatida) vybraných lokalít Turskej doliny (Lúˇcanská Malá Fatra). Folia faun. Slov. 2012, 17, 261–272. 52. Mangová, B.; Hulejová Sládkoviˇcová, V.; Krumpál, M.; Kozánek, M. The impact of different urban conditions on structural characteristics of oribatid mite communities. Biologia 2019, 74, 153–168. [CrossRef] 53. Salminen, J.; Korkama, T.; Stroömmer, R. Interaction modification among impairs ecosystem processes in lead- polluted soil. Environ. Toxicol. Chem. 2002, 21, 2301–2309. [CrossRef][PubMed] 54. Seniczak, A.; Seniczak, S. The effect of cadmium on Archegozetes longisetosus (Acari, Oribatida) in laboratory conditions. Europ. J. Soil Biol. 2002, 38, 315–317. [CrossRef] 55. Seniczak, S.; Kaczmarek, S.; Klimek, A.; Seniczak, A. The effect of some air pollutants on the vertical distribution of mites (Acari) in soils of young Scots pine forests in Poland. Europ. J. Soil Biol. 2002, 38, 311–314. [CrossRef] 56. Zaitsev, A.S.; Van Straalen, N.M. Species diversity and metal accumulation in oribatid mites (Acari, Oribatida) of forests affected by a metallurgical plant. Pedobiologia 2001, 45, 467–479. [CrossRef] 57. Skubala, P.; Zaleski, M. Heavy metal sensitivity and bioconcentration in Oribatid mites (Acari, Oribatida) Gradient study in meadow ecosystems. Sci. Total Environ. 2012, 414, 364–372. [CrossRef][PubMed] 58. Haimi, J.; Mätäsniemi, L. Soil community in heavy-metal contaminated coniferous forest with and without liming. Europ. J. Soil Biol. 2002, 38, 131–136. [CrossRef] 59. Migliorini, M.; Pigino, G.; Bianchi, N.; Bernini, F.; Leonzio, C. The effect of heavy metal contamination on the soil arthropod community of a shooting range. Environ. Pollut. 2004, 129, 331–340. [CrossRef][PubMed] 60. Zvereva, E.L.; Kozlov, M.V. Responses of terrestrial arthropods to air pollution: A meta-analysis. Environ. Sci. Pollut. Res. 2010, 17, 297–311. [CrossRef][PubMed] Appl. Sci. 2021, 11, 3537 14 of 14

61. Feketeová, Z.; Mangová, B.; Hulejová Sládkoviˇcová, V.; Šimkovic, I. Biological activity of the metal-rich post-flotation tailings at an abandoned mine tailings pond (four decades after experimental afforestation). Environ. Sci. Pollut. Res. 2015, 22, 12174–12181. [CrossRef] 62. Feketeová, Z.; Mangová, B.; Hulejová Sládkoviˇcová, V.; Pogányová, A.; Šimkovic, I.; Krumpál, M. Biological properties of extremely acidic cyanide-laced mining waste. Ecotoxicology 2016, 25, 202–212. [CrossRef] 63. Siepel, H. Applications of microarthropod life-history tactics in nature management and ecotoxicology. Biol. Fertil. Soils 1995, 19, 75–83. [CrossRef] 64. Błoszyk, J.; Adamski, Z.; Napierała, A.; Dylewska, M. Parthenogenesis as a life strategy among mites of the suborder Uropodina (Acari: Mesostigmata). Can. J. Zool. 2004, 82, 1503–1511. [CrossRef] 65. Honciuc, V. Edaphic mite populations (Acari, Oribatida). In Species Monitoring in the Central Parks of Bucharest; Onete, M., Ed.; Editura Ars Docenti: Bucharest, Romania, 2008; pp. 53–67. 66. Honciuc, V.; Manu, M. Ecological study on the edaphic mites populations (Acari: Mesostigmata-Gamasina, Oribatida) in urban areas from Romania. Rom. J. Biol. Zool. 2010, 55, 3–18. 67. Honciuc, V.; Manu, M. Ecological research on soil mites populations (Acari: Mesostigmata-Gamasina, Oribatida) from forest ecosystems near Bucharest city. Rom. J. Biol. Zool. 2010, 55, 19–30. 68. Niedbała, W.; Błoszyk, J.; Kaliszewski, M.; Ka´zmierski,A.; Olszanowski, Z. Structure of soil mite (Acari) communities in urban green of Warsaw. Fragm. Faun. 1990, 33, 21–44. [CrossRef] 69. Gulvik, M.E. Mites (Acari) as indicators of soil biodiversity and land use monitoring: A review. Polish J. Ecol. 2007, 55, 415–440. 70. Verhoef, H.A. Water balance in Collembola and its relation to habitat selection: Water content, haemolymph osmotic pressure and transpiration during an instar. J. Physiol. 1981, 27, 755–760. [CrossRef] 71. Frampton, G.; Van den Brink, P.; Gould, P. Effect of spring drought and irrigation on farmland arthropods in southern Britain. J. Appl. Ecol. 2001, 37, 865–883. [CrossRef] 72. Holmstrup, M.; Maraldo, K.; Krogh, P.H. Combined effect of copper and prolonged summer drought on soil Microarthropods in the field. Environ. Pollut. 2007, 146, 525–533. [CrossRef] 73. Hågvar, S.; Abrahamsen, G. Microarthropods and Enchytraeidae (Oligochaeta) in naturally lead-contaminated soils: A gradient study. Environ. Entomol. 1990, 19, 1263–1277. [CrossRef] 74. Liiri, M.; Haimi, J.; Setälä, H. Community composition of soil microarthropods of acid forest soils as affected by wood ash application. Pedobiologia 2002, 46, 108–124. [CrossRef]