InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Environmental Sciences » "Emerging Pollutants in the Environment - Current and Further Implications", book edited by Marcelo L. Larramendy and Sonia Soloneski, ISBN 978-953-51-2160-2, Published: September 2, 2015 under CC BY 3.0 license. © The Author(s).

Chapter 2

Heavy Metals — Soil Communities at Ecological Risk

By P. Šalamún, V. Hanzelová and T. Brázová
DOI: 10.5772/60404

Article top


Maps of the study areas with indication of the emission source and sampling points (white arrows - K1-K4 for Krompachy, J1-J4 for Jelšava); prevailing winds are indicated in the compass rose.
Figure 1. Maps of the study areas with indication of the emission source and sampling points (white arrows - K1-K4 for Krompachy, J1-J4 for Jelšava); prevailing winds are indicated in the compass rose.
Redundancy analysis (RDA) performed on physicochemical soil properties, trace elements and selected indices in relation to sampling sites with data explained 73.1 % of the variation in the first two axes; F=10.9; P=0.002 (Krompachy) [1].
Figure 2. Redundancy analysis (RDA) performed on physicochemical soil properties, trace elements and selected indices in relation to sampling sites with data explained 73.1 % of the variation in the first two axes; F=10.9; P=0.002 (Krompachy) [1].
Structure and Enrichment Index graphic representation (Krompachy) [1]. A - disturbed food web, bacterial decomposition; B - maturing food web, balanced decomposition; C - structured food web, fungal decomposition; D - degraded food web, fungal decomposition [11].
Figure 3. Structure and Enrichment Index graphic representation (Krompachy) [1]. A - disturbed food web, bacterial decomposition; B - maturing food web, balanced decomposition; C - structured food web, fungal decomposition; D - degraded food web, fungal decomposition [11].
Redundancy analysis (RDA) performed on physicochemical soil properties, trace elements and selected indices in relation to sampling sites with data explained 64.0 % of the variation in the first two axes; F=7.1; P=0.002 (Jelšava) [53].
Figure 4. Redundancy analysis (RDA) performed on physicochemical soil properties, trace elements and selected indices in relation to sampling sites with data explained 64.0 % of the variation in the first two axes; F=7.1; P=0.002 (Jelšava) [53].
Structure and Enrichment Index graphic representation (Jelšava) [53]. A - disturbed food web, bacterial decomposition; B - maturing food web, balanced decomposition; C - structured food web, fungal decomposition; D - degraded food web, fungal decomposition [11].
Figure 5. Structure and Enrichment Index graphic representation (Jelšava) [53]. A - disturbed food web, bacterial decomposition; B - maturing food web, balanced decomposition; C - structured food web, fungal decomposition; D - degraded food web, fungal decomposition [11].

Heavy Metals — Soil Communities at Ecological Risk

P. Šalamún1, V. Hanzelová1 and T. Brázová1

1. Introduction

The quality of human life depends on the chemical composition of the environment, including food, water, and many other factors necessary for life. The development of human society, however, has brought substantial increases in ecological problems and contamination of the environment [1]. Some parts of the environment (e.g. atmosphere and hydrosphere) are able to partially remove the contaminants or to disperse them over wide areas and eliminate their potentially ecotoxic effects. The pedosphere, however, cannot, and important “hot-spots” with extremely high levels of a variety of pollutants have developed concurrently with human society in recent decades [2]. Soil is one of the most important factors for terrestrial life through the production of food, so its degradation could reduce the diversity or abundance of various valuable organisms necessary for maintaining balanced ecological functions [3].

Understanding the impact of pollution on important soil processes, where soil organisms play key roles, requires observing not only the quantity of the pollutants in a system, but also their qualitative effects and ability to change the environmental conditions for soil biota [4]. Various approaches using different soil communities were thus developed for biomonitoring potentially endangered ecosystems. The wide variety of communities, either directly in the soil or above it, tightly linked to soil processes makes the selection of an appropriate bioindicator challenging. Soil mesofauna, composed of various invertebrate communities, are the most suitable [5]. Soil invertebrate communities are able to buffer short-term resource shortages without significant losses in abundance or diversity [6], but they also have relatively rapid and identifiable reactions to contamination, because they are integrators of the physical, chemical, and biological properties related to their food resources [5]. Within the mesofauna, soil nematodes [7], springtails [8], and mites [9] have been considered for use as reliable bioindicators. Nematodes may be the most promising candidates for systems of bioindication, because more information on their taxonomy and ecology [10] is available than for other animal groups.

1.1. Nematodes as bioindicators

Nematodes are amongst the most abundant multicellular organisms in nearly all soil ecosystems, occupying key positions at most trophic levels in the soil food web [11]. They are involved in organic matter mineralisation, plant growth, and crop productivity by their regulation of detritivore populations and thus contribute to the stability of the food web [12, 13, 14]. Their life cycles have been widely studied, partly because some taxa are pests [15]. The nematode Caenorhabditis elegans is a key model for genetic studies, which could significantly improve the bioindicating potential of nematodes on molecular and gene-expression levels.

The introduction of the maturity index by [7], one of the best-known indicators of soil health [15], was a crucial breakthrough in nematode bioindication. Nematodes are classified in this system along a coloniser-persister (c-p) continuum and assigned to c-p groups ranging from 1 (colonisers) to 5 (persisters). Similar to the r/K dichotomy, colonisers are characteristic of early successional stages, and persisters are found mainly in more developed and complex ecosystems. The advantage of the c-p system is that nematodes only need to be identified to the generic level, which does not affect the resolution of the indication.

Several other derived or new ecological indices have been developed from the maturity index. Ferris et al. [11] provided a framework for determining the enrichment (enrichment index, EI) and structure (structure index, SI) of food webs based on the relative weighted abundance of the functional guilds of nematodes. A matrix of feeding habits, with life history characteristics embodied in the c-p classification, is used to enhance the resolution of nematode analyses and to connect the functional groups with the enrichment and structure of the food web and channels of decomposition.

We applied these approaches at two locations, Kovohuty JSC Krompachy and Slovak Magnesite Works JSC Jelšava, with soils contaminated by dust fallout containing high amounts of potential ecotoxic elements. The levels of contamination at these sites have been analytically determined, but little is known about the soil communities and their interactions or about the general health of the soil. The aims of this study were thus to improve our knowledge of the polluted sites and to better understand the biotic processes under environmental pressure and the consequences of industrial pollution for soil food chains.

The specific objectives were

  1. To clarify the impact of contamination on the quality of soil abiotic and biotic components.

  2. To evaluate the ability of soil nematode communities to detect ecosystemic disturbances.

  3. To define the suitability of this system of bioindication for the detection of soil degradation.


Figure 1.

Maps of the study areas with indication of the emission source and sampling points (white arrows - K1-K4 for Krompachy, J1-J4 for Jelšava); prevailing winds are indicated in the compass rose.

2. Material and methods

Study areas

The study was performed in the model towns of Krompachy and Jelšava that have important emission sources, smelters Kovohuty JSC Krompachy and Slovak Magnesite Works Jelšava (SMZ JSC). Krompachy has been known for its processing of iron and copper ore since the 14th century, and copper has been distributed all over Europe. The magnesite industry in Jelšava has flourished since the end of the 19th century. Soil contamination in both areas is caused by the fallout of solid particles containing various elements during the mining and processing of the ores [16]. More than an estimated 1500 tonnes of solid particles yearly containing accessorial elements (Cd, Cr, Cu, Pb, and Ni) from Kovohuty JSC Krompachy, and as many as 165, 000 tonnes of emissions containing magnesium in total from SMZ JSC were emitted to the atmosphere between the 1960s and the 1990s [17].

Sampling and data analysis

Soil samples were collected in April 2010 for Krompachy and in May 2009 for Jelšava at four sampling sites each along transects downwind of the emission sources (Figure 1). Four 1 kg replicates (each consisting of four subsamples) were collected from the surface horizons (0-20 cm) of permanent grasslands at each sampling site by quadrat sampling. Each replicate was stored in a plastic bag at 4 °C until analysed. The soil samples were analysed for geochemical properties and nematode community structure.

2.1. Geochemical properties and heavy metals

The soil samples were processed as follows:

  1. Soil moisture was measured gravimetrically by drying the replicates to a constant weight at 105 °C [18].

  2. Organic matter content (Corg) was determined by titration with K2Cr2O7/H2SO4 [18].

  3. NH4+ was determined spectrophotometrically using Nessler reagent [18].

  4. NO3- was determined by an ion-selective electrode [18].

  5. pH was determined in a solution of CaCl2 [18].

  6. The total and mobilisable concentrations of heavy metals (As, Cd, Cu, Cr, Mg, Ni, Pb, and Zn) were determined by inductively coupled plasma mass spectrometry on an Agilent 7500 C (Agilent Technologies, USA) following the manufacturer’s instructions. The replicate samples were air-dried before analysis in an oven at 30 °C for 48 h or until constant weight, ground, and sifted through a 0.2-mm sieve. Metals were extracted by 2M HNO3 (10 g of soil in 50 mL of HNO3) for 6 h for total content and by 0.05M Na2EDTA neutral solution (5 g of soil in 50 mL of Na2EDTA) for mobilisable content [19].

Steps 1-5 were performed by the Laboratory of the Central and Testing Institute in Agriculture in Košice, Slovakia according to certified methods [18].

Krompachy Sampling site
K1 K2 K3 K4
GPS position 48°55'01"N, 20°53'07"E48°55'17"N, 20°54'13"E48°55'52"N, 20°54'35"E48°54'56"N, 20°55'55"E
Distance (m) 750166022503820
Altitude (m) 460380407404
pH 4.90±0.125.75±0.215.70±0.375.15±0.25
C org (%)2.52±0.873.92±0.403.56±0.412.90±0.38
NO 3 - (±1.775.83±0.535.63±1.814.53±1.94
NH 4 + (±1.3712.10±2.018.78±2.7619.73±8.12
Ninorg. (±2.9317.93±2.1314.40±4.5624.25±9.87
Moisture (%) 18.87±3.1728.43±1.3519.47±0.7322.77±3.11
Jelšava J1 J2 J3 J4
GPS position 48°38'39''N, 20°13'40''E48°38'05''N, 20°14'46''E48°36'44''N, 20°15'26''E48°35'49''N, 20°15'13''E
Distance (m) 487 m2196 m4630 m6096 m
Altitude (m) 338 m326 m260 m261 m
pH > 9.07.4±0.17.5±0.16.9±0.1
C org (%)2.4±0.12.0±0.42.3±0.41.9±0.2
NO 3 - (±0.94.8±3.05.5±1.02.3±0.2
NH 4 + (±1.515.0±2.312.2±3.014.0±2.2
Ninorg. (±2.319.8±5.017.7±3.516.3±2.3
Moisture (%) 22.0±5.523.3±1.018.7±1.219.9±0.9

Table 1.

Study sites characteristics [1, 53].

[i] - * - by Slovak Soil Science and Conservation Research Institute

Loamy Cambisols with intense stony skeletons dominated in both areas. Despite some variations in soil properties, the basic characteristics of the soils were similar (Table 1), with soil moistures ca. 20% and inorganic nitrogen contents (Ninorg.) between 10 and 20 mg kg-1. The Krompachy soils were richer in Corg and were more acidic than the Jelšava soils (soil pH >9).

2.2. Analysis of nematode community structure

Nematode communities were isolated by extraction from 100 g of each replicate soil sample using a modified Baermann procedure [20]. The extracted nematodes were fixed in Ditlevsen’s solution (formalin/acetic acid/alcohol) [21] and were counted using a microscope. All nematodes in the samples were classified to order, family, and genus. The nematode communities were analysed for (i) total abundance (number of individuals per 100 g of soil), (ii) trophic structure (bacterial feeders; fungal feeders; plant feeders; omnivores; and predators) [22], and (iii) the ecological indices (1) MI2-5 [23] with Nemaplex c-p values, (2) the Shannon-Weaver index (H´) for generic diversity [24], (3) generic richness, expressed as mean number of genera at a site, (4) EI based on the proportion of opportunistic bacterivorous and fungivorous nematodes, (5) SI based on the relative weighted abundance of guilds representing the structure of the food web, and (6) the channel index (CI) [11].

Trace elements Sampling site
K1 K2 K3 K4 Limit
Total concentration (HNO3 extraction)
Arsenic (As) 78.654±15.287a 19.574±5.402b 24.377±5.685b 5.609±1.306c 5
Cadmium (Cd) 1.494±0.181a 1.061±0.478b 0.339±0.062c 0.161±0.013c 0.3
Chrome (Cr) 5.781±7.019a 1.458±0.159a 1.475±0.206a 1.012±0.021a 10
Copper (Cu) 257.643±38.452a 175.398±80.216b 65.796±14.021c 18.322±2.129c 20
Nickel (Ni) 2.617±0.262b 3.598±0.349a 0.627±0.314c 1.021±0.027c 10
Lead (Pb) 25.140±8.380a 29.038±11.551a 8.216±3.761b 5.195±0.089b 30
Zinc (Zn) 80.62±18.019a 60.253±29.127a 8.067±4.449b 3.201±0.290b 40
Mobilisable fraction (Na2EDTA extraction)
Arsenic (As) 4.779±0.688a 2.870±0.690b 1.074±0.494c 0.105±0.021d n.a.
Cadmium (Cd) 0.407±0.101a 0.266±0.097b 0.066±0.016c 0.023±0.003c n.a.
Chrome (Cr) 0.081±0.02a 0.065±0.031a 0.073±0.036a 0.010±0.003b n.a.
Copper (Cu) 76.660±23.554a 46.945±18.585b 14.320±3.788c 4.380±0.011c n.a.
Nickel (Ni) 0.54±0.115a 0.568±0.112a 0.167±0.197b 0.101±0.008b n.a.
Lead (Pb) 18.565±2.744a 16.293±7.951a 5.413±1.681b 4.055±0.337b n.a.
Zinc (Zn) 63.285±17.193a 26.945±12.502b 3.600±1.377c 0.643±0.085c n.a.

Table 2.

Total and mobilisable concentration of trace elements in sampling sites from Krompachy ( [1].

[i] - abcd Means followed by the same letters on the same rows are not statistically different according to Least Significant Difference Test (P<0.05)

[ii] - Limit - limits posted by The Decree of the Ministry of Land Management of the Slovak Republic No. 531/1994-540 on the admissible values of harmful substances in uncontaminated soil. n.a. – not available

2.3. Data analysis

Spearman’s nonparametric correlation coefficient (rs) was calculated using STATISTICA v. 9.0 to test the relationships between the characteristics of a nematode community and the concentrations of mobilisable heavy metals at a site [25]. Correlations at P<0.05 and P<0.01 were considered significant. Differences in the mean traits and indices of a community amongst sites were tested by Duncan’s tests. We used a constrained ordination redundancy analysis (RDA) in CANOCO 5 to analyse the ecological distances between sites (nematode community and soil parameters). The significance of an axis was tested by Monte Carlo permutation [25].

The effects of contamination on soil ecosystems can be categorised as direct or indirect. Alterations in the soil communities near Kovohuty JSC Krompachy were likely due to direct toxicity from the high levels of heavy metals in the soil samples. The contamination acted mostly indirectly near SMZ JSC, altering the basic soil properties. We will present and discuss the results from these two areas separately.


Figure 2.

Redundancy analysis (RDA) performed on physicochemical soil properties, trace elements and selected indices in relation to sampling sites with data explained 73.1 % of the variation in the first two axes; F=10.9; P=0.002 (Krompachy) [1].

3. Direct toxic effects of industrial emissions – A case study at Kovohuty JSC Krompachy

3.1. Results

3.1.1. Heavy metals

The total heavy metal content (except Pb and Ni) along the Krompachy transect was highest at the K1 site and decreased downwind (Table 2). The patterns of the decreases in concentration, however, varied amongst the metals. As, Cd, and Cu decreased significantly (P<0.05) and relatively continuously towards K4. The concentrations of Ni, Pb, and Zn were relatively high at sites K1 and K2 but were significantly lower at K3 (P<0.05). Cr content did not vary significantly (P>0.05) along the transect. The mobile fractions of all metals had similar decreasing trends (Table 2). As, Cd, Cu, and Zn had high total concentrations, exceeding the limits for uncontaminated soils in Slovakia [26] (Table 2). Corg and pH were likely the most important factors influencing the mobility of the metals, and the RDA (Figure 2) indicated that they were negatively correlated with the concentrations of all heavy metals except Cr.

3.1.2. Composition of the nematode fauna: trophic and c-p groups

A total of 58 genera were identified, including 20 bacterivores, five fungivores, 10 omnivores, eight predators, and 15 plant feeders (Table 3). The most common genera in the communities were Acrobeloides, Aphelenchoides, Geocenamus, Helicotylenchus, Paratylenchus, Pratylenchus, and Rhabditis (Table 3). The mean number of genera (17.75) in the replicates was lowest near the industrial plant (K1), which also had the lowest mean population density of 386.5 per 100 g of soil (Table 4). Both parameters increased with distance from the plant to 35.75 and 833.25, respectively, at K4. Ninorg. was significantly correlated (P<0.01) with generic richness, and soil pH (P<0.01) and Ninorg. (P<0.05) were significantly correlated with nematode abundance.

The industrial complex affected the densities of the trophic groups. The trophic structure differed at each sampling site (Table 4). Plant feeders were the most abundant trophic group, with a proportion in the community near 66% at the most contaminated site, K1. Bacterivores were the second most common group, representing more than 50% of the total nematode populations at K2 and K3. This group was positively correlated with Corg and pH (P<0.01) (Table 5). Fungivores, considered relatively insensitive to disturbances, including heavy metal pollution, occurred in relatively low densities, except at K2 with 8.98% of the total population. Omnivores and predators were most affected by the pollutants, with very low densities at K1-K3. The low resistance of these two groups to disturbances is in agreement with their negative correlation (P<0.01) with Cr (predators) or with all heavy metals (omnivores). Both trophic groups were also negatively correlated with Ninorg., and soil moisture (Table 5). The distribution of the c-p groups near the plant was consistent with their life strategies and survival abilities. General opportunistic nematodes, c-p 2, were most common at the most polluted sites (K1-K3), followed by c-p 1. The other three c-p groups represented only ca. 10% of the identified nematodes. The numbers of c-p 4 and 5 nematodes were significantly higher at the most remote site, K4, where representatives of c-p 4 were the most abundant within the nematode community (Table 4).

3.1.3. Ecological indices

Mean H′ and generic richness were highest at K4, indicating the richest composition of nematode taxa (Table 4). MI2-5 and SI, measures of ecosystemic maturity, showed analogous development (Figure 2), with stable values from K1 towards K3 and steep increases at K4 (Table 4). EI values indicated that K2 stored the highest amounts of organic matter available for degradation, but K1 had the highest rate of degradation, as indicated by the lowest CI value (P<0.05) (Table 4).

Nematode genera K1 K2 K3 K4
c-p Abundance Abundance Abundance Abundance
Bacterial feeders
Acrobeloides 244.00±26.1584.25±25.13202.3±39.830.00±13.04
Alaimus 41.00±2.000.50±0.580.75±1.5015.25±12.74
Amphidelus 4---7.75±15.50
Anaplectus 2---0.75±1.50
Aulolaimus 3-0.50±1.00--
Cephalobus 27.50±4.0444.75±6.6582.75±36.7932.00±26.77
Diplogaster 1---2.50±4.36
Eucephalobus 23.25±4.7216.00±18.49-9.50±10.66
Eumonhystera 1---1.50±1.73
Geomonhystera 1---0.75±0.96
Heterocephalobus 22.00±2.8331.25±29.6236.00±25.343.00±2.58
Mesorhabditis 17.25±13.20--1.25±1.50
Monhystera 11.50±3.0022.50±11.598.25±7.54-
Panagrolaimus 11.00±0.82--9.25±13.84
Plectus 25.25±3.3020.00±16.1923.00±10.6827.75±3.40
Prismatolaimus 31.50±2.387.00±10.920.25±0.509.25±6.34
Protorhabditis 14.75±9.50---
Rhabditis 132.75±18.30286.5±63.241.75±24.5155.00±16.25
Teratocephalus 310.00±9.9011.00±12.685.25±3.770.25±0.50
Wilsonema 2-2.25±1.710.50±0.580.75±1.50
Fungal feeders
Aphelenchus 20.25±0.5083.75±43.663.75±4.2715.00±7.79
Diphtherophora 3--3.50±3.00-
Dorylaimoides 4---4.50±9.00
Filenchus 21.25±1.893.00±6.00-11.25±10.18
Tylencholaimus 4---4.25±7.85
Aporcelaimellus 5-0.75±1.504.00±4.8337.00±9.31
Aporcelaimus 51.00±2.00-0.75±1.5021.25±25.99
Campydora 4---2.75±5.50
Dorylaimus 42.00±2.834.50±6.4511.25±9.6045.75±30.61
Enchodelus 4---0.75±1.50
Epidorylaimus 4--4.00±4.69-
Eudorylaimus 4-1.25±2.5013.50±19.3634.50±7.33
Mesodorylaimus 5--1.50±3.002.00±1.63
Microdorylaimus 4---3.75±2.63
Prodorylaimus 5---1.00±2.00
Anatonchus 4-0.50±1.00-2.25±3.30
Clarkus 4-4.50±9.000.25±0.5014.25±12.50
Mononchus 41.75±2.3630.25±23.773.25±4.030.75±1.50
Mylonchulus 4-1.75±2.87-8.75±5.32
Nygolaimus 5---1.00±2.00
Oxydirus 5---17.25±16.09
Prionchulus 4-5.25±8.54-0.25±0.50
Tripyla 3-2.00±2.450.75±1.5015.75±7.41
Plant feeders
Aglenchus 25.50±1.7327.50±21.556.00±3.7431.50±21.58
Aphelenchoides 28.75±7.80147.5±74.968.75±23.4310.50±13.10
Axonchium 51.00±2.00-0.75±1.50-
Boleodorus 210.00±9.0620.75±9.6415.00±10.1017.25±3.30
Coslenchus 2--0.50±1.001.25±1.89
Criconema 3-0.75±0.9616.25±30.511.25±2.50
Geocenamus 3149.8±129.82.00±2.835.00±3.922.75±1.26
Helicotylenchus 31.25±0.5073.75±100.134.00±34.81218.8±208.6
Heterodera 31.25±2.50--20.75±25.18
Macrotrophurus 3---0.50±1.00
Malenchus 20.50±1.002.25±4.502.25±3.867.50±11.21
Paratylenchus 214.50±14.1531.25±30.61159.0±177.921.25±19.10
Pratylenchus 366.00±69.3394.00±54.3089.00±94.7053.00±52.35
Rotylenchus 3---0.50±1.00
Tylenchus 2--0.50±1.000.50±1.00

Table 3.

C-p values, the average abundance (±SD), and dominance of nematode genera in individual sampling sites from Krompachy [1].

Several significant relationships between the heavy metal contents in the soil and the ecological indices were observed across all sampled sites (Table 5). As an indicator of diversity, generic richness was negatively correlated (P<0.05 or <0.01), with nearly all heavy metals except Ni, similar to H′ (except Ni and Pb). SI and MI2-5 were negatively correlated (P<0.01) with Cr content.

Indices, trophic & c-p groups Sampling sites
K1 K2 K3 K4
Trophic & c-p groups
Bacterial feeders 31.87±7.32a 50.25±5.60b 51.83±14.57b 26.12±8.59a
Fungal feeders 0.36±0.44b 8.98±5.10a 0.98±0.48b 4.21±0.87b
Omnivores 1.19±1.30c 0.79±1.00c 3.35±2.43c 17.90±5.71b
Predators 0.71±1.21c 4.00±2.20bc 0.47±0.67c 7.50±2.99ab
Plant feeders 65.86±8.88b 35.99±9.23a 43.37±12.62a 44.27±12.60a
c-p 1 33.86±10.47c 46.80±6.56d 10.42±4.68a 15.63±4.84ab
c-p 2 51.93±7.77cd 42.67±6.24bc 79.76±11.88e 29.73±10.60ab
c-p 3 8.39±5.75b 3.00±2.06ab 2.04±1.49a 5.78±2.73ab
c-p 4 5.39±7.82b 7.42±4.01b 6.51±5.40b 31.19±12.30a
c-p 5 0.43±0.86c 0.12±0.23c 1.26±1.11c 17.67±5.19b
Ecological indices
Genera richness 17.75±2.06b 22.25±4.27b 23.25±3.86b 35.75±2.36a
Abundance 386.5±228.8a 1063.75±290.2b 844.25±385.0ab 833.25±261.0ab
Enrichment Index 71.10±10.57a 82.21±4.1a 34.30±13.04b 69.96±3.49a
Structure Index 42.86±15.39b 44.81±13.53b 30.20±20.45b 88.48±5.97a
Channel Index 1.11±1.30c 6.67±2.85ab 3.05±3.3bc 9.24±2.52a
Maturity Index 2-5 2.31±0.17c 2.34±0.14c 2.22±0.176c 3.43±0.291b
Shannon-Weaver Index 2.01±0.22b 2.34±0.15b 2.25±0.148b 2.85±0.366a

Table 4.

Percentage of individual nematode trophic groups and average of ecological indices values calculated for Krompachy sampling sites [1].

[i] - abc Means followed by the same letters on the same rows are not statistically different according to Least Significant Difference Test (P<0.05)

3.2. Discussion

Heavy metals can both directly and indirectly affect soil environments, directly by modulating the physiology and behaviour of the soil fauna and flora, or indirectly by altering environmental conditions such as pH or resource availability [27, 28]. The solubility of heavy metals in soil, for which substrate pH is the main driver, determines their toxicity to soil biota [3]. The soil conditions (mostly pH) in the Krompachy area suggest that the solubility of heavy metals could be relatively low, but only the solubilities of As and Cr were low. The solubility and thus the accessibility of the other heavy metals with >20% potential mobility at all sites (K1-K4) could be considered a consequence of chronic contamination of the soil profile and/or a low ability of the soil to bind these elements. The values of the soil parameters and the trend of heavy metal contamination along the transect ruled out the impact of natural geological and/or physicochemical properties. Kovohuty JSC Krompachy was thus the most likely source of the higher heavy metal levels. This conclusion is in agreement with previous data collected from this area [17, 29, 30].

The influence of polymetallic pollution is evident in all aspects of soil ecosystems, including communities of soil nematodes and their trophic structures. Bacterivores and plant feeders were the most resistant to the adverse effects of pollution, consistent with the findings of a previous study in this area [31]. On the other hand, acidic soil combined with heavy metal pollution can positively influence mostly fungi and decrease bacterial densities [32, 33], but our results do not support these effects. We found relatively low proportions of fungivores and high proportions of bacterivores in the acidic soils at most sites. Numerous studies [34, 35, 36, 37] have described differences in the survival of soil nematodes with similar ecological and life strategies under harsh abiotic circumstances, including heavy metal contamination. Fungivores, such as representatives of the genus Aphelenchus tolerant to contamination [38, 39], should thus be more abundant under such environmental contamination. The density of Aphelenchus, however, was significantly higher only at K2, which could indicate that interacting environmental factors (e.g. moisture and texture) might have an important impact on the population density of this trophic group [40].

Indices, trophic & c-p groups As Cd Cr Cu Ni Pb Zn pH Moisture Ninorg. Corg
Bacterial feeders 0.1620.1350.3380.1410.1380.1240.0350.664**0.1-0.0880.699**
Fungal feeders -0.384-0.321-0.215-0.313-0.022-0.19-0.4090.340.802**0.611*0.352
Omnivores -0.753**-0.846**-0.693**-0.857**-0.748**-0.831**-0.796**-0.178-0.0690.475-0.035
Predators -0.469-0.41-0.723**-0.44-0.074-0.28-0.4750.030.522*0.628**0.276
Plant feeders 0.3030.3290.0030.350.2090.2970.444-0.54*-0.529*-0.479-0.600*
c-p 1 0.576*0.629**0.1850.612*0.768**0.662**0.588*0.0120.4790.0740.102
c-p 2 0.3290.3180.612*0.276-0.10.1910.2970.275-0.509*-0.609*0.169
c-p 3 -0.0150.112-0.2320.0910.2260.0790.206-0.494-0.212-0.176-0.727**
c-p 4 -0.615*-0.679**-0.674**-0.671**-0.506*-0.532*-0.656**-0.1430.2240.594*0.116
c-p 5 -0.804**-0.746**-0.854**-0.777**-0.507*-0.737**-0.712**-0.153-0.0920.433-0.221
Genera richness -0.836**-0.736**-0.845**-0.771**-0.409-0.601*-0.765**0.1440.2670.656**0.085
Abundance -0.329-0.182-0.162-0.168-0.029-0.018-0.1760.680**0.4440.529*0.297
Enrichment Index 0.3820.4680.0180.4680.753**0.521*0.4440.0590.662**0.2880.104
Structural Index -0.447-0.474-0.750**-0.459-0.25-0.341-0.453-0.3030.20.553*-0.081
Channel Index -0.564*-0.542*-0.338-0.514*-0.281-0.416-0.612*0.1030.614*0.581*0.118
Maturity Index 2-5 -0.456-0.485-0.741**-0.479-0.259-0.359-0.468-0.2850.2590.597*-0.046
Shannon Weaver Index -0.632**-0.568*-0.688**-0.635**-0.347-0.432-0.638**0.1690.3820.579*0.338

Table 5.

Correlations among heavy metal concentrations, soil characteristics, nematode community structure, and selected ecological indices from Krompachy [1].

[i] - * and ** represent significant level at P<0.05 and P<0.01, respectively

Omnivores and predators are generally considered the most sensitive to any disturbances and stresses [41, 42, 43]. Our data supported the high sensitivity of these two trophic groups and their preference for more eco-friendly conditions over polluted environments. In addition, only these two groups were strongly negatively correlated (P<0.01) with heavy metal pollution. Specifically, omnivores correlated negatively with all heavy metals (P<0.01), and predators strongly correlated negatively (P<0.01) with Cr. Nagy et al. [44] reported similar findings from a field experiment in Hungary, where the lowest concentration of the Cr mobile fraction that produced an observable effect was approximately 0.5 mg kg-1. The presence of the more toxic and mobile Cr6+, which our assay could not detect, may have been responsible for the significant influence of Cr on these two trophic groups, despite the lower Cr concentrations in our study. A synergistic contribution of other heavy metals may also have been responsible for the stress response of these trophic groups.


Figure 3.

Structure and Enrichment Index graphic representation (Krompachy) [1]. A - disturbed food web, bacterial decomposition; B - maturing food web, balanced decomposition; C - structured food web, fungal decomposition; D - degraded food web, fungal decomposition [11].

Community characteristics and ecological indices (H′, MI2-5, and SI) indicated that the soil conditions at our sites were significantly stressed by the contaminants. Generic richness and H′ were lower under higher loads of heavy metals. The lightly contaminated K4 site had the highest diversity of nematode genera, belonging to the entire c-p scale, including opportunists (e.g. Rhabditis, Acrobeloides, and Aphelenchus) and persisters (e.g. Dorylaimus, Oxydirus, and Nygolaimus). The proportion of persisters with higher c-p values rapidly decreased with increasing pollution, which consequently led to lower nematode diversity and a simplification of the food chain. Diversity was mostly negatively correlated with the degree of pollution, but nematode abundance increased at moderately contaminated K2 and K3, but rapidly decreased under a heavy load at K1. Sánchez-Moreno et al. [39] reported similar findings in southern Spain after the Aznalcollar mining spill. A higher generic diversity reflects a higher resilience of ecosystems and complexity of relationships at all trophic levels [45, 46]. A decrease in generic diversity may cause an ecosystem to regress to an earlier successional stage or even to collapse. On the other hand, an increase in total density under a low or moderate degree of pollution does not necessarily indicate an uncontaminated environment but rather better survival conditions, as in our study (the growth in the population of opportunistic nematodes). Contamination can in some cases act as a direct or indirect stimulus to a particular component of the ecosystem and lead to inaccurate assessments of the actual environmental state [4]. Nevertheless, the generic composition of the nematode communities and several of the ecological indices (EI, SI, and MI2-5) suggest that a certain amount of regression is occurring in the study area closer to the pollution source.

Georgieva et al. [42] have studied the combined effect of heavy metals. The maturity of an ecosystem, measured by MI2-5, was negatively affected more under significantly lower polymetallic contamination than if the heavy metals were applied separately. A nematode community was highly insensitive to Cd, showing no significant changes, at soil concentrations up to 160 mg kg-1 [47]. Similar insensitivities in nematode communities were also described for increased levels of As [48]. The toxicity of some elements, though, can increase in combination with other heavy metals and may have negative impacts on nematodes [42]. This synergistic effect could thus be responsible for the decrease in MI2-5 and SI towards the pollution source in our study, despite their nonsignificant correlation with the heavy metals.

The analysis of SI and EI at Krompachy indicated significant differences amongst the sites in structural complexity and resource availability of the ecosystem. The graphic representation of the SI and EI parameters (Figure 3) mapped K4 in quadrat B (maturing ecosystem with low or moderate disturbance), K1 and K2 in quadrat A (ecosystem under high disturbance and with a disturbed food web), and K3 in quadrat D (stressed ecosystem with a degraded food web). These results indicated differences in the composition of the structural component (c-p groups 3-5) rather than the enrichment component (c-p 1) [49]. The shift in the ecosystem in our study under heavy metal pollution corresponded well with the conditions reported by [50, 51, 52].

The decomposition of organic matter can also have a large impact on the development of a soil ecosystem and its maturity. The rate of breakdown of organic matter depends on soil conditions and the participation of a variety of decomposers [11]. We studied decomposition using CI and found that the breakdown of organic matter followed distinct degradation pathways at the different sites. EI indicated high levels of energy resources at all sites except K3, but CI indicated that the quality of these resources (C:N ratio) differed: K1 and K3 had mainly bacterial pathways of degradation, dominant in N-enriched conditions, but the other sites had a slower breakdown, with an increasing importance of fungi breaking down organic matter with a higher C:N ratio. The change from fast to slow decomposition was also apparent in the nematode community structure and the substitution of bacterial enrichment opportunists (mainly Rhabditis) by more efficient general opportunists (Cephalobus and Acrobeloides).

4. Indirect effects of industrial emissions – A case study of SMZ JSC Jelšava

4.1. Results

4.1.1. Heavy metals

The total concentrations of heavy metals in the soil were relatively low and did not exceed the limits for uncontaminated soils (Table 6) set by [26]. Only the concentration of Mg increased, where the total and mobilisable concentrations were more than ten-fold higher at J1 (568.2 and 423.9 mg kg-1, respectively) than at J4 (46.03 and 16.4 mg kg-1, respectively). The increasing Mg gradient, but also for other heavy metals (Cr, Cu, Ni, and Zn) towards the pollution source (J1) is shown in Figure 4.

4.1.2. Composition of the nematode fauna: trophic and c-p groups

A total of 52 nematode genera were identified (Table 7). The abundance of most genera, e.g. Helicotylenchus, Aphelenchoides, Axonchium, and Oxydirus, increased with distance from the pollution source, but some genera, e.g. Acrobeloides and Rhabditis, dominated only at the most polluted site, J1. The generic richness was approximately two-fold higher at J2-J4 than at J1 (Table 8) and, as with abundance, was significantly (P<0.01) limited by pH (Table 9).

Trace elements Sampling site
J1 J2 J3 J4 Limit
Total concentration (HNO3 extraction)
Arsenic (As) 0.01±0.07a 0.06±0.04a 0.03±0.02a 0.07±0.05a 5
Cadmium (Cd) 0.03±0.01a 0.031±0.01a 0.03±0.01a 0.02±0.01a 0.3
Chrome (Cr) 0.19±0.05b 0.04±0.01a 0.03±0.004a 0.12±0.15ab 10
Copper (Cu) 7.15±1.34a 3.73±1.17a 2.73±0.29a 1.85±0.27a 20
Magnesium (Mg) 568.18±98.13c 68.25±8.06a 105.40±2.83ab 46.03±3.46a n.a
Nickel (Ni) 0.95±0.36b 0.35±0.02a 0.66±0.08ab 0.50±0.22a 10
Lead (Pb) 2.34±1.06ab 5.78±1.15c 1.56±0.06a 1.73±0.20a 30
Zinc (Zn) 0.37±0.10d 0.14±0.03ab 0.22±0.03bc 0.12±0.06a 40
Mobilisable fraction (Na2EDTA extraction)
Arsenic (As) 0.02±0.01a 0.05±0.03a 0.02±0.02a 0.02±0.002a n.a.
Cadmium (Cd) 0.014±0.005bc 0.023±0.002d 0.008±0.001a 0.013±0.002b n.a.
Chrome (Cr) 0.10±0.03c 0.02±0.004a 0.01±0.01a 0.02±0.004ab n.a.
Copper (Cu) 5.13±2.27c 2.32±1.29ab 1.02±0.08a 1.64±0.26a n.a.
Magnesium (Mg) 423.88±183.43c 60.13±6.12ab 47.75±14.50a 16.43±8.58a n.a.
Nickel (Ni) 0.52±0.09d 0.30±0.02b 0.20±0.02a 0.34±0.02bc n.a.
Lead (Pb) 0.48±0.14a 0.91±0.04b 0.75±0.11ab 0.72±0.42ab n.a.
Zinc (Zn) 0.21±0.04d 0.12±0.03bc 0.06±0.01a 0.09±0.01ab n.a.

Table 6.

Total and mobilisable concentration of trace elements in sampling sites from Jelšava ( [53].

[i] - abcd Means followed by the same letters on the same rows are not statistically different according to Least Significant Difference Test (P<0.05)

[ii] - Limit - limits posted by The Decree of the Ministry of Land Management of the Slovak Republic No. 531/1994-540 on the admissible values of harmful substances in uncontaminated soil. n.a. – not available


Figure 4.

Redundancy analysis (RDA) performed on physicochemical soil properties, trace elements and selected indices in relation to sampling sites with data explained 64.0 % of the variation in the first two axes; F=7.1; P=0.002 (Jelšava) [53].

The soil pH and concentrations of Mg, Cr, Cu, Ni, and Zn strongly (P<0.01) influenced all trophic groups, predators and omnivores in particular (Tables 8 and 9). Bacterivores, with more than 90% proportion in community at J1, were the only trophic group with a positive correlation (P<0.01) with pH, Cu, Mg, and Zn (Table 9). The proportion of plant feeders, the second most abundant trophic group, was relatively low (8.1%) near the pollution source, while fungivores reached a proportion of only 0.41%, omnivores 0.22%, and predators were absent in community (Table 8).

The distribution of the c-p groups amongst the sites supported their mapping according to r/K characteristics. The group most tolerant to pollution, c-p 1, dominated under strong contamination (Table 8) and was able to survive even at the higher concentrations of Mg, Zn, and Cu (Table 9). The abundance of c-p 3 was positively correlated only with Cr and Zn (Table 9) and was highest at J1 (Figure 4), mostly due to a high frequency of bacterivores (Prismatolaimus and Teratocephalus) and plant feeders (Pratylenchus) (Table 7). In contrast, c-p 4 and 5 were sensitive to Cu, Mg, and Zn contamination and alkaline soil (Table 9), and occurred predominantly at the sites farthest from the pollution source. Representatives of c-p 2 were positively influenced mainly by Corg.

4.1.3. Ecological indices

The ecological indices indicated degradation of the soil ecosystem at J1, most likely due to the exceptionally high alkalinity and Mg content (Figure 4). J3 had the highest diversity, as shown by H′, which decreased towards J1 (Table 8). MI2-5 and SI supported this trend, with significant changes (P<0.01) in values between J1 and J2 (Table 8).

The amount of available food resources in the soil system identified by EI was similar at all sites, but J3 had slightly less favourable conditions (Table 8). A graphic representation of EI and SI suggested a higher level of ecological disruption at J1 relative to the other sites (Figure 5). The trend in CI values reflected the increasing importance of fungi in the decomposition of organic matter with distance from the pollution source. The high concentration of Mg; higher levels of Cu, Zn, and Corg; and soil alkalinity were the most important factors influencing the structure and maturity of the nematode communities (Table 9, Figure 4).

Nematode genera J1 J2 J3 J4
c-p Abundance Abundance Abundance Abundance
Bacterial feeders
Acrobeles 2-1.00±2.00--
Acrobeloides 239.75±9.9124.50±23.1627.00±13.1411.50±7.85
Alaimus 40.25±0.502.75±2.503.00±3.566.75±4.99
Anaplectus 2---2.75±3.59
Aulolaimus 3--0.50±1.002.00±4.00
Cephalobus 24.00±3.8321.00±14.4921.75±14.7311.50±3.11
Cervidellus 26.50±4.51---
Diplogaster 1-1.25±2.502.25±2.87-
Eucephalobus 28.75±4.118.25±5.508.25±9.541.25±1.89
Heterocephalobus 23.25±3.77--0.50±1.00
Chiloplacus 24.00±5.66---
Mesorhabditis 111.50±7.724.75±9.50--
Monhystera 2-1.00±2.001.00±0.821.50±1.73
Panagrolaimus 15.00±5.48---
Plectus 20.50±0.585.00±4.089.50±5.807.25±1.89
Prismatolaimus 35.25±3.406.25±7.59--
Rhabditis 1176±76.6286±132.469.75±58.01144.3±114.4
Teratocephalus 311.75±6.55---
Fungal feeders
Aphelenchus 21.25±1.899.00±3.5619±8.0419.75±9.67
Diphtherophora 3-1.50±1.732±1.830.5±0.58
Dorylaimoides 4-0.5±1--
Filenchus 2-2.75±5.5-3.75±3.3
Aporcelaimellus 50.5±118±5.8911.25±9.917.50±3.11
Dorylaimus 4-47±36.5218.5±26.837±12.03
Epidorylaimus 4-11.25±8.060.5±0.580.5±0.58
Eudorylaimus 40.25±0.52.25±2.634±1.833.5±1.29
Mesodorylaimus 5-1.00±0.82-11.25±4.92
Prodorylaimus 5-6.75±6.89.75±9.8412.75±11.67
Anatonchus 4-0.25±0.506.25±2.752±1.41
Clarkus 4-0.5±1-15.50±9.61
Discolaimus 5--1.25±1.26-
Miconchus 4-2.00±2.312.00±2.160.25±0.50
Mononchus 4-0.50±0.580.50±1.00-
Mylonchulus 4-4.50±2.897.25±4.113.00±4.24
Nygolaimus 5-0.25±0.50--
Oxydirus 5-16.50±22.2217.25±9.1893.75±88.28
Prionchulus 4-0.25±0.50--
Tripyla 3-12.25±15.520.75±1.50-
Plant feeders
Aglenchus 21.00±2.0012.00±14.240.25±0.500.50±1.00
Aphelenchoides 27.00±6.7822.50±11.7038.50±40.346.50±1.73
Axonchium 5-22.75±29.0311.50±8.70210±64.10
Boleodorus 2-10.75±13.792.75±4.197.00±5.60
Coslenchus 2---0.50±1.00
Ditylenchus 20.25±0.50---
Geocenamus 30.25±0.500.50±1.00--
Helicotylenchus 31.00±1.4188.25±109.230.75±32.9464.00±35.80
Malenchus 2-14.75±16.76-3.75±5.68
Paratylenchus 20.25±0.5025.00±38.4311.25±13.403.75±1.71
Pratylenchus 315.50±17.906.75±5.680.25±0.502.75±4.19
Rotylenchus 3-21.50±32.1313.25±16.6824.25±13.87
Tylenchus 2-16.50±14.067.00±6.4818.75±22.97
Xiphinema 5-24.00±28.32--

Table 7.

C-p values and the average abundance (±SD) of nematode genera in individual sampling sites from Jelšava [53].

4.2. Discussion

In this case study, the long-term industrial emission of magnesite has affected the soil ecosystem, mostly by changing the basic soil conditions. Such high loads of Mg in the soils as recorded in our study are rare, so the impact of Mg emissions on soil fauna is not often studied and is not yet understood [53]. Similar steep increases in soil alkalinity as a direct consequence of industrial inputs of Mg as was found in our assay described also Machín and Navas [54]. Most studies have investigated the importance of Mg deficiency in higher plants as a trace element important in photosynthesis. Hronec [17] described a decrease in soil fertility and the replacement of native plant species by tolerant halophytes in an area with high Mg input. Moreover, the long-term fallout of Mg dust (>600 t km-2 y-1) may degrade natural ecosystems and create corrosive crusts that could affect natural interactions in the soil even more [17].

The absence of sensitive nematode groups, likely caused by an ineffective regulation of osmotic pressure by their cuticles due to changes in the concentrations of various ions in the soil pore water, can be a direct effect of altered soil conditions [55]. The concentrations of ions in soil solutions under extreme changes in pH may negatively affect the regulation of water by organisms [56], thereby disrupting homeostasis. A change in the concentrations of soil ions may thus have been an important factor responsible for the decrease in the higher c-p groups (mainly predators and omnivores with highly permeable cuticles) at the sites with high pH and Mg concentration. The indirect influence of contamination was likely a restriction in the availability of resources and the interactions necessary for survival or an effect on the abiotic attributes of the environment [32].

Indices, trophic & c-p groups Sampling sites
J1 J2 J3 J4
Trophic & c-p groups
Bacterial feeders 91.31±6.28a 50.47±12.29bc 38.98±7.86cd 25.20±13.23d
Fungal feeders 0.41±0.52a 3.02±3.49ab 6.16±2.38b 3.19±0.77ab
Omnivores 0.22±0.26a 11.32±3.05a 11.24±9.64a 10.09±1.87a
Predators 0±0a 4.38±3.22ab 12.67±9.17ab 14.94±8.21b
Plant feeders 8.06±6.18a 30.80±14.86ab 30.95±14.16ab 46.59±5.97b
c-p 1 67.96±5.10a 58.85±7.00ab 25.66±16.27c 33.42±20.39bc
c-p 2 25.57±5.57a 15.14±4.42a 38.96±19.95a 15.35±3.55a
c-p 3 6.11±1.54a 3.81±3.81a 1.04±0.90a 0.94±1.75a
c-p 4 0.22±0.44a 14.06±4.71b 16.58±7.91b 17.64±3.58b
c-p 5 0.13±0.25a 8.14±3.10a 17.75±8.25ab 32.65±19.04b
Ecological indices
Genera richness 14.50±0.58a 28.75±2.22b 25.25±3.50b 27.50±1.00b
Abundance 303.75±104.6a 764±376.06a 358.75±169.6a 742.00±224.8a
Enrichment Index 91.33±2.30a 93.94±2.01a 68.96±27.63a 87.50±8.31a
Structure Index 38.12±7.22a 88.25±4.03b 80.73±11.18b 93.24±4.87b
Channel Index 0.16±0.20a 1.47±1.70a 20.42±32.01a 5.35±3.10a
Maturity Index 2-5 2.22±0.050a 3.39±0.25b 3.21±0.47b 3.93±0.44b
Shannon-Weaver Index 1.56±0.08a 2.33±0.24b 2.57±0.13b 2.29±0.11b

Table 8.

Percentage of individual nematode trophic groups and average of ecological indices values calculated for Jelšava sampling sites [53].

[i] - abcd Means followed by the same letters on the same rows are not statistically different according to Least Significant Difference Test (P<0.01).

In areas with neutral soil pHs (as at J4), strong shifts in soil composition are usually expected after an increase in pH. Alterations in edaphic diversity after pH manipulation have been confirmed experimentally and by field studies. Earthworms, enchytraeids, and nematodes were more abundant after liming [57, 58], but Acari preferred acidic conditions [59, 60]. The altered soil conditions in our study as a consequence of Mg contamination and an extremely high pH have led to direct and indirect effects on the nematode communities and most likely influenced the composition of the entire soil fauna. The plant feeders showed the first signs of the indirect effects of high Mg levels. A fluctuation in plant feeder density is generally considered a sign of changes in primary production [61]. Excessive concentrations of Mg under alkaline conditions induce a deficiency of essential nutrients (e.g. Ca and P) and impede the development of plant root systems [62]. Relatively impoverished phytocoenoses composed mainly of halophytes such as Puccinellia distans or Agrostis stolonifera are usually found near pollution sources [63]. A shift towards bacterivorous nematodes in the community was another indication of an indirect influence of pollution in our study, because bacteria prosper more than fungi at alkaline pHs [63].

Indices, trophic & c-p groups As Cd Cr Cu Mg Ni Pb Zn pH Moisture Ninorg. Corg
Bacterial feeders 0.1160.1680.3490.779**0.888**0.471-0.1620.642**0.750**0.253-0.4380.433
Fungal feeders -0.266-0.479-0.653**-0.467-0.472-0.696**0.24-0.794**-0.361-0.3960.542*-0.052
Omnivores 0.1080.303-0.513*-0.306-0.437-0.3810.620*-0.462-0.711**-0.0320.414-0.432
Predators -0.056-0.301-0.535*-0.641**-0.839**-0.545*0.225-0.681**-0.603*-0.3840.446-0.539*
Plant feeders 0.0190.078-0.259-0.726**-0.862**-0.4060.2-0.499-0.789**-0.0740.465-0.471
c-p 1 0.2210.410.4910.788**0.738**0.606*-0.1260.793**0.555*0.482-0.4030.188
c-p 2 -0.361-0.576*0.013-0.0470.406-0.091-0.279-0.1410.550*-0.321-0.2210.727**
c-p 3 0.1180.3510.653**0.3190.619*0.481-0.2740.700**0.526*0.425-0.3220.329
c-p 4 -0.16-0.126-0.495-0.596*-0.681**-0.4570.265-0.726**-0.735**-0.2630.375-0.271
c-p 5 -0.081-0.244-0.526*-0.646**-0.791**-0.563*0.257-0.651**-0.617*-0.330.419-0.526*
Genera richness 0.1750.406-0.500*-0.543*-0.527*-0.4910.618*-0.396-0.725**0.0180.359-0.632**
Abundance 0.0430.394-0.247-0.313-0.475-0.130.215-0.18-0.715**0.171-0.012-0.613*
Enrichment Index 0.4580.661**0.2280.4760.2880.3060.1940.553*0.0550.597*0.012-0.303
Structural Index 0.080.121-0.422-0.429-0.732**-0.3380.353-0.421-0.753**-0.1210.365-0.753**
Channel Index -0.241-0.445-0.501*-0.581*-0.720**-0.593*0.104-0.753**-0.478-0.3130.583*-0.148
Maturity Index 2-5 0.10.105-0.451-0.415-0.729**-0.350.359-0.439-0.744**-0.1440.359-0.736**
Shannon Weaver Index -0.162-0.292-0.679**-0.829**-0.456-0.779**0.362-0.777**-0.426-0.3790.421-0.295

Table 9.

Correlations among heavy metal concentrations, soil characteristics, nematode community structure, and selected ecological indices from Jelšava [53].

[i] - * and ** represent significant level at P<0.05 and P<0.01, respectively.

An abundance of food for bacteria often increases the proportion of c-p 1 bacterivores, as confirmed experimentally by [58]. Representatives of the genus Rhabditis were the most abundant nematodes along the entire transect in our study. The survival and dominance of Rhabditis and other bacterivores in this hostile environment could, however, be a complex phenomenon, comprising not only food resources, but also behavioural and physiological adaptations to the hostile conditions [64]. For example, C. elegans can withstand high pH and salinity and survive under extreme osmotic pressure [55, 65]. Moreover, c-p 1 nematodes can produce inactive dauer larvae, which increase the chances of withstanding harmful or limited conditions. Other studies [1, 4] have shown that c-p 2 nematodes are also able to withstand severe ecological conditions and comprise the dominant c-p group in highly disturbed ecosystems [66]. Our results only partially supported these findings. More than 25% of the nematodes at J1 belonged to the c-p 2 group, and of these, species of Cervidellus and Chiloplacus were the most tolerant. The proportion of c-p 2 nematodes within the community, however, was significantly lower than that of c-p 1 nematodes. The discrepancy between the hypothetical dominance of stress-tolerant c-p 2 nematodes and our results was likely due to various physiological and ecological aspects of c-p 1 and 2 nematodes and their population dynamics (e.g. reproduction rate, length of life cycle, and resource availability). The abundance of c-p 1 nematodes in the soil may also be due to, aside from soil fertility reported by [11], the changes in soil pH. Further study is necessary to confirm these possibilities.


Figure 5.

Structure and Enrichment Index graphic representation (Jelšava) [53]. A - disturbed food web, bacterial decomposition; B - maturing food web, balanced decomposition; C - structured food web, fungal decomposition; D - degraded food web, fungal decomposition [11].

This study has demonstrated the negative ecological impact of contamination by quantifying different ecological indices (e.g. SI, MI2-5, and H′) and measuring key parameters of nematode communities (generic richness and abundance). Not all communities traits showed a clear response to stress in the ecosystem, and some responded inconsistently, e.g. abundance similar as in the Krompachy case study did not reflect the stress of the ecosystem. Total abundance should thus not be used as an indicator for either heavy metals, as suggested by [42], or physical or chemical changes in soil properties such as pH. Generic richness and H′, though, indicated substantial differences in species diversity in the ecosystems with diverse contamination and soil salinity, in accordance with diversity shifts observed in other nematode communities [58] and in other edaphic groups after pH manipulation [59, 67]. SI and MI2-5 were the best tools in our study for showing the impact of high Mg concentrations and soil salinity on nematode community structure and the maturity of the ecosystem. SI combined with EI identified J1 as severely disturbed (quadrat A), and the other sites showed only a low degree of stress. A decrease in the intensity of stress and an increase in the complexity of soil community structure with distance from a pollution source have been reported in several other areas of intensive industry [1, 44]. Based on low CI values in our study, indicating predominant bacterial pathway of decomposition, c-p 1 enrichment opportunists were eliminated from maturity assessment (MI2-5), as they respond more to the food availability than pollution [47]. Eliminating the influence of enriched conditions provided us with clear gradual increase in environment stability and maturity from the site J1 onwards.

5. Conclusions

Heavy metal contamination of soil is long-term ecological problem. High levels of pollutants may lead to losses in production, breaches in the otherwise strict sanitary limits for agricultural products, and health-related problems for local inhabitants. This study demonstrated that the background geochemical heavy metal content was significantly exceeded, with measurable consequences for the local environments. The sensitive groups of soil nematodes, mainly omnivorous and predacious nematodes of the higher c-p groups, were significantly suppressed in the ecosystems, and tolerant groups, with better physiological and behavioural adaptations to harsh conditions, survived. Such alterations in community structure may have far-reaching consequences, because some of the ecological regulation or environmental feedbacks provided by nematodes may be suppressed in, or eliminated from, the ecosystem.

Despite the recent decrease in emissions and the application of various treatments to prevent their leakage at both metallurgical plants, the potential threat of further environmental deterioration remains relatively high. The high concentrations of possible toxic elements bound to the soil matrix and temporally unavailable to the soil biota are the main threat in the Krompachy area. The release of the pollutants to the soil solution by sudden shifts in local soil conditions, e.g. an increase in soil acidity by acid rain, could lead to their increased toxicity to soil communities or to groundwater contamination. The threat from pollution to the ecosystem around Jelšava is less imminent, because the levels of toxic elements were not significantly higher. The high Mg concentrations and soil salinity, and consequently the extreme soil alkalinity, however, are the most important problems in this area. Under these conditions, only a few organisms with adaptations for surviving in such an environment are able to prosper, as our results clearly showed. In both case studies, the nematode community structures collapsed near the pollution sources, with the domination of nematodes able to survive inhospitable soil conditions.


The authors acknowledge the support of the Slovak Research Development Agency (project LPP-0085-09) (0.6), VEGA (Grant No. 2/0193/14) (0.3), and the project “Application Centre to protect humans, animals and plants against parasites” (Code ITMS: 26220220018) with the support of the Operational Program “Research and Development” funded by the European Regional Development Fund (0.1).


1 - Šalamún P., Renčo M., Kucanová E., Brázová T., Papajová I., Miklisová D., Hanzelová V. Nematodes as bioindicators of soil degradation due to heavy metals. Ecotoxicology 2012; 21 2319–2330.
2 - Háněl L. Colonization of chemical factory wastes by soil nematodes. Pedobiologia 2004; 48 373–381.
3 - Kabata-Pendias A. Trace Elements in Soils and Plants. CRC Press, Taylor and Francis Group, Boca Raton, USA; 2011.
4 - Wilson MJ., Kakouli–Duarte T. Nematodes as Environmental Indicators. Wallingford, UK: CAB International; 2009.
5 - Neher DA. Role of nematodes in soil health and their use as indicators. Journal of Nematology 2001; 33 161–168.
6 - Nannipieri P., Grego S., Ceccanti B. Ecological significance of the biological activity in soil. In: Bollag JM., Stozky G. (eds.) Soil Biochemistry. Marcel Dekker, New York; 1990. p293–355.
7 - Bongers T. The maturity index: an ecological measure of environmental disturbance based on nematode species composition. Oecologia 1990; 83 14–19.
8 - Frampton GK. The potential of Collembola as indicators of pesticide usage: evidence and methods from the UK arable ecosystem. Pedobiologia 1997; 41 179–184.
9 - Ruf A. A maturity index for predatory soil mites (Mesostigmata: Gamasina) as an indicator of environmental impacts of pollution on forest soil. Appl. Soil Ecol. 1998; 9 447–452.
10 - Gupta VVSR., Yeates GW. Soil microfauna as indicators of soil health. In: Pankhurst CE., Doube BM., Gupta VVSR. (eds.) Biological Indicators of Soil Health. CAB International; 1997. p201–233.
11 - Ferris H., Bongers T., De Goede RGM. A framework for soil food web diagnostics: extension of the nematode faunal analysis concept. Applied Soil Ecology 2001; 18 13–29.
12 - Verhoef HA., Brussaard L. Decomposition and nitrogen mineralization in natural and agro-ecosystems: The contribution of soil animals. Biogeochemistry 1990; 11 175–211.
13 - Brussaard L. Soil fauna, guilds, functional groups and ecosystem processes. Applied Soil Ecology 1998; 9 123–35.
14 - De Deyn GB., Raaljmakers CE., Zoomer HR., Berg MP., De Rulter PC., Verhoef HA., Bezemer TM., Van der Putten WH. Soil invertebrate fauna enhances grassland succession and diversity. Nature 2003; 422 711–713.
15 - Van Straalen NM. The use of soil invertebrates in ecological survey of contaminated soils. In: Doelman P., Eijsackers HJP. (eds.) Vital Soil: Function, Value and Properties. Elsevier, Amsterd; 2004. p159–194.
16 - Fargašová A. Metal contamination in Slovakia. 2009. (accessed 15 May 2011).
17 - Hronec O. Exhaláty – Pôda – Vegetácia. Prešov: TOP s.r.o. Bratislava: Slovenská poľnohospodárska a potravinárska komora; 1996.
18 - MP SR, 2005. Vyhláška Ministerstva pôdohospodárstva Slovenskej Republiky č. 338/2005 Z.z. (accessed 13 September 2014).
19 - Sabienë N., Brazauskienë DM., Rimmer D. Determination of heavy metal mobile forms by different extraction methods. Ekologija 2004; 1 36–41.
20 - Southey JF. Laboratory methods for work with plant and soil nematodes. Ministry of Agriculture, Fisheries and Food, Reference Book 402, London, H.M. Stationery Off.; 1986.
21 - Van Bezooijen J. Methods and techniques for nematology. 2006. p112 (accessed 4 May 2012).
22 - Yeates GW., Bongers T., De Goede RGM., Freckman D.W., Georgieva S.S. Feeding habits in nematode families and genera - an outline for soil ecologists. J. Nematol 1993; 25 315–331.
23 - Bongers T., Korthals GW. The Maturity Index, an instrument to monitor changes in the nematode community structure. 1993 conference proceedings, Summaries of the 45th International Symposium on Crop Protection. Ghent, Belgium, p. 80.
24 - Shannon CE., Weaver W. The mathematical theory of communication. University of Illinois, Urbana, Illinois, USA; 1949.
25 - Ter Braak C.J.F., Šmilauer P. CANOCO reference manual and user's guide: software for ordination, version 5.0. Microcomputer Power, Ithaca, USA: Biometris; 2012.
26 - MP SR 1994. Správa o stave životného prostredia Slovenskej republiky: Pôda. Rozhodnutie MP SR o najvyšších prípustných hodnotách škodlivých látok v pôde a o určení organizácií oprávnených zisťovať skutočné hodnoty týchto látok č. 531/1994 – 540. (accessed 20 May 2011).
27 - Impellitteri CA., Allen HE., Yin Y., You SJ., Saxe JK. Soil properties controlling metal partitioning. In: Selim HM., Sparks DL. (eds.) Heavy Metal Release in Soils. Lewis Publ., Boca Raton, FL; 2001. p149–165.
28 - Carrillo-Gonzáles R., Šimunek J., Sauvé S., Adriano D. Mechanisms and pathways of trace elements mobility in soils. Advances in Agronomy 2006; 91 112–178.
29 - Banášová V., Lackovičová A. The disturbance of grasslands in the vicinity of copper plant in the town of Krompachy (Slovenské rudohorie Mts., NE Slovakia). Bull. Slov. Bot. Spoločn. 2004; 26 153–161.
30 - Takáč P., Kozáková Ľ., Vaľková M., Zeleňák F. Ťažké kovy v pôdach stredného Spiša. Acta Montanistica Slovaca 2008; 13 82–86.
31 - Valocká B., Sabová M. Structure of soil nematode communities in some regions of Slovakia with different pollution levels. Biológia 1997; 52 609–613.
32 - Korthals GW., Alexiev AD., Lexmond TM., Kammenga JE., Bongers T. Long–term effects of copper and pH on the nematode community in an agroecosystem. Environmental Toxicology and Chemistry 1996a; 15 979–985.
33 - Bakonyi G., Nagy P., Kádár I. Long–term effects of heavy metals and microelements on nematode assemblage. Toxicology Letters 2003; 140–141 391–401.
34 - Korthals GW., Popovici J., Van Megen HHB. Soil nematodes in heathland around a zinc-smelter near Budel, The Netherlands. In: De Goede RGM., Bongers T. (eds), Nematode Communities of Northern Temperate Grassland Ecosystems. Giessen, Germany, Focus Verlag; 1998. p155–158.
35 - Zhang WD., Wang XF., Li Q., Jiang Y., Liang WJ. Soil nematode responses to heavy metal stress. Helminthologia 2007; 44 87–91.
36 - Pen-Mouratov S., Shururov N., Steinberger Y. Influence of industrial heavy metal pollution on soil free-living nematode population. Environ. Pollut. 2008; 152 172–183.
37 - Chen G., Qin J., Shi D., Zhang Y., Ji W. Diversity of soil nematodes in area polluted with heavy metals and polycyclic aromatic hydrocarbons (PAHs) in Lanzhou, China. Environmental Management 2009; 44 163–172.
38 - Li F., Neher DA., Darby BJ., Weicht TR. Observed differences in life history characteristics of nematodes Aphelenchus and Acrobeloides upon exposure to copper and benzo(a)pyrene. Ecotoxicology 2005; 14 419–429.
39 - Sánchez-Moreno S., Camargo JA., Navas A. Ecotoxicological assessment of the impact of residual heavy metals on soil nematodes in the Guadiamar river basin (Southern Spain). Environmental Monitoring and Assessment 2006; 116 245–262.
40 - Ettema CHH. Soil nematode diversity: Species coexistence and ecosystem function. Journal of Nematology 1998; 30 159–169.
41 - Weiss B., Larink O. Influence of sewage sludge and heavy metals on nematodes in arable soil. Biology and Fertility of Soils 1991; 12 5–9.
42 - Georgieva SS., Mcgrath SP., Hooper DJ., Chambers BS. Nematode communities under stress: the long–term effects of heavy metals in soil treated with sewage sludge. Applied Soil Ecology 2002; 20 27–42.
43 - Pen-Mouratov S., Shukurov N., Steinberger Y. Soil free-living nematodes as indicators of both industrial pollution and livestock activity in Central Asia. Ecological Indicators 2010; 10 955–967.
44 - Nagy P., Bakónyi G., Bongers T., Kádár I., Fábián M., Kiss I. Effects of microelements on soil nematode assemblages seven years after contaminating an agricultural field. Science of the Total Environment 2004; 320 131–143.
45 - Odum PE. Trends expected in stressed ecosystems. Bioscience 1985; 35 419–422.
46 - Peterson G., Allen CR., Holling CS. Ecological resilience, biodiversity, and scale. Ecosystems 1998; 1 16–18.
47 - Korthals GW., Van de Ende A., Van Megen H., Lexmond TM., Kammenga JE., Bongers T. Short–term effects of cadmium, copper, nickel and zinc on soil nematodes from different feeding and life–history strategy groups. Applied Soil Ecology 1996b; 4 107–117.
48 - Nagy P. Effects of an artificial metal pollution on nematode assemblage of a calcareous loamy chernozem soil. Plant and Soil 1999; 212 35–43.
49 - Park BY., Lee JK., Ro HM., Kim YH. Effects of heavy metal contamination from an abandoned mine on nematode community structure as an indicator of soil ecosystem health. Applied Soil Ecology 2011; 51 17–24.
50 - Sánchez–Moreno S., Navas A. (2007) Nematode diversity and food web condition in heavy metal pollutes soils in river basin in southern Spain. European Journal of Soil Biology 2007; 43 166–179.
51 - Shao Y., Zhang W., Shen J., Zhou L., Xia H., Shu W., Ferris H., Fu S. Nematodes as indicators of soil recovery in tailings of a lead/zinc mine. Soil Biology and Biochemistry 2008; 40 2040–2046.
52 - Li Q., Zhong S., Li F., Lou Y., Liang W. Nematode community structure as bioindicator of soil heavy metal pollution along an urban-rural gradient in southern Shenyang, China. International Journal of Environment and Pollution 2011; 45 297–309.
53 - Šalamún P., Kucanová E., Brázová T., Miklisová D., Renčo M., Hanzelová V. Diversity and food web structure of nematode communities under high soil salinity and alkaline pH. Ecotoxicology 2014; 23 1367–1376.
54 - Machín J., Navas A. Soil pH changes induced by contamination by magnesium oxides dust. Land Degradation and Development 1999; 11 37–50.
55 - Van Gundy SD. Factors in survival of nematodes. Annual Review of Phytopathology 1965; 3 43–68.
56 - Bååth E., Berg B., Lohm U., Lundgren B., Lundkvist H., Rosswall T. Effects of experimental acidification and liming on soil organisms and decomposition in a Scots pine forest. Pedobiologia 1980; 20 85–100.
57 - Potthoff M., Asche N. Stein B. Muhs A. Beese F. Earthworm communities in temperate beech wood forest soils affected by liming. European Journal of Soil Biology 2008; 44 247–254.
58 - Räty M., Huhta V. Earthworms and pH affect communities of nematodes and enchytraeids in forest soil. Biology and Fertility of Soils 2003; 38 52–58.
59 - Hågvar S., Abrahamsen G. Colonisation by Enchytraeidae, Collembola and Acari in sterile soil samples with adjusted pH levels. Oikos 1980; 34 245–258.
60 - Van Straalen NM., Verhoef HA. The development of a bioindicator system for soil acidity based on Arthropod pH preferences. Journal of Applied Ecology 1997; 34 217–232.
61 - Yeates GW. How plants affect nematodes. Advances in Ecological Research 1987; 17 61–113.
62 - Li YZ. et al. The effects of excessive magnesium in soil on physiological indicators of 19 soybean. In Agro-Environmental Protection 1990; 9 41–42.
63 - Kautz G., Zimmer M., Zach P., Kulfan J., Topp W. Suppression of soil microorganisms by emissions of a magnesite plant in the Slovak republic. Water, Air and Soil Pollution 2001; 125 121–132.
64 - McSorley R. Adaptations of nematodes to environmental extremes. Florida Entomologist 2003; 86 138–142.
65 - Khanna N., Cressman CP., Tatara CP., Williams PL. Tolerance of the nematode Caenorhabditis elegans to pH, salinity, and hardness in aquatic media. Archives of Environmental Contamination and Toxicology 1997; 32 110–114.
66 - Bongers T. The maturity index, the evolution of nematode life history traits, adaptive radiation and cp-scaling. Plant and Soil 1999; 212 13–22.
67 - Hågvar S. Reactions to soil acidification in microarthropods: Is competition a key factor? Biology and Fertility of Soils 1990; 9 178–181.