Arsenic Speciation in Mekong Delta Sediments Depends on Their

Feb 16, 2018 - Faculty of Environment and Natural Resources, Ho Chi Minh City University of Technology, VNU-HCM, 268 Ly Thuong Kiet, Ho Chi Minh City,...
0 downloads 4 Views 3MB Size
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