Trace element and radiological characterisation of ash and soil at a legacy site in the former Raša coal-mining area

The Štrmac legacy site is situated on the Istrian peninsula to the NW of the former Adriatic carbonate platform (30) and is a part the Outer Dinarides (31). In the late Carboniferous, the Adriatic platform was located in the north of Gondwana and subsequently formed a distinct carbonate platform in the early Jurassic only to disintegrate in the Cretaceous (32). Considering its paleogeology, the majority of the Istrian bedrock consists of carbonates, most of which are Mezozoic and early Cenozoic limestone. The Istrian peninsula consists three distinct geological units: 1) the Jurassic- Cretaceous-Eocene carbonate plain in the south and west, 2) the Cretaceous-Eocene carbonate-clastic zone across the thrust structure in east and north-east, and 3) the Eocene flysch sediments in the central part (32).

The study site (Figure 1) is located in the Raša River basin containing deposits of bituminous Paleocene Kozina limestone (33, 34). The settlement of Štrmac is situated in the southeastern part of the Istrian peninsula on “terra rossa” cambic soil formed by insoluble residues of Eocene limestone with flysch sediments and loess (30). It is close to the Plomin bay, which is a part of the northern Adriatic Sea (about 5 km of air distance).

Two field campaigns were carried out, in May and July 2019. In the first field campaign, ash and soil samples were collected from 13 locations in the Štrmac area, while in the second field campaign only ash samples were collected from 10 sampling locations (Table 1), with subsamples of different horizons along the vertical profiles (Figure 1d and 1e).

The sampling locations were chosen randomly, depending on accessibility, with the aim to investigate the potential influence of atmospheric events on the surrounding environment. Previous investigations at the site showed potential cytotoxicity of soil (33). Additionally, during sampling, we noticed private vegetable gardens in the vicinity, owned by the local community.

Soil and ash samples were air-dried, crushed, sieved through a 1 mm sieve, and homogenised in an agate mortar for further analysis.

Particle induced X-ray emission (PIXE) elemental analysis was performed at the Ruđer Bošković Institute Tandem Accelerator Facility with 2 MeV proton beam obtained from 1 MV Tandetron accelerator (HVE Tandetron 4110, High Voltage Engineering Europa B.V., Amersfoort, Netherlands). The beam spot size on the sample was 3 mm in diameter. Emitted characteristic X-ray lines were measured simultaneously with two X-ray detectors, while backscattered ions were measured with a surface barrier detector using Rutherford backscattering spectrometry (RBS). For low-energy X-rays we used a silicon drift detector (SDD) (Vitus H20, KETEK GmbH, Munich, Germany) with the active area of 10 mm² and 8 μm thick Be window. It is placed at the distance of 110 mm from the sample at an angle of 150 ° to the incident beam. For energies above 3.5 keV we used a Si(Li) detector (SSL80165, Canberra, Meriden, CT, USA) with 30 mm² active area and 25 μm thick Be window, positioned at the angle of 135 ° and covered with a Mylar filter (275 μm). For backscattered ions we used a surface-barrier detector (Ultra, Ortec, Oak Ridge, TN, USA), placed at the angle of 165 ° from the incident beam. The schematic of the standard chamber for PIXE/RBS experimental setup is shown in Figure 2.

Figure 2

Standard IBA experimental chamber for PIXE and RBS measurements with marked positions of detectors, sample holder, suppression, Faraday cup, and incoming beam

The energy spectrum was determined using a multichannel analyser coupled with a homemade data acquisition system SPECTOR. Each sample was irradiated with 0.6 nA with a total charge of 0.6 μC. The PIXE experimental setup was calibrated with two standards: Standard Reference Material (SRM) 2710 Montana Soil and PTXRFIAEA08 Natural Soil, provided as a reference material from the International Atomic Energy Agency in 2014 (35). Samples were ground, and their powder pressed into pellets of 1 cm in diameter using a standard press with a pressure of 6.4 t/cm² and then mounted on the sample holder with carbon tape.

Proton backscattered spectra were analysed with the SIMNRA simulation software (SIMNRA 7.0, Max Planck Institute for Plasma Physics, Garching, Germany) as described elsewhere (36) to determine the concentrations of major low Z elements (Z<10). For quantitative analysis of the obtained PIXE spectra we used the GUPIXWIN software (GUPIXWIN 2.2.0, University of Guelph, Guelph, ON, Canada) (37) with a fixed matrix solution approach on thick targets and input matrix composition taken from SIMNRA results.

The accuracy of PIXE analysis was checked using the Standard Reference Material (SRM) 2710 Montana Soil and PTXRFIAEA08 Natural Soil.

For radioactivity analyses, soil and ash samples were prepared as described elsewhere (38, 39). Briefly, the samples were sieved (maximum grain size of 2 mm), dried at 105 °C for three days, ashed at 400 °C, and packed in 200 mL sealed cylindrical containers.

Radionuclide activity in the obtained ash was determined from the activity of decay product with a shorter half-life (T_½) under the assumption of a secular equilibrium between radionuclides. In undisturbed soil, the secular equilibrium between ²³⁸U and ²³⁴Th and between ²³²Th and ²²⁸Ac was established naturally. However, the loss of gaseous ²²²Rn from the surface layer of soil and during sample preparation resulted in a disequilibrium between ²²⁶Ra and ²¹⁴Pb. In order to restore the equilibrium, sealed samples were left to rest for more than 30 days.

Radionuclide activity concentrations were measured using highresolution gamma-ray spectrometer with a high-purity germanium coaxial detector (Ortec GMX, Oak Ridge, TN, USA) at relative efficiency of 74.3 % and energy resolution of 2.23 keV, all at ⁶⁰Co 1.33 MeV. The counting time was 80,000–250,000 seconds. Efficiency and energy were calibrated using certified standards (CBSS2 MIX, Eurostandard CZ, Czech Republic). The participating laboratory is accredited according to the ISO/IEC 17025 standard and the quality assurance was carried out by participation in proficiency testing provided by the IAEA and the European Commission’s Joint Research Centre (40).

The detailed procedure for the analysis of soil samples, peak analysis, and self-attenuation corrections has been described in detail elsewhere (41–44). The typical values of the detection limit (for measurements of 80,000 s) were 1 Bq/kg for ²³²Th and ²²⁶Ra, 2 Bq/kg for ⁴⁰K, 3 Bq/kg for ²¹⁰Pb, and 4 Bq/kg for ²³⁸U.

The radiological effects of the legacy site and the surrounding soil on the local population were assessed by calculating the radium equivalent index (Ra_eq), external absorbed dose rate (Ḋ), and annual effective dose (E). The Ra_eq is used to define a uniform value with respect to radiation exposure in case the ash from the waste site were to be used as building material and is calculated using Equation 1. 1Raeq = ARa + 1.43 ATh+ 0.077 AKR{a_{eq}}{\rm{ = }}{A_{Ra}}{\rm{ + 1}}{\rm{.43 }}{A_{Th}}{\rm{ + 0}}{\rm{.077 }}{A_K} where A_Ra, A_Th and A_K are the activity concentrations of ²²⁶Ra, ²³²Th and ⁴⁰K, respectively.

The weights were based on the estimation that 370 Bq/kg of ²²⁶Ra, 259 Bq/kg of ²³²Th, and 4810 Bq/kg of ⁴⁰K produce the same gamma-ray dose. If the Ra_eq values exceed the threshold value of 370 Bq/kg, the material will produce exposure higher than 1.5 mSv/year to local residents (45–47).

The UNSCEAR guidelines (9) provide the absorbed dose rates (Ḋ) in nGy/h due to gamma radiation in air at 1 m above the ground for uniform distribution of naturally occurring radionuclides (²²⁶Ra, ²³²Th, and ⁴⁰K). The absorbed dose rates were calculated as follows: 2Ḋ=0.462ARa+0.621ATh+0.0417AK\dot D = 0.462{A_{Ra}} + 0.621{A_{Th}} + 0.0417{A_K}

The annual effective dose (E) in mSv was also calculated according to the UNSCEAR guidelines (9): 3E=Ḋ(nGyh−1)×8760(hyr−1)×0.2×0.7(SvGy−1)E = \dot D\left( {nGy{h^{ - 1}}} \right) \times 8760\left( {hy{r^{ - 1}}} \right) \times 0.2 \times 0.7\left( {SvG{y^{ - 1}}} \right) where 0.7 Sv/Gy is the conversion coefficient from the absorbed dose in air to the effective dose received by adults, while 0.2 is the outdoor occupancy factor assuming that adults spend 20 % of their time outdoors.

Correlation between trace elements and Bonferroni correction were run on the PAST statistical software (PAST version 4.17, University of Oslo, Norway) assuming the significance level of P<0.05. For correlations between radionuclides ²³⁸U, ²³⁵U, ²²⁶Ra, ²³²Th, ⁴⁰K, and trace elements we used Spearman’s rank correlation ran on IBM SPSS software (IBM SPSS Statistics for Windows, version 23.0, Armonk, NY, USA) also assuming the significance level of P<0.05.

Table 2 compares the values of major elements (given in %) and of minor and trace elements (given as mg/kg) measured in 12 soil samples with the SRM 2710 Montana Soil. Relative errors for most major elements are up to 10 %, while for some trace elements, errors are higher due to very low elemental concentrations that are around the PIXE limits of the detection.

Figure 3 shows typical PIXE spectra collected with the SDD detector optimised for low energy X-rays (Figure 3a) and Si(Li) detector for higher X-ray energy (Figure 3b).

Figure 3

PIXE spectra of the soil sample taken with (a) SDD and (b) Si(Li) detector induced with 2 MeV protons

Mean concentrations and standard deviations of major, minor, and trace elements of the 12 studied soil samples are listed in Tables 3–5. Table 3 shows the concentrations of major elements in percentages, while Tables 4 and 5 show concentrations of minor and trace elements as mg/kg. The last row of each table shows average limits of detection (LOD) for each element.

Earlier, Petrović (48) reported soil elemental composition at Štrmac measured with high-resolution mass spectrometry with inductively coupled plasma (HR-ICP-MS) (Tables 4 and 5). Only one soil sample was taken near our sampling location 1. The reported concentration ranges mostly correspond to ours, save for Cu, Sr, and Pb, which are higher (Table 5). Differences in Pb concentrations are probably owed to differences in the measurement methods employed, as sample preparation for ICP-MS involves dissolving samples in acids, which may render Pb and Rb extraction from the soil incomplete. Even so, this comparison confirms that our PIXE method is satisfactorily accurate.

Compared to the values from the literature for world soils (Table 6), our measurements show much higher levels for Ni, Cr, P, and S. Higher P levels apply only to samples 7–12. Pb, Zn, and Ni levels are also elevated in almost all samples, but to a lesser extent. Similar findings were reported by previous studies (26–29). Moreover, the concentration of Se is extremely elevated, but since these measurements are around the PIXE’s LOD, it is difficult to draw definitive conclusions about soil contamination with Se at Štrmac. Previous measurement with ICP-MS (16, 26–29) clearly confirm that high Se levels are owed to its high content in Raša coal, which is also true for coal in general (49). According to the Geochemical Atlas of the Republic of Croatia (50), coastal Croatia has the highest concentrations of most potentially toxic elements (PTE) in Croatia and higher than the world/European average.

S and V were chosen for spatial comparison because both elements are known to have elevated concentrations in Raša coal and ash (27), and their elevated concentrations in soil samples indicate soil contamination originating from the landfill. Although the comparison shows a rather weak correlation between S and V, their distribution by locations above and below the median mostly coincides (Table 7). The exceptions are location 12, with concentration above the median for S and below the median for V, and location 3, with concentration below the median for S and above the median for V. Generally, locations further away from the landfill have lower concentrations of either element, with the exception of locations 3 for S and 12 for V. This suggests that the source of soil pollution with S and V in Štrmac is ash from the landfill, but this connection is not unequivocal, as there are spatial variations in the distribution of elements.

Previous research (26–28) showed positive and statistically significant (p<0.05) Kendall-Tau correlation between S, Se, and V in Raša coal and Raša coal ash samples. In this study we determined a positive (r>0.5) and statistically significant (p<0.05) correlation between S and As (r=0.58; p=0.013), and a negative (r<-0.5) and a statistically significant correlation between S and Al (r=-0.587; p=0.012) and S and Fe (r=-0.514 and p=0.0278).

The R value for the S-Se correlation is 0.0237 (p=0.919) and for S-V 0.443 (p=0.0581). We believe that the absence of significant correlation between S and Se is owed to the unreliability of the PIXE method in characterising Se in these samples. The S-V correlation is close to significance, and it is possible that the rounding of measurement figures had an influence on this outcome. However, the Bonferroni correction showed that only Al-Fe (p=0.021), Rb-Fe (p=0.049), and Na-Ti (p=0.049) correlations were statistically significant. This indicates that some of the other correlation coefficients could be falsely significant, probably due to heterogeneous composition of soil samples, and that the correlation between these elements is the most reliable.

Due to the limitations of the PIXE method, the levels of rare earth elements – Ce, Be, Sc, Ri, and Li – were not determined in this study. However, we did find a positive and statistically significant correlation between Al and Ti, Fe and Ni, and Fe and Mn. The Cu- Zn correlation was positive but not statistically significant (r=0.112; p=0.631). The results for the positive S-As correlation, which would normally indicate that As contamination at the landfill originated from Raša coal, are not reliable, as As was mostly below the detection limit. For the same reason, we found no correlations between Co and other elements.

Correlation between Al and Fe (Fe also correlates with Rb and K, Ni and Zn) stands out due to high correlation coefficients and small p-values (r=0.927, p=7.18*10^-5). Considering the importance of Al, K, and Fe in the lithosphere and the formation of minerals, it is probably a natural correlation, but the additional correlation with Ni and Zn indicates a certain influence of the local foundry and/or heating plant in Štrmac.

Pollution indices (PIs) are a simple way to assess the degree of pollution of soil, but they should be taken with reserve (55). The simplest of PIs, the single pollution index (SPI), is calculated by dividing a single heavy metal or PTE concentration in soil with “background geochemical values” or reference values such as those found in the Geochemical Atlas of Croatia (50). Other PIs for individual elements require additional data from different soil horizons or from a greater number of samples (56), which was out of our study’s scope.

If we apply pollution categorisation described by Kowalska et al. (56), the soil pollution at the Štrmac landfill is very strong for S, Se, Co, and Cu (Table 8). This is in line with reports of high soil pollution indices in the wider area of Labin, including Štrmac for Hg, Cd, V, Se, Pb, Cr, Zn, Cu, U, and S (26–28), in which S had the highest and V the lowest index. Those indices were based on corresponding background levels reported in the Geochemical Atlas of Croatia (50). Our study shows a similar result regarding the relationship of S and V with other elements.

The lowest PIs are those calculated with the reference background levels taken from Croatian legislation for agricultural land (51). This legislation focuses on a small number of PTEs that the legislator considers relevant and sets their pollution limits. Even in these terms, the pollution is high for Cr and moderate for Cu and Ni. In contrast, the PIs for As, Pb, and Zn are below these limits for the soil to be considered polluted.

Tables 9 and 10 show activity concentrations on naturally occurring radionuclides ²³²Th, ²³⁸U, ²²⁶Ra, ²¹⁰Pb, and ⁴⁰K in soil and ash samples. Activity concentrations of ²³⁸U, ²²⁶Ra and ²¹⁰Pb were higher in ash samples, while ²³²Th and ⁴⁰K activity concentrations were higher in soil samples. These results are in line with previous reports (9, 10, 19). Furthermore, the ranges of ²³²Th, ²³⁸U, ²²⁶Ra, ²¹⁰Pb, and ⁴⁰K in soil samples (Table 10) are in good agreement with previous measurements of soil samples in this part of Croatia (43, 44).

Table 11 summarises radium equivalent indices (Ra_eq), absorbed dose rates (D), and annual effective doses (E) in ash and soil calculated from activity concentrations of ²²⁶Ra, ²³²Th, and ⁴⁰K. The Ra_eq in most ash samples exceeded the threshold of 370 Bq/kg. Its wide range (122–751 Bq/kg) is probably owed to the inhomogeneity of the mined coal used for heating. To make use of this ash, it should be mixed it with other non-radioactive materials in ratios that comply with Croatian legislation (57).

The calculated absorbed dose rates originating from ash are higher than the background values for the Istrian region (70– 80 nGy/h) (44), considering that the mean rate is 184 nGy/h, and the maximum is 346 nGy/h, about four times higher than the background rates.

The mean annual effective doses range from 0.05 mSv in soil to 0.42 mSv in ash. For comparison, the worldwide mean dose in soil is 0.07 mSv and ranges from 0.3 to 0.6 mSv across countries (9).

Spearman’s rank correlation shows positive (r>0.9) and statistically significant (P<0.05) correlation between radionuclides ²³⁸U, ²³⁵U ²²⁶Ra, ²³²Th, ⁴⁰K and trace elements Na, P, and Ti. Only Sr correlated negatively (r=-0.98; P=0.005) with ²³⁸U, ²³⁵U, ²²⁶Ra, and ⁴⁰K. We did not find any correlations between ²¹⁰Pb and any trace elements. However, due to a small number of soil samples, we cannot make solid conclusions on these correlations. In case of possible future remediation efforts, a more detailed investigation with a larger number of samples is required.

The future of the Štrmac legacy site is not known. Our findings, however, can inform remediation operations, should such decision be made. We already have positive precedents of remediation of coal fly ash deposition sites in Croatia (10, 58). Moreover, fly ash can be used in construction industry. In European countries (France, Denmark, Italy, Germany, and the Netherlands) 85–100 % of produced coal fly ash is used for the production of cement, concrete, and ceramics (59, 60). There is also potential to extract rare earth elements, which are highly concentrated in coal fly ash, for use in batteries, lightweight alloys, and medical equipment, amongst several other applications, as reported elsewhere (61).