Comprehensive Study of Reclaimed Soil, Plant, and Water Chemistry Relationships in Highly S-Contaminated Post Sulfur Mine Site Jezi ó rko (Southern Poland)

: The aim of the work was a comprehensive study of the soils (pH, EC, SOC, N T , S T ), surface waters (pH, EC, Ca 2+ Mg 2+ , Na + , NO 3 − , SO 42 − , Cl − , HCO 3 − ), and reactions of trees and herbaceous plants in the restored forest ecosystem of a former sulfur mine. Common birch and Scots pine growth reaction, vitality (according to IUFRO standards- International Union of Forest Research Organizations), nutrient supply (Na, K, P, Ca, Mg, K), and Calamagrostis epigejos (L.) Roth chemical composition (Na, K, P, Ca, Mg, K) were assayed. The chemistry dynamics (pH, EC, DOC, N T , Ca, Mg, and S at the beginning and end of the experiment) of soil leaching and the sulfur load leached from the sulfur-contaminated soil substrates were evaluated. The remediation effects of birch and pine litter were assayed in an experiment under controlled conditions. It was found that reclamation was effective in the majority of the post-mining site; however, hotspots with sulfur contamination reaching even 45,000 mg kg − 1 , pH < 2.0 and electrical conductivity (EC) of 6500 µ S cm − 1 were reported. Surface waters typically displayed elevated concentrations of sulfate ions (average 935.13 mg L − 1 ), calcium ions (up to 434 mg L − 1 ), and high EC (average 1797 µ S cm − 1 ), which was related both to sulfur contamination and the sludge lime that was used in neutralization. Calamagrostis epigejos was found to be a species that adapted well to the conditions of elevated soil salinity and sulfur concentration. It was observed that the application of organic matter had a signiﬁcant beneﬁcial impact on the chemistry of soil solutions, but did not show a remediation effect by increased sulfur leaching in a short-term study.


Introduction
Currently, Poland is the only country where native sulfur is extracted on a large scale. In 2015, 627.6 thousand Mg sulfur was mined [1,2]. Sulfur deposits occur in southern Poland, in the areas near Tarnobrzeg (Osiek, Baranów, Machów, Jeziórko deposits), Staszów (Solec and Grzybów deposits), and Lubaczów (Basznia field) [2]. Sulfur deposits were first extracted in the 1960s using the open strip mining method (Piaseczno Mine and Machów Mine), and the borehole Frasch method [3,4]. This method of melting an underground sulfur deposit was patented by Herman Frasch in 1894 in the United States. It involves an injection of superheated water at a temperature of 140-160 • C into a sulfur field. Subsequently, sulfur is melted underground and pumped to the surface in liquid form [5]. This method was used in Poland in the following mines: Grzybów (1966Grzybów ( -1996, Basznia (1977Basznia ( -1993, Machów II (1985-1993, Jeziórko , and around the world in mines in the United States Reclamation treatment of the FSMJ site started in 1993. Of the over 2000 ha previously occupied by the mine, over 700 ha in total were reforested, of which 216.5 ha-where this research was conducted-are currently managed by State Forest National Forest Holding (Figure 1b). Reclamation treatment included removing the mine infrastructure (such as mine wells, pipelines, and access roads), improving the hydrographic conditions, landscaping, neutralizing the pH of sulfurous and excessively acidic soils by liming (average 400-500 Mg ha −1 of sludge lime), fertilizing (70 kg ha −1 P2O5, 60 kg ha −1 K2O), and sowing grass seeds [4]. These treatments were followed by reforestation, mainly with one-year-old Scots pine seedlings (Pinus sylvestris L.) and an admixture of silver birch (Betula pendula Roth) and boreal oak (Quercus rubra L.) [3].

Field and Laboratory Study
As part of field work, a total of 88 sampling points were located in the FSMJ site in a regular grid of squares with sides measuring 150 m ( Figure S1). The sampling points were first identified on the map using ArcGIS software (ESRI Poland (Environmental Systems Research Institute), Warszawa, Poland) [14], and then targeted and stabilized using the Garmin GPS receiver (GPSMAP ® 60CSx, Garmin Ltd., Olathe, KS, USA).
Soil samples were collected from the top layer (0-20 cm) of the monitoring points [15]. In the laboratory, the collected soil samples were analyzed for grain size using the Fritsch GmbH Laser Particle Sizer ANALYSETTE 22 (Fritsch GmbH, Idar-Oberstein, Germany); pH was measured potentiometrically in 1 mol L −1 KCl (maintaining a ratio of soil:solution 1:2.5), and electrical conductivity (EC) was measured at a ratio of soil:solution 1:5. Soil organic Carbon (SOC), total Nitrogen (NT), and total Sulfur (ST), were measured with the TruMacLeco CNS analyzer (Leco Corporation, Saint Joseph, MI, USA) [16]. The degree of soil sulfur contamination was assessed on the basis of the guidelines of the Institute of Soil Cultivation and Fertilization Institute of Crop and Soil Fertilization (IUNG, Poland) [17].
Based on GPS measurements, the exact position of the research sites was determined, and data was collected to make a photo-interpretation key (geotagging of photographs). After determining land-cover classes (orthophoto screen vectorization in the ArcMap ArcGIS 10.2, Esri program [14]), patches of microhabitats were selected in three categories (Figure 3), i.e., an effectively reclaimed and forested site category (pine P and birch B stands, four plots each), and a degraded D stand with little herbaceous vegetation cover and occasional trees from natural succession (four plots). Detailed research was conducted on 0.01-ha (r = 5.64 m) circular study plots, from which composite soil samples were collected from a 0-5 cm and 5-40 cm horizon (sampling at five points arranged in a Reclamation treatment of the FSMJ site started in 1993. Of the over 2000 ha previously occupied by the mine, over 700 ha in total were reforested, of which 216.5 ha-where this research was conducted-are currently managed by State Forest National Forest Holding (Figure 1b). Reclamation treatment included removing the mine infrastructure (such as mine wells, pipelines, and access roads), improving the hydrographic conditions, landscaping, neutralizing the pH of sulfurous and excessively acidic soils by liming (average 400-500 Mg ha −1 of sludge lime), fertilizing (70 kg ha −1 P 2 O 5 , 60 kg ha −1 K 2 O), and sowing grass seeds [4]. These treatments were followed by reforestation, mainly with one-year-old Scots pine seedlings (Pinus sylvestris L.) and an admixture of silver birch (Betula pendula Roth) and boreal oak (Quercus rubra L.) [3].

Field and Laboratory Study
As part of field work, a total of 88 sampling points were located in the FSMJ site in a regular grid of squares with sides measuring 150 m ( Figure S1). The sampling points were first identified on the map using ArcGIS software (ESRI Poland (Environmental Systems Research Institute), Warszawa, Poland) [14], and then targeted and stabilized using the Garmin GPS receiver (GPSMAP ® 60CSx, Garmin Ltd., Olathe, KS, USA).
Soil samples were collected from the top layer (0-20 cm) of the monitoring points [15]. In the laboratory, the collected soil samples were analyzed for grain size using the Fritsch GmbH Laser Particle Sizer ANALYSETTE 22 (Fritsch GmbH, Idar-Oberstein, Germany); pH was measured potentiometrically in 1 mol L −1 KCl (maintaining a ratio of soil:solution 1:2.5), and electrical conductivity (EC) was measured at a ratio of soil:solution 1:5. Soil organic Carbon (SOC), total Nitrogen (N T ), and total Sulfur (S T ), were measured with the TruMacLeco CNS analyzer (Leco Corporation, Saint Joseph, MI, USA) [16]. The degree of soil sulfur contamination was assessed on the basis of the guidelines of the Institute of Soil Cultivation and Fertilization Institute of Crop and Soil Fertilization (IUNG, Poland) [17].
Samples of surface water were collected from 15 monitoring points on reservoirs and watercourses ( Figure 2) in summer, autumn, winter, and spring. The water samples were analyzed for pH, EC, and ionic content (Ca 2+ Mg 2+ , Na + , NO 3− , SO 4 2− , Cl − , HCO 3 − ) on a Dionex Ion Chromatography 5000 apparatus [18,19]. The spatial variability of ionic content and hydrogeochemical classification were assessed according to Szczukariew-Prikłoński classification [20]. Based on GPS measurements, the exact position of the research sites was determined, and data was collected to make a photo-interpretation key (geotagging of photographs). After determining land-cover classes (orthophoto screen vectorization in the ArcMap ArcGIS 10.2, Esri program [14]), patches of microhabitats were selected in three categories (Figure 3), i.e., an effectively reclaimed and forested site category (pine P and birch B stands, four plots each), and a degraded D stand with little herbaceous vegetation cover and occasional trees from natural succession (four plots). Detailed research was conducted on 0.01-ha (r = 5.64 m) circular study plots, from which composite soil samples were collected from a 0-5 cm and 5-40 cm horizon (sampling at five points arranged in a regular pattern: one point in the middle, and four in the corners of the square, as inscribed in the circular study plots). Next, pH, EC, S T , SOC, and N T were analyzed in the samples, according to standard methodology (described previously). The degree of herbaceous plant cover on study plots was determined using the Braun-Blanquet scale, and an inventory of the existing trees was made. DBH and height Ht were measured in case of trees.
Sustainability 2018, 10, x FOR PEER REVIEW 4 of 17 regular pattern: one point in the middle, and four in the corners of the square, as inscribed in the circular study plots). Next, pH, EC, ST, SOC, and NT were analyzed in the samples, according to standard methodology (described previously). The degree of herbaceous plant cover on study plots was determined using the Braun-Blanquet scale, and an inventory of the existing trees was made. DBH and height Ht were measured in case of trees.   Sustainability 2018, 10, x FOR PEER REVIEW 4 of 17 regular pattern: one point in the middle, and four in the corners of the square, as inscribed in the circular study plots). Next, pH, EC, ST, SOC, and NT were analyzed in the samples, according to standard methodology (described previously). The degree of herbaceous plant cover on study plots was determined using the Braun-Blanquet scale, and an inventory of the existing trees was made. DBH and height Ht were measured in case of trees.   Samples for the chemistry of foliage and the nutrient supply assessment of plants were collected from sites in previously defined categories (D, P, B), including one-year-old (current year P c ) and two-year-old (P c+1 ) Scots pine needles and silver birch leaves. Foliage samples were collected from tree tops (southwest exposition) growing on reclaimed and reforested sites (pine and birch stands, i.e., P and B category plots). Also, mixed samples (one per each plot) of wood small-reed leaves (Calamagrostis epigejos (L.) Roth) were collected. This was because it was an undergrowth species that occurred in every research plot. The content of N and S in the plant tissue samples was determined using LecoTruMac, while the other macroelements (Na, K, P, Ca, Mg, and K) were determined after mineralization in a mixture of HNO 3 and HClO 4 (in ratio of 3:1) with an ICP OES ICAP 6000 Series spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) [16].
Spatial imaging (map compositions, Figures S1 and S2) of variations in the studied soil features, i.e., sulfur concentration in the soil, pH, and EC, and herbaceous and woody vegetation surface cover, was obtained with the IDW (inverse distance weighted) interpolation algorithm in ArcMap ArcGIS 10.2 [14].
Using the Statistica 12 software (StatSoft Polska, Kraków, Poland) [21], correlations between the concentration of S T in topsoil horizons and pH, EC, and vegetation surface cover were investigated. The correlation between soil chemistry and plant macronutrient supply was also studied using the ANOVA test and Pearson correlation.
As part of the studies on the biogeochemical transformations of mineral sulfur and soil dissolution chemistry, changes in the highly contaminated environment of the experiment under controlled conditions were planned. In the experiment, the possible inhibitory impact of organic matter in the form of litterfall collected from birch and pine stands was taken into account. We planned two different variants of soil substrate sulfur contamination: LS (low sulfur substrate)-average 5090 mg kg −1 S, and HS (high sulfur substrate)-42,500 mg kg −1 S. At the beginning of the experiment, the input (starting) properties of the substrates prepared and used in the experiment were determined (soils: grain size, pH, EC, S T , SOC, N T , Ca, and Mg content, litter: pH, S, N, C, Ca, and Mg content using the methodology described previously, and the C:N ratio was calculated [22]). Soil material was collected in the post-mining site (substrates from plots in the previously described D category-degraded [3]). Litterfall was collected from managed stands to eliminate the impact of sulfur from the mine. In the course of the experiment, the soil-litter composites were rinsed with 100 mL of distilled water twice a week, and once a week, the filtrate was collected to determine the pH, EC, dissolved organic carbon (DOC), N T , Ca, Mg, and S content [22]. After the experiment, the properties of the used substrates and litter were analyzed again (substrate: pH, EC, S T , SOC, N T , Ca, and Mg content; litter: pH, S, N, C, Ca, and Mg content, and the C:N ratio was calculated) [22].

Soil pH and Contamination Spatial Variability
Soil pH KCl varied from 2.6 to 7.6 (Table 1, Figure S2a). Sampling points with acidic pH (pH KCl from 2.7 to 6.5) constituted 49%, while neutral pH (pH KCl from 6.6 to 7.2) and alkaline pH (pH KCl >7.3) were found on 38% and 13% of the study area, respectively. The soil EC ranged from 15 µS cm −1 to 2080 µS cm −1 (Table 1, Figure S2b). Only at two monitoring points was the EC value exceeded at 2000 µS cm −1 , which according to the FAO (Food and Agriculture Organization of the United Nations) guidelines [23] means that the growth of sensitive plants may be impeded. The S T content in the topsoil layers was from 10 mg kg −1 to 45,740 mg kg −1 (Table 1). S T content, which is similar to the natural content according to IUNG guidelines [17] was found on 35% of the analyzed site. According to these guidelines, elevated S T content was reported on 28%, high S T content (slight contamination) was reported on 12%, and very high content (heavy S contamination) was reported on 26% of the study area ( Figure S2c). Extreme sulfur concentrations exceeded the highest IUNG scale by more than 10-fold. S T correlated positively with EC (r = 0.80) and pH (r = 0.42). In some spots that were extremely Sustainability 2018, 10, 2442 6 of 16 contaminated with sulfur, there was more than 30% herbaceous plant cover (three and four on the Braun-Blanquet scale). The high S T concentration in the soil did not result in heavy acidification due to neutralization with high doses of sludge lime. However, high EC and limited vegetation were reported alongside a high sulfur concentration in soil samples, which may be considered a good bioindicator of ineffective reclamation.
Based on the conducted research, it may be concluded that the reclamation to forest was not fully effective due to significant sulfur contamination. However, locations with high sulfur contamination were considered hotspots (about 2% of the monitored site), and the distribution of topsoil contamination was irregular and dispersed. The most contaminated regions were concentrated mainly in the northern and central parts of the studied site.

Surface Water Chemistry Variability
The average water pH range, which was calculated on the whole database, was from pH 2.6 to 7.9 (Table 2). In the monitored site, hotspots were reported with extremely low pH values ranging from 2.6 to 3.1. Differences in pH between the studied seasons were insignificant. The average EC for the entire investigated period was 1797 µS cm −1 , and ranged from 1083 µS cm −1 to 3430 µS cm −1 , as shown in Table 2. According to standards [24], high EC indicating the salinity of waters is EC >1500 µS cm −1 , and we found that 82% of the analyzed samples came in this category. SO 4 2− content ranged from 472 mg L −1 to 1503 mg L −1 ( Table 2). The concentration of sulfate ions in seasonal variability was four-fold to 12-fold higher than defined in the current standards [24] for good quality waters (i.e., SO 4 2− < 138.5 mg L −1 ). Thus, generally, water mineralization in the studied site was high. The average total mineralization for the entire monitoring period was 1489 mg L −1 . Sulfate ions (63%) and calcium ions (22%) accounted for the largest share. High SO 4 2− ion content was related to high soil sulfation, while increased Ca 2+ ion concentration was due to applying doses of sludge lime for neutralization in the course of site reclaim treatments. To sum up, the problem in the case of the investigated surface waters was mainly the outcome of extremely low pH values recorded at some points, as well as high EC and high SO 4 2− concentration. The concentration of the remaining ions in the surface waters of the investigated site was relatively low (Table 2), and had no impact on the overall water chemical quality assessment. SO 4 2− and Ca 2+ ions determined the chemistry and hydrogeochemical classification of water.
According to Szczukariew-Prikłoński hydrogeochemical classification [20], the investigated water was mainly of sulfate-calcium class (58%). The investigated water was not natural mineral water, as indicated in Szczukariew-Prikłoński hydrogeochemical classification, but rather extensively anthropogenically transformed surface waters. Water chemistry was directly related to the specific geochemical system of the anthropogenically-transformed environment due to sulfation caused by borehole sulfur extraction and liming in the course of technical reclamation.

Soil and Vegetation Relations
It was confirmed that areas classified as degraded (D) typically displayed high sulfur contamination, i.e., from 50 mg kg −1 to 28,900 mg kg −1 , acidic pH (pH KCl from 1.9 to 6.5), and high salinity (EC from 149 µS cm −1 to 1160 µS cm −1 , as shown in Table 3). The soils also displayed low SOC and N T content (Table 3). Such unfavorable soil properties led to significant forest growth inhibition or even dieback. The density, i.e., the average number of trees ha −1 (N) was 612, the average DBH (diameter at breast height) was 7.7 cm, and the average height H t was 5.5 m (Table 4). In these conditions, there is a significant inhibition of vascular plant succession (herbaceous vegetation surface cover on average 26%, as shown in Table 4). Category D surfaces are an example of inefficient reclamation, especially in terms of neutralization and the isolation of highly S-contaminated horizons. Surfaces that were classified as successfully reclaimed and reforested (categories P and B) had a lower soil sulfur concentration. In the P category, S T concentration in the soil ranged from 10 mg kg −1 to 1200 mg kg −1 (Table 3), while in the B category, it ranged from 90 mg kg −1 to even 5300 mg kg −1 ( Table 3). Soil pH in these categories was also higher than in the category D, although it also varied widely from 3.9 to 7.1 and from 4.5 to 6.6 in category P and category B, respectively (Table 4). In the P and B site categories, significantly more trees with apparently better growth parameters were reported as opposed to category D, where trees were rare and deformed (Table 3).
However, it should be emphasized that category D covered a small total area of 4.39 ha in the entire facility. This, of course, compared with hundreds of effectively reclaimed hectares of land may seem little, but on the other hand, these areas represent extreme site conditions for plants, microorganisms, and animals, and the uniqueness of the reconstructed biochemical system being studied.

Foliage Chemistry and Macronutrient Supply
Macronutrient supply to trees in the reclaimed post-mining sites were generally adequate, i.e., within the ranges determined by ICP-Forest (the International Co-operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests) [25]. The mean P content in birch foliage was 2320.4 mg kg −1 ; in pine needles (in P c ) the mean was 1348.1 mg kg −1 , whereas in older needles (P c+1 ), the mean was 1045.4 mg kg −1 . The mean K content in birch leaves was 0.80%; in pine needles (in P c ), the mean was 0.50%, and in P c+1 , the mean was 0.50%, respectively ( Table 5). The mean Ca content in birch leaves was 1.35%; in pine needles (in P c ), the mean was 0.33%, and in P c+1 , the mean was 0.75%, respectively ( Table 5). The mean Mg content in birch leaves was 1526 mg kg −1 ; in pine needles (in P c ), the mean was 705 mg kg −1 , and in P c+1 , the mean was 715 mg kg −1 , respectively. The mean Na content (Na) in birch leaves was 34.9 mg kg −1 ; in pine needles (in P c ), the mean was 10.3 mg kg −1 , and in P c+1 , the mean was 10.0 mg kg −1 , respectively (Table 5). Significant differences were found in P, K, Ca, Mg, and Na content between the birch and the pine, while differences between needle ages were found only in the case of Ca ( Table 5).
The highest S content was found in birch leaves (mean 1954 mg kg −1 ). The S content in pine needles in P c was 1272 mg kg −1 , and in P c+1 , it was 1221 mg kg −1 . Significant differences were found in the S content between the species (Table 5).
The N content, as a deficient nutrient in these conditions, amounted to 1.71% in birch leaves and to 1.12% and 1.11% in pine needles for P c and P c+1 , respectively. In the above-ground parts of wood small-reed, a deficit of N and P was found (mean N 0.62%, mean P 904 mg kg −1 in category D; mean N 0.93%, mean P 1351 mg kg −1 in category P; and mean N 0.75 mg kg −1 , mean P 1132 mg kg −1 in category B, as shown in Table 6). In category D, no elevated S content was found in wood small-reed tissue (1779 mg kg −1 ) compared with the other categories (2249 mg kg −1 for P and 1717 mg kg −1 for B, as shown in Table 6). The S content in tree foliage correlated positively with S T content in the soil (r = 0.61) and EC (r = 0.68), and negatively with N T content in the soil (r = −0.66). In the case of wood small-reed on all of the plots in the site categories, linear correlations between the mentioned features did not occur. It was found that wood small-reed probably applied a strategy based on the mechanism of blocking the uptake of sulfur from highly contaminated soil. It is definitely a species that is resistant to the existing conditions of environmental stress and sulfurous soils. Thus, we found that wood small-reed, as a succession species, has a beneficial phytostabilizing impact, and may be used in the remediation of such types of soils.

Chemistry of Soil Solution in the Leaching Experiment under Controlled Conditions
It was predictable that the rate and amount of leached elements depended firstly on the amount of sulfur in the substrate (LS or HS) after 12 weeks of the experiment. In the highly sulfated HS composites, there was a significant reduction in sulfur concentration solely due to rinsing them (HS 42.521 mg kg −1 , HS-c 35,634 mg kg −1 , HS-B 33,247 mg kg −1 , and HS-P 34,157 mg kg −1 , as shown in Table 7). It was observed that despite a significant reduction in S t content in the HS substrate composites, the level of sulfation was still very high, and exceeded the values provided by IUNG standards (the Institute of Crop and Soil Fertilization). The S content is very high (heavy contamination), i.e., above 1000 mg kg −1 [17] up to 35-fold. The results of the experiment in a less sulfated variant (LS) indicate that the rinsing of substrates has a significant impact on a big reduction of EC in the substrate (LS 1.87 mS cm −1 , LS-c 0.39 mS cm −1 , LS-B 0.59 mS cm −1 , and LS-P 0.23 mS cm −1 , as shown in Table 7). It was shown that pine litter has a significant impact on raising the pH of the soil substrate (LS-B 2.1 and LS-P 2.9, as shown in Table 7); however, it still remained phytotoxic (pH below 3.0). Low soil substrate pH during litter decay is a normal phenomenon in natural conditions, and it is due to the organic acids that are released, causing a drop in soil pH [26].
It was found that birch litter resulted in elevated N t leaching from HS and Mg leaching from both substrates (Table 8). At the beginning of the experiment, birch litter was richer in N t than pine litter; therefore, there was more leaching of this element in the case of composites with an addition of birch litter.
The obtained experiment results indicated that the application of litter, especially in the case of the most sulfur-contaminated soil, is not sufficient for detoxification in the short term (12 weeks). Probably in hotspots where previous neutralization was ineffective or not thorough, it will have to be repeated with higher doses of the neutralizer in the form of sludge lime, as was done on most of the former Jeziórko mine site [4].  Table 8. Mean values and standard deviations for pH, EC, content of dissolved organic carbon (DOC), and the concentrations of St, N t , Ca, Mg, and Al in leachates at the beginning and after 12 weeks of leaching through soil substrates and composites in controlled conditions.

Discussion
The total S content in uncontaminated soil usually ranges from 100 mg kg −1 to 1000 mg kg −1 [17,28]. The allowable level of sulfur concentration in agricultural soils is 500 mg kg −1 [29]. This level has been exceeded at 33 of the 86 points in the FSMJ. In areas susceptible to elevated sulfur content, such as in the vicinity of sulfur mines, processing plants, or in industrial areas, the recommended levels are frequently significantly exceeded. An example of a sulfur-contaminated site is the industrial areas of Alberta, Canada (a sour gas processing plant). There, the level of sulfur concentration in the soil ranged from 500 mg kg −1 to 6000 mg kg −1 [30]. In our study of the FSMJ, sulfur concentrations in the topsoil horizons reached up to 45,740 mg kg −1 , exceeding by 45-fold the typical soil levels provided by Kabata-Pendias et al. [17], Stevenson and Cole [28], and by 90-fold the standards provided by the Canadians [29]. This demonstrates that the extensive remediation efforts were not sufficient at many spots of the sampled area.
High sulfur concentrations in the soil are usually associated with acidification and low pH [31]. In areas affected by lignite mining, with Miocene deposits and in areas with iron and gold mines, high sulfur contamination is often connected with iron, aluminum, or ferric sulfate [11,32,33]. High concentrations of Al and Fe, which are common in post-mining soils, have a strong negative impact on soil biota and vegetation [34]. In these studied soils in the reclaimed FSMJ area, Al and Fe concentrations were not excessive when compared with values for natural soils, which are considered to be 10,000 mg kg −1 to 35,000 mg kg −1 for Al, and 300 mg kg −1 to 50,000 mg kg −1 for Fe [35]. The sulfur in these studied soils was found to be native (elemental) sulfur, and its origin is connected with the transformation of gypsum and anhydrite deposits in a specific underground environment. Thus, the nature of soil sulfur contamination here is different than in the case of sulfur derived from the weathering of pyrites and marcasites, so Al and Fe soil concentrations in the FJSM area are not a limiting factor for vegetation. In the FSMJ, there were no significant correlations between soil sulfur concentration and soil acidification because of the impact of neutralization with flotation lime in various doses during the initial reclamation process in this area. Sulfur concentration in the topsoil had no relation to vegetation cover and abundance, and it proved impossible to determine the direct toxic effect of sulfur concentration on vegetation. In the FSMJ, many points were found where sulfur content in the topsoil horizons exceeded 2500 mg kg −1 and where pH was over 7.0 [4]. In some spots that were heavily contaminated with sulfur, there was up to 70% herbaceous vegetation surface cover. This could be, as mentioned above, a result of using high doses of flotation lime to neutralize the sulfur during the initial reclamation, which clearly could increase the pH of the soil if properly applied. Obviously, other soil characteristics, including available water, soil compaction, and nutrient availability, may be limiting factors for vegetation cover. However, healthy tree stands (points with pine and birch stands) were observed in some of the investigated soils, where high soil sulfur concentration was a significant differing factor in comparison to the degraded sites (points without vegetation cover).
Based on data from an analysis of dry needle weight, it may be concluded that tree supply with nutrients Ca, Mg, and K is suitable, while N is a deficit nutrient for the trees occurring in the investigated areas. The content of biogenic elements in the foliage of the investigated tree species is within the range provided by ICP-Forest [25] as excellent element content [7]. With regard to the nutrition of pine and birch trees on reclaimed land from the Lusatian Mine District (Eastern Germany) [36], tree supply with N in the investigated former sulfur mine sites was exceedingly low for the pine (supply grade 1), and low for the birch (supply grade 2). P, K, and Mg content in foliage was sufficient or good [7]. In assessing forest tree stand nutrition, it is also important to examine the proportions between macronutrients, especially the N:P ratio [37]. The N:P ratio in pine needles was 7.6 (in birch leaves), 8.3 (P c ), and 10.6 (P c+1 ) [7]. The N:P ratio in pine needles obtained by Pietrzykowski et al. [38] in the former Piaseczno sulfur mine was 10. Pietrzykowski et al. [38] found no significant differences between the N:P in mine land and the control (managed stands in natural habitats). The results obtained in the FJSM yielded similar values as in the cited work. It may be thought that supply in biogenic element is rather good.
Calamagrostis epigejos is a plant that Holub found to be highly tolerant to abiotic stress factors [39]. Calamagrostis epigejos is a dominant species in central European post-mining lands at the initial successional stage [40]. We indicated that Calamagrostis epigejos is tolerant to high S concentrations in the soil. In degraded areas (elevated sulfur concentration), higher amounts of this element were not reported in the above-ground tissue parts of the plant compared with afforested areas, which displayed a lower sulfur content in the soil. In areas with pine stands (category P), N content in above-ground parts of Calamagrostis epigejos was similar to the content that was observed in natural habitats in the Podyjí National Park (0.84-1.03% N, where they amounted to between 1600-1700 mg kg −1 P) [39].
The condition of the surface waters was rated as good, but there were some points with low pH, high EC, and high SO 4 2− concentration (hotspots). High concentrations of Ca 2+ ions were found (in all of the sampling zones), which was associated with the release of this ion into the water from lime sludge for the neutralization and isolation of sulfur soils. Local hotspots with a low pH and high sulfate concentration indicated an AMD threat, which should be monitored. The results after 12 weeks of the experiment conducted under controlled conditions show that the application of organic matter, especially in the case of the soils that were the most contaminated with sulfur, is not sufficient for fast detoxification. The results indicate that the rate and amount of the leached elements depended primarily on the amount of sulfur concentration in the substrate. Birch litter intensified N t leaching from the substrate with HS, while Mg did likewise in LS and HS. In the experiment with rinsing technosol substrates and litter application in controlled conditions, Woś and Pietrzykowski [21] showed a significant impact of birch litter on the intensification of soil-forming processes and increased leaching of DOC, N t , and Mg. It is obvious that the application of organic matter, irrespective of its type, has a significant impact on the rate of DOC leaching. Menyailo et al. [41] found that some species affected soil pH, DOC, and Mg content, but did not observe the impact of individual species on Ca and N t content in soils. Also, Chodak and Niklińska [42] found that tree species had a significant impact on the chemical and microbial properties of mine soils, mainly through pH changes with SOC and N t content. Apparently, the load of base cations contained in the litter was not enough to neutralize the toxic acidity of the investigated substrates. To mitigate the negative effects of excessive sulfur concentration in soil associated with the sulfur industry, the use of lime (limestone, hydrated lime, quick lime) is recommended [43]. McTee [44] reported that CaCO 3 application was the most efficient way to improve the ecological properties of acidic soils contaminated with elemental sulfur. Although the pH of the substrates increased over 12 weeks of the experiment, it was still within the range of toxic values, and soil pH was found to be the major factor affecting soil microorganisms [45].
Probably in hotspots where the neutralization was ineffective or not thorough enough, it would be necessary to re-apply the process, taking into account higher doses of flotation lime, as has been done so far in most areas of the former Jeziórko sulfur mine [4]. In the next stage, it would be possible to introduce afforestation, which will gradually introduce litterfall and organic matter, which play a key role in soil-forming processes. Soil organic matter, even in its initial phase of accumulation, plays an important role in the nutrient balance in newly developed ecosystems on reclaimed and afforested post-mine sites [18]. Decaying organic matter and biochemical processes produce more intensive leaching and displacement of macronutrients into the soil profile, and soil development will occur, but without neutralization, it will be a very long-term process.

Summary and Conclusions
Based on comprehensive research, it was concluded that on reclaimed and afforested areas of the Jeziórko mine: • Hotspots that are highly contaminated with sulfur still occur. The hotspots cover a relatively small area of about 2% of the monitored site (216.5 ha). However, they typically display a significant disturbance of soil and surface water chemistry, affecting the reaction of vegetation.

•
The investigated surface waters mostly display relatively high pH (6.5 on average). Only in several points (three out 15), pH ranged from 2.6 to 3.1, indicating heavy acidification and AMD risk, so the area must be monitored. • It was found that at sulfur concentration in soils of up to 3% after suitable neutralization with sludge lime (at a minimum dose of 500 Mg kg −1 ), the introduced birch and pine stands displayed good growth parameters and viability. The remaining ones occur in spots constituting about 2% of heavily contaminated sites, where neutralization was not done properly. In order to introduce vegetation and achieve a satisfactory yield, thorough neutralization treatment should be repeated. • Wood small-reed proved extremely resistant to sulfur contamination and the relatively high salinity of soils, and it occurred even on extremely sulfated soils. Therefore, there is a perspective and an indication of the use of this species in the natural succession of similar sites as an alternative temporary biological stabilization. • Soil detoxification will not be achieved without neutralization. Changes in the chemistry of soil solution after the application of organic matter were due to the process of litter mineralization and the leaching of nutrients. This effect was more pronounced in the case of birch litter application, where more intensive nitrogen and magnesium leaching took place, which is important in the perspective of soil-forming processes and plant nutrition.