Subscriber access provided by Hong Kong University of Science and Technology Library
Article
Arsenic speciation in Mekong Delta sediments depends on their depositional environment Yuheng Wang, Pierre Le Pape, Guillaume Morin, Maria Asta, Georgina King, Barbora Bartova, Elena Suvorova, Manon Frutschi, Maya Ikogou, Vu Hoai Cong Pham, Phu Le Vo, Frédéric Herman, Laurent Charlet, and Rizlan Bernier-Latmani Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b05177 • Publication Date (Web): 16 Feb 2018 Downloaded from http://pubs.acs.org on February 16, 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.
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 25
Environmental Science & Technology
1 2
Arsenic speciation in Mekong Delta sediments depends on their
3
depositional environment
4 5
Yuheng Wang1, Pierre Le Pape2, Guillaume Morin2, Maria P. Asta1, Georgina King3, Barbora
6
Bártová1, Elena Suvorova1, Manon Frutschi1, Maya Ikogou2, Vu Hoai Cong Pham4, Phu Le
7
Vo4, Frédéric Herman5, Laurent Charlet6 and Rizlan Bernier-Latmani1*
8 9 10
1
11 12
Ecole Polytechnique Fédérale de Lausanne (EPFL), Environmental Microbiology Laboratory (EML), EPFLENAC-IIE-EML, Station 6, CH-1015 Lausanne, Switzerland
2
Centre National de la Recherche Scientifique (CNRS) – Université Pierre et Marie Curie (UPMC Paris 6),
13
Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC, CNRS-UPMC-IRD-MNHN
14
UMR 7590), Campus Jussieu, 4 place Jussieu, 75005 Paris, France
15
3
University of Bern, Institute of Geological Sciences, Baltzerstrasse 1+3, CH-3012, Bern
16
4
Ho Chi Minh City University of Technology, Faculty of Environment and Natural Resources, VNU-HCM, 268
17 18
Ly Thuong Kiet, Ho Chi Minh City, Vietnam 5
19 20 21
University of Lausanne, Institute of Earth Surface Dynamics, Géopolis Building, CH-1015 Lausanne, Switzerland
6
Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, IRD, IFSTTAR, ISTerre, 38000 Grenoble, France
22 23 24
*correspondence to:
[email protected] 25 26 27
ACS Paragon Plus Environment
Environmental Science & Technology
28
Page 2 of 25
TOC Art
29
Himalaya
Sea
Mekong Delta
Before 7200 BP Floodplain
transgression
Mangrove Sea level
O-As(V)
7200 - 7000 BP
transgression
O-As NOM-As(III)
7000 - 5000 BP
regression
O-As NOM-As(III)
A.er 5000 BP
Moderate NOM
30 31
2
O-As NOM-As(III) Pyrite-As
O-As(V)
O-As NOM-As(III)
High NOM
O-As NOM-As(III) Pyrite-As
ACS Paragon Plus Environment
O-As NOM-As(III) Pyrite-As
Page 3 of 25
32
Environmental Science & Technology
Abstract
33
Arsenic (As) contamination in groundwater is pervasive throughout deltaic regions of
34
Southeast Asia and threatens the health of millions. The speciation of As in sediments
35
overlying contaminated aquifers is poorly constrained. Here, we investigate the chemical and
36
mineralogical compositions of sediment cores collected from the Mekong Delta in Vietnam,
37
elucidate the speciation of iron and arsenic, and relate them to the sediment depositional
38
environment. Gradual dissolution of ferric (oxyhydr)oxides with depth is observed down to 7
39
m, corresponding to the establishment of reducing conditions. Within the reduced sediment,
40
layers originating from marine, coastal or alluvial depositional environments are identified
41
and their age is consistent with a late Holocene transgression in the Mekong Delta. In the
42
organic matter- and sulfur-rich layers, arsenic is present in association with organic matter
43
through thiol-bonding and in the form of arsenian pyrite. The highest arsenic concentration
44
(34 – 69 ppm) is found in the peat layer at 16 m and suggests the accumulation of arsenic due
45
to the formation of thiol-bound trivalent arsenic (40% – 55%) and arsenian pyrite (15% –
46
30%) in a paleo-mangrove depositional environment (~7200 yr BP). Where sulfur is limited,
47
siderite is identified, and oxygen- and thiol-bound trivalent arsenic are the predominant
48
forms. It is also worth noting that pentavalent arsenic coordinated to oxygen is ubiquitous in
49
the sediment profile, even in reduced sediment layers. But the identity of the oxygen-bound
50
As species remains unknown. This work shows direct evidence of thiol-bound trivalent
51
arsenic in the Mekong Delta sediments and provides insight to refine the current model of the
52
origin, deposition, and release of arsenic in the alluvial aquifers of the Mekong Delta.
3
ACS Paragon Plus Environment
Environmental Science & Technology
53
Introduction
54
Arsenic (As) contamination in aquifers of the Mekong Delta, as well as in other deltaic
55
areas of South and Southeast Asia1-4 impacts the health of millions. It is widely accepted that
56
arsenate, typically associated with oxidized environments, is bound to iron (oxyhydr)oxides
57
that are transported in river water from the Himalayas and downstream watersheds to those
58
deltaic areas. Arsenate, together with ferric (oxyhydr)oxides, is then deposited by river water
59
at the floodplain sediment surface. Bioreduction of Fe(III) and As(V) during diagenetic burial
60
leads to the subsequent release of As5, 6, and natural organic matter (NOM) released from the
61
margins of peat layers7 fuels this microbial process by serving as an electron donor6, 8.
62
Additionally, other As-bearing phases were uncovered in deltaic sediments. For instance,
63
sequential extractions from sediments in Bangladesh have suggested that arsenic could also
64
bind to NOM9. Micro- and bulk-scale As K-edge X-ray absorption spectroscopy (XAS) has
65
provided evidence for the presence of arsenate, arsenite, and sulfide-bound arsenic in the form
66
of arsenian pyrite, in sediments in the Bengal Basin10 and in the Cambodian Mekong Delta11,
67
12
68
these previous studies10, 11, 12, the presence of additional sulfide-bound As species could not be
69
excluded. Consequently, both the nature and abundance of the sulfide-bound As pool in
70
deltaic sediment layers rich in clay remain unknown7,11,13.
. However, although arsenian pyrite has been evidenced in organic-rich sediment layers in
71
The geological formation history of the Mekong Delta is recorded in the stratigraphic
72
structure of the sediments, which is dominated by alluvial deposition, tidal events and marine
73
deposition14-16. The latter is mainly due to the Holocene transgression, by which the
74
postglacial sea-level rise turned alluvial depositional systems into estuaries and shallow
75
marine sediments17. Mangrove forests followed the migration of the coastline and intruded
76
into more elevated areas in the Mekong Delta18. At about 6,000 – 5,000 yr BP, the sea water
77
level reached its maximum, and marine deposition prevailed throughout the Mekong Delta17,
78
18
. Later, due to the slow regression of the sea level and high fluvial sediment load over the 4
ACS Paragon Plus Environment
Page 4 of 25
Page 5 of 25
Environmental Science & Technology
79
last 4,550 yr BP, the Mekong delta extensively prograded to its current form14, 16. While the
80
stratigraphic structure of Mekong Delta sediments and its relation to the depositional
81
environment are well documented, systematic studies of As speciation along sediment
82
stratigraphic profiles are scarce in the literature. Nonetheless, previous studies have
83
documented that geogenically-derived As in aquifers appears to be associated with sediments
84
deposited during the Holocene marine transgression.19, 20
85
Here, we retrieved sediment cores from the Mekong Delta in Vietnam in an attempt to
86
obtain a comprehensive view of the As-bearing phases in the sediments, to investigate arsenic
87
speciation in more detail and to link it to the sediment depositional environment. We probed
88
the redox-preserved sediments at a range of spatial scales, including bulk mineralogical and
89
chemical analyses, electron microscopy, and XAS at Fe and As K-edges. We uncovered the
90
complex speciation of Fe and As in the sediment relative to its stratigraphic structure and
91
dated the sediment to relate its composition to its depositional environment. More specifically,
92
As is present in a combination of four species: sulfur-bound As(III) associated with NOM,
93
arsenic in pyrite, oxygen-bound As(III) and oxygen-bound As(V). Their fractions in the total
94
As pool varied with depth and thus with depositional environment. This work shows direct
95
evidence of thiol-bound trivalent arsenic in the Mekong Delta sediments and is the first
96
integrated and systematic study of As speciation along a sediment stratigraphic profile in a
97
deltaic environment.
98 99
Materials and Methods
100
Sampling. Sampling trips to the Quoc Thai Commune, An Giang Province, Vietnam took
101
place in January 2014 and January 2015 (Figure S1). The cores (QTC2 and QTC3) were
102
sampled and cut into sections as a function of depth (Figure S2) under an Argon stream or in
103
an N2-atmosphere field anoxic chamber, and the sections were sealed in individual MYLAR®
104
bags. A portion of each segment was dried under vacuum using a desiccator in an anoxic 5
ACS Paragon Plus Environment
Environmental Science & Technology
105
chamber for a few days, ground and homogenized to a fine powder, and stored in sealed
106
serum bottles until analysis. The sediment samples were analyzed in a similar way to that
107
described by Wang et al.21, 22.
108 109
Analysis of sediment chemistry. Water content of the sediment samples was obtained by
110
weighing the samples before and after drying. Sediment total organic carbon (TOC) content
111
was measured in dried homogenized samples with a Shimadzu® TOC-V CPH/CPN analyzer.
112
To study the sediment chemical composition and As concentration, sediment samples were
113
analyzed with a PANalytical® Axios-mAX X-ray fluorescence (XRF) spectrometer using
114
wax-fused pellets for both major and trace elements. The detection limits are 0.01% for major
115
elements and 1 to 5 ppm for trace elements. Accuracy of the analyses was assessed by
116
analysis of standard reference materials.
117 118
Electron microscopy. Scanning electron microscopy (SEM) observation on gold-coated
119
finely ground dry sediment powder was either conducted using a Carl Zeiss® Merlin
120
microscope with a GEMINI II column at 1 kV or using a Gemini Zeiss® Ultra55 microscope
121
at 15 kV. For scanning transmission electron microscopy and X-ray energy dispersive
122
spectroscopy (STEM-EDS), finely ground dry sediment powder was loaded onto C-coated Cu
123
or Ni TEM grids and transferred under anoxic conditions to the microscope (an FEI TecnaiTM
124
Osiris microscope at 200 kV equipped with four windowless Super-X SDD EDX detectors).
125
Quantitative EDS results were obtained with the Bruker® Esprit software. The morphology
126
were studied by high-resolution TEM (HRTEM).
127 128
X-ray diffraction (XRD). XRD analysis was conducted under anoxic conditions. Finely
129
ground dry sediment samples preserved in serum bottles were loaded onto a Si single crystal
130
“Zero Diffraction Plate” in a Jacomex® anaerobic chamber and then inserted in an anaerobic 6
ACS Paragon Plus Environment
Page 6 of 25
Page 7 of 25
Environmental Science & Technology
131
sample chamber equipped with a Kapton® window. XRD measurements were performed
132
under anoxic conditions using CoKa (wavelength = 1.79 Å) radiation on a Panalytical®
133
X’Pert Pro MPD diffractometer mounted in Bragg-Brentano configuration and equipped with
134
an X'Celerator® detector. Data were recorded in the continuous-scan mode in the 3 - 80° 2θ
135
range with a step of 0.033°.
136 137
Iron K-edge X-ray absorption spectroscopy. Iron K-edge extended X-ray absorption fine
138
structure (EXAFS) spectra were collected in transmission detection mode at 80 K using liquid
139
N2 cryostat on the XAFS beamline at ELETTRA (Trieste, Italy) and on the 4-1 beamline at
140
the Stanford Synchrotron Radiation Lightsource (SSRL, CA, USA) using a Si(111) double
141
crystal monochromator. The raw spectra were normalized and background subtracted using
142
the Athena Code45. Linear combination fitting (LCF) of EXAFS data was performed on k3-
143
weighted EXAFS data over a k-range of 3-14 Å-1 and E0 was set to 7,125 eV. Respective
144
signals of the model compounds used for the LCF are presented in Figure S3. Details of the
145
data collection and processing are provided in the Supplementary Information.
146 147
Arsenic K-edge X-ray absorption spectroscopy. Arsenic K-edge X-ray near edge
148
spectrosopcy (XANES) and EXAFS spectra were collected in fluorescence detection mode at
149
10 - 15 K using a liquid He cryostat at the 11-2 wiggler beamline at SSRL (California, USA)
150
with a Si(220) double crystal monochromator. The raw spectra were normalized and
151
background subtracted using the Athena code. LCF of XANES data was performed from -20
152
to 60 eV with E0 set to 11,868 eV. LCF of EXAFS data were performed on k3-weighted
153
EXAFS data over a k-range of 3-12 A-1 and E0 was set up at 11,868 eV. Model compound
154
spectra are presented in Figure S4. Details of the data collection and processing are provided
155
in the Supplementary Information. 7
ACS Paragon Plus Environment
Environmental Science & Technology
Page 8 of 25
156 157
Additional analyses. Sediment porosity and hydraulic conductivity analyses, and sediment
158
dating using either 14C for the NOM-rich layer at 16 m or optically stimulated luminescence
159
(OSL) measurements for selected silicate-rich layers were carried out. The experimental
160
details of these analyses are reported in the Supplementary Information.
161 162
Results
163
Chemical and mineralogical composition, iron speciation and age of the sediments. The
164
sampling area is located in the upper Mekong Delta in Vietnam (Figure S1), in a rice paddy
165
200 m away from the Bassac River, one of the main branches of the Mekong River, in an area
166
where heavy As groundwater contamination was reported
167
collected in January 2014, covers the depth of 0-20 m. A second core (QTC3), collected in
168
January 2015, extends down to 40 m. The two cores, collected from immediately adjacent
169
locations, possess very similar lithologies and chemical compositions as a function of depth
170
(Figure 1 and S5, Tables S1 and S2). XRD analyses reveal that the major minerals are quartz,
171
mica/illite, chlorite, and albite, and that they are distributed largely uniformly in the sediment
172
cores (Figure S6). Additionally, Fe(II)-bearing minerals—pyrite and siderite—were identified
173
at specific depths below 7 m using XRD (Figure S6) and scanning electron microscopy
174
(SEM) (Figure S7).
2-4
. A sediment core (QTC2),
175
In order to determine Fe speciation in the sediment cores, bulk Fe K-edge XAS
176
analysis was carried out for QTC2 samples and the results are shown in Figure S8. LCF of the
177
EXAFS data was performed using relevant model compounds described in the Supplementary
178
Material and Methods section, and the results are reported in Figure 2a and Table S3. Ferric
179
oxyhydroxides—goethite and ferrihydrite—were detected in the depth range of 0-6 m with a
180
decreasing abundance with depth. Phyllosilicate minerals, e.g., illite and chlorite, also host a
8
ACS Paragon Plus Environment
Page 9 of 25
Environmental Science & Technology
181
significant fraction of total Fe. Below 7 m, Fe-reducing conditions begin to prevail and an
182
alternation of pyrite and siderite is observed. Biotite is the main Fe-bearing phyllosilicate
183
mineral in the reduced sediment layers, where about 9%-19% of Fe is also found to be
184
associated with NOM.
185
Based on Fe K-edge XAS data, sediment chemical and mineralogical composition,
186
and visual observation of the cores, we describe five distinct sediment types henceforth
187
referred to as type I through V (Figure 1, Table S4). Type I (orange layer in Figure 1a) is a
188
clay layer in the depth range of 0-7 m, with a color transitioning from brown-orange, to light
189
green and to light gray. It is NOM-poor (TOC = 0.2 – 0.6% w/w) (Figure 1d, Table S4). In
190
this sediment layer, the Fe-bearing phases are phyllosilicates (illite and chlorite) and ferric
191
oxyhydroxides (goethite and ferrihydrite). The latter two account for ca. half of the Fe pool at
192
0 - 6 m, but disappear at 7 m where traces of pyrite are detected (Figure 2a). Thus, this layer
193
(type I) represents oxic sediments undergoing a redox transition from 6 to 7 m. The layer at 2
194
m was dated to 2,560 ± 360 yr (Table 1). Type II (dark gray layers in Figure 1a) is
195
represented by two gray clay layers at 8 – 11 m and 14 m, respectively, and includes moderate
196
to high NOM (TOC = 1.4 - 7.5% w/w) (Figure 1d, Table S4), and high pyrite content (1.2 -
197
1.9% w/w) (Figure 2a, Table S4). The 8 m layer is at the shallowest depth at which pyrite
198
precipitation is observed with XRD, suggesting that Fe(III)- and sulfate-reducing conditions
199
were extant at the time of formation of this mineral. This sediment type contains pyrite and
200
biotite as the main Fe-bearing components and Fe in NOM as a minor component. The 10 m
201
layer was dated to 7,830 ± 890 yr (Table 1). Type III (light gray layer in Figure 1a) is a gray
202
clay layer at 12 m and is moderate in NOM (TOC = 1.0-1.4% w/w) (Figure 1d, Table S4) but
203
rich in siderite (1.6% w/w) (Figure 2a, Table S4). The sulfur content is significantly lower
204
than in adjacent sediment layers (Figure 1c). This sediment contains biotite and siderite as the
205
main Fe components, and pyrite as a minor phase (Figure 2a, Table S4). This suggests that
206
sulfate was less abundant during the deposition than for the adjacent type II sediment. This 9
ACS Paragon Plus Environment
Environmental Science & Technology
207
layer (at 12 m) was dated to 6,480 ± 790 (Table 1). Type IV (black layer in Figure 1a) is at 16
208
m and is a peat layer with a very high NOM (TOC = up to 33.9% w/w) (Figure 1d, Table S4)
209
and pyrite contents (2.6% w/w) (Figure 2a, Table S4). Consistent with this composition, Fe
210
speciation is dominated by pyrite and NOM-associated Fe (Figure 2a). The composition also
211
underscores an abundant sulfate supply during sediment deposition. This layer was dated
212
using 14C to 8,079 ± 708 cal. yr BP (Table 1). Type V (green layer in Figure 1a) is a gray clay
213
layer at 17 - 33 m, that is relatively poor in NOM (TOC = 0.1 – 0.9% w/w, except for the
214
layer at 17 m where TOC = 2.7 – 4.8% w/w, as it lies immediately beneath type IV NOM-rich
215
layer) and pyrite-free (Table S4), and in which siderite content varies with depth (Figure 2a,
216
Figure S6). Fe is mainly present in biotite and siderite and, to a much lesser extent, associated
217
with NOM (Figure 2a). This suggests that sulfate was scarce during the deposition of this
218
layer. This layer was dated at two depths: to 8,840 ± 1,670 at 18m (near the top of the layer)
219
and at 16,870 ± 2,190 yr at 33 m (near the bottom of the layer). Taking into account the 2σ
220
uncertainties of the determined ages, we confirm a younger-to-older chronological sequence
221
from the top to the bottom of the sediment profile (Table 1).
222 223
Arsenic redox states in the sediment. The As K-edge XANES spectra (Figure 3a) of type II
224
(8 and 14 m), type III (12 m), type IV (16 m), and type V (18 m, 20 m and 35 m) sediments
225
display two absorption maxima at about 11,869 eV and 11,875 eV. The first absorption peak
226
(11,869 eV) corresponds to a mixture of As in arsenian pyrite (at 11,868.7 eV), S-bound
227
As(III) (at 11,869.5 eV) and O-coordinated As(III) (arsenite) (at 11,871.4 eV)23-25. The
228
second absorption maximum corresponds to O-coordinated As(V) (arsenate) (at 11,875.0
229
eV)24-26 (Figure S4). The XANES spectra of the type I (2 m and 7 m) sediment exhibit a main
230
absorption maximum at 11,875 eV, corresponding to O-bound As(V), along with a
231
contribution from O-bound As(III) at 11,871.4 eV. Thus, in sediment types I to V, mixed As
232
redox states prevail and their composition is quantified using LCF analysis (Figure 3a, Table 10
ACS Paragon Plus Environment
Page 10 of 25
Page 11 of 25
Environmental Science & Technology
233
S5, Figure S9). Even in the oxidized sediment (type I), we detect a contribution from As(III),
234
presumably from As(V) reduction. In the sediment types that include NOM (types II, III, IV),
235
the XANES spectra are dominated by As in arsenian pyrite, and S-bound As(III) with minor
236
O-bound As species. In type V sediments, XANES data indicate contributions from S-bound
237
As(III) and O-bound As(III) and As(V) species. LCF analysis of the As K-edge EXAFS data
238
shows the same results and thus reinforces the reliability of the fitting components (Figure 2b,
239
Table S5). It is worth noting that the LCF analysis of both the As K-edge XANES and
240
EXAFS data must invoke As(V) oxyanions to obtain acceptable fits for all samples including
241
reduced sediments of type II to V (Figure 2b, Figure 3, Table S5).
242 243
Arsenian pyrite in the sediment. LCF results for As EXAFS and XANES data (Figure 2b,
244
Figure 3, Table S5) underscore the significant contribution of arsenian pyrite to sediment
245
types II to IV, while evidencing its absence in types I and V. Scanning Electron Microscopy
246
(SEM) observation of type II, III and IV samples revealed the presence of both framboidal
247
and euhedral pyrite particles of varying size (Figure S7). STEM coupled to EDS and SAED
248
analysis of sediment types II (8, 10 and 14 m), type III (12 m), and type IV (16 m) confirms
249
that both As-bearing and As-free pyrite particles occur and that the former group contains As
250
atomic concentrations of 0.06 - 0.4%. (Figure S10-S14, Table S6). Pyrite particles in the
251
sediment samples were readily identified during the STEM observation due to their
252
distinctive morphology. Moreover, using a positive correlation between the abundance of
253
pyrite-associated arsenic and pyrite abundance in the sediment, we estimated an average As
254
concentration of ~280 ppm wt. in pyrite in these sediments (Figure S15).
255 256
Thiol-bound arsenic(III) in the sediment. Arsenian pyrite alone does not account for the
257
entire pool of S-bound As in sediment types II (8 and 14 m), III (12 m), IV (16 m), and V (18
258
m) according to LCF results for As K-edge XAS data (Figure 3b). A sulfur-bound As(III) 11
ACS Paragon Plus Environment
Environmental Science & Technology
259
species lacking second shell coordination is needed to fit the EXAFS data. This As species
260
corresponds rather to As(III) bound to thiol groups27,28 than to amorphous As2S3
261
precipitates29,30, as the former exhibits a shorter As-S distance and yields a better fit (Figure
262
S16). Neither amorphous As2S3 nor NOM-associated As could be observed by STEM-EDS,
263
likely because NOM and S are 3-4 orders of magnitude more abundant than As, diluting the
264
pool of As associated with either. Thiol-bound As(III) is particularly abundant in type IV
265
organic-rich sediment (16 m, peat layer) as shown in Figure S17, which supports the finding
266
that As(III) is bound to NOM thiol groups.
267 268
Discussion
269
Iron speciation in relation to sediment stratigraphic structure and depositional
270
environment. The sediment core structure evidences the alternation of pyrite-rich (type II (8
271
– 11 m, 14 m) and type IV (16 m)) and pyrite-poor but siderite-rich (type III (12 m) and type
272
V (17 - 37 m)) sediments (Figure 2a, Figure S6). We hypothesize that this originates from
273
varying sulfate abundance in the depositional environment, which resulted from a late
274
Holocene transgression of the Mekong delta. The Mekong delta has prograded from the
275
Vietnam-Cambodia border to the South China Sea in the last 20,000 yrs17. Thus, sediment
276
layers at the bottom of the cores (type V) represent an alluvial depositional environment
277
where limited sulfate was available, and siderite, rather than pyrite, was formed during
278
diagenetic burial. The dating of the sediment layer at 33 m to 16,870 ± 2,190 yr (Table 1) is
279
broadly consistent with this progradation phase. However, a transgressive phase that lasted
280
from ~15,000 to ~5,000 yr BP14, 17, 18, and during which the sea water level rose by ~75 m17, is
281
consistent with the pyrite-bearing middle of the core (type II to IV). The sediment layer at 18
282
m (type V) dated at 8,840 ± 1,670 yr does not include pyrite but rather siderite, and represents
283
a time at which sea level was rising but had not reached this particular location. The peat layer
284
at 16 m (type IV) is dated to 8,079 ± 708 cal. yr BP (Table 1) and corresponds to the short 12
ACS Paragon Plus Environment
Page 12 of 25
Page 13 of 25
Environmental Science & Technology
285
period during which a mangrove forest covered the current sampling location18. Mangrove
286
forests usually thrive at marine coastlines and promote Fe(III) and sulfate reduction and pyrite
287
precipitation31 as well as NOM accumulation in sediments32. The common presence of
288
framboidal and euhedral pyrite (Figure S7) accompanied by high NOM content (Figure 1,
289
Table S4) observed in these layers is reminiscent of present-day mangrove sediments33. At
290
about 6,000 – 5,000 yr BP, the sea water level reached its maximum17, 18 and the depositional
291
environment turned to a shallow marine depositional environment implying lower NOM input
292
under sulfate-reducing conditions. Such conditions explain the presence of pyrite in the
293
sediment layers at 12 m (type III) and at 10 and 14 m (type II) above the type IV peat layer.
294
Consistent with these findings, Stuckey et al.10 observed a peat layer at a shallower depth (6
295
m) at a location upstream from our site along the Mekong River11 that was dated at ~5,500 yr
296
BP, commensurate with continued transgression until about 5,000 BP. Finally, type I
297
sediment represents the most recent alluvial deposition in the latest 5,000 yrs18, and formed
298
above sea level after marine regression. The layer at 2 m is dated to 2,560 ± 360 yr (Table 1),
299
which is consistent with that interpretation.
300
Therefore, we propose that type IV sediments represent deposition in a paleo-
301
mangrove environment and type II and III represent shallow marine environments (Figure 4).
302
The lower pyrite and S content in type III sediment and the presence of siderite resulting from
303
Fe (oxyhydr)oxides reduction could be due to a transient high input of alluvial sediment or to
304
a lower sea level between two transgression events. A previous study in the Ganges-
305
Brahmaputra Delta in Bangladesh has evidenced the lack of pyrite in a 20-100 m depth
306
interval and similarly surmised a low sulfur flux during deposition10. Overall, the dating
307
results correspond well to the iron speciation, providing strong support to the interpretation of
308
the depositional environment for each sediment type.
309
13
ACS Paragon Plus Environment
Environmental Science & Technology
310
Evidence of thiol-bound arsenic associated with NOM. Both NOM- and pyrite-associated
311
As are predominant S-bound As species found in the sediment cores investigated here. Thiol-
312
bound As(III) associated with NOM is particularly abundant in the organic-rich peat layer
313
(type IV) but also occurs in lower amounts in most of the other reduced layers (Figure 2b.
314
Table S4). This species identified here with As K-edge XAS analysis likely accounts for a
315
significant part of the As-sulfide pool reported in a peat layer in the Cambodian Mekong
316
Delta by Stuckey et al.11, 12, in addition to arsenian pyrite. Arsenic binding to high affinity
317
NOM-thiol groups was recently reported in As-rich peat from a wetland27 but was never
318
demonstrated in deltaic environments. Here, accumulation of thiol-bound As(III) in the
319
organic-rich layer at 16 m could have occurred in the paleo-mangrove that was extant around
320
7,200 yr BP and could continue to be active in trapping dissolved arsenic released by As(V)
321
bioreduction6 at the sediment redox transition boundary. Alternatively, it could serve as a
322
local source of organic carbon and As to the aquifer in addition to the reductive dissolution of
323
iron (oxyhydr)oxides and this possibility remains to be ascertained.
324 325
Authigenic origin of pyrite-associated arsenic. Arsenian pyrite was previously identified in
326
a peat layer in sediments from the Mekong delta in Cambodia11. Here, we show that arsenian
327
pyrite is present in the NOM-rich peat layer (type IV) interpreted as paleo-mangrove deposits,
328
as well as in less NOM-rich clay layers (types II and III) interpreted as shallow marine
329
deposits. In these environments, both the formation of authigenic arsenian pyrite and the
330
preservation of detrital arsenian pyrite are possible. The latter could originate from arsenic-
331
rich coal deposits in the upper Mekong regions34, where it is the dominant arsenic host35,36.
332
However, it is generally expected to undergo oxidation during erosional processes37. Thus,
333
arsenian pyrite is likely to have formed authigenically in the mangrove sediments,
334
incorporating As into the pyrite structure38, 39, 40. Indeed, both the framboidal and euhedral
14
ACS Paragon Plus Environment
Page 14 of 25
Page 15 of 25
Environmental Science & Technology
335
pyrite particles found in the sediments are reported to form in low-temperature sulfate-
336
reducing environments40, such as marine41 and mangrove sediments33.
337 338
Oxygen-bound arsenic in the sediment. In the present study, LCF analysis of both the As
339
K-edge XANES and EXAFS data must invoke sorbed As(III) and As(V) oxyanions to obtain
340
acceptable fits (Figure 2b, Figure 3, Table S5). This is consistent with previous studies, which
341
have shown that bulk As XAS analyses of reduced deltaic sediments exhibited As(III) and
342
As(V) along with arsenian pyrite10,
343
(oxyhydr)oxides that were not bioavailable in the sediment. However, this interpretation has
344
been controversial due to the need to implicate the persistence of ferric minerals in reduced
345
sediments42. In the present cores, at depths below 7 m (i.e., in sediment types II-V), ferric
346
oxyhydroxides are under the detection limit of EXAFS (~10 % of total Fe) and were not
347
observed by STEM, suggesting that As(V) could be associated with another phase in these
348
reduced layers. Moreover, STEM-EDS analysis shows that siderite43 or phyllosilicate
349
minerals44, 45 known to potentially bear As(III) and As(V) did not show detectable amount of
350
As in any of the reduced sediment samples (data not shown). Therefore, the identities of the
351
oxygen bound As(III) and As(V) pools in the reduced sediment remain unknown but are
352
unlikely to be ferric (oxyhydr)oxides.
12
and invoked As(V) adsorbed onto ferric
353 354
Environmental implications. The correlation of As speciation with both Fe speciation and
355
the chronostratigraphic record is consistent with the distribution of As species being governed
356
by depositional history. These processes occurring during diagenesis include NOM
357
accumulation and early diagenetic pyritization resulting from Holocene transgression in the
358
Mekong Delta. Consequently, the distribution of As species in these sediments, shown to
359
result from sedimentary processes, helps constrain the conceptual model of the origin of As
360
species. The model for As release from sediment to groundwater in South and Southeast Asia 15
ACS Paragon Plus Environment
Environmental Science & Technology
361
that has emerged from previous work invokes the reduction of As-bearing Fe(III)
362
(oxyhydr)oxides in the near surface environment of flooded soils6, 8, 46. Here, we confirm the
363
presence of Fe(III) (oxyhydr)oxides in the first 7 m of the sediment from the Mekong Delta in
364
Vietnam and their gradual disappearance at depth. Further, we show that arsenian pyrite and
365
thiol-bound trivalent arsenic in NOM are the predominant As species in the reduced sediment.
366
Accumulation of both pyrite-associated and NOM-bound As in the reduced sediment layers
367
are proposed to have occurred in the paleo-mangrove environment and in shallow marine
368
sediments. Complexation of As by NOM could also be a relevant process for the Bengal
369
Delta as arsenic enrichment by NOM during marine transgression has also been proposed for
370
sediments there47.
371
Immobilization of As by pyrite and NOM in the reduced sediment layers could remain
372
an active mechanism to date, trapping dissolved arsenic9, 11, 27 released by As(V) bioreduction6
373
at the sediment redox boundary, and advected down the sediment profile48 (which is
374
supported by hydraulic conductivity data for the sediments (Table S7)). On the other hand,
375
pyrite-associated and NOM-bound As could serve as local sources of As. To date, no study
376
has specifically tackled this potential and, at this time, we cannot exclude that possibility.
377
Further investigations, potentially laboratory-based experiments on these sediments and
378
sediment porewaters analyses, are needed to evaluate the reactivity of the As-rich sediment
379
layers inherited from depositional processes with respect to their ability to trap or release As
380
to advected porewaters.
381 382 383 384 385 386 387 388 389 390
ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI: xxx AUTHOR INFORMATION Corresponding author *E-mail:
[email protected] 16
ACS Paragon Plus Environment
Page 16 of 25
Page 17 of 25
391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441
Environmental Science & Technology
Notes The authors declare no competing financial interest. ACKNOWLEDGMENTS This study is funded by the Swiss National Science Foundation under grant number 200021_157007. It was conducted under the framework of a CARE-RESCIF initiative. We acknowledge funding from the CODEV-EPFL Seed Money Fund and from RESCIF-EPFL. We also acknowledge funding for the OSL dating from the Swiss National Science Foundation under grant number PZ00P2_167960. We are grateful to CARE for access to the laboratory facility at HCMUT. We thank the following individuals in Vietnam for help with fieldwork: Mr. Nguyen Huu Loc (DONRE of An Giang), Mr. Pham Van Si (Division of Natural Resources & Environment, An Phu district, An Giang province), Mr. Hieu (staff of Quoc Thai Commune), Mr. Huynh (Land owner at Quoc Thai), and Vu Viet Anh, Le Hoang Anh, Nguyen Thi Bao Tu, Pham Kim Bao Ngan, Hua Bao Anh and Le Khanh (students of Faculty of Environment & Natural Resources, HCMUT - VNU). We thank Elena Rossel from EPFL for help with TOC analyses, CEAL of EPFL for analytical support, Jean-Claude Lavanchy from University of Lausanne for assistance during XRF analysis, and CIME at EPFL for access to electron microscopes. We also thank Ludovic Delbes and Benoit Baptiste for their help in anaerobic XRD measurements as well as Imène Esteve for her help in SEM data acquisition at IMPMC. The Fe K-edge XAS experiments were performed on the XAFS beamline at ELETTRA, Trieste, Italy and on the 4-1 beamline at the Stanford Synchrotron Radiation Lightsource (SSRL). The As K-edge XAS experiments were performed on the 11-2 beamline at SSRL. We are grateful for the technical assistance received during the synchrotron-based analyses from Ryan Davis and John Bargar at SSRL, and from Luca Olivi, Clara Guglieri Rodriguez and Giuliana Aquilanti at ELETTRA. Use of the SSRL, SLAC National Accelerator Laboratory, is supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC02-76SF00515.
References 1. Buschmann, J.; Berg, M.; Stengel, C.; Winkel, L.; Sampson, M. L.; Trang, P. T. K.; Viet, P. H. Contamination of drinking water resources in the Mekong delta floodplains: Arsenic and other trace metals pose serious health risks to population. Environment International 2008, 34, 756-764. 2. Nguyen, K. P.; Itoi, R. Source and release mechanism of arsenic in aquifers of the Mekong Delta, Vietnam. Journal of Contaminant Hydrology 2009, 103, 58-69. 3. Hoang, T. H.; Bang, S.; Kim, K.-W.; Nguyen, M. H.; Dang, D. M. Arsenic in groundwater and sediment in the Mekong River delta, Vietnam. Environ. Pollut. 2010, 158, 2648-2658. 4. Erban, L. E.; Gorelick, S. M.; Zebker, H. A.; Fendorf, S. Release of arsenic to deep groundwater in the Mekong Delta, Vietnam, linked to pumping-induced land subsidence. Proc. Natl. Acad. Sci. U.S.A. 2013, 110, 13751-13756. 5. Nickson, R.; McArthur, J.; Burgess, W.; Ahmed, K. M.; Ravenscroft, P.; Rahmann, M. Arsenic poisoning of Bangladesh groundwater. Nature 1998, 395, 338-338. 6. Stuckey, J. W.; Schaefer, M. V.; Kocar, B. D.; Benner, S. G.; Fendorf, S. Arsenic release metabolically limited to permanently water-saturated soil in Mekong Delta. Nature Geosci 2016, 9, 70-76. 7. McArthur, J. M.; Banerjee, D. M.; Hudson-Edwards, K. A.; Mishra, R.; Purohit, R.; Ravenscroft, P.; Cronin, A.; Howarth, R. J.; Chatterjee, A.; Talukder, T.; Lowry, D.; Houghton, S.; Chadha, D. K. Natural organic matter in sedimentary basins and its relation to 17
ACS Paragon Plus Environment
Environmental Science & Technology
442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491
arsenic in anoxic ground water: the example of West Bengal and its worldwide implications. Applied Geochemistry 2004, 19, 1255-1293. 8. Kocar, B. D.; Polizzotto, M. L.; Benner, S. G.; Ying, S. C.; Ung, M.; Ouch, K.; Samreth, S.; Suy, B.; Phan, K.; Sampson, M.; Fendorf, S. Integrated biogeochemical and hydrologic processes driving arsenic release from shallow sediments to groundwaters of the Mekong delta. Applied Geochemistry 2008, 23, 3059-3071. 9. Anawar, H. M.; Akai, J.; Komaki, K.; Terao, H.; Yoshioka, T.; Ishizuka, T.; Safiullah, S.; Kato, K. Geochemical occurrence of arsenic in groundwater of Bangladesh: sources and mobilization processes. Journal of Geochemical Exploration 2003, 77, 109-131. 10. Lowers, H. A.; Breit, G. N.; Foster, A. L.; Whitney, J.; Yount, J.; Uddin, M. N.; Muneem, A. A. Arsenic incorporation into authigenic pyrite, Bengal Basin sediment, Bangladesh. Geochimica et Cosmochimica Acta 2007, 71, 2699-2717. 11. Stuckey, J. W.; Schaefer, M. V.; Kocar, B. D.; Dittmar, J.; Pacheco, J. L.; Benner, S. G.; Fendorf, S. Peat formation concentrates arsenic within sediment deposits of the Mekong Delta. Geochimica et Cosmochimica Acta 2015, 149, 190-205. 12. Stuckey, J. W.; Schaefer, M. V.; Benner, S. G.; Fendorf, S. Reactivity and speciation of mineral-associated arsenic in seasonal and permanent wetlands of the Mekong Delta. Geochimica et Cosmochimica Acta 2015, 171, 143-155. 13. Polizzotto, M. L.; Harvey, C. F.; Sutton, S. R.; Fendorf, S. Processes conducive to the release and transport of arsenic into aquifers of Bangladesh. Proceedings of the National Academy of Sciences of the United States of America 2005, 102, 18819-18823. 14. Tjallingii, R.; Stattegger, K.; Wetzel, A.; Van Phach, P. Infilling and flooding of the Mekong River incised valley during deglacial sea-level rise. Quaternary Science Reviews 2010, 29, 1432-1444. 15. Tamura, T.; Saito, Y.; Sieng, S.; Ben, B.; Kong, M.; Choup, S.; Tsukawaki, S. Depositional facies and radiocarbon ages of a drill core from the Mekong River lowland near Phnom Penh, Cambodia: Evidence for tidal sedimentation at the time of Holocene maximum flooding. Journal of Asian Earth Sciences 2007, 29, 585-592. 16. Tamura, T.; Saito, Y.; Sieng, S.; Ben, B.; Kong, M.; Sim, I.; Choup, S.; Akiba, F. Initiation of the Mekong River delta at 8 ka: evidence from the sedimentary succession in the Cambodian lowland. Quaternary Science Reviews 2009, 28, 327-344. 17. Xue, Z.; Liu, J. P.; DeMaster, D.; Van Nguyen, L.; Ta, T. K. O. Late Holocene Evolution of the Mekong Subaqueous Delta, Southern Vietnam. Marine Geology 2010, 269, 46-60. 18. Lap Nguyen, V.; Ta, T. K. O.; Tateishi, M. Late Holocene depositional environments and coastal evolution of the Mekong River Delta, Southern Vietnam. Journal of Asian Earth Sciences 2000, 18, 427-439. 19. Ravenscroft, P.; Burgess, W. G.; Ahmed, K. M.; Burren, M.; Perrin, J. Arsenic in groundwater of the Bengal Basin, Bangladesh: Distribution, field relations, and hydrogeological setting. Hydrogeology Journal 2005, 13, 727-751. 20. Winkel, L. H. E.; Trang, P. T. K.; Lan, V. M.; Stengel, C.; Amini, M.; Ha, N. T.; Viet, P. H.; Berg, M. Arsenic pollution of groundwater in Vietnam exacerbated by deep aquifer exploitation for more than a century. Proceedings of the National Academy of Sciences 2011, 108, 1246-1251. 21. Wang, Y.; Frutschi, M.; Suvorova, E.; Phrommavanh, V.; Descostes, M.; Osman, A. A. A.; Geipel, G.; Bernier-Latmani, R. Mobile uranium(IV)-bearing colloids in a miningimpacted wetland. Nature Communications 2013, 4, 2942. 22. Wang, Y.; Bagnoud, A.; Suvorova, E.; McGivney, E.; Chesaux, L.; Phrommavanh, V.; Descostes, M.; Bernier-Latmani, R. Geochemical Control on Uranium(IV) Mobility in a Mining-Impacted Wetland. Environmental Science & Technology 2014, 48, 10062-10070.
18
ACS Paragon Plus Environment
Page 18 of 25
Page 19 of 25
492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542
Environmental Science & Technology
23. Wang, Y.; Morin, G.; Ona-Nguema, G.; Menguy, N.; Juillot, F.; Aubry, E.; Guyot, F.; Calas, G.; Brown Jr., G. E. Arsenite sorption at the magnetite-water interface during aqueous precipitation of magnetite: EXAFS evidence for a new arsenite surface complex. Geochim. Cosmochim. Acta 2008, 72, 2573-2586. 24. Wang, Y.; Morin, G.; Ona-Nguema, G.; Juillot, F.; Guyot, F.; Calas, G.; Brown Jr., G. E. Evidence for Different Surface Speciation of Arsenite and Arsenate on Green Rust: An EXAFS and XANES Study. Environ. Sci. Technol. 2010, 44, 109-115. 25. Wang, Y.; Morin, G.; Ona-Nguema, G.; Brown, G. E. Arsenic(III) and Arsenic(V) Speciation during Transformation of Lepidocrocite to Magnetite. Environmental Science & Technology 2014, 48, 14282-14290. 26. Wang, Y.; Morin, G.; Ona-Nguema, G.; Juillot, F.; Calas, G.; Brown Jr., G. E. Distinctive Arsenic(V) Trapping Modes by Magnetite Nanoparticles Induced by Different Sorption Processes. Environ. Sci. Technol. 2011, 45, 7258-7266. 27. Langner, P.; Mikutta, C.; Kretzschmar, R. Arsenic sequestration by organic sulphur in peat. Nature Geosci 2012, 5, 66-73. 28. Miot, J.; Morin, G.; Skouri-Panet, F.; Férard, C.; Aubry, E.; Briand, J.; Wang, Y.; Ona-Nguema, G.; Guyot, F.; Brown, G. E. XAS Study of Arsenic Coordination in Euglena gracilis Exposed to Arsenite. Environ. Sci. Technol. 2008, 42, 5342-5347. 29. Newman, D. K.; Beveridge, T. J.; Morel, F. M. M. Precipitation of arsenic trisulfide by Desulfotomaculum auripigmentum. Applied and Environmental Microbiology 1997, 63, 2022-2028. 30. Burton, E. D.; Johnston, S. G.; Kocar, B. D. Arsenic Mobility during Flooding of Contaminated Soil: The Effect of Microbial Sulfate Reduction. Environmental Science & Technology 2014, 48, 13660-13667. 31. Otero, X. L.; Ferreira, T. O.; Huerta-Díaz, M. A.; Partiti, C. S. M.; Souza Jr, V.; VidalTorrado, P.; Macías, F. Geochemistry of iron and manganese in soils and sediments of a mangrove system, Island of Pai Matos (Cananeia — SP, Brazil). Geoderma 2009, 148, 318335. 32. Chmura, G. L.; Anisfeld, S. C.; Cahoon, D. R.; Lynch, J. C. Global carbon sequestration in tidal, saline wetland soils. Global Biogeochemical Cycles 2003, 17, 1111. 33. Noël, V.; Marchand, C.; Juillot, F.; Ona-Nguema, G.; Viollier, E.; Marakovic, G.; Olivi, L.; Delbes, L.; Gelebart, F.; Morin, G. EXAFS analysis of iron cycling in mangrove sediments downstream a lateritized ultramafic watershed (Vavouto Bay, New Caledonia). Geochimica et Cosmochimica Acta 2014, 136, 211-228. 34. Zhou, Y.; Ren, Y. Distribution of arsenic in coals of Yunnan Province, China, and its controlling factors. International Journal of Coal Geology 1992, 20, 85-98. 35. Huggins, F. E.; Shah, N.; Zhao, J.; Lu, F.; Huffman, G. P. Nondestructive determination of trace element speciation in coal and coal ash by XAFS spectroscopy. Energy & Fuels 1993, 7, 482-489. 36. Yudovich, Y. E.; Ketris, M. P. Arsenic in coal: a review. International Journal of Coal Geology 2005, 61, 141-196. 37. Nickson, R. T.; McArthur, J. M.; Ravenscroft, P.; Burgess, W. G.; Ahmed, K. M. Mechanism of arsenic release to groundwater, Bangladesh and West Bengal. Appl. Geochem. 2000, 15, 403-413. 38. Savage, K. S.; Tingle, T. N.; O’Day, P. A.; Waychunas, G. A.; Bird, D. K. Arsenic speciation in pyrite and secondary weathering phases, Mother Lode Gold District, Tuolumne County, California. Applied Geochemistry 2000, 15, 1219-1244. 39. Le Pape, P.; Blanchard, M.; Brest, J.; Boulliard, J.-C.; Ikogou, M.; Stetten, L.; Wang, S.; Landrot, G.; Morin, G. Arsenic Incorporation in Pyrite at Ambient Temperature at Both Tetrahedral S–I and Octahedral FeII Sites: Evidence from EXAFS–DFT Analysis. Environmental Science & Technology 2017, 51, 150-158. 19
ACS Paragon Plus Environment
Environmental Science & Technology
543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567
40. Blanchard, M.; Alfredsson, M.; Brodholt, J.; Wright, K.; Catlow, C. R. A. Arsenic incorporation into FeS2 pyrite and its influence on dissolution: A DFT study. Geochimica Et Cosmochimica Acta 2007, 71, 624-630. 41. Wilkin, R. T.; Barnes, H. L. Formation processes of framboidal pyrite. Geochimica et Cosmochimica Acta 1997, 61, 323-339. 42. Polizzotto, M. L.; Harvey, C. F.; Li, G.; Badruzzman, B.; Ali, A.; Newville, M.; Sutton, S.; Fendorf, S. Solid-phases and desorption processes of arsenic within Bangladesh sediments. Chem. Geol. 2006, 228, 97-111. 43. Jönsson, J.; Sherman, D. M. Sorption of As(III) and As(V) to siderite, green rust (fougerite) and magnetite: Implications for arsenic release in anoxic groundwaters. Chem. Geol. 2008, 255, 173-181. 44. Goldberg, S. Competitive adsorption of arsenate and arsenite on oxides and clay minerals. Soil Science Society of America Journal 2002, 66, 413-421. 45. Chakraborty, S.; Wolthers, M.; Chatterjee, D.; Charlet, L. Adsorption of arsenite and arsenate onto muscovite and biotite mica. Journal of Colloid and Interface Science 2007, 309, 392-401. 46. Islam, F. S.; Gault, A. G.; Boothman, C.; Polya, D. A.; Charnock, J. M.; Chatterjee, D.; Lloyd, J. R. Role of metal-reducing bacteria in arsenic release from Bengal delta sediments. Nature 2004, 430, 68-71. 47. Chatterjee, D.; Roy, R. K.; Basu, B. B. Riddle of arsenic in groundwater of Bengal Delta Plain—role of non-inland source and redox traps. Environmental Geology 2005, 49, 188-206. 48. Polizzotto, M. L.; Kocar, B. D.; Benner, S. G.; Sampson, M.; Fendorf, S. Near-surface wetland sediments as a source of arsenic release to ground water in Asia. Nature 2008, 454, 505-508.
568
20
ACS Paragon Plus Environment
Page 20 of 25
Page 21 of 25
569 570 571
Environmental Science & Technology
Table 1. Dating results for representative depths of all layer types in core QTC3 determined using optically stimulated luminescence (OSL) dating or 14C dating (cal. yr BP). Note: σ represents standard deviation of sediment dating results. Depth (m)
Sediment type Dating method
Age ± 2σ (yr)
2
I
OSL
2,560 ± 360
10
II
OSL
7,830 ± 890
12
III
OSL
6,480 ± 790
16
IV
18
V
OSL
8,840 ± 1,670
33
V
OSL
16, 870 ± 2,190
14
C
8,079 ± 708
572
21
ACS Paragon Plus Environment
Environmental Science & Technology
Page 22 of 25
573 0 5 10
Type I: clay + Fe(III) oxides
Type II: clay + OM + pyrite
25 30
c
d
10
Type III: clay + siderite + pyrite Type II Type IV: OM + pyrite
20
b
5
QTC3
Depth (m)
15
0
Depth (m)
a
QTC2
574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598
15 20 25
Type V: clay + siderite
30 35
0 20 40 60 80 0 1 2 3 4 5 Sediment As content Sediment S content (ppm) (%)
0
10 20 30 40 Sediment TOC (%)
Figure 1. (a) Scheme showing the lithology of the QTC2 and QTC3 sediment cores based on on-site observation and powder XRD analysis of sediment samples (Figure S6). (b) arsenic, (c) sulfur, and (d) total organic carbon content profiles of the QTC2 (red circles, Table S1) and QTC3 (open blue squares, Table S2) sediment samples as a function of depth. The term “clay” in the figure refers to the mineralogy.
22
ACS Paragon Plus Environment
Page 23 of 25
Environmental Science & Technology
599 600 601 602 603 604 605 606 607 608 609 610 611
2m
type I: 2,560 ± 360 yr BP
3m
type I
4m
type I
6m
type I
7m
type I
7m
type I
8m
type II
8m
type II
10 m
type II
2m
b
type I As in pyrite Thiol-bound As(III) Arsenite Arsenate
Fe in clays Fe in oxides Fe in pyrite Fe in siderite Fe in NOM
Depth
Depth
a
10 m
type II: 7,820 ± 890 yr
12 m
type III: 6,480 ± 790 yr
12 m
type III
14 m
type II
14 m
type II
16 m
type IV: 8,079 ± 708 cal. yr BP
16 m
17 m
type V
18 m
type V: 8,840 ± 1,670 yr
18 m
type V
20 m
type V
20 m
type V
type IV
0 20 40 60 80 As Concentration by Speciation (ppm)
0 1 2 3 4 5 Fe Concentration by Speciation (wt. %)
Figure 2. (a) Fe quantitative speciation in the depth range of 0-20 m and age of selected sediment layers. Quantitative speciation is calculated by multiplying the Fe concentration in QTC2 sediment samples (Table S1) and the fraction of Fe as each of the five species (Table S3) obtained by LCF of Fe K-edge EXAFS data (Figure S8). The Fe-rich clay fraction combines Fe in various reference clay minerals that are reported in Table S3. The age of selected sediment layers were determined using 14C or OSL methods and the analytical details were provided in the Supporting Information. (b) Arsenic quantitative speciation in the depth range 0-20 m. The As quantitative speciation is calculated by multiplying As concentration in QTC2 sediment samples (Table S1 and the fraction of As as each of the four species (Table S5) obtained by LCF of As K-edge EXAFS data (Figure 3).
23
ACS Paragon Plus Environment
Environmental Science & Technology
Page 24 of 25
612 1
a
10
b
c
Type I (2m) Type I (7m)
Type II (10m)
Type II (14m)
3
Type III (12m)
k χ(k)
Normalized absorption
Type II (8m)
Type IV (16m)
Type V (18m)
Type V (20m)
Type V (35m)
11850
613 614 615 616 617 618 619
11880 11910 Energy (eV)
4
6
8 10 -1 k (Å )
12
Figure 3. As K-edge XAS data. As K-edge XAS data collected at 10 K for selected QTC2 sediment samples and results of linear combination fit (LCF). Experimental and calculated curves are plotted as solid and dotted lines, respectively. (a) XANES data. (b) Unfiltered EXAFS k3χ(k) functions. (c) Fourier transform (FT) of unfiltered k3χ(k) functions in the k range of 3-12 Å-1. The LCF results are reported in Table S5.
24
ACS Paragon Plus Environment
Page 25 of 25
Environmental Science & Technology
Fe(III)
Fe oxides
flood water
O-As(V)
O-As(V)
reduc3ve dissolu3on
aqueous As
Py + NOM
S-As(III)
As
-O-As(V) + -O-As(III)
Sd + Py
S-As(III)
As
-O-As(V) + -O-As(III)
Py + NOM
S-As(III)
As
-O-As(V) + -O-As(III)
Py + NOM
S-As(III)
As
-O-As(V) + -O-As(III)
Sd + NOM
S-As(III)
-O-As(V) + -O-As(III)
Fe(II)
620 621 622 623 624 625 626 627 628 629 630 631 632 633
Py : pyrite Sd: siderite
: Fe oxides
: pyrite
: NOM
Alluvial deposi-on layer (type I)
Shallow-marine layer (types II, 7830 yr)
Shallow-marine layer (types III, 6480 yr)
Shallow-marine layer (types II)
Paleo-mangrove layer (types IV, 8079 cal. yr BP)
Pre-transgression layer (type V)
-S-As(III): sulfur-bound As(III) -O-As(V): oxygen-bound As(V) -O-As(III): oxygen-bound As(III)
Figure 4. Conceptual model for As speciation in Mekong Delta sediment. Flood water from the Mekong River brings As-bearing Fe(III) (oxyhydr)oxides. Before marine transgression, Fe (oxyhydr)oxides are buried during sediment deposition. Marine transgression at ∼8,079 yr yields a paleo-mangrove layer in which As is trapped in the form of pyrite- and NOMassociated As during diagenesis. Further extent of the marine transgression leads to shallow marine deposits in which As is trapped in arsenian pyrite and, to a lesser extent, by NOM. Arsenate and arsenite persist in these reduced paleo-layers. Marine regression after 5,000 BP resets the deposition of As-bearing Fe(III) (oxyhydr)oxides. Fe(III) (oxyhydr)oxides begin to be reduced after burial, releasing As into sediment porewater. The size of the pyrite symbol indicates its relative abundance in different sediment layers; the grayscale of the type II – IV layers represents the relative abundance of NOM. Note that the thickness of the sediment layers is not to scale.
25
ACS Paragon Plus Environment