Subscriber access provided by ORTA DOGU TEKNIK UNIVERSITESI KUTUPHANESI
Article
Zn isotope fractionation during sorption onto kaolinite Damien Guinoiseau, Alexandre Gelabert, Julien Moureau, Pascale Louvat, and Marc F. Benedetti Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b05347 • Publication Date (Web): 13 Jan 2016 Downloaded from http://pubs.acs.org on January 13, 2016
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 free 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 accessible to all readers and 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.
Environmental Science & Technology 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
1
Zn isotope fractionation during sorption onto
2
kaolinite
3
Damien Guinoiseau, Alexandre Gélabert, Julien Moureau, Pascale Louvat and Marc F.
4
Benedetti
*
5
INSTITUT DE PHYSIQUE DU GLOBE DE PARIS – SORBONNE PARIS CITE –
6
UNIVERSITE PARIS DIDEROT – CNRS UMR 7154, PARIS, FRANCE
7 8
* corresponding author :
[email protected] Abstract
9
In this study, we quantify zinc isotope fractionation during its sorption onto kaolinite, by
10
performing experiments under various pH, ionic strength and total Zn concentrations. A
11
systematic enrichment in heavy Zn isotopes on the surface of kaolinite was measured, with
12
∆66Znadsorbed-solution ranging from 0.11 ‰ at low pH and low ionic strength to 0.49 ‰ at high pH
13
and high ionic strength. Both the measured Zn concentration and its isotopic ratio are correctly
14
described using a thermodynamic sorption model that considers two binding sites: external basal
15
surfaces and edge sites. Based on this modelling approach, two distinct Zn isotopic fractionation
16
factors were calculated: ∆66Znadsorbed-solution = 0.18 ± 0.06 ‰ for ion exchange onto basal sites, and
17
∆66Znadsorbed-solution = 0.49 ± 0.06 ‰ for specific complexation onto edge sites. These two distinct
1 ACS Paragon Plus Environment
Environmental Science & Technology
Page 2 of 33
18
factors indicate that Zn isotope fractionation is dominantly controlled by the chemical
19
composition of the solution (pH, ionic strength).
20
Introduction
21
As bioessential element, Zinc is necessary for growth and reproduction of living organisms1
22
and Zn deficiency in oceanic or continental environments as well as in human diet is of a
23
growing concern. By contrast, toxic levels of Zn are frequently observed in highly polluted sites
24
affected by mining or smelting activities. An accurate understanding of Zn biogeochemical cycle
25
is thus required at the Earth’s surface. Over the last two decades, Zn stable isotope signatures2
26
have proved effective in tracking Zn sources in polluted environments3-8. But in ecosystems like
27
soils or rivers, interactions between Zn in solution and in minerals can lead to Zn isotope
28
fractionation that may hamper Zn source identification. Besides its environmental toxicity, Zn
29
availability or mobility is clearly dependent on Zn speciation, which is difficult to determine. Zn
30
isotope measurements can help identify the ways in which Zn interacts between solution and
31
solids. In the environment, Zn isotope signatures reflect both Zn sources and/or in-situ processes.
32
In order to distinguish between these two factors, isotopic fractionation during Zn complexation
33
by reactive phases within the aquatic environment has to be addressed. Among the three main
34
groups of reactive phases, metal oxides, organic matter, and phyllosilicates, Zn isotopic
35
fractionation factors during adsorption are only known for the first two.
36
Enrichment in heavy Zn isotopes at the surface of iron oxides with respect to coexisting solution
37
was reported, with ∆66Znsolid-solution of 0.29 ± 0.07 ‰9 for goethite and 0.53 ± 0.07 ‰9,
38
ferrihydrite. This systematic enrichment is due to stronger Zn-O bonds at the mineral surface
39
than in solution, as evidenced by EXAFS spectroscopy9. For Zn adsorption on birnessite
10
for
2 ACS Paragon Plus Environment
Page 3 of 33
Environmental Science & Technology
40
(manganese oxides), Zn isotopes show only limited fractionation at low ionic strength, while
41
heavy Zn isotopes are sorbed preferentially at high ionic strength (∆66Znsolid-solution from 0.52 to
42
0.77 ‰)11. The impact of surface loading was also noted, with substantial fractionation occurring
43
at low surface coverage and high ionic strength (∆66Znsolid-solution up to 2.74 ‰)11. Such variable
44
Zn fractionation factors had previously been reported for sorption on Fe, Al and Mn oxides12,
45
with preferential sorption of light Zn isotopes onto goethite and of heavy ones onto birnessite, a
46
result inconsistent with other studies9, 11.
47
Regarding organic matter, Zn binds to the phenolic sites of purified humic acid with a ∆66Znsolid-
48
solution
49
heavy Zn isotope enrichment onto diatom surfaces (∆66Zndiatom-solution = 0.35 ± 0.10 ‰)14, but
50
preferential light Zn isotope internalization by cells (∆66Zncell-solution -0.2 to -0.8 ‰)15. Zn isotope
51
fractionation in higher plants is complex: despite isotope equilibrium between free, organic and
52
inorganic Zn complexes in soil solutions, heavy Zn isotopes rather sorb onto the root surface,
53
while lighter δ66Zn values in the aerial parts result from preferential light Zn isotope transport
54
into stems and leaves16-18.
55
Given the large range of Zn isotopic fractionation associated with these mineral and biological
56
sorption effects, most of the observed δ66Zn range in geological records on Earth (-0.4 to 1.4
57
‰)19 could be explained by such processes. However, Zn isotope fractionation during Zn
58
sorption onto one of the most important mineralogical groups, the phyllosilicates, is still missing
59
and needs to be investigated.
60
In this study, and for the first time, we experimentally determine Zn isotope fractionation during
61
its sorption onto kaolinite, the dominant clay mineral in highly weathered soils, e.g. in tropical
62
and subtropical regions. Cu isotope fractionation during sorption on kaolinite was recently
= 0.24 ± 0.06 ‰13. Additionally, biological activity impacts Zn cycling, with preferential
3 ACS Paragon Plus Environment
Environmental Science & Technology
63 64
reported (∆65Cusolid-solution = -0.29 ‰)20;
65
Page 4 of 33
Cu enrichment being rather attributed to isotopic
fractionation between Cu(H 0) and Cu(H 0) in solution and the preferential sorption of the
65
latter than to water-mineral surface processes. These results contradict previous observations of a
66
systematic enrichment of heavy Cu isotopes on Al and Fe oxide surfaces (∆65Cusolid-solution from
67
0.6 to 1.3 ‰)10,
68
solution.
69
Zn sorption on kaolinite has been extensively studied22-25 and is invariably described by a two-
70
site model involving; (1) an ion exchange reaction at low pH on permanent negatively charged
71
sites located on basal sheets and (2) a specific binding at higher pH on pH-dependent-charge
72
sites located on the clay edges.
73
The objective of this study is twofold: (1) define and correctly model zinc sorption onto kaolinite
74
under multiple pH, ionic strength and zinc concentration conditions, (2) measure the
75
corresponding Zn isotope signatures in the same experimental conditions in both dissolved and
76
particulate phases, to determine Zn isotope fractionation factors during its sorption onto
77
kaolinite.
21
, due to stronger Cu-O bonds formed at the mineral surface compared to
78
Materials and methods
79
Starting material
80
20 grams of Source Clay KGa-2 kaolinite (Warren County, Georgia) from the Clay Minerals
81
Society were decanted in a sodium hexametaphosphate (dispersant) - milli-Q water (mQ)
82
mixture for 16 h, and the fraction < 2 µm was recovered. The residual metal oxides were
83
removed by DCB extraction using Holmgren’s protocol26. The clay was washed several times
84
with mQ water, conditioned under Na form (in NaNO3 solution) and washed again three times to 4 ACS Paragon Plus Environment
Page 5 of 33
Environmental Science & Technology
85
remove the salt excess. The clay was then freeze-dried, and stored at room temperature. Based on
86
XRD measurements treated with Rietveld refinement, the starting sample is composed of around
87
99 % kaolinite and 1 % anatase (TiO2). The specific surface area (SSA) determined by N2-BET27
88
measurement is equal to 20.7 m2/g, similar to a previous study23.
89
All reagents were prepared from distilled acids, Titrinorm bases and mQ water to avoid Zn
90
contamination. Zn stock solution used in all experiments was a 1000 ppm AAS specpure solution
91
(Alfa Aesar) in 5 % HNO3.
92
Sorption experiments
93
Three sets of experiments were performed:
94
Experiment set 1: acid-base titrations of kaolinite at 20 g/L in 0.1 M, 0.01 M and 0.005 M
95
NaNO3 (background electrolyte) were conducted to determine protonation parameters and
96
reactive site densities. Forward and back titrations were performed twice for each experiment
97
with a Titrando unit (Metrohm) under continuous N2 flux (details in SI.1)
98
Experiment set 2: sorption edge experiments at constant initial [Zn] and variable pH were
99
conducted in batch in a glove box (Jacomex) under a N2 atmosphere to avoid the presence of
100
carbonate species and minimize any possible Zn isotopic fractionation between Zn2+ and ZnCO3
101
aqueous species. A suspension of kaolinite at 5 g/L was equilibrated overnight in a N2-saturated-
102
NaNO3 solution (0.1 M in a first series and 0.01 M in a second one) prior to the addition of a Zn
103
solution (final [Zn] = 50 ± 3 µM). In each tube, the fixed pH value (between 3 and 9) was
104
adjusted daily by addition of small volumes of 1N HNO3 and/or 1N NaOH. Experiment
105
equilibration time was 48 h, enough to reach sorption equilibrium and short enough to avoid clay
106
dissolution and/or precipitation of Zn-Al-Layer-Double-Hydroxide (LDH)28. After the final pH
5 ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 33
107
was measured, the suspensions were centrifuged, supernatants filtered at 0.22 µm and decanted
108
kaolinite samples freeze-dried.
109
Experiment set 3: Sorption isotherms at constant pH and variable total [Zn] were performed
110
under conditions similar to set 2. A suspension of kaolinite at 5 g/L was equilibrated overnight in
111
0.01 M N2-saturated-NaNO3. Total Zn concentration in each batch ranged from 5 to 1500 µM.
112
pH was fixed at 4.00 ± 0.10 for a first series, and at 6.00 ± 0.10 for a second one. After two days,
113
samples were treated as in set 2.
114
In sets 2 and 3, sample solutions were acidified to 0.5 M HNO3 and analyzed for Zn and Al
115
concentrations with ICP-OES (Thermo, ICAP 6000 series). As Al concentration remained lower
116
than 4 µM, kaolinite dissolution was therefore negligible, even at acidic pH. Suspended fractions
117
were digested with HNO3-HF at 100°C and Zn concentrations were determined with ICP-OES.
118
Zn isotopes analysis
119
Prior to Zn isotope analysis, sample solutions were evaporated and kaolinite residues were
120
digested. In a HEPA 13 air-filtered and over-pressurized cleanroom, samples were then dissolved
121
in HCl 6 M for chromatographic elution. Zn was isolated following a protocol adapted from
122
Marechal et al.2 using polypropylene columns (0.8 I.D. x 4 cm length) and AG-MP1 200-400
123
mesh resin (Biorad). Yields were measured after Zn extraction and were always in the range 100
124
± 5 %. Total procedural blanks were 7 ± 4 ng Zn (n=13), which represented less than 1 % of the
125
sample Zn amount. This low blank and the equilibrated isotopic mass balance between Zn in
126
solution and in kaolinite for almost all samples (Fig. S7) means that the procedural blank had a
127
negligible impact on measured Zn isotope ratio measurement, contrary to other experiments11. A
128
Cu Alfa Aesar stock solution (65/63Cu = 0.44574)29 was added to each sample as an internal
129
spike, with a Zn/Cu ratio of 2, in order to correct for instrumental mass bias (details in SI.2).
6 ACS Paragon Plus Environment
Page 7 of 33
Environmental Science & Technology
130
Samples were analyzed on a MC-ICP-MS Neptune Plus (ThermoFinnigan) using a sample-
131
standard bracketing technique (SSB)30 with the Zn and Cu Alfa Aesar solutions (66/64Zn =
132
0.565055 ± 0.000004 (2σ) and
133
referred to in the text as ZnIPG or CuIPG, respectively. Zn isotope composition was expressed as
65/63
Cu = 0.44574 ± 0.000004 (2σ))29 as in-house standards,
Zn Zn % δ
Zn =
− 1$ ∗ 1000 (1) Zn Zn #
134 135
As samples obeyed a mass dependent fractionation law (Fig. S2), only δ66Zn is reported in this
136
study. The in-house ZnIPG solution was calibrated with respect to IRMM 3702 and JMC Lyon2
137
international standards, giving δ66ZnIPG/IRMM = -0.23 ± 0.04 ‰ (2σ, n=82) and δ66ZnIPG/JMC =
138
0.04 ± 0.03 ‰ (2σ, n=55), respectively. Additionally, independent cross calibration of ZnIRMM
139
over ZnJMC gave δ66ZnIRMM/JMC = 0.28 ± 0.03 ‰ (2σ, n=14) in the range of previously reported
140
signatures4, 31-33. Repeated Zn chemical purification from the certified BCR-2 basalt gave δ66Zn
141
of 0.24 ± 0.03 ‰ (2σ, n=11) relative to JMC, in excellent agreement with former published
142
values34,
143
solution. For kaolinite samples, measured δ66Znsolid was corrected from the structural Zn initially
144
present in the mineral lattice ([Zn]struct..= 37 ± 5 ppm (n=2) and δ66Znstruct. = 0.47 ± 0.02 ‰ (n=2)
145
determined after acid digestion) following:
35
. All sample δ66Zn values are expressed relative to our in-house Zn IPG standard
δ Zn'()*' =
+Zn,(-' ∗ δ
Zn(-' − +Zn,.)/0.. ∗ δ
Zn.)/0.. +Zn,'()*'
(2)
146
The isotopic fractionation between Zn adsorbed on kaolinite and Zn in solution, ∆66Znadsorbed-
147
solution was
defined as:
7 ACS Paragon Plus Environment
Environmental Science & Technology
Page 8 of 33
∆
Zn'()*'4(/.-( = δ
Zn'()*' − δ
Zn(/.-( (3) 148
Two standard deviations (2σ) reported in table S1 and in figures were calculated from the three
149
replicates of each sample δ66Zn measurement. If the resulting 2σ was lower than 0.04 ‰, the
150
long-term reproducibility of IRMM 3702 (0.04 ‰, 2σ) was reported instead.
151
Sorption modeling
152
Previous studies22,
23, 25, 36
revealed that trace metal sorption onto kaolinite occurs at two
153
different types of sites. The first are permanent negatively charged sites, X- sites, allowing ion
154
exchange with the solution through electrostatic interactions. They are located at the surface of
155
clay sheets and result from substitution in Al or Si sheets. The second ones correspond to pH-
156
dependent edge sites, S-OH0.5-sites that favor specific Zn-solid complexes.
157
The sorption of cations (H+, Na+ and Zn2+) was modeled with ECOSAT v.4.7 (2001)37 software
158
using Gaines-Thomas ion exchange formalism for basal site binding. According to previous
159
models22, 24, 25, the X- sites are fully saturated with Na+, and Zn2+ is exchanged with a Zn:Na ratio
160
of 0.5. The input equations are:
XNa + H ⇋ XH + Na logK ?@ (4)
2XNa + Zn ⇋ X Zn + 2Na logK ?B (5) 161
1-pK Triple Layer complexation Model (1-pK TLM) was chosen for edge site sorption with a
162
single deprotonation step rather than the 2-pK TLM model because of the reduced number of
163
adjustable parameters38. Heidmann et al23 also used a unique deprotonation step but with a 1-pK
164
Basic Stern Model. Edge sites are expressed as SOH0.5- in their deprotonated form due to the
165
incomplete neutralization of O2- charge by the octahedrally coordinated Al3+ comprising the
166
dominant aluminol edge sites. AFM studies39,
167
represents only 30 % of the KGa-2 kaolinite. Furthermore, Zn complexation is expected to occur
40
reported that the edge site surface area
8 ACS Paragon Plus Environment
Page 9 of 33
Environmental Science & Technology
168
mainly as an inner-sphere complex (ISC), in agreement with EXAFS study28 results. The binding
169
model also takes into account the possible effect of Na+ and NO3- competition on these sites: S − OH F.4 + H ⇋ S − OH F. logK GHIB J.KL (6) S − OH F.4 + Na ⇋ S − OH − NaF. logK GHI4@J.KL (7) S − OH F.4 + H + N0O 4 ⇋ S − OH − N0O F.4 logK
GHIB 4@FP J.KQ
(8)
S − OH F.4 + Zn ⇋ S − OHZnS. logK GHIT.KL (9) 170
In the 1-pK TLM model, ISC compensations of H+ and Zn2+ charges are assumed onto the 0-
171
plane at the kaolinite surface. The non-specific ions, Na+ and NO3-, are compensated on the
172
inner-diffuse-layer-plane instead.
173
For the model adjustment, first, the number of sites, capacitance, protonation/deprotonation
174
constants and influence of NaNO3 electrolyte were fitted, using acid/base titrations (experiment
175
set 1, fig. S1). Then, Zn binding constants were estimated from sorption edge and isotherm
176
experiments (sets 2 and 3, Figure 1 and Table 1). Model parameters were adjusted to obtain the
177
best R2 value between measured and modeled [Zn] (Fig. S3, R2 = 0.98 and p < 0.00001).
178
Results
179
Titration and sorption experiments
180
Potentiometric titrations: Acid-base titrations of KGa-2 in the absence of Zn were performed at
181
0.005, 0.01 and 0.1 M NaNO3. The adsorbed H+ density is lower at high ionic strength, due to
9 ACS Paragon Plus Environment
Environmental Science & Technology
Page 10 of 33
182
stronger competition between Na+ and H+. A common intersection point around pH = 5.3 defines
183
the point of zero charge and agrees with published data22.
184
Sorption edges: The two Zn sorption edges performed at pH from 3 to 9 in 0.01 and 0.1 M
185
NaNO3 (Figure 1A and 1B respectively) differ from each other. At low ionic strength (Figure
186
1A), with lower competition between Zn2+ and Na+, Zn sorption begins below pH 3 and is
187
complete by pH 7.5. At high ionic strength (Figure 1B), due to strong competition with Na+ on
188
negatively charged basal sites, Zn sorption is lower than 0.05 µmol Zn/m2 below pH 6 and
189
increases to reach complete sorption with 0.49 µmol Zn/m2 at pH 7.5.
190
Sorption isotherms: At pH 4, Zn sorption increases from 0.02 µmol Zn/m2, to a plateau at 0.48
191
µmol Zn/m2 when free [Zn2+] increases from 3 to 1400 µM (Figure 1C). At pH 6, (Figure 1D),
192
[Zn]adsorbed increases from 0.04 to 0.96 µmol Zn/m2, with free [Zn2+] raising from 3 to 1500 µM.
193
Conversely to the pH 4 isotherm, no plateau is reached at pH 6.
194 195
Zn isotope ratio changes during sorption
196
Sorption edges: Five and six samples from the 0.1 M and 0.01 M sorption edge experiments,
197
respectively, were analyzed (supernatant and solid fraction) (Figure 2). Zn isotopes behave in a
198
similar manner at both ionic strengths: continual enrichment in heavy isotopes in the solid phase
199
and in light isotopes in the solution when pH increases. Mass balance calculations for almost all
200
samples yield an isotope composition indistinguishable from that of the starting solution (ie. 0.00
201
± 0.04 ‰, fig. S7).
202
At high ionic strength (Figure 2A), δ66Znsolution evolves towards a lighter isotopic signature, from
203
0.00 ‰ when Zn sorption is minimal to -0.49 ‰ at 95 % of adsorbed Zn. Simultaneously,
10 ACS Paragon Plus Environment
Page 11 of 33
Environmental Science & Technology
204
positive δ66Znadsorbed on kaolinite varies from 0.31 ‰ at low coverage to 0.04 ‰ at 95 %
205
sorption, with a maximum of 0.35 ‰ for approximately 25 % of sorption.
206
At low ionic strength (Figure 2B), δ66Znadsorbed decreases almost continuously from 0.18 to 0.07
207
‰. δ66Znsolution is slightly negative (0.00 to -0.07 ‰) up to 50 % sorption then sharply decreases
208
up to 100 % sorption (-0.07 to -0.41 ‰). At 30 % of sorption, the low δ66Znadsorbed (0.05 ‰) is
209
considered as an outlier, as will be discussed later.
210
For each experiment, the corresponding ∆66Znadsorbed-solution is high and evolves with Zn surface
211
coverage from 0.32 to 0.52 ‰ at 0.1 M NaNO3, and from 0.18 to 0.49 ‰ at 0.01 M NaNO3
212
(Figure 2 and table S1). The difference of ∆66Znadsorbed-solution between the two ionic strengths
213
indicates that at least two different sorption processes implying that two different mechanisms
214
driving Zn isotope fractionation are operating.
215
Sorption isotherms: Five and six samples were analyzed from the experiments at pH 4 and 6,
216
respectively (both the supernatant and the solid fraction). As for sorption edges, heavy Zn
217
isotopes are preferentially adsorbed on kaolinite surfaces compared to solution (Figure 3).
218
Isotopic mass balances for solid and solution give the Zn isotope signature of the stock solution
219
used (0.00 ± 0.04 ‰, Fig. S7), except for the first point of each isotherm where δ66Znadsorbed is
220
the highest.
221
At pH 4 (Figure 3A), δ66Znsolution (0.02 ± 0.02 ‰) remains close to the initial δ66Znstock-solution,
222
while δ66Znadsorbed increases from 0.10 to 0.34 ‰. At 500 µM Zn2+eq., as for sorption edge, an
223
outlier is observed with unexpectedly low δ66Znadsorbed. Excluding the high ∆66Znadsorbed-solution
224
(0.32 ‰) determined for the first isotherm point at low [Zn2+], the ∆66Znadsorbed-solution remains
225
stable at 0.11 ± 0.01 ‰ for [Zn2+] between 160 to 980 µM. 11 ACS Paragon Plus Environment
Environmental Science & Technology
Page 12 of 33
226
At pH 6 (Figure 3B), δ66Znsolution is enriched in light isotopes for low [Zn]eq. compared to
227
δ66Znstock-solution. Adsorbed Zn is enriched in heavy isotopes (from 0.12 to 0.22 ‰); δ66Znadsorbed
228
first decreases for [Zn] between 0 and 120 µM down to 0.12 ‰ and then increases up to 0.2 ‰.
229
The corresponding ∆66Znadsorbed-solution ranges from 0.17 to 0.28 ‰ with an almost constant value
230
of 0.18 ‰ for Zn addition between 96 and 394 µM (Table S1).
231
Discussion
232
Evolution of Zn-kaolinite binding with pH: a two-site model
233
Results from sorption of H+ and Zn2+ were well fitted using our two-site model. Adjusted
234
parameters for the model are given in Table 1. The concentration of SOH0.5- sites (171.9
235
mmol/kg), representing 16.6 sites/nm2 of edge surface area, is nine times that of X- sites (19
236
mmol/kg). The former value is slightly higher than the 12.2 sites/nm2 previously reported for the
237
same KGa-2 clay23. The high uncertainties concerning edge surface area determinations for
238
KGa-2 and other poorly ordered kaolinites (between 12 and over 30 %)39-41 may explain this
239
small discrepancy between the two studies. The low amount of X- sites, which arise from rare
240
substitutions in the kaolinite lattice, is slightly higher than in previous estimations (13.6
241
mmol/kg)42. The complexation constants logK ?B (Eq. (5)) and logK GHIT.KL (Eq. (9)), are
242
fitted at 0.8 and 1.2 (Table 1). For non-specific ionic exchange, this logK ?B value agrees with
243
previous determination using a similar formalism22,
244
protonation for edge sites22,
245
single protonation step for SOH0.5- sites. However, Heidmann et al.23 assumed a bidendate
246
complex involving the hydrolyzed-metal forms of Cu and Pb. Such Zn hydrolyzed species were
247
implemented in our model for Zn sorption but did not improve the fit.
24, 25
23
. For ISC, all studies used a two-step-
, except for Heidmann et al.23 and in this study who adopt a
12 ACS Paragon Plus Environment
Page 13 of 33
Environmental Science & Technology
248
According to our model, at low pH, Zn is bound exclusively onto pH-independent basal sites
249
(Figure 1) and the strong influence of ionic strength during sorption edge experiments at this pH
250
(3 to 5) indicates an outer-sphere complexation (OSC). At higher pH, binding to edge sites
251
becomes dominant (Figure 1B), and the negligible effect of Na+ and NO3- concentrations on Zn
252
sorption suggests that zinc is adsorbed as an ISC. This corroborates the findings of Nachtegaal et
253
al.28 who find a Zn monodendate ISC on kaolinite at pH 5 and 1 mM [Zn]tot. During Zn sorption,
254
precipitation of Zn-Al LDH was observed for kaolinite28, pyrophyllite43 or montmorillonite44 by
255
EXAFS spectroscopy. The degree of precipitation is proportional to [Zn], pH and aging time.
256
Our low Zn concentrations (50µM) and equilibration time (2 days) for sorption edges mean that
257
LDH precipitations is not widespread. However, higher [Zn] were introduced for the final points
258
of the sorption isotherms. The local environment of Zn once sorbed on kaolinite at pH 4, 6 and 8,
259
with a [Zn] of 800 µM but an aging time of 29 days has been checked by X-ray absorption
260
spectroscopy. Preliminary data treatment (not shown here) seems to indicate an absence of Zn-
261
Al-LDH precipitates at pH 4, traces of LDH at pH 6 and a significant amount of Zn into the
262
newly-formed Zn-Al-LDH phase at pH 8. The presence of Zn-Al-LDH after only 2 days of
263
reaction, in our pH 4 and 6 sorption isotherms is thus highly improbable.
264
As in the sorption edge experiments, our model faithfully reproduces the Zn sorption isotherms.
265
All Zn is complexed onto basal sites at pH 4 because edge sites are protonated. At pH 6, this
266
ionic exchange also dominates at low Zn concentrations until saturation. The specific binding
267
onto edge sites prevails at high [Zn] (>0.7 mM). The successive involvement of the different
268
binding sites was also previously reported for Cu sorption onto kaolinite23.
269
Zinc isotope fractionation during binding
13 ACS Paragon Plus Environment
Environmental Science & Technology
Page 14 of 33
270
δ66Zn measurements in the sorption edge experiments (Figure 2A and 2B) reveal distinct Zn
271
isotope behavior according to the ionic strength. In both cases, enrichment of heavy Zn isotopes
272
onto the kaolinite surface is measured, with ∆66Znadsorbed-solution between 0.18 and 0.52 ‰. The
273
two patterns observed (Figure 2A and 2B) are explained by an equilibrium fractionation (a
274
Rayleigh fractionation process is rejected on the basis of details in Fig. S4) with two distinct
275
fractionation factors: one for low coverage when Zn sorption occurs predominantly onto
276
exchange sites, and one for higher coverage during Zn sorption onto edge sites. For each
277
coverage value, the adsorbed δ66Zn signature results from the relative proportions of each
278
binding site type. The whole δ66Znadsorbed vs. Zn coverage curves allow for the determination of
279
the Zn isotope fractionation factor of each site.
280
In the sorption edge experiment at high pH (pH > 6.5) and high ionic strength (Figure 1B), 50
281
to 100 % of Zn is adsorbed, an average ∆66Znadsorbed-solution of 0.49 ‰ is observed (Figure 2A) and
282
only SOH0.5- sites contribute to the binding. Therefore, the isotope fractionation factor for Zn
283
sorption onto edge sites is ∆GHIT.KL = 0.49 ± 0.06 ‰. For sorption edge at low pH (pH < 5.5)
284
and low ionic strength (0.01 M NaNO3, Figure 1A), 20 to 40 % of zinc is adsorbed, a
285
∆66Znadsorbed-solution of 0.18 ‰ is observed for the first point (Figure 2B), and the model predicts
286
an ionic exchange binding exclusively on basal X- sites (Figure 1A). The isotope fractionation
287
factor for Zn sorption onto exchange sites is thus ∆?B = 0.18 ± 0.06 ‰.
288
To confirm that these fractionation factors correctly explain the δ66Znadsorbed and δ66Znsolution
289
when both binding sites are involved during sorption, we calculate the theoretical δ
Zn.X( (/.-( .
290
The calculation is based on the proportion of Zn adsorbed (f), the proportion of each site (pX2Zn
291
and pSOHZn1.5+), and their associated fractionation coefficient (Y?B and YGHIT.KL ):
14 ACS Paragon Plus Environment
Page 15 of 33
Environmental Science & Technology
δ
Zn.X( (/.-( = − 1000f
292
p? ]α?B − 1_ pGHIS. (αGHIT.KL − 1) + (10) 1 − f + fα?B 1 − f + fαGHIT.KL
The theoretical δ
Zn.X( '()*' is obtained by:
.X( δ
Zn.X( '()*' = δ Zn'-. + p?B ∗ ∆?B + pGHIT.KL ∗ ∆GHIT.KL (11)
293
Theoretical calculations are reported in Figure 2 for sorption edges and show the δ
Zn.X( (/.-(
294
(red dotted line) and the δ
Zn.X( '()*' (red solid line). Additionally, the blue dotted line
295
separates the proportion of Zn linked to basal sites (green area) from those linked to edge sites
296
.X( (blue area). δ
Zn.X( (/.-( and δ Zn'()*' values are in good agreement with measured
297
δ66Znsolution and δ66Znadsorbed in sorption edges, validating: (1) the dominance of edge sites at high
298
pH or high ionic strength (Figure 2A), (2) the control by basal sites at low pH and low ionic
299
strength (Figure 2B), and (3) intermediate proportions in intermediate conditions. This good
300
agreement between our two-site complexation model used to describe [Zn] evolution and its
301
robustness to describe zinc isotope fractionation also validates the two zinc isotope fractionation
302
factors we propose for Zn sorption on kaolinite edge and basal sites.
303
A similar calculation is done for sorption isotherms at pH 4 and 6 (Figure 3). Even if the model
304
validates the enrichment of heavy isotopes onto the kaolinite surface, the δ
Zn.X( '()*' (red
305
solid line in Fig. 3) slightly overestimates measured δ66Znadsorbed, except for the lowest [Zn]eq.
306
where δ66Znadsorbed is underestimated. For these low [Zn] experiments and in the first point of
307
Figure 2A, the isotopic mass balances are not fulfilled (Fig.S7). The amount of Zn adsorbed is
308
very low and close to that initially present into the crystal lattice of kaolinite. The uncertainty on
309
[Zn]struct. within the kaolinite lattice (37 ± 5 ppm) combined with the correction of this crystalline
310
Zn pool (δ66Znstruct. = 0.47 ± 0.02 ‰) in the δ66Znsolid measured can explain this difference
15 ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 33
311
(details in Fig. S5). For higher [Zn], the uncertainty in crystalline Zn cannot be responsible for
312
the observed offset, however other issues may arise:
313
1) Suspended kaolinite was separated from the solution by centrifugation, however some
314
solution may remain in the kaolinite pellet porosity and contribute to lower the measured
315
δ66Znadsorbed. This hypothesis is applicable for the isotherm at pH 4 (Figure 3A), where dissolved
316
[Zn] is much higher than adsorbed [Zn]. For instance, 200 µL of solution would be enough to
317
shift the adsorbed signature by 0.1 ‰. But at pH 6, the volume of interstitial solution needed to
318
decrease δ66Znadsorbed is too large (> 2 mL) to account for the discrepancy. The same issue was
319
reported by Balistrieri et al.10 and Juillot et al.
320
dissolved Zn contamination. However, this may overprint previous fractionation by attempting to
321
establish a new equilibrium with the washing solution.
322
2) For the isotherm at pH 6 (Figure 3B), a slight overestimation of the SOH0.5- proportion (blue
323
9
who tested washing their oxides to remove
area in fig 3B) may be the cause for a δ
Zn.X( '()*' higher than measured. A 10 %
324
overestimation of the SOH0.5- proportion would shift the δ
Zn.X( '()*' by 0.05 ‰. Thus,
325
uncertainties on edge and basal sites proportions may explain most of the discrepancy reported
326
between theoretical and measured δ66Znadsorbed.
327
3) A change in sorbed Zn speciation at high surface coverage (high [Zn]) cannot be excluded.
328
For birnessite, a progressive change from tetrahedral to tetrahedral/octahedral coordination of
329
adsorbed complexes occurred at high [Zn]11. For kaolinite, such a modification from octahedral
330
to tetrahedral coordination at high [Zn] could have happened and would result, by strengthening
331
of Zn complex bonds, in an increase of δ66Znadsorbed values, which is not observed here. However,
332
Zn could alternatively be adsorbed on a third type of binding sites at high [Zn]. 16 ACS Paragon Plus Environment
Page 17 of 33
Environmental Science & Technology
333
Zinc speciation impact on Zn isotope fractionation processes
334
The preferential sorption of heavier Zn isotopes on kaolinite relative to aqueous solution can be
335
reconciled with molecular complexation mechanisms. Two main processes can cause such an
336
equilibrium fractionation: (1) isotopic partitioning between Zn aqueous species in solution, and
337
(2) fractionation between dissolved and adsorbed Zn species at the mineral-solution interface.
338
The first step is to constrain the multiplicity of Zn species in solution (1). To this end, aqueous
339
Zn speciation calculated with ECOSAT v.4.7 (2001)37 for sorption edge solutions with pH < 8
340 341 342 343 344
shows that Zn(H 0)
is the major species, with a very small contribution of ZnNOO (H 0) and
of ZnOH(H 0) , whereas Zn(OH) (H 0) is negligible below pH 7 (Fig. S6). The formation of ZnNOO (H 0) and ZnOH(H 0) may induce isotope fractionation among Zn species in
solution. Several studies reported ab initio isotopic partition coefficients between Zn aqueous
species in solution45, 46, but those for ZnNOO (H 0) and ZnOH(H 0) are currently unknown.
345
Extension of such calculations to these species would be useful to see if such processes could
346
account for the measured δ66Znadsorbed.
347 348
66 Independent of the proportion of ZnNOO (H 0) , both [Zn] and δ Zn are well modeled by
sorption of Zn(H 0)
. ZnNOO (H 0) is not likely to be adsorbed on kaolinite edge sites.
349
ZnOH(H 0) is not needed to describe sorption experiments (see section Sorption modeling).
350
However, its appearance at around pH 6.0 (Fig. S6A & B) coincides with the steep increase of
351
Zn sorption on edge sites in Figure 1A and 1B. Under our experimental conditions, we cannot
352
completely rule out a possible sorption of this species on kaolinite. But its impact on Zn isotope
353
fractionation in solution is uncertain since its partition coefficient is unknown. Spectroscopic
17 ACS Paragon Plus Environment
Environmental Science & Technology
Page 18 of 33
354
studies of the kaolinite surface have thus far been unable to identify the coordination of Zn
355
during edge site binding.
356
Adsorption of aqueous species onto the surface of minerals is frequently associated with bond
357
stiffening owing to the different complexing environment, and should thus also engender an
358
isotopic fractionation. Stable isotope theory47,
359
found in stronger bonding environments. All else being equal, low coordination number and/or
360
stronger-field ligands concentrate heavy isotopes. During sorption onto iron oxides (goethite and
361
ferrihydrite) enrichment in Zn heavy isotopes at the mineral surface results from stronger surface
362
complexes compared to dissolved Zn9. For Zn binding by organic matter, Jouvin et al.13 assumed
363
that carboxylate and phenolate complexes act as analogues of the low and high affinity sites of
364
humic acid, respectively. The length difference (2.00 Å for carboxylate and 1.91 Å for
365
phenolate) explains the lack of isotope fractionation onto low affinity sites, and the significant
366
fractionation for high affinity sites. Previous sorption studies14-18 onto diatoms or roots also
367
underlined this trend, with heavy isotope enrichment at diatom and root surfaces due to a change
368
in Zn coordination; from octahedral in solution to tetrahedral on the surfaces.
369
Our results confirm a higher intensity of isotopic fractionation during Zn sorption onto specific
370
edge sites as ISC (∆GHIT.KL = 0.49 ‰), compared to ionic exchange (∆?B = 0.18 ‰) where
371
OSC is foreseen. First, Zn sorption on edge sites (Figure 2A) induces an isotope fractionation
372
close to that resulting from zinc sorption on ferrihydrite surface (∆ = 0.53 ± 0.07 ‰)9. For
373
ferrihydrite, EXAFS data show a Zn-O bond length of 1.96 Å and a tetrahedral complex. A
374
longer Zn-O bond (2.06 Å) and an octahedral conformation for aqueous zinc in solution is
375
observed49, inducing a stronger Zn-O bond for ferrihydrite compared to Zn in solution9. In
376
kaolinite, Zn sorption on edge sites was previously observed on a unique sample as a
48
predicts that heavy isotopes are preferentially
18 ACS Paragon Plus Environment
Page 19 of 33
Environmental Science & Technology
377
monodendate inner-sphere complex by X-ray absorption spectroscopy28. An octahedral complex
378
at the surface is suggested by a 6 fold coordination geometry28, and the measured Zn-O distance
379
(2.07 Å) is similar to that reported for Zn(H 0)
. Thus, bond length and coordination number
380
of Zn in solution seem to be preserved during sorption, however a significant zinc fractionation
381
occurs in our kaolinite experiment (∆GHIT.KL = 0.49 ± 0.06 ‰). Since Nachtegaal et al.28 study
382
reported only one measurement and at a much higher [Zn] than in our study, additional EXAFS
383
spectroscopic studies should be performed to constrain the molecular mechanism responsible for
384
this ∆GHIT.KL fractionation.
385
Moreover, as shown by the observed competition between Na+ and Zn2+ at high ionic strength,
386
the basal site complexation of Zn would correspond to an OSC, where Zn remains in a six-fold
387
coordination geometry with its hydration shell and would exhibit bond lengths similar to the
388
aqueous species, as was recorded for Cu50, 51. Since the mechanism leading to heavy Zn isotope
389
enrichment during OSC complexation formation (∆?B = 0.18 ± 0.06 ‰) remains unclear,
390
additional spectroscopic studies should again be performed to identify it. Nevertheless, several
391
possibilities may be proposed, such as a symmetry modification or a distortion of Zn hydration
392
sphere during its sorption, as was reported for Se complexing with maghemite52.
393
Environmental implications
394
This study adds a new piece to the complex puzzle of Zn isotopic cycle at the Earth’s surface.
395
Even if simplified conditions were used in our experiments (absence of carbonate and a single
396
mineral phase), their results may be applied to natural systems such as soils where Zn availability
397
is high and directly linked to the clay and/or organic matter content given their high CEC
398
value53. Except in certain Zn deposits like smelter slags, most of the pedological δ66Zn values of 19 ACS Paragon Plus Environment
Environmental Science & Technology
Page 20 of 33
399
both natural and anthropogenic origins span a range (-0.33 to 0.49 ‰)54 as large as that reported
400
in this study for Zn sorbed onto kaolinite surface (0.00 to 0.35 ‰). The identification of Zn
401
sources can be efficient, and particularly strengthened by combining Zn isotope tracer with
402
spectroscopic observations5. However, some δ66Zn fluctuations observed in mineral horizons55, 56
403
or sub-surficial layers5 of soils still remain misunderstood and may be attributed to in-situ
404
fractionating processes, such as Zn sorption on kaolinite clay. This observation is particularly
405
relevant in systems like in tropical soils, where the weathering reactions are intense and clays are
406
abundant.
407
Associated content: Additional information is presented as detailed in the text. This material is
408
available free of charge via the Internet at http://pubs.acs.org.
409
Author information:
410
Corresponding author: e-mail:
[email protected] Equipe de Géochimie des Eaux – IPGP-
411
SORBONNE PARIS CITE – UNIVERSITE PARIS DIDEROT – CNRS UMR 7154 – 1 rue
412
Jussieu, 75238 Paris Cedex 05
413
Acknowledgements: Authors greatly thank Jean-Yves Piquemal for BET measurements, Laure
414
Cordier for technical assistance with ICP-OES analyses and Sophie Nowak for XRD acquisition.
415
Farid Juillot and Marc Blanchard are thanked for fruitful discussions. The manuscript was deeply
416
improved by comments of P. Sossi, four anonymous reviewers and D. Giammar. This work was
417
partly funded by EC2CO FIZCAMO project (CNRS-INSU) and is part of CAPES-COFECUB
418
project n°713-2011.
20 ACS Paragon Plus Environment
Page 21 of 33
Environmental Science & Technology
419 420
Figure 1: Zinc sorption evolution with pH for sorption edge experiments at 0.01 M NaNO3 (A)
421
and 0.1 M NaNO3 (B). Zn sorption evolution with Zn concentration in solution at 0.01 M NaNO3
422
for sorption isotherms at pH 4 (C) and pH 6 (D). Experimental points (red dots) are well fitted
423
(black straight line) with various proportions of ionic exchange (green dotted line), and specific
424
edge binding sites (blue dotted line). Samples analyzed for their Zn isotope composition are
425
marked by black squares.
21 ACS Paragon Plus Environment
Environmental Science & Technology
Page 22 of 33
426 427
Figure 2: Zinc isotope compositions for sorption edge experiments ([Zn]ini.=50±3µM and varying
428
pH) plotted against the amount of Zn adsorbed onto kaolinite at high (A) and low (B) ionic
429
strength. White and black dots refer to dissolved and adsorbed δ66Zn, respectively, with an initial
430
δ66Znstock-solution at 0.00 ± 0.04 ‰. Red solid and dotted lines depict theoretical δ66Znadsorbed and
431
δ66Znsolution obtained from our two-site model. The blue dotted lines separate the proportion of 22 ACS Paragon Plus Environment
Page 23 of 33
Environmental Science & Technology
432
∆66Znadsorbed-solution explained by a Zn sorption onto exchange basal sites (green area) or onto edge
433
sites (blue area).
434 435
Figure 3: Zinc isotope composition as a function of [Zn]aq. added at pH 4 (A) and pH 6 (B).
436
White and black dots refer to dissolved and adsorbed δ66Zn, respectively, with an initial 23 ACS Paragon Plus Environment
Environmental Science & Technology
Page 24 of 33
437
δ66Znstock-solution at 0.00 ± 0.04 ‰. Red solid and dotted lines depict theoretical δ66Znadsorbed and
438
δ66Znsolution obtained from our two-site model. The blue dotted lines separate the proportion of
439
∆66Znadsorbed-solution explained by a Zn sorption onto exchange basal sites (green area) or onto edge
440
sites (blue area).
441
Table 1: Fitted parameters of the two-site model used to describe Zn evolution. Site density and
442
protonation constant for edge sites were derived from acid base titrations of kaolinite. Zn
443
complexation constants in both binding sites are fitted according to sorption experiments. General features [kaol.] (g/L) Specific Surface Area (m2/g) Site
-
SOH0.5
5 or 20 20.7 Parameter on Triple Layer Model Site density (sites/nm2) Inner capacitance (κ) (F/m2) Outer capacitance (κ) (F/m2) log KSOH20.5+ log KSOH-Zn1.5+
16.6 0.4 5.0 5.3 1.2
log KSOH-Na0.5+
-1.1
0.5-
log KSOH2-NO3
XNa
444
4.1
Parameter on Gaines-Thomas ionic exchange Site concentration (mmol/kg) 18.0 log KXH 1.3 log KX2-Zn
0.8
445
446
447
448
24 ACS Paragon Plus Environment
Page 25 of 33
449
Environmental Science & Technology
TOC/Abstract Art
450 451 452 453 454 455 456 457 458 459 460 461
25 ACS Paragon Plus Environment
Environmental Science & Technology
462 463 464
1.
Page 26 of 33
Frassinetti, S.; Bronzetti, G. L.; Caltavuturo, L.; Cini, M.; Croce, C. D., The role of zinc
in life: A review. 2006, 25, (3), 597-610. 2.
Maréchal, C. N.; Télouk, P.; Albarède, F., Precise analysis of copper and zinc isotopic
465
compositions by plasma-source mass spectrometry. Chemical Geology 1999, 156, (1–4), 251-
466
273.
467
3.
Borrok, D. M.; Wanty, R. B.; Ian Ridley, W.; Lamothe, P. J.; Kimball, B. A.; Verplanck,
468
P. L.; Runkel, R. L., Application of iron and zinc isotopes to track the sources and mechanisms
469
of metal loading in a mountain watershed. Applied Geochemistry 2009, 24, (7), 1270-1277.
470
4.
Cloquet, C.; Carignan, J.; Libourel, G., Isotopic composition of Zn and Pb atmospheric
471
depositions in an urban/periurban area of Northeastern France. Environmental Science &
472
Technology 2006, 40, (21), 6594-6600.
473
5.
Juillot, F.; Maréchal, C.; Morin, G.; Jouvin, D.; Cacaly, S.; Telouk, P.; Benedetti, M. F.;
474
Ildefonse, P.; Sutton, S.; Guyot, F.; Brown Jr, G. E., Contrasting isotopic signatures between
475
anthropogenic and geogenic Zn and evidence for post-depositional fractionation processes in
476
smelter-impacted soils from Northern France. Geochimica et Cosmochimica Acta 2011, 75, (9),
477
2295-2308.
478
6.
Sivry, Y.; Riotte, J.; Sonke, J. E.; Audry, S.; Schäfer, J.; Viers, J.; Blanc, G.; Freydier, R.;
479
Dupré, B., Zn isotopes as tracers of anthropogenic pollution from Zn-ore smelters The Riou
480
Mort–Lot River system. Chemical Geology 2008, 255, (3–4), 295-304.
26 ACS Paragon Plus Environment
Page 27 of 33
481
Environmental Science & Technology
7.
Chen, J.; Gaillardet, J.; Louvat, P.; Huon, S., Zn isotopes in the suspended load of the
482
Seine River, France: Isotopic variations and source determination. Geochimica et Cosmochimica
483
Acta 2009, 73, (14), 4060-4076.
484
8.
Chen, J.; Gaillardet, J. r.; Louvat, P., Zinc isotopes in the Seine River waters, France: A
485
probe of anthropogenic contamination. Environmental Science & Technology 2008, 42, (17),
486
6494-6501.
487
9.
Juillot, F.; Maréchal, C.; Ponthieu, M.; Cacaly, S.; Morin, G.; Benedetti, M.; Hazemann,
488
J. L.; Proux, O.; Guyot, F., Zn isotopic fractionation caused by sorption on goethite and 2-Lines
489
ferrihydrite. Geochimica et Cosmochimica Acta 2008, 72, (19), 4886-4900.
490
10. Balistrieri, L. S.; Borrok, D. M.; Wanty, R. B.; Ridley, W. I., Fractionation of Cu and Zn
491
isotopes during adsorption onto amorphous Fe(III) oxyhydroxide: Experimental mixing of acid
492
rock drainage and ambient river water. Geochimica et Cosmochimica Acta 2008, 72, (2), 311-
493
328.
494
11. Bryan, A. L.; Dong, S.; Wilkes, E. B.; Wasylenki, L. E., Zinc isotope fractionation during
495
adsorption onto Mn oxyhydroxide at low and high ionic strength. Geochimica et Cosmochimica
496
Acta 2015, 157, (0), 182-197.
497
12. Pokrovsky, O. S.; Viers, J.; Freydier, R., Zinc stable isotope fractionation during its
498
adsorption on oxides and hydroxides. Journal of Colloid and Interface Science 2005, 291, (1),
499
192-200.
27 ACS Paragon Plus Environment
Environmental Science & Technology
Page 28 of 33
500
13. Jouvin, D.; Louvat, P.; Juillot, F.; Marechal, C. N.; Benedetti, M. F., Zinc isotopic
501
fractionation: Why organic matters. Environmental Science & Technology 2009, 43, (15), 5747-
502
5754.
503
14. Gélabert, A.; Pokrovsky, O. S.; Viers, J.; Schott, J.; Boudou, A.; Feurtet-Mazel, A.,
504
Interaction between zinc and freshwater and marine diatom species: Surface complexation and
505
Zn isotope fractionation. Geochimica et Cosmochimica Acta 2006, 70, (4), 839-857.
506
15. John, S. G.; Geis, R. W.; Saito, M. A.; Boyle, E. A., Zinc isotope fractionation during
507
high-affinity and low-affinity zinc transport by the marine diatom Thalassiosira oceanica.
508
Limnology and Oceanography 2007, 52, (6), 2710-2714.
509
16. Jouvin, D.; Weiss, D. J.; Mason, T. F. M.; Bravin, M. N.; Louvat, P.; Zhao, F.; Ferec, F.;
510
Hinsinger, P.; Benedetti, M. F., Stable isotopes of Cu and Zn in higher plants: Evidence for Cu
511
reduction at the root surface and two conceptual models for isotopic fractionation processes.
512
Environmental Science & Technology 2012, 46, (5), 2652-2660.
513
17. Moynier, F.; Pichat, S.; Pons, M. L.; Fike, D.; Balter, V.; Albarede, F., Isotopic
514
fractionation and transport mechanisms of Zn in plants. Chemical Geology 2009, 267, (3-4), 125-
515
130.
516 517
18. Weiss, D. J.; Mason, T. F. D.; Zhao, F. J.; Kirk, G. J. D.; Coles, B. J.; Horstwood, M. S. A., Isotopic discrimination of zinc in higher plants. New Phytologist 2005, 165, (3), 703-710.
518
19. Cloquet, C.; Carignan, J.; Lehmann, M.; Vanhaecke, F., Variation in the isotopic
519
composition of zinc in the natural environment and the use of zinc isotopes in biogeosciences: a
520
review. Anal Bioanal Chem 2008, 390, (2), 451-463. 28 ACS Paragon Plus Environment
Page 29 of 33
521 522
Environmental Science & Technology
20. Li, D.; Liu, S.-A.; Li, S., Copper isotope fractionation during adsorption onto kaolinite: Experimental approach and applications. Chemical Geology 2015, 396, (0), 74-82.
523
21. Pokrovsky, O. S.; Viers, J.; Emnova, E. E.; Kompantseva, E. I.; Freydier, R., Copper
524
isotope fractionation during its interaction with soil and aquatic microorganisms and metal
525
oxy(hydr)oxides: Possible structural control. Geochimica et Cosmochimica Acta 2008, 72, (7),
526
1742-1757.
527 528
22. Gu, X.; Evans, L. J., Surface complexation modelling of Cd(II), Cu(II), Ni(II), Pb(II) and Zn(II) adsorption onto kaolinite. Geochimica et Cosmochimica Acta 2008, 72, (2), 267-276.
529
23. Heidmann, I.; Christl, I.; Leu, C.; Kretzschmar, R., Competitive sorption of protons and
530
metal cations onto kaolinite: experiments and modeling. Journal of Colloid and Interface
531
Science 2005, 282, (2), 270-282.
532 533 534 535 536 537 538 539
24. Ikhsan, J.; Johnson, B. B.; Wells, J. D., A comparative study of the adsorption of transition metals on kaolinite. Journal of Colloid and Interface Science 1999, 217, (2), 403-410. 25. Srivastava, P.; Singh, B.; Angove, M., Competitive adsorption behavior of heavy metals on kaolinite. Journal of Colloid and Interface Science 2005, 290, (1), 28-38. 26. Holmgren, G. G. S., A rapid citrate-dithionite extractable iron procedure. Soil Sci. Soc. Am. J. 1967, 31, (2), 210-211. 27. Brunauer, S.; Emmett, P. H.; Teller, E., Adsorption of gases in multimolecular layers. Journal of the American Chemical Society 1938, 60, (2), 309-319.
29 ACS Paragon Plus Environment
Environmental Science & Technology
Page 30 of 33
540
28. Nachtegaal, M.; Sparks, D. L., Effect of iron oxide coatings on zinc sorption mechanisms
541
at the clay-mineral/water interface. Journal of Colloid and Interface Science 2004, 276, (1), 13-
542
23.
543
29. Chen, J.-B.; Louvat, P.; Gaillardet, J.; Birck, J.-L., Direct separation of Zn from dilute
544
aqueous solutions for isotope composition determination using multi-collector ICP-MS.
545
Chemical Geology 2009, 259, (3–4), 120-130.
546
30. Louvat, P.; Moureau, J.; Paris, G.; Bouchez, J.; Noireaux, J.; Gaillardet, J., A fully
547
automated direct injection nebulizer (d-DIHEN) for MC-ICP-MS isotope analysis: application to
548
boron isotope ratio measurements. Journal of Analytical Atomic Spectrometry 2014, 29, (9),
549
1698-1707.
550
31. Borrok, D. M.; Gieré, R.; Ren, M.; Landa, E. R., Zinc isotopic composition of particulate
551
matter generated during the combustion of coal and coal + tire-derived fuels. Environmental
552
Science & Technology 2010, 44, (23), 9219-9224.
553
32. Moeller, K.; Schoenberg, R.; Pedersen, R.-B.; Weiss, D.; Dong, S., Calibration of the
554
new certified reference materials ERM-AE633 and ERM-AE647 for copper and IRMM-3702 for
555
zinc isotope amount ratio determinations. Geostandards and Geoanalytical Research 2012, 36,
556
(2), 177-199.
557
33. Petit, J. C. J.; De Jong, J.; Chou, L.; Mattielli, N., Development of Cu and Zn isotope
558
MC-ICP-MS measurements: Application to suspended particulate matter and sediments from the
559
Scheldt estuary. Geostandards and Geoanalytical Research 2008, 32, (2), 149-166.
30 ACS Paragon Plus Environment
Page 31 of 33
Environmental Science & Technology
560
34. Bigalke, M.; Weyer, S.; Kobza, J.; Wilcke, W., Stable Cu and Zn isotope ratios as tracers
561
of sources and transport of Cu and Zn in contaminated soil. Geochimica et Cosmochimica Acta
562
2010, 74, (23), 6801-6813.
563
35. Sossi, P. A.; Halverson, G. P.; Nebel, O.; Eggins, S. M., Combined separation of Cu, Fe
564
and Zn from rock matrices and improved analytical protocols for stable isotope determination.
565
Geostandards and Geoanalytical Research 2014, n/a-n/a.
566 567
36. Angove, M. J.; Johnson, B. B.; Wells, J. D., Adsorption of cadmium(II) on kaolinite. Colloids and Surfaces A: Physicochemical and Engineering Aspects 1997, 126, (2–3), 137-147.
568
37. Keizer, M. G.; van Riemsdijk, W. H., ECOSAT: Equilibrium calculation of speciation
569
and transport, user manual, Version 4.7. Wageningen Agricultural University, The Ntherlands
570
1999.
571
38. Lützenkirchen, J., Comparison of 1-pK and 2-pK Versions of Surface Complexation
572
Theory by the Goodness of Fit in Describing Surface Charge Data of (Hydr)oxides.
573
Environmental Science & Technology 1998, 32, (20), 3149-3154.
574 575 576 577
39. Bickmore, B. R.; Nagy, K. L.; Sandlin, P. E.; Crater, T. S., Quantifying surface areas of clays by atomic force microscopy. American Mineralogist 2002, 87, (5-6), 780-783. 40. Zbik, M.; Smart, R. S. C., Nanomorphology of kaolinites; comparative SEM and AFM studies. Clays and Clay Minerals 1998, 46, (2), 153-160.
578
41. Sutheimer, S. H.; Maurice, P. A.; Zhou, Q., Dissolution of well and poorly crystallized
579
kaolinites; Al speciation and effects of surface characteristics. American Mineralogist 1999, 84,
580
(4), 620-628. 31 ACS Paragon Plus Environment
Environmental Science & Technology
581 582 583 584
Page 32 of 33
42. Schroth, B. K.; Sposito, G., Surface charge properties of kaolinite. Clays and Clay Minerals 1997, 45, (1), 85-91. 43. Ford, R. G.; Sparks, D. L., The nature of Zn precipitates formed in the presence of pyrophyllite. Environmental Science & Technology 2000, 34, (12), 2479-2483.
585
44. Lee, S.; Anderson, P. R.; Bunker, G. B.; Karanfil, C., EXAFS study of Zn sorption
586
mechanisms on montmorillonite. Environmental Science & Technology 2004, 38, (20), 5426-
587
5432.
588
45. Black, J. R.; Kavner, A.; Schauble, E. A., Calculation of equilibrium stable isotope
589
partition function ratios for aqueous zinc complexes and metallic zinc. Geochimica et
590
Cosmochimica Acta 2011, 75, (3), 769-783.
591
46. Fujii, T.; Moynier, F.; Blichert-Toft, J.; Albarède, F., Density functional theory
592
estimation of isotope fractionation of Fe, Ni, Cu, and Zn among species relevant to geochemical
593
and biological environments. Geochimica et Cosmochimica Acta 2014, 140, 553-576.
594 595 596 597
47. Bigeleisen, J.; Mayer, M. G., Calculation of equilibrium constants for isotopic exchange reactions. The Journal of Chemical Physics 1947, 15, (5), 261-267. 48. Schauble, E. A., Applying stable isotope fractionation theory to new systems. Reviews in Mineralogy and Geochemistry 2004, 55, (1), 65-111.
598
49. Kuzmin, A.; Obst, S.; Purans, J., X-ray absorption spectroscopy and molecular dynamics
599
studies of Zn2+ hydration in aqueous solutions. Journal of Physics: Condensed Matter 1997, 9,
600
(46), 10065.
32 ACS Paragon Plus Environment
Page 33 of 33
Environmental Science & Technology
601
50. Morton, J. D.; Semrau, J. D.; Hayes, K. F., An X-ray absorption spectroscopy study of
602
the structure and reversibility of copper adsorbed to montmorillonite clay. Geochimica et
603
Cosmochimica Acta 2001, 65, (16), 2709-2722.
604 605
51. Strawn, D. G.; Palmer, N. E.; Furnare, L. J.; Goodell, C.; Amonette, J. E.; Kukkadapu, R. K., Copper sorption mechanisms on smectite. Clays and Clay Minerals 2004, 52, (3), 321-333.
606
52. Jordan, N.; Ritter, A.; Foerstendorf, H.; Scheinost, A. C.; Weiß, S.; Heim, K.; Grenzer, J.;
607
Mücklich, A.; Reuther, H., Adsorption mechanism of selenium(VI) onto maghemite. Geochimica
608
et Cosmochimica Acta 2013, 103, (0), 63-75.
609 610
53. Appelo, C. A. J.; Postma, D., Geochemistry, groundwater and pollution. A.A. Balkema Publishers: 1996; p 537.
611
54. Yin, N.-H.; Sivry, Y.; Benedetti, M. F.; Lens, P. N. L.; van Hullebusch, E. D.,
612
Application of Zn isotopes in environmental impact assessment of Zn–Pb metallurgical
613
industries: A mini review. Applied Geochemistry 2016, 64, 128-135.
614
55. Viers, J.; Oliva, P.; Nonell, A.; Gélabert, A.; Sonke, J. E.; Freydier, R.; Gainville, R.;
615
Dupré, B., Evidence of Zn isotopic fractionation in a soil–plant system of a pristine tropical
616
watershed (Nsimi, Cameroon). Chemical Geology 2007, 239, (1–2), 124-137.
617
56. Weiss, D. J.; Rausch, N.; Mason, T. F. D.; Coles, B. J.; Wilkinson, J. J.; Ukonmaanaho,
618
L.; Arnold, T.; Nieminen, T. M., Atmospheric deposition and isotope biogeochemistry of zinc in
619
ombrotrophic peat. Geochimica et Cosmochimica Acta 2007, 71, (14), 3498-3517.
620
33 ACS Paragon Plus Environment