Subscriber access provided by UNIV OF WESTERN ONTARIO
Environmental Processes
Mercury isotopic fractionation during pedogenesis in a tropical forest soil catena (French Guiana): deciphering the impact of historical goldmining Stéphane Guédron, David Amouroux, Emmanuel Tessier, Catherine Grimaldi, Julien Pierre, Gilbert Barre, Sylvain Berail, Vincent Perrot, and Michel Grimaldi Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b02186 • Publication Date (Web): 17 Sep 2018 Downloaded from http://pubs.acs.org on September 21, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 33
Environmental Science & Technology
202 Hg (‰) -4
-3
-2
-1
0
anthropogenic
goldmined pristine
Litterfall Palm tree leaves
Litter
Mineral horizon
Hg0 (reduction / evaporation)
Hg(II)
Hg0liq
mixing (litter – mineral)
Hg(II)
Complexation (OM & oxide)
Oxic (-) MDF (+) Anoxic (Fe oxides)
Hg (II)
(Fe reduction)
ACS Paragon Plus Environment
Hg(II)-DOM
Soil
Organo -mineral horizon
Hg0
Atmosphere
background
Foliar uptake Hg(II)
1
Environmental Science & Technology
Page 2 of 33
Mercury isotopic fractionation during pedogenesis in a tropical forest soil catena (French Guiana): deciphering the impact of historical goldmining S. Guédron*1, D. Amouroux2, E. Tessier2, C. Grimaldi3, J. Barre2, S. Berail2, V. Perrot1 and M. Grimaldi4,5.
1
Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, IFSTTAR, ISTerre, 38000
Grenoble, France. 2
CNRS / Univ. Pau & Pays Adour / E2S UPPA, Institut des Sciences Analytiques et de
Physico-chimie pour l’Environnement et les Matériaux – IPREM, UMR5254, 64000, PAU, France. 3
UMR SAS, INRA, Agrocampus Ouest, 35000 Rennes, France
4
Sorbonne Universities, Science Faculty, Paris 06, IRD, CNRS, INRA, UPEC, Univ Paris
Diderot, Institute of Ecology and Environmental Sciences, iEES Paris, 75005 Paris, France; 5
Institut de Recherche pour le Développement, Centre IRD de Cayenne, 97323 Cayenne
cedex, France
*corresponding author: S. Guédron Email:
[email protected] Keywords: Mercury, stable isotopes, amazon region, goldmining, pedogenesis, soil catena
1 ACS Paragon Plus Environment
Page 3 of 33
Environmental Science & Technology
1
Abstract:
2
We used mercury (Hg) stable isotopes to investigate the Hg cycle in a rainforest soil catena
3
(French Guiana) partially goldmined during the early 1950´s. Litterfall showed homogeneous
4
∆199Hg values (-0.18 ± 0.05 ‰, i.e. modern gaseous elemental Hg (GEM) isotopic signature).
5
After litter decomposition, Hg bound to organic matter (OM) is mixed with Hg from pristine
6
(-0.55 ± 0.22 ‰) or goldmined (-0.09 ± 0.16 ‰) mineral materials. Negative ∆199Hg values in
7
deep pristine mineral horizons (-0.60 ± 0.16 ‰) suggest the transfer of Hg bound to dissolved
8
OM depleted in odd isotopes due to mass independent fractionation during Hg abiotic
9
reduction. Perennial palm tree leaves collected above goldmined and pristine soil recorded
10
contrasted ∆199Hg signatures likely resulting from GEM re-emission processes from soils and
11
leaf surfaces. Upslope, soil δ202Hg signatures showed a negative shift (ε= ~ -1 ‰) with depth
12
attributed to mass dependent fractionation during Hg sorption and complexation onto iron
13
oxides and dissolved OM. Downslope, higher δ202Hg values in soils resulted from
14
hydromorphy (lower humification, higher Hg(II) reduction…). The unique Hg isotopic
15
signatures of Amazonian soils probably results in multistep fractionation processes during
16
pedogenesis (millions of years) and in potentially different Hg isotopic signature of pre-
17
anthropogenic background GEM.
2 ACS Paragon Plus Environment
Environmental Science & Technology
Page 4 of 33
18
Introduction
19
Tropical soils of the Amazon region have accumulated atmospheric mercury (Hg) during their
20
development that lasted millions of years 1. Their inherited Hg geochemical background is
21
higher than those reported in boreal and temperate climate soils 2-6. In French Guiana, highest
22
Hg concentrations were reported in upslope oxic ferralsols associated to the clay-size iron
23
oxides and sulfur-bearing functional groups of the organic matter (OM)
24
associated into aggregates 7-9. In the lowlands, where the water table comes up to the top soils,
25
gleysols are depleted in Hg because of prevailing reducing conditions that lead to the
26
dissolution of iron oxides. To this natural geochemical Hg background adds past and current
27
artisanal small-scale goldmining (ASGM) Hg contamination. Indeed, ASGM activities in the
28
Amazonian region use Hg to extract eluvial gold concentrated in soil and sediment, and emit
29
large amounts of Hg in the atmosphere (i.e., burning of amalgams) or in soils and rivers
30
through direct discharges of elemental and particulate Hg 10-12.
31
Natural Hg stable isotopes show large mass dependent (MDF) and mass independent (MIF)
32
fractionation ranges (~ 7 ‰ for δ202Hg and ~ 10 ‰ for ∆199Hg) which allow tracing natural or
33
anthropogenic Hg sources and identifying biogeochemical Hg pathways 13. In particular, Hg
34
MIF signature is reported to be efficient tracer of specific Hg sources since it is not altered
35
during MDF processes, but can be changed significantly by photochemical reduction, the
36
abiotic reduction of inorganic Hg(II) by dissolved organic matter (DOM)14 and the
37
evaporation of liquid elemental Hg15 causing anomalies of odd-mass relative to even-mass
38
isotopes 16, 17. For example, a recent study of ASGM impacted soils in Amapá (Brazil) showed
39
that ∆199Hg values were higher in mine tailings (~ 0 to −0.01 ‰) than in pristine soils (~ -
40
0.55 ‰)18. Besides, MDF signatures can be used as tracer of biogeochemical reactions since
41
they result from numerous processes including redox reactions, volatilization, complexation
42
and OM binding … 16, 19-21. Fractionation of Hg isotopes in soil during pedogenesis has been
3
which are closely
3 ACS Paragon Plus Environment
Page 5 of 33
Environmental Science & Technology
43
increasingly investigated during the last decade. Uptake of atmospheric Hg by foliage was
44
reported to result in a large shift of δ202Hg values (up to -2.9 ‰) suggestively caused by
45
kinetic MDF during foliar uptake of gaseous elemental Hg (GEM) while little fractionation
46
was observed during the incorporation of deciduous material in the soil 22-25. Downward in the
47
soil profile, Hg isotopic compositions can be linked to geogenic Hg mixed with atmospheric
48
Hg and biogeochemical processes controlled by hydro-pedological conditions (redox)
49
affecting Hg isotopic fractionation mostly by reduction and complexation processes 22, 25.
50
In this work, we studied Hg isotopic fractionation during pedogenesis in a soil catena located
51
in a small catchment area of French Guiana covered by tropical rain forest combining upslope
52
undisturbed soils and downstream reworked soils contaminated with Hg during ancient
53
goldmining activities dating the early 1950´s. Perennial palm tree (PT) leaves collected at
54
various distances from the goldmined site were used as bio-indicators of Hg atmospheric
55
signature below the canopy. We used Hg isotopes to track the dissemination of anthropogenic
56
Hg in the catchment and studied the Hg isotopic composition along with the downward Hg
57
incorporation in soil (from litter to deep mineral horizons) with consideration of geochemical
58
and hydro-pedological processes.
59
Material and methods
60
Research areas
61
Soils were sampled in October 2012 on the Combat Creek watershed (52°23’W, 4°35’N), a
62
small catchment of ~1 km² covered by tropical rain forest in French Guiana. Details on the
63
geology and morphology of the watershed are given in Guédron et al. 3. Four soil profiles
64
were collected along a toposequence (Fig. 1) including: (i) a ferralsol (IT4) located upslope
65
and having a high clay (< 2 µm size fraction) content and ferruginous nodules throughout the
66
entire profile with a micro-aggregated structure extending over 1 m depth allowing a good
67
vertical water drainage, (ii) a midslope acrisol (IT1) with a massive weathered schist (C) 4 ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 33
68
horizon having high fine silt content at shallow depth (< 1 m), (iii) a downslope acrisol (IT2)
69
with an hydromorphic weathered schist (Ch) horizon and (iv) a gleysol (IT3) typical of
70
hydromorphic conditions of the low lands with a permanent shallow groundwater table
71
imposing reducing conditions. In these soils, very common in French Guiana3, 9, quartz,
72
kaolinite, and Fe-Al-(oxyhydr)oxides dominate mineralogy, quartz being the sole remnant
73
from the primary rock-forming minerals.
74
Samples from the 3 soil profiles of the ancient gold mined “flat” (i.e.; III-0, SL3 and SL6)
75
were collected in 20053. Gold from the former flat was extracted using Hg for amalgamation
76
according to the ancient ‘Long Tom’ sluices process dating from the early 1950´s. The
77
resulting soils (initially gleysols) are strongly disorganized with gravel, sand and silt as
78
dominant grain-size fractions, the finest fraction being lost during mining operations. Hg
79
droplets and Au-Hg amalgams were identified in these gold-mined gleysols 3.The area where
80
SL6 and III-0 were collected has been deforested and re-exploited (2008-2010) without Hg
81
use (gravimetrical gold extraction).
82
Sampling collection and analysis
83
Soil and leaves collection. Soil samples were collected every 10 cm with an auger down to
84
120 cm depth when possible, homogenized and stored in clean polyethylene bags. Litter
85
samples were collected according to different layers of the organic horizon (O) with Oi the
86
recent or slightly humified litter and Oe the humified organic matter. Perennial palm tree
87
(Oenocarpus Bataua Mart.) leaves were collected in the vicinity of each soil profile along the
88
toposequence (IT4, IT1 and IT3), on a upslope ferralsol of the Creek Mouche (CM) located
89
upstream in the catchment, on the former goldmined flat (SL6) and at its edge (FP) (Fig. 1
90
and S.I.1). All samples were kept in a fridge at 4°C before being freeze-dried and crushed to
91
obtain a powder of d < 63 µm. More details for sample collection and conditioning are given
92
in S.I.1. 5 ACS Paragon Plus Environment
Page 7 of 33
Environmental Science & Technology
93
Determination of soil, litter, and leaves chemical composition. Total Hg concentration
94
([THg]) in solids was determined by atomic absorption spectroscopy after dry mineralization
95
and gold amalgamation (Altec, Model AMA 254). Triplicates were performed for each
96
analysis (S.I.2). Concentrations obtained for certified reference material (CRM MESS-3 –
97
National Research Council Canada) were in the published range of concentration (0.091 ±
98
0.008 µg g−1).
99
Total carbon ([C]), sulfur ([S]), and nitrogen ([N]) concentrations in samples were determined
100
from the dry combustion of soil sample aliquots using a Fisons® (model 1500CHNS) infra-
101
red analyzer. Granulometry was determined with standard sieving and sedimentation
102
procedures
103
sample aliquot in aqua regia solution (HCl/HNO3; 3:1 v/v ratio) for 10 h at 70 °C.
104
Mercury isotopes analysis (see details in S.I.2). Hg isotope ratios were determined after
105
digestion of 0.2 to 0.5 gram of the soil sample in 50 mL PP vials with addition 4 mL of
106
HNO3/HCl/H2O2 mixture and heated at 85°C for 4 hours according to the procedure
107
established by Foucher and Hintelmann
108
leaves, high pressure acid digestion has been achieved following the protocol validated by
109
Estrade et al.
110
according to the total Hg content) was digested in 4 mL suprapur HNO3 at 300°C and 130
111
bars for 3 hours (S.I.2). Then, all sample extracts were centrifuged at 1400 rpm for 10 min
112
and the supernatant was removed and diluted for isotopic measurement to about 10% of acid
113
mixture. Hg isotope analyses were performed according previous works
114
vapor generation with SnCl2 reduction coupled to a MC-ICP-MS (Nu Plasma, Nu
115
Instruments). To correct the instrumental mass bias, internal standard of Tl (NIST 997,
116
205
26, 27
. Fe and Al contents were determined by AAS after the digestion of a soil
30
28
and used previously in our lab29. For palm tree
for lichen specimen. Homogenized sample matrix (between 0.5 to 1.2 g,
16, 31, 32
using cold
Tl/203Tl = 2.38714) and sample standard bracketing with NIST 3133 (prepared in the same
6 ACS Paragon Plus Environment
Environmental Science & Technology
Page 8 of 33
117
acid matrix than samples) standard solution were used. MDF of Hg are reported relative to the
118
NIST 3133 Hg solution16:
119
δXXXHg = (XXX/198Hgsample / XXX/198HgNIST3133 -1) * 1000 in ‰ 199
Hg
200
Hg,
201
Hg and
204
120
where xxx is the studied isotopes. For
Hg isotopes, MIF of Hg is
121
reported as the difference between the theoretical value predicted by MDF of δXXXHg and the
122
measured values noted ∆XXXHg in ‰ 16:
123
∆XXXHg = δXXXHg - δ202Hg * βXXX
124
where βXXX are equal to 0.2520, 0.5024, 0.7520 and 1.493 for the studied isotopes (xxx) 199,
125
200, 201 and 204 isotopes, respectively. Enrichment factor ε is used to compare fractionation
126
or mixing of Hg isotopes for MDF between two pools (i.e., a sink or product (B) relative to a
127
source or reactant (A)) and expressed as follow: ε202HgB-A = δ202Hg B − δ202Hg A.
128
Analytical uncertainty was evaluated by multiple measurements of certified reference material
129
for Hg concentrations. All details for the method, Quality Assurance/Quality Control and
130
tables of Hg analysis data are given in S.I. 2.
131
Statistical treatment and data presentation. Linear regressions were performed only when
132
normality test, constant variance test and alpha power of the performed regression (α test, p >
133
0.8) passed. Correlation coefficient (R) and p values are reported for the computed linear
134
regressions. In the manuscript, mean concentrations (arithmetic means) are always presented
135
with the associated standard deviation (mean ± SD) and isotopic compositions with their
136
associated 2SD. When the range (2SD) of measured values of a group of samples is higher
137
than the 2SD values of the analytical measurement, mean values are presented with this 2SD
138
range and the number of corresponding measurements (N).
139
Results and discussion
140
Discrimination of isotopic Hg pools in a tropical forest soil catena.
7 ACS Paragon Plus Environment
Page 9 of 33
Environmental Science & Technology
141
The composition of Hg even isotopes generally follows the mass dependent fractionation
142
(MDF) theoretical line and exhibits a large range of delta values for both δ202Hg (~ 3 ‰ - Fig.
143
2) and δ200Hg (~ 1.5 ‰ - S.I.3.a). Mass independent fractionation of even isotopes (MIFeven;
144
∆200Hg or ∆204Hg) was weak and generally not significant in the samples investigated (S.I.3),
145
even if larger dispersion for ∆204Hg was observed resulting probably from larger analytical
146
uncertainty. Significant Hg mass independent fractionation of odd isotopes (MIFodd) is
147
observed within a 0.7 ‰ range for both ∆199Hg and ∆201Hg (Fig. 2 and S.I.3.b). They both
148
plot on the 1:1 empirically derived relationship resulting of fractionation during
149
photochemical reduction of inorganic mercury from aqueous solutions in the presence of
150
DOM 16 (S.I. 3). Owing to the contrasted signatures of both odd and even isotopes, our plot of
151
∆199Hg and δ202Hg allows identifying three major pools: (i) former goldmining soils with the
152
highest δ202Hg and ∆199Hg values of our dataset; (ii) foliage and soil litter that have
153
intermediate values and (iii) undisturbed soil horizons (A, B and C) enriched in lighter Hg
154
isotopes together with lower ∆199Hg values (Fig. 2).
155
Hg contaminated soils of the former goldmined flat ([THg] = 2360 ± 2533 ng g-1) have
156
isotopic composition for both MIFodd and MDF (∆199Hg = -0.09 ± 0.16 ‰; δ 202Hg = -0.65 ±
157
0.65‰, N=12) in agreement with other Amazonian goldmining tailings18 and cinnabar or Hg
158
ores
159
isolated elemental Hg (Hg0) droplets or Hg0 amalgamated with micrometric gold particles
160
inherited from goldmining operations dating from the early 1950´s in the Combat watershed 4.
161
During ancient mining operations, gleysols of the lowland were suspended with creek water
162
forming a slurry to remove fine particles and concentrate heavy gold-rich particles before
163
their onsite amalgamation with Hg 3. Thus, in opposition to undisturbed gleysols that belong
164
to the same lowland area, mined soils have ∆199Hg values negatively correlated with coarse
23, 33-35
. These near-zero ∆199Hg values are also consistent with previous observations of 3
8 ACS Paragon Plus Environment
Environmental Science & Technology
Page 10 of 33
165
sand fraction (R = -0.88, p 2mm) contents (bottom right panels). In upslope and midslope soils of the undisturbed toposequence (IT1, IT2 and IT3), the coarse fraction is essentially composed of ferruginous nodules. Vertical profiles of total Hg in soils impacted by former goldmining activities (right corner panel) are from Guédron et al. 3. Total carbon is equivalent to organic carbon since no carbonates are present in this rain forest setting. Soil of the toposequence are defined undisturbed because of their neighboring position to the flat. The recently reworked goldmined flat (grey area) is an open field area (deforested without organic horizon). The 3 soil profiles of this ancient gold mined “flat” (i.e.; III-0, SL3 and SL6) were collected in 2005 before its rework while palm tree leaves were collected in 2012. Figure 2. Mass-independent fractionation (Δ199Hg) vs mass-dependent fractionation (δ202Hg) in palm tree (PT) leaves and soil horizons and their comparison with literature data. The shaded areas correspond approximately to i.e., precipitation (blue)22, 24 ,47, 61, ores and mine tailings (orange) 18, 23, 33, 35, atmospheric TGM and GEM (grey)24, 46, 47, 60, and litterfall, foliage and lichens (green) 22, 24, 25, 36, 47, 48 data from the literature. A bidirectional error bar represent 2σ of replicated isotopic measurements of NIST RM 8610 (UM Almaden). Figure 3. Mass-independent fractionation (Δ199Hg) vs mass-dependent fractionation (δ202Hg) in palm tree (PT) leaves collected over forested ancient goldmined soils (SL6), at the edge of the reworked goldmined flat (FP), over undisturbed soils (IT4, IT1 and IT3) of the toposequence IT and over pristine remote soils (CM). Are also presented the mean value with 2SD of THg isotopic composition for the four litter samples of the undisturbed toposequence IT (dark green diamond), the soil horizons (A, B and C) below litter (red diamond, N = 22) collected in the oxic undisturbed toposequence soil profiles and the soil horizons (A and B) below the litter of the goldmined soils (blue diamond, N = 12). The colored arrows illustrate the isotopic trends in the products of the reaction for abiotic reduction (abiotic red.)55, 56 (red arrow), the foliar uptake processes22, 25, 47 (grey arrow) and the Hg photochemical reduction on foliage and forest floor 16,25 (green arrow). The blue arrow illustrates the isotopic trends in the substrate of the reaction for Hg(0) evaporation (Hg(0) evapo.)15. Blue and gray shaded rectangles refer to available literature data for the mean with 2SD isotopic extent of atmospheric GEM in pristine background areas (δ202Hg= 0.5 ± 0.5 ‰ and ∆199Hg= -0.25 ± 0.15 ‰) 22, 25, 46 , 47 and anthropogenic contaminated areas (δ202Hg= -0.65 ± 0.37 ‰ and ∆199Hg= -0.05 ± 0.15 ‰)24, 46 which is likely similar to re-emission over goldmined soils. Figure 4. Depth profiles for Δ199Hg (left panel) and δ202Hg (right panel) isotopes in pristine and contaminated goldmined soils. Blue dotted lines refer to the mean position of the water table recorded along the year (S.I.6) for acrisol IT1 (~ 65cm depth) and IT2 (~ 75cm depth).
ACS Paragon Plus Environment
Environmental Science & Technology
Undisturbed soil profiles
60
Palm tree leaves
Undisturbed soils Ferralsol IT4 Acrisol IT1 Acrisol IT2 Nodule acrisol IT2 Gleysol IT3 Goldmined soils III-0 SL-3 SL6
40
30
20
Combat creek watershed 10
IT4
400
sandy clayey hydromorphic weathered schist horizon (Ch) IT2
300
IT3
200
100
0
Distance from the flat (m) [THg] (ng g-1)
IT1 0
IT3
200
400
600 0
-20 Depth (cm)
CM
SL6
-40 -60 -80
III-0
-100
0 100m
-120
ACS Paragon Plus Environment
3
6
Log[THg] (ng g-1)
Coarse fraction > 2 mm (%)
[C] (%)
0
IT2
FP
bright-brown clayey compact mineral horizon (B)
Soil horizons Fresh Litter (Oi) Humified Litter (Oe) Organo-mineral (A) Mineral (B) Weathered schist (C) Weathered schist hydromorphic (Ch)
I2
SL3
Reddish-brown macroporous mineral horizon (B)
reddish brown micro-aggregated horizon (B) IT1
50
Elevation (m)
N
litter, root mat (Oi, Oe) and organo-mineral horizon (A)
IT4
Goldmined soil profiles Goldmined flat (Open field)
red compact weathered schist horizon (C)
Page 30 of 33
40
80
0
20
40
60
80
100
1000
10000
Page 31 of 33
Environmental Science & Technology
0.6
Soils
0.4 0.2
Undisturbed soils Ferralsol IT4 Acrisol IT1 Acrisol IT2 Nodule acrisol IT2 Gleysol IT3 Goldmined soils III-0 SL-3 SL6
Precipitation
Ores and
0.0
199 Hg (‰)
Soil litter & palm tree leaves
mine tailings Goldmined soils atmospheric TGM & GEM
-0.2 Young palm tree leaves
-0.4
Soil horizons Fresh Litter (Oi) Humified Litter (Oe) Organo-mineral (A) Mineral (B) Weathered schist (C) Weathered schist hydromorphic (Ch)
Litterfall, foliage and lichens
Old palm tree leaves Undisturbed soil horizons
Palm tree leaves 2SD
SL6 forested ancient mine flat CM over Ferralsol (young) FP edge of the mine flat (young) IT4 over Ferralsol IT1 over Acrisol IT3 over Gleysol CM over Ferralsol FP edge of the mine flat
-0.6
-0.8 -3
-2
-1
0
202 Hg (‰)ACS Paragon Plus Environment
1
Environmental Science & Technology
Page 32 of 33
Palm tree leaves SL6 forested ancient mine flat CM over Ferralsol (young) FP edge of the mine flat (young) IT4 over Ferralsol IT1 over Acrisol IT3 over Gleysol CM over Ferralsol FP edge of the mine flat
'199 Hg (‰) 0.4
R
=
0.9
, p
0.0