Fabrizio
Bardelli
*ab,
Valentina
Rimondi
c,
Pierfranco
Lattanzi
d,
Mauro
Rovezzi
e,
Marie-Pierre
Isaure
f,
Andrea
Giaccherini
c and
Pilario
Costagliola
c
aUniversità di Roma La Sapienza, Roma, Italy. E-mail: fabrizio.bardelli@gmail.com
bCNR-Nanotec, Italy
cUniversità di Firenze, Italy. E-mail: pilario.costagliola@unifi.it; andrea.giaccherini@gmail.com; valentina.rimondi@unifi.it
dCNR, IGG, Italy. E-mail: pierfrancolattanzi@gmail.com
eUniv. Grenoble Alpes, CNRS, IRD, Irstea, OSUG, FAME, Météo France, Grenoble, France. E-mail: mauro.rovezzi@esrf.eu
fUniversité de Pau et des Pays de l'Adour, E2S UPPA, CNRS, IPREM, Pau, France. E-mail: marie-pierre.isaure@univ-pau.fr
First published on 5th August 2022
Tree bark near former mercury (Hg) mines and roasting plants is known to have exceptionally high (up to several mg kg−1) Hg concentrations. This study explores the change of Hg speciation with depth (down to 25–30 mm from the outermost surface) in black pine (Pinus nigra) bark by means of high-resolution X-ray absorption near edge structure (HR-XANES) spectroscopy at the Hg LIII-edge. Principal component analysis and linear combination fitting applied to the HR-XANES spectra suggested that in the outermost layer (∼0–2 mm from the surface), roughly 50% of Hg is in the form of nanoparticulate metacinnabar (nano-β-HgS). A progressive increase in Hg-organic species (Hg bound to thiol groups) is found in deeper bark layers, while nano-β-HgS may decrease below the detection limit in the deepest layers. Notably, bark layers did not contain cinnabar (α-HgS), which was found in the nearby soils along with β-HgS (bulk), nor Hg0, which is the main Hg species in the atmosphere surrounding the sampled trees. These observations suggested that nano-β-HgS, at least in part, does not originate from mechanically trapped wind-blown particulates from the surrounding soil, but may be the product of biochemical reactions between gaseous elemental Hg and the bark tissue.
Environmental significanceA detailed knowledge of Hg speciation in environmental samples is critical to correctly predict the potential impacts. However, performing X-ray near-edge absorption spectroscopy (XANES) in samples with Hg amounts <100 mg kg−1 is challenging, because of the low signal-to-noise ratio of the resulting spectra. Moreover, the large core-hole lifetime broadening at the Hg LIII-edge results in poor spectral features, preventing a robust analysis of the data. The use of high-resolution XANES (HR-XANES) allows overcoming both limitations. In this work, HR-XANES was exploited for the first time to study Hg-contaminated tree bark samples. The results suggest that the bark, which was considered inert, may instead actively control the Hg biogeochemistry by transforming gaseous Hg into an insoluble, and thus scarcely bioavailable, Hg(II)-sulfide. |
The largest part of Hg present on or in the vegetation originates from gaseous (atmospheric), elemental Hg.9,10 Forests, in particular, are essential to regulate the Hg biogeochemistry cycle, due to the strong affinity of Hg to organic matter in soils and vegetation. For this reason, vegetation has been proposed as a biomonitor to evaluate the extent of Hg deposition in terrestrial ecosystems.11–14 On the other hand, forests may recycle Hg back to the atmosphere, because Hg uptake may involve both reversible and irreversible mechanisms.15 Wildfires are another variable to consider; indeed, in the 1997–2006-decade, Hg release due to biomass combustion accounted for 8% of the global (natural + anthropogenic) Hg emissions.16,17 While this contribution is minor on a global scale, local or regional impacts could be much more pronounced,16 and biomass burning is believed to be a major driver of the inter-annual Hg variations in the troposphere.18
Forest vegetation removes atmospheric gaseous Hg0 through stomatal and non-stomatal uptake.19–21 Similarly, tree bark may concentrate significant amounts of total Hg, and is proposed to be efficient monitors of airborne Hg.12,14,22–26 In particular, a study investigating the potential of Pinus pinaster as a biosorbent of Hg2+ suggested that the metal is adsorbed from aqueous solution by an ion exchange mechanism.27 However, the mechanism of Hg uptake at the bark/atmosphere interface, where Hg is essentially present as Hg0,28 is only partially known.29 Hanson et al. (1997)30 indicated that Hg0 is also re-emitted from bark, although at rates probably minor compared to those of soil or foliage.
After plant tissue senescence, the plant detritus (i.e., litterfall), composed also of tree bark, brings Hg back to the soil reservoir.31–33 Here, Hg may eventually be sequestered by forming insoluble Hg-phases,34 evade back to the atmosphere, or even methylate,35 possibly increasing its geo/bio availability.36 It is therefore important to know the exact speciation of Hg contained in tree bark, which, in turn, would contribute to control the Hg exchange rates with the environment.
Black pine (Pinus nigra) bark from Monte Amiata (Italy) is known to contain exceptionally high (up to several mg kg−1) Hg content, because of the proximity to former Hg mines and roasting plants of the third largest Hg district in the world.12,37 Over a century of mining and roasting caused Hg dispersion in the atmosphere and in the drainage network. Moreover, the area hosts a geothermal field, which is exploited to produce energy. These geothermal power plants are an additional source of Hg in the environment.38
In our previous study, conventional X-ray absorption near edge structure spectroscopy (XANES) data suggested that Hg speciation in superficial bark (<16 mm) included metacinnabar (β-HgS), possibly Hg0, and organic Hg-species (namely, Hg bound to cysteine or tannic acid).29 However, the comparatively low energy resolution of the data and the limited set of Hg references prevented a reliable determination of the trend of Hg speciation with depth. The study of the Hg content and speciation with depth is important to evaluate if bark can be used as a biomonitor for atmospheric Hg pollution. For example, Arnold et al., 201839 claimed that Pinus bark is not a good passive sampler for Hg air concentrations, because they found no correlation between Hg bark and total gaseous Hg concentrations. However, it should be noted that the total gaseous Hg concentration that the authors used to expose saplings was much lower than that detected in the area surrounding the trees sampled in this work. Moreover, these authors did not differentiate between inner and outer bark.
In this study, the change of Hg speciation with depth (down to ∼30 mm from the outermost surface) in Pinus nigra outer bark from Monte Amiata is explored by exploiting, for the first time for this subject, high resolution XANES (HR-XANES). The Hg speciation in bark is compared with that in the corresponding soils collected nearby the sampled trees, and the implications for the Hg biogeochemical cycle are discussed.
Fig. 1 (a) Study area and sampling sites considered in the present work; (b) sampling sites considered in Rimondi et al. (2020).37 (c) Zoomed view of the investigated area. The mining and smelting area near to the Abbadia San Salvatore town (darker grey area) is indicated by the dotted line contour. The trees selected for the HR-XANES study are reported (green tree symbols). |
Sample ID | Tree age (years) | Depth (mm) | Hg in bark (mg kg−1) | Hg in soil37 (mg kg−1) | Hg0 in air40 (ng m−3) |
---|---|---|---|---|---|
A11 | 100 | 0–2 | 8.6 ± 1.9 | 480 | 9.86–15.7 |
2–4 | 20.4 ± 6.4 | ||||
4–8 | 5.3 ± 0.7 | ||||
A13 | 150 | 0–2 | 16.5 ± 2.4 | 186 | 16.0–17.9 |
2–8 | 9.8 ± 0.3 | ||||
8–15 | 8.9 ± 1.0 | ||||
15–25 | 1.4 ± 0.08 | ||||
A19 | 100 | 0–5 | 13.7 ± 1.3 | 66 | 3.82–4.24 |
5–10 | 11.3 ± 0.3 | ||||
10–19 | 7.6 ± 0.3 | ||||
19–29 | 2.6 ± 0.03 | ||||
A22 | 150 | 0–2 | 8.3 ± 0.05 | 66 | 7.48–14.8 |
2–9 | 4.3 ± 1.3 | ||||
9–15 | 3.8 ± 0.01 | ||||
15–25 | 3.0 ± 0.1 | ||||
A49 | 100 | 0–2 | 9.0 ± 2.3 | 163 | 25–116 |
2–6 | 7.7 ± 0.02 | ||||
6–13 | 5.9 ± 0.03 | ||||
13–27 | 1.5 ± 0.01 |
The Pinus nigra bark typically consists of easily split layers, having a thickness of some millimeters. Bark blocks of about 30 mm from the bark surface were sampled. The phloem (living tissue) was not reached, so that the sampled bark layers were composed of tissue with no active translocation of soluble organic compounds. More details about the sampling procedure can be found in Rimondi et al., 2020.37 Soils were collected as composite samples: three splits from the first 10 cm ground layer from random locations around the tree, outside the tree canopy, and after removal of the tree litter.
To achieve the best possible signal-to-noise ratio in absorption spectra, nineteen bark samples from five trees with the highest Hg total content (>1 mg kg−1) were chosen. The Hg content of the sample portions (about 0.5 g) brought to the synchrotron facility are reported in Table 1. Despite the relatively high amount of Hg in the samples, HR-XANES measurements remained challenging. In particular, the acquisition of spectra of the samples containing less than 10 mg kg−1 of Hg with statistics comparable to those of the higher Hg-concentration samples would have taken too long compared to the available beamtime. Therefore, the deepest bark layer spectra, which correspond to the lowest Hg-concentrations, have a lower signal to noise ratio.
Soils were air dried under a chemical hood, sieved at 2 mm and ground for geochemical investigations. Total Hg contents, determined with a Milestone DMA-80 instrument using the same protocol described for bark samples to prevent the loss of volatile Hg species, exhibited high analytical variability. Therefore, a larger mass of sample (1 g) was digested in a sand bath with aqua regia (HCl/HNO3 3:1) and analyzed by inductively coupled plasma optical emission spectrometry (ICP-OES, PerkinElmer Optima 8000 instrument) equipped with a hydride generator (PerkinElmer B0507957). Two of the samples were processed in duplicate to evaluate method precision, which was within 15%. The analytical accuracy, checked against CCRMP STSD-1 (ICP-OES) international standards, was within 10%. The mercury concentration in soil samples ranged between 66 and 480 mg kg−1 (Table 1).
Portions of the same samples used for the chemical characterization were newly ground for synchrotron analysis in an agate mortar, then mixed and homogenized with cellulose and pressed into a pellet. This procedure should have ensured a more homogeneous distribution of Hg in the analyzed samples. Extreme care was taken to avoid cross contamination among the samples having different concentrations. Depending on the concentration of Hg, five to twenty consecutive scans were acquired and then averaged to improve statistics.
Reference Hg species measured by HR-XANES included: (i) inorganic Hg compounds, in particular, sulphides (α-HgS, β-HgS), which may derive from entrapment of windblown particles from the surrounding soil, and (ii) organic-bound Hg species, which may derive from the binding of Hg with organic matter in soil and/or in the bark. The Hg references HgCl2, HgO, CH3HgCl, and Hg-acetate (Hg(CH3COO)2) were purchased from Merck or Sigma Aldrich (purity > 95%). α-HgS (cinnabar) and β-HgS (metacinnabar) were natural specimens,42 which were evaluated for purity by X-ray powder diffraction. The elemental Hg reference (Hg0) was prepared by thoroughly mixing liquid Hg with epoxy resin, which was then cut with a microtome to make a film as thin as possible, suitable for acquisition in the transmission mode. The Hg-tannic acid reference (Hg-phenol) was prepared by dissolving 500 mg of tannic acid in Milli-Q® water, adding 2.5 mL of Hg(NO3)2 1000 mg L−1 solution, adjusting the pH to 5, and then freeze drying the solution. Liquid references, Hg-cysteine ([Hg] = 18 mM, M/L = 1/2, pH = 2.2, and Hg(Cys)2), CH3HgOH ([Hg] = 43 mM; pH = 9.4), CH3Hg-cysteine ([Hg] = 21 mM, M/L = 1/2, pH = 7.2, and CH3Hg-(Cys)2), and CH3Hg-glycine ([Hg] = 10 mM, M/L = 1/1, and pH = 7.1), were synthetized as described in Isaure et al. (2020), and measured as frozen liquids at 20 K. In Hg(SR)2 and Hg(SR)3, Hg is bonded with two or three thiolate groups, respectively (SR = cyclohexanethiolate).44,45 Hg(SR)2 has a polymeric structure consisting of bent S–Hg–S bonds (160.4°) linked together by weaker and farther Hg–S bonds. The Hg atoms have distorted tetrahedral coordination, and each thiolate acts as an asymmetric bridge between two Hg atoms. Hg(SR)3 contains Hg atoms in a distorted trigonal-planar coordination environment, with one angle larger than and one lower than 120°.46 The details of the synthesis of the above references are reported in Isaure et al. (2020),43 except for the Hg(SR)2 and Hg(SR)3 references, whose preparations are reported in Alsina et al. (1992).46
Further references are available in the ESI† of some published papers.45,47,48 These included nanoparticulate β-HgS (nano-β-HgS) and Hg bound to cysteine groups or to L-glutathione groups, both synthetized at pH 7.5 (Hg(Cys)2 and Hg(GSH)2, with GSH = γ-Glu–Cys–Gly),47 bulk and nanoparticulate HgSe,48 and Hg bound to thiolate and NO3 groups (Hg(SR)2 + (NO3)2).45 The HR-XANES spectra of all reference compounds are reported in Fig. S3† in the ESI.†
Principal component analysis (PCA) and linear combination fitting (LCF) were applied to the nineteen HR-XANES spectra to obtain quantitative information on the sample speciation. PCA allows determining the number of abstract components responsible for the main spectral features, which, in an actual system, can be related to the number of Hg species. It can also suggest the best candidates among the set of available references through a procedure named target transformation (TT49–51).
LCF consists of a weighted linear combination of the reference standards able to best fit a given experimental spectrum. PCA and LCF were performed using the SIXPACK and IFEFFIT packages, respectively.52,53
Notably, thanks to the higher energy resolution of the data, the presence of elemental Hg or organic-bound Hg (at least in amounts larger than ∼5–10%, which was estimated to be the maximum sensitivity of the HR-XANES data), which was a possibility in our previous study,42 can be ruled out. The different spectral features arising from the contribution of different relative amounts of α-HgS and β-HgS are clear, as highlighted in Fig. 2, where the spectra of the two representative soil samples, displaying the highest (lowest) and lowest (highest) β-HgS (α-HgS) contents (i.e. A19 and A49, respectively) and the corresponding LCF, are shown.
Fig. 2 (a) Comparison of β-HgS (metacinnabar) and α-HgS (cinnabar) and of two representative soil spectra (A19 and A49) highlighting the different spectral features. (b) and (c) Linear combination fitting (LCF) of the soil samples A49 (b) and A19 (c) performed using the β-HgS and α-HgS reference spectra, showing how different relative amounts of β-HgS and α-HgS reflect in different spectral features in the mixed phase. The numerical values of the LCF are reported in Table S2.† |
Principal component analysis, performed on the nineteen bark samples' spectra, suggested that all the spectra can be satisfactorily described with three principal components (see Fig. S7†). The target transformation procedure indicated β-HgS, nano-β-HgS, Hg(SR)2, Hg(GSH)2, CH3Hg(Cys)2, Hg(Cys)2, CH3HgCl, and Hg(SR)3 as the acceptable candidates among the set of reference compounds (see Table S3† for a list of the most relevant statistical indicators for the TT). Despite the information from PCA, all combinations of three (or more)-component LCF were attempted. However, in the final refinements, fits with more than two-components were discarded, because of the negligible changes in the χν2 and R2 statistical indicators with respect to two-component fits. More details on the LCF analyses are reported in a dedicated section at the end of the ESI† section.
Quantitative two-component LCF of the surface layer samples resulted in good matches using a combination of the nano-β-HgS and Hg(Cys)2 references (Fig. 3). The nano-β-HgS spectrum has a lower amplitude structure compared with its bulk counterpart, as is expected from a nanoparticulate material, and it was found to significantly improve the fits in the region ∼25 eV above the edge with respect to the bulk β-HgS one. When the nano-β-HgS fraction diminished below the sensitivity of the LCF (i.e. 5–10%, as estimated considering negligible changes in the χν2 and R2 statistical indicators), which usually occurred in the case of the deepest layers (Table S1†), the best fits were obtained using a combination of the Hg(Cys)2 and Hg(SR)2 references. This can be explained by the fact that the spectral features of Hg(Cys)2 and Hg(SR)2 are similar, while those of nano-β-HgS are significantly different. Therefore, when nano-β-HgS is present, the sensitivity required to discern between Hg(Cys)2 or Hg(SR)2 is highly reduced, and only one among these Hg organic references suffices to obtain good matches with the experimental data.
Fig. 3 (a) Comparison of the inorganic nano-β-HgS (nanoparticulate metacinnabar) and Hg(SR)2 and Hg(Cys)2 organic references, and of bark spectra representative of the outermost (0–2 mm) and inner (8–15 mm) layers of sample A13. (b) Linear combination fitting (LCF) of the bark samples A13 (layer 0–2 mm) and A13 (layer 8–15 mm) (c) using references Hg(SR)2 and Hg(Cys)2. The numerical values of the LCF are reported in Table S2.† |
It has to be noted that, for deeper layer samples, good results were also obtained with two-component fits performed using Hg(SR)2 or Hg(Cys)2 and the CH3HgCl reference.
With only one exception out of nineteen spectra, the distribution of the Hg species with depth showed the same general trend: the inorganic β-HgS fraction decreases from the outermost to the innermost layer, with a corresponding increase of the organic-bound Hg species. This trend is highlighted in a bar plot in Fig. 4, where the average total (bulk) Hg content is also reported. The inorganic Hg fraction is dominant in the outermost layer (46 ± 10% on average) and decreases down to 16 ± 15% on average in the innermost layer, while the organic Hg fraction increases accordingly. Likewise, the total Hg content decreases monotonically from the outermost to the innermost layer. In agreement with the above results, SEM-EDS observations detected the rare presence of discernible HgS particles (about one micrometer size) in the outermost bark layers (Fig. S5†), while no such particles could be observed in the inner portions.
Fig. 4 Bar plot highlighting the trend of the inorganic (nano-βHgS, red color) and Hg-organic (blue color) fractions (%) as a function of the depth of the bark layer. Total Hg content (mg kg−1) is also reported (gray color). The fractions shown are the averages of the values reported in Table S1,† for each layer. The bark layer thicknesses have been averaged, because they were slightly different for each bark sample (see the Materials and methods). The data show that the organic Hg fraction increases with depth at the expense of the inorganic one, and the total Hg content decreases with depth. The error bars represent the standard deviations for each layer (n = 5). |
In a previous study,29 it was suggested that β-HgS in the bark originated from mechanically trapped soil particles, likely as resuspension of the surrounding soil, where Hg sulphides were detected. This is supported by SEM-EDS observations, which revealed the presence of micrometer-size HgS particles at the bark surface. However, in some samples (Table S1†), β-HgS also occurs in significant amounts (>30%) in the deepest layers (down to 29 mm), where micron-sized HgS particles were not detected by SEM-EDS, suggesting that in deeper layers, the origin of β-HgS may be different. In addition, bulk β-HgS was detected in the soil surrounding the samples, while nanoparticulate β-HgS was detected in bark samples, suggesting a different origin for the Hg-sulphide species in soils and bark. Indeed, β-HgS may occur as a biogenic nanoparticulate, in agreement with recent studies,21,44,57 which show that β-HgS nanoparticulate clusters may be produced in low temperature systems, through formation of complexes with thiolates or cysteine. Specifically, Manceau et al.21 reported the formation of nanoparticulate β-HgS at the surface of leaves exposed to a Hg-bearing atmosphere by the reaction of gaseous Hg species via intermediate compounds, such as Hg(SR)2, which was detected in the deeper bark layers along with Hg-cysteine. At present, we cannot discriminate beyond any doubt between β-HgS of inorganic (metallurgical waste via soil) or biochemical origin. However, the results of the present work indirectly favor the second option because α-HgS was identified in the soil in almost the same proportion as that of β-HgS, but was notably absent from the bark matrix, suggesting that the origin of the two Hg-sulphides may be different.
Manceau and coauthors21,58 proposed two possible mechanisms for the formation of the Hg sulphide species in plant leaves or organic matter, both based on a biomineralization process: enzymatic sulphidization or chemical dealkylation. The authors, however, note that both mechanisms require Hg atoms to face each other, a possibility that increases with the total Hg content. In agreement with this scenario, in this work, higher nanoparticulate β-HgS amounts correlate with the total Hg content at the bark surface (Fig. 4).
The present HR-XANES data confirm that Hg is bonded with two S atoms, but the match with the Hg(Cys)2 or Hg(SR)2 references, or with a linear combination of them, is not exact. The actual local Hg coordination in the organic fraction is probably a more complex one, with the two primary thiol groups at an intermediate angle (>160° and <175°) and possibly, with one or more secondary bonds with amine, carboxyl, or hydroxyl groups, as proposed by Manceau et al.21
The possible formation of methylated species would be obviously of paramount importance for its ecotoxicological consequences. Indeed, LCF of some samples gives permissive evidence of the possible presence of methylated species, such as CH3HgCl. However, methylated Hg species in general are expected to be present at very low levels in the environment compared to other Hg species, μg kg−1vs. the mg kg−1 levels measured in the bark samples,59i.e. at much lower levels than it is required to introduce a significant contribution in the HR-XANES spectra (>5–10% of the total Hg content). Therefore, it is likely that the improvement of LCF observed when including the CH3HgCl reference is an artefact, probably arising from the poor ability of HR-XANES to distinguish adjacent elements, such as S and Cl. Substantial evidence against the presence of CH3HgCl comes from SEM-EDS observations, which did not show Hg bound to Cl.
The absence of Hg0 in bark layers, present as the dominant airborne Hg form in the surrounding atmosphere of the Monte Amiata mining district60 and which may have been present in the liquid form in the bark samples (detectable as solid (frozen) Hg0 in this work), is noteworthy. The results of this study suggest that gaseous elemental Hg is transformed at the bark surface by the reaction with thiol-containing proteins; although we cannot totally exclude a loss back to the atmosphere, despite the great care taken to prevent it (see the Materials and methods section). Possible transformations into comparatively scarcely mobile forms in the bark layers represent a sink that may slow down the diffusion of airborne Hg in the environment. On the other hand, bark may contribute to the occurrence of Hg in the litterfall through bark foliage and ultimately to soils. It has to be noted that no organic Hg compounds have been observed in soils surrounding trees, possibly because their concentration is lower than the detection limit of the HR-XANES technique (5–10% of the total Hg). This indicates that Hg sulphidation in bark may provide Hg to litterfall mainly as a sulphide, a theoretically poorly geo-available form of Hg, as found to occur in tree leaves.21 This conclusion supports the concept that the bark, which was considered inert, may actively contribute to the Hg biogeochemistry and eventually to the Hg cycle in the environment. However, it should be noted that the reactivity of nanoparticulate Hg species is not known at present, and more research on this subject would be important.
Footnote |
† Electronic supplementary information (ESI) available. See https://doi.org/10.1039/d2em00239f |
This journal is © The Royal Society of Chemistry 2022 |