Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES
Article
Fate of arsenic during Red River water infiltration into aquifers beneath Hanoi, Vietnam. Dieke Postma, Nguyen Thi Hoa Mai, Vi Mai Lan, Pham Thi Kim Trang, Helle Ugilt Sø, Pham Quy Nhan, Flemming Larsen, Pham Hung Viet, and Rasmus Jakobsen Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b05065 • Publication Date (Web): 13 Dec 2016 Downloaded from http://pubs.acs.org on December 21, 2016
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 30
1
Environmental Science & Technology
Revision submitted to EST
v3
2
3
Fate of arsenic during Red River water
4
infiltration into aquifers beneath Hanoi,
5
Vietnam.
6
7
8
Dieke Postma*1, Nguyen Thi Hoa Mai2 , Vi Mai Lan2, Pham Thi Kim Trang2, Helle Ugilt Sø1,
9
Pham Quy Nhan3, Flemming Larsen1, Pham Hung Viet2, Rasmus Jakobsen1
10 11
1
Geological Survey of Denmark and Greenland, Øster Voldgade 10, DK-1350 Copenhagen, Denmark.
12
2
Research Centre for Environmental Technology and Sustainable Development (CETASD), Hanoi
13
University of Science (VNU), Hanoi, Vietnam.
14
3
15
* corresponding author
[email protected] phone +45 30273702
Dept. of Hydrogeology, Hanoi University of Mining and Geology (HUMG), Hanoi, Vietnam.
16 17 1 ACS Paragon Plus Environment
Environmental Science & Technology
Page 2 of 30
18
ABSTRACT
19
Recharge of Red River water into arsenic contaminated aquifers below Hanoi was investigated. The
20
groundwater age at 40 m depth in the aquifer underlying the river was 1.3 ± 0.8 year, determined by
21
tritium-helium dating. This corresponds to a vertical flow rate into the aquifer of 19 m/yr. EC and
22
PCO2 indicate that water recharged from the river is present in both the sandy Holocene and gravelly
23
Pleistocene aquifers and is also abstracted by the pumping station. Infiltrating river water becomes
24
anoxic in the uppermost aquifer due to the oxidation of dissolved organic carbon. Further
25
downward, sedimentary carbon oxidation causes the reduction of As-containing Fe-oxides. Because
26
the release of arsenic by reduction of Fe-oxides is controlled by the reaction rate, arsenic entering
27
the solution becomes highly diluted in the high water flux and contributes little to the groundwater
28
arsenic concentration. Instead the As concentration in the groundwater of up to 1 µM is due to
29
equilibrium controlled desorption of arsenic, adsorbed to the sediment before river water started to
30
infiltrate due to municipal pumping. Calculations indicate that it will take several decades of river
31
water infiltration to leach arsenic from the Holocene aquifer to below the WHO limit of 10 µg/l.
32
33
2 ACS Paragon Plus Environment
Page 3 of 30
Environmental Science & Technology
34
INTRODUCTION
35
Many large Asian cities rely on groundwater for their water supply and the abstraction of
36
groundwater has commonly caused the development of large depression cones1-5. In addition, some
37
of these cities, exemplified by Kolkata6,7, and Hanoi8,9 are situated on top of aquifers that are
38
seriously contaminated by geogenic arsenic. Generally the uppermost Holocene aquifers are most
39
contaminated by arsenic while underlying Pleistocene aquifers in most cases are low in arsenic10,11.
40
Because groundwater abstraction for municipal water supply occurs from the deeper Pleistocene
41
aquifers, with groundwater low in arsenic, the change in groundwater flow patterns may cause the
42
migration of arsenic rich groundwater from upper aquifers into the underlying low arsenic
43
aquifers8,9,12,13. The effect of municipal groundwater abstraction pumping also influences the water
44
quality on a regional scale because depression cones extend far beyond the city limits5,14 .
45
Both Hanoi and Dhaka are located along major rivers and hydrogeological studies
46
have shown1,2,15,16 groundwater abstraction to cause a very significant infiltration of river water into
47
the underlying aquifers. Recharge of aquifers by losing rivers, often called bank infiltration,
48
enhances the available water resource and managed and monitored bank infiltration is a strategy
49
that has been used extensively particularly in Europe17-20. Apart from increasing the quantity of the
50
groundwater resource, bank infiltration also may improve the quality of the infiltrating surface
51
water. During bank infiltration, processes taking place at the surface water-groundwater interface
52
can remove nutrients, organic carbon, including organic contaminants, and pathogens such as
53
microbes and viruses21,22. However, the effect of Red River bank infiltration on the water quality of
54
the aquifers below Hanoi is completely unknown even though it is considered as the main source of
55
groundwater recharge below the city2.
3 ACS Paragon Plus Environment
Environmental Science & Technology
56
Page 4 of 30
The city of Hanoi has a leaky sewerage system as well as open channels carrying
57
wastewaters and untreated sewage water through the city, both of which are leaking into the
58
underlying aquifers12,23,24. In addition the aquifers beneath Hanoi are contaminated by arsenic with
59
concentrations exceeding 300 µg/L8,9. In contrast the water of the Red River running along the city
60
of Hanoi is very low in dissolved arsenic12,25,26. Therefore one could hope that the large scale
61
infiltration of river water2 can improve the water quality in the aquifers and also flush out the
62
arsenic. The time scale needed for flushing arsenic from the aquifer is currently unknown but must
63
depend on infiltration rates, groundwater flow patterns and residence times as well as the
64
geochemical processes controlling the groundwater arsenic content. The latter includes the
65
mobilization of arsenic by reductive dissolution of Fe-oxides and the adsorption and desorption of
66
the mobilized arsenite12,25,27. The objective of the present study is to investigate the rate of
67
infiltration of Red River water into the underlying aquifer, clarify the geochemical processes
68
occurring between the infiltrating water and the aquifer sediment and predict the effect of river
69
water infiltration on the groundwater arsenic content.
70
71
MATERIALS AND METHODS
72
Field Area Description. Our field locality at Nam Du (20o58´25.43”N, 105o54’09.73”E) is located
73
on the eastern boundary of the city of Hanoi, Vietnam, directly on the shore of the Red River. The
74
geology of the area has been investigated extensively in relation with the establishment of the Nam
75
Du well field in 200412. Generally, the geological sequence in the area consist of a top clay layer
76
underlain by a up to 40 m thick Holocene aquifer consisting of fine to medium grained sand,
77
comprising some more fine or coarse grained layers. The sandy aquifer is underlain by an aquitard
78
of variable thickness (1-10 m) followed by a coarse grained sand to gravel aquifer of Pleistocene
4 ACS Paragon Plus Environment
Page 5 of 30
Environmental Science & Technology
79
age12. Remote sensing shows the Nam Du area to belong to the meander belts of the Red River28. A
80
detailed geological profile of our study area, based on natural gamma logging of boreholes and
81
sediment description is shown in Fig.1. Uppermost there is an up to 11 m thick layer of clay and
82
silty clay which is absent close to the river. Below follows an up to 40 m layer of fine to medium
83
sized sand. Underneath there is an up to 20 m thick layer of coarse sand and gravel. Between the
84
two sand formations a clay layer is found in some borings but not in all. Finally at the bottom there
85
is a layer of claystone-siltstone. In accordance with the geological description of the area12, the
86
upper aquifer with fine grained sand is considered to be of Holocene age while the underlying
87
coarse sand-gravel layer is interpreted as being of Pleistocene age. The intermediate clay lenses are
88
probably remnants of the aquitard separating the Holocene and the Pleistocene aquifers which
89
elsewhere is found as a more continuous layer12. The claystone-siltstone corresponds to the
90
lithology of the Neogene in the area.
91
The Red River at this location is about 400 m wide. The depth of the river channel at
92
Nam Du varies between 13.5 and 17.0 m below the ground level at the nest of boreholes on the
93
shore (Fig.1). The mid channel depth is 15.3 m and is situated within the Holocene aquifer. The
94
groundwater below Hanoi has been exploited since 1894 and groundwater abstraction presently
95
amounts to > 900,000 m3/day27. Groundwater is being pumped mainly from six major well fields
96
distributed throughout the city2. As a result, large coalescing depression cones have developed
97
below the city of Hanoi in both the Holocene and Pleistocene aquifers which are partially
98
connected8,12. The well field at Nam Du comprises a gallery of pumping stations close to the river
99
(Norrman 2008). The pumping station included in our field site (Fig 1) is part of this gallery of
100
pumping stations but the flow towards the pumping station is the result of the drawdown caused by
101
the entire pumping activity below Hanoi12. The hydrostatic pressure in our wells screened in the
102
Pleistocene aquifer is about 6 m lower than those in the nest screened in the Holocene aquifer (Fig. 5 ACS Paragon Plus Environment
Environmental Science & Technology
103
1), indicating a strong downward gradient induced by the pumping. The pumping combined with
104
the less permeable Holocene deposit on top of the high permeable Pleistocene causes
105
predominantly vertically downward groundwater flow, providing the shortest passage through the
106
Holocene. The pressure difference between the Holocene and Pleistocene aquifers is also found
107
elsewhere in the city but its magnitude varies from locality to locality2,12.
108
Groundwater sampling and analysis. To investigate the hydrogeological and geochemical
109
processes we have collected three sets of data (Fig. 1); A) Detailed water sampling of the sediment
110
porewater performed on a sandbank in the river. B) Water sampling in a nest of boreholes located
111
on the river shore with screens installed at different depths in the Holocene aquifer, C) Water
112
sampling in boreholes and the nearest pumping station screened in the Pleistocene aquifer.
113
Page 6 of 30
To sample the interface between river and groundwater we used a 1 cm diameter
114
drive-point piezometer, and a 50 mL syringe for extracting the water. Sampling was carried out in
115
shallow water, approximately 25 cm deep, on a sandbank in the river located about 50 m from the
116
shore, with a depth resolution of 0.2-0.5 m until 3 m below the surface (Fig. 1). At each depth the
117
internal volume of the piezometer was emptied at least three times before taking the sample. On the
118
shore a nest of 11 wells (Fig 1) was installed in the depth range 4.9-40 m. Additionally, four wells
119
ranging in depth from 40 to 60 m were installed in the Pleistocene aquifer between the pumping
120
station and the nest at the shore (Fig. 1). The wells were equipped with Ø60 mm PVC-casings, a 1
121
m long screen and 1 m sand trap. A quartz sand filter pack was installed, and the well was sealed
122
using bentonite. The top end of the PVC-casing was sealed to prevent the entrance of surface water
123
during flooding. To further avoid leakage along the casing, the top of the wells were embedded in a
124
concrete platform, each well within a protective steel casing and steel screw cap. Directly after
125
completion the wells were pumped to remove the water affected by the drilling operation. Sampling
6 ACS Paragon Plus Environment
Page 7 of 30
Environmental Science & Technology
126
of the wells was done more than three months after their installation in 2010-2011. The pumping
127
well is screened 48-65 m below the surface.
128
Groundwater was sampled from the boreholes using a Grundfos MP1 submersible
129
pump. Five borehole volumes were flushed before taking the sample. O2, pH, and electrical
130
conductivity (EC) were measured with probes in a flow cell mounted directly on the sampling tube.
131
Samples for other parameters were collected in 50 mL syringes and filtered through 0.2 µm
132
Sartorius Minisart cellulose acetate filters. Aqueous As(V) and As(III) were separated by passing
133
the water through a disposable anion exchange cartridge29. The cartridge was together with the
134
filter, carefully flushed with N2 before use. Arsenite was determined as the As concentration in the
135
water filtered through a cartridge, and As(V) was calculated as the difference between the total As
136
and As(III) concentrations. In the field Fe(II) was measured spectrophotometrically and alkalinity
137
determined by Gran-titration. Samples for cations were preserved with 2% of a 7 M HNO3 solution
138
and refrigerated until analysis in the laboratory. Samples for NH4+, Cl- and SO42- were collected in
139
20-mL polypropylene vials and frozen immediately after sampling. In the laboratory cations were
140
analyzed by flame atomic absorption spectrophotometry. Arsenic was determined on the same
141
instrument using a HVG hydride generator and a graphite furnace. Anions were analyzed by ion
142
chromatography using a Shimadzu LC20AD/HIC-20A Super. Detection limits were; As 0.013 µM,
143
Mn 0.91 µM, Ca 0.50 µM, NH4 5.6 µM, PO4 1.1 µM, NO3 3.2 µM, SO4 2.1 µM. For further details
144
on procedures and methods consult Postma et al.25.
145
Samples for tritium-helium dating were taken in duplicate in clamped copper tubes in
146
2015 and analyzed by the Institut für Umweltphysik, Universität Bremen. 3H concentrations were
147
determined by the 3He in growth method. The accuracy of the dates is given by the standard
148
deviation as ± 0.8 years (Jürgen Sültenfuß, pers. comm.)
7 ACS Paragon Plus Environment
Environmental Science & Technology
Page 8 of 30
149
150
RESULTS AND DISCUSSION
151
River Water Infiltration. Tritium-helium dating of water sampled from the boreholes in the nest
152
screened in the Holocene aquifer (Fig. 1) is shown in Fig.2. The groundwater is extremely young in
153
screens at all depths indicating fast infiltration of river water into the underlying aquifer. The oldest
154
groundwater age of 8 years is found at 30 m depth. This piezometer is located in one corner of the
155
well cluster and here a clay layer (not shown in Fig. 1) was encountered during drilling from 26 to
156
28 m depth which was not present in the other drillings just a few meters away. Below this clay
157
layer, the groundwater apparently has a larger residence time yielding a higher groundwater age.
158
The 30 m depth screen also shows a maximum in the dissolved As(III) concentration (Fig. 3) and in
159
the NH4 concentration (not shown) reflecting that the longer groundwater residence time allows
160
more As(III) and NH4 to build up. The remaining screens yield groundwater ages ranging from less
161
than one year to 1.3 year in the screen at 40 m depth. Duplicate measurements are, except for the 30
162
m screen (Fig. 2), within 0.05 year indicating a high precision. However, the estimated accuracy of
163
± 0.8 year is high compared to the measured ages. In addition there is uncertainty as to whether
164
infiltration occurs from the deeper or the more shallow parts of the channel. Using the vertical
165
distance between the mean channel depth of 15.3 m and the groundwater age of 1.3 ± 0.8 years at
166
40 m yields a vertical transport rate of 19 m/yr with an uncertainty range of 12-49 m/yr.
167
A further indicator for river water intruding into the aquifer is the electrical
168
conductivity (EC) of the groundwater (Fig. 2). Shallow groundwater on the shores of the Red River
169
30 km NW of Hanoi, unaffected by pumping, has an EC in the range 750-1000 µS/cm25. River
170
water, however, has an EC of around 200 µS/cm and the difference is due to mineral weathering
171
mainly occurring in the soil zone. The strong contrast in EC between river water and normal
8 ACS Paragon Plus Environment
Page 9 of 30
Environmental Science & Technology
172
groundwater makes the EC a good indicator for intrusion of river water into aquifers. As shown in
173
Fig. 1 and 2, river water, groundwater in the Holocene and Pleistocene aquifers, as well as water of
174
the pumping station all have an EC in the range 200-320 µS/cm indicating rapid recharge from the
175
river to be the dominant source of water in the aquifer.
176
An additional useful tracer parameter is the partial pressure of CO2. River water has a
177
low PCO2 = 10-3.3, close to that in the atmosphere (PCO2 = 10-3.4). However water infiltrating through
178
soils, or through organic rich sediment, is exposed to root respiration and organic decay that
179
generate CO2. At our field site 30 km upstream along the Red River, the groundwater therefore has
180
a PCO2 in the range 10-1.0 - 10-1.4 ref 25. Fig 2 shows how all the water in the aquifer below the river
181
has a PCO2 in the range 10-2.3 to 10-2.8. These low partial pressures of CO2 also indicate the recharge
182
of river water into the aquifer without being exposed to extensive organic matter decay or
183
respiration. This is consistent with fast recharge of water through the bottom of the river.
184
Redox Processes. The river water has an O2 concentration of 0.24 mmol/L, only slightly below the
185
concentration of 0.27 mmol/L in equilibrium with the atmosphere30. In addition river water contains
186
28 µmol/L of nitrate. At a depth of a few meters below the surface, both oxygen and nitrate have
187
disappeared from the groundwater (Fig. 3). At the same time the river water, and one sample at 1 m
188
depth, contain about 0.6 mmol/L dissolved organic carbon (DOC), while the remaining
189
groundwater has a DOC concentration near 0.1 mmol/L. As the river water infiltrates, the DOC
190
appears to be filtered by the uppermost sediment. Stoichiometric calculation suggest that the
191
oxidation of 0.27 mmol/L DOC would be sufficient for the complete removal of dissolved oxygen
192
and nitrate from the infiltrating river water. There is no evidence for Fe-oxides functioning as
193
electron acceptor for DOC oxidation in the top sediment layers because Fe(II) appears in the
194
groundwater at greater depth. This pattern of DOC filtration by the uppermost sediment layer and
195
concomitant removal of dissolved O2 and nitrate, was also observed at other bank infiltration 9 ACS Paragon Plus Environment
Environmental Science & Technology
196
sites18,19,21,31. Below 10 m depth the Fe(II) concentration starts to increase in the infiltrating water
197
(Fig. 3), probably because the degradation of sedimentary organic carbon causes reductive
198
dissolution of Fe-oxides. Similar Fe(II) increases attributed to reductive dissolution of Fe-oxides
199
have been observed at other sites with ongoing bank infiltration18.21.
Page 10 of 30
200
The distribution of arsenic in the river and groundwater is displayed in Fig. 1 and in
201
detail for the Holocene aquifer in Fig. 3. The river water and the shallow groundwater are free of
202
arsenic, while it is present in the deeper groundwater of the Holocene, the Pleistocene aquifers and
203
the water of the pumping station. This distribution indicates that the infiltrating river water is
204
mobilizing arsenic from the aquifer sediment. Postma et al.32,33 found the rate of arsenic
205
mobilization from the sediment into the groundwater to depend on the burial age of the aquifer
206
sediment with the youngest sediments being most reactive. The geochemical processes controlling
207
the groundwater arsenic content were quantified in a 1-D reactive transport model, up to a burial
208
age of 6000 years, using the code PHREEQC-334. In the model the oxidation of organic matter
209
causes the reduction of arsenic-containing Fe-oxides, resulting in the release of Fe(II) and As(III) to
210
the groundwater. Both organic carbon and Fe-oxides become consumed over time and less reactive
211
as well. The oxidation of organic matter and the reduction of Fe-oxides are controlled by separate
212
rate equations that quantify their decreasing rates over time (Supporting Information Eqn. 1 and 2).
213
The rate equations were calibrated using the groundwater chemistry of four sites on the Red River
214
floodplain 30 km NW of Hanoi33, where both the burial age of the Holocene aquifer sediment and
215
the vertical flow rate of the groundwater are known.
216
The application of the model developed for the field area 30 km N of Hanoi to the
217
Nam Du field site assumes that floodplain sediment is similar at the two sites. This appears a
218
reasonable assumption given that remote sensing28 shows both sites to belong to the Red River
219
meander belts and they also have comparable Holocene sediments. Accordingly we have applied 10 ACS Paragon Plus Environment
Page 11 of 30
Environmental Science & Technology
220
the model to the Nam Du data (see Supporting Information for the PHREEQC input file). The main
221
difference between the two field sites is the much faster vertical flow rate at Nam Du of 19 m/yr
222
compared to 0.5 m/yr in the undisturbed aquifers north of Hanoi32. The burial age of aquifer
223
sediments adjacent to the river at Nam Du is not known but as end-member case we assume the
224
aquifer to be very young (< 100 yr). With that assumption the organic carbon and iron oxides have
225
the highest possible reactivity and calculated release rates of Fe(II) and As(III) will accordingly be
226
maximum rates. Any older sediment would consequently show lower Fe(II) and As(III) release
227
rates. The composition of the solution entering into the 1-D model is that of river water. The results
228
calculated by the model for 70 years of infiltration are shown as the solid lines in Fig 3 but would
229
be near identical for other times < 100 yr. To simulate in addition the consumption of DOC in the
230
uppermost sediment we have added a constant rate of organic carbon degradation, of 5.1 mM C/yr
231
in the upper layer of the model, matching the removal of O2 and NO3 in the incoming water by
232
DOC oxidation (Fig. 3). Below the top layer the oxidation rate of sedimentary organic carbon, is
233
calculated from the rate equation and the main electron accepting process is the reductive
234
dissolution of Fe-oxides. Fig. 3 shows the modelled release of Fe(II) to the groundwater which
235
corresponds to a rate of 35 µM Fe/yr. In the bottom 4 m of the profile a small amount of Fe(II)
236
becomes precipitated as siderite but higher up in the profile all Fe(II) stays in solution. Comparison
237
with field data shows the modelled values to be slightly higher than the field data but the
238
assumption that the Holocene aquifer at Nam Du is young appears reasonable. Otherwise the
239
modelled Fe(II) concentration would be much higher than the measured Fe(II) values because the
240
Fe(II) release of an old low reactive sediment would be significantly lower.
241
Fe-oxides being reduced in the model contain 1.2 mmol As(V) per mol of Fe, a value
242
based on analysis of Fe-oxides by reductive dissolution in sediment samples from four sites with a
243
different burial age of the aquifer sediment on the Red River floodplain33. During the reduction of
11 ACS Paragon Plus Environment
Environmental Science & Technology
Page 12 of 30
244
Fe-oxides in the aquifer, As(V) becomes reduced to As(III) and is released to the groundwater.
245
Using the model we can evaluate how much of the As measured in the groundwater could possibly
246
be derived from reductive dissolution of As-containing Fe-oxides from sediment that is considered
247
to be < 100 yr old. For any older sediment the arsenic release rate would be lower. To keep all
248
As(III) released from reduced Fe-oxides in solution we have disabled the adsorption of As(III) in
249
the model. Figure 3 compares the measured As(III) in the groundwater with the model prediction
250
assuming reductive dissolution of As-containing Fe-oxide being the only source of As(III) in the
251
water. The modelled As profile in Fig. 3 corresponds to an As(III) release rate of 0.052 µM
252
As(III)/yr. For the 19 m/yr vertical flow rate, the modelled As(III) concentration at 40 m depth is
253
0.11 µM while the uncertainty range of 12-49 m/yr in the vertical flow rate corresponds to a As(III)
254
concentration range of 0.04-0.17 µM.
255
Fig. 3 clearly shows the modelled As concentration, released by reductive dissolution
256
of Fe-oxides in 70 year old sediment, to be much lower than the As concentration measured in the
257
aquifer. Any older sediment would have a lower reactivity and therefore a lower As(III) release
258
rate. Therefore most of the As(III) present in the groundwater cannot be released by reductive
259
dissolution of As-containing Fe-oxides. Note that the same release rate of 0.052 µM As(III)/yr
260
yields a modelled groundwater arsenic concentration of up to 4 µmol As(III)/L at the Dan Phuong
261
locality upstream along the Red River33. Because the reduction of Fe-oxides, and thereby the release
262
of As, is kinetically controlled the resulting groundwater concentration will depend on the
263
groundwater flow rate35. At the Dan Phuong site the downward groundwater velocity is 0.5 m/yr25
264
while at Nam Du it is 19 m/yr. The same As(III) release rate therefore produces groundwater with a
265
As(III) concentration that is almost 40 times higher at Dan Phuong than at Nam Du! The
266
importance of reaction kinetics on arsenic mobilization is also evident from the much higher arsenic
267
concentration found in porewater in sediment along the shores of the Red River, 6 km further
12 ACS Paragon Plus Environment
Page 13 of 30
Environmental Science & Technology
268
downstream at Van Phuc26. Vertical flow rates here are 0.1-3.1 m/yr26, which is an order of
269
magnitude lower than the vertical flow rate we derive for the bottom of the river channel, and
270
therefore result in a higher arsenic concentration. In fact, the main role of redox processes for
271
arsenic mobility in the Nam Du aquifer is to maintain anoxic conditions in the aquifer, thereby
272
keeping arsenic in its As(III) form. If the oxidation of arsenic to As(V) should occur it would cause
273
a much stronger arsenic retardation because As(V) adsorbs much stronger to the aquifer sediment36.
274
As(III) retardation. Groundwater must have been high in As(III) before water supply pumping by
275
the city of Hanoi induced river water infiltration into the aquifer, because high As(III)
276
concentrations are seen everywhere else on the Red River floodplain in aquifers contained in young
277
sediments25,32. Dissolved As(III) in the groundwater is in equilibrium with As(III) adsorbed onto the
278
sediment and the results of a forced gradient experiment37 and push pull tests38 indicate that arsenic
279
sorption rapidly responds to changes in the water chemistry. The infiltration of low arsenic river
280
water is therefore expected to cause desorption of arsenic from the aquifer sediment. The
281
distribution of As(III) in the groundwater over depth, with low As(III) in the upper 15 m and higher
282
As(III) further down (Fig. 3), is consequently interpreted as a leaching front and its position will be
283
determined by the number of pore volumes of groundwater flushed through the aquifer and the
284
adsorption isotherm.
285
In a previous study36 we experimentally investigated the adsorption of As(III) onto the
286
Holocene aquifer sediment at the Nam Du field site using a sediment sample from 5.5-6 m depth.
287
The results showed a linear adsorption isotherm up to an aqueous concentration of at least 10 µM
288
As(III) corresponding to a K’d of 12.6 L/kg. The adsorption of As(III) was also found to be
289
completely reversible. For the measured aqueous concentration of 1 µM As(III) (Fig. 3) and a
290
porosity of 0.439 the isotherm predicts an adsorbed concentration of 59 µmol As(III) per liter of
291
contacting groundwater corresponding to a retardation of 60. Accordingly it will take the flushing 13 ACS Paragon Plus Environment
Environmental Science & Technology
Page 14 of 30
292
of 60 pore volumes of groundwater through the Holocene aquifer to desorb all adsorbed As(III).
293
Starting with an initial concentration of 1 µM As(III) in the groundwater, close to what is measured
294
in the groundwater, and the corresponding adsorbed As(III) concentration of 59 µmol/L, we have
295
calculated how the As(III) versus depth profile would look like at different times and compared it to
296
the measured As(III) concentration (Fig. 4). As the water moves downward the As(III)
297
concentration increases along the isotherm until it reaches the concentration that was present in the
298
aquifer before river water infiltration began. Because the adsorption isotherm is linear, the shape of
299
the leaching front will be independent of the As(III) concentration. In the model, the vertical flow
300
rate is 19 m/yr and the leaching front moves downward with 0.32 m/yr and after 70 years of
301
modelled river water infiltration there is a reasonable correspondence with the field data (Fig. 4).
302
Using the uncertainty in the infiltration rate from 12 to 49 m/yr the range in front progression is
303
0.20 to 0.82 m/yr and the time for the leaching front to reach the same depth from 112 to 27 years.
304
Although groundwater abstraction in Hanoi started 100 years ago8 we do not know when river
305
water started to infiltrate into the aquifer at Nam Du. Van Geen et al.27 constructed a flow model for
306
Van Phuc, on the next bend of the river 6 km further south. Here the pumping of the Hanoi water
307
supply has caused a flow reversal from a gaining to a losing river. The timing of the flow reversal
308
was estimated to be between 1951 and 1971. Our estimate of 70 ± 42 years of river water
309
infiltration is within the time space of groundwater abstraction by the Hanoi water supply but
310
clearly associated with considerable uncertainty.
311
Environmental Perspective. Groundwater flow modeling2 suggest that infiltrating Red River water
312
is the main source of groundwater recharge into the aquifers beneath Hanoi which are heavily
313
pumped by the municipal water supply. Infiltrating river water may be beneficial for the
314
groundwater quality in at least two ways. First of all, the river water contains no arsenic and one
14 ACS Paragon Plus Environment
Page 15 of 30
Environmental Science & Technology
315
may hope that its infiltration will lower the groundwater arsenic content over time. The second
316
effect will be dilution of the sewage leaking into the aquifer23,24.
317
It is of major importance to understand the hydrogeological and geochemical
318
processes occurring at the river-groundwater interface to assess the future water quality in the
319
aquifers below Hanoi. Our data from Nam Du suggest that very fast recharge is taking place from
320
the river channel. However another recent study at the Van Phuc site26 located 6 km downstream
321
from our site at Nam Du suggests much slower infiltration to occur along the shores of the river. In
322
future work it would be important to delineate the relative importance of the two modes of recharge,
323
particularly because they seem to result in a different groundwater arsenic content. While we find at
324
Nam Du that due to the high infiltration the kinetically controlled reductive dissolution of As-
325
containing Fe-oxides only contributes insignificantly to the groundwater arsenic content, the study
326
at Van Phuc suggest that the slow infiltration rates allow much higher arsenic concentrations to
327
build-up in the groundwater26. An additional reason for the much higher arsenic concentrations
328
found in porewater of the shoreface sediment could also be that Van Phuc is located just
329
downstream of the channel at Yen My which drains most of the untreated urban wastewater from
330
the city of Hanoi into the Red River24, providing an ample source of reactive carbon.
331
While the infiltrating river water may readily displace the original groundwater as the
332
result of pumping, our results indicate that a decrease in the arsenic concentration in the pumped
333
water will take much longer. Our data indicates that it probably will take another 60 ± 36 years to
334
flush the adsorbed As(III) from the Holocene aquifer. To reach the screens of the groundwater
335
abstraction wells the groundwater thereafter flows through Pleistocene deposits. Van Geen et al.27
336
determined the retardation of As(III) in Pleistocene sands to be 16-20 which is a factor 3-4, lower
337
than we find in the Holocene. In addition many Pleistocene sediment layers consist of highly
338
permeable gravels that probably have an even lower retardation. Flushing of arsenic from the 15 ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 30
339
Pleistocene aquifers is therefore expected to proceed much faster than from the Holocene aquifers.
340
Thus bank infiltration of Red River water into the Holocene and Pleistocene aquifers will likely
341
cause the groundwater arsenic concentration to decrease but it will take several decades.
342
343
SUPPORTING INFORMATION
344
Description of the reactive transport model and model input file
345
346
ACKNOWLEDGEMENT
347
This research has been funded by the European Research Council under the ERC Advanced Grant
348
ERG-2013-ADG. Grant Agreement Number 338972.
349
350
REFERENCES
351
(1) Hoque, M. A.; Hoque, M. M.; Ahmed, K. M. Declining groundwater level and aquifer
352
dewatering in Dhaka metropolitan area, Bangladesh: causes and quantification. Hydrogeol. J.
353
2007,15, 1523–1534.
354 355
(2) Jusseret, S.; Tam, V. T.; Dassargues, A. Groundwater flow modeling in the central zone of
356
Hanoi, Vietnam. Hydrogeol. J. 2009, 17, 915-934.
357
16 ACS Paragon Plus Environment
Page 17 of 30
Environmental Science & Technology
358
(3) Haque, S. J.; Onodera, S.; Shimizu, Y. An overview of the effects of urbanization on the
359
quantity and quality of groundwater in South Asian megacities. Limnology, 2013, 14, 135-145.
360 361
(4) Zhou, Y.; Dong, D.; Liu, J.; Li, W. Upgrading a regional groundwater level monitoring network
362
for Beijing Plain, China. Geoscience Frontiers 2013, 4, 127-138.
363 364
(5) Khan, M. R.; Koneshloo, M.; Knappett, P. S. K.; Ahmed, K. M., Bostick, B. C.; Mailloux, B. J.;
365
Mozumder, R. H.; Zahid, A.; Harvey, C. F.; van Geen, A.; Michael, H. A. Megacity pumping and
366
preferential flow threaten groundwater quality. Nature Comm. 2016, DOI: 10.1038/ncomms12833
367 368
(6) Chakraborti, D.; Das, B.; Rahman, M. M.; Chowdhury, U. K.; Biswas, B.; Goswami, A. B.;
369
Nayak, B.; Pal, A.; Sengupta, M. K.; Ahamed, S.; Hossain, A.; Basu, G.; Roychowdhury, T.; Das,
370
D. Status of groundwater arsenic contamination in the state of West Bengal, India: A 20-year study
371
report. Mol. Nutr. Food Res. 2009, 53, 542 – 551.
372 373
(7) Malakar, A.; Islam, S.; Ali, M. A.; Ray, S. Rapid decadal evolution in the groundwater arsenic
374
content of Kolkata, India and its correlation with the practices of her dwellers. Environ. Monit.
375
Assess. 2016, 188: 584, DOI 10.1007/s10661-016-5592-9.
376 377
(8) Berg, M.; Trang, P. T. K.; Stengel, C.; Buschmann, J.; Viet, P. H.; Dan, N. V.; Giger, W.;
378
Stuben, D. Hydrological and sedimentary controls leading to arsenic contamination of groundwater
379
in the Hanoi area, Vietnam: the impact of iron arsenic ratios, peat, river bank deposits, and
380
excessive groundwater abstraction. Chem. Geol. 2008. 249, 91–112.
381
17 ACS Paragon Plus Environment
Environmental Science & Technology
Page 18 of 30
382
(9) Winkel L. H. E.; Pham, T. K. T.; Vi, M. L.; Stengel, C.; Aminia, M.; Nguyen, T. H.; Viet, P. H.;
383
Berg, M. Arsenic pollution of groundwater in Vietnam exacerbated by deep aquifer exploitation for
384
more than a century. P. Natl. Acad. Sci. USA (PNAS) 2011, 108, 1246-1251.
385 386
(10) DPHE. Arsenic contamination of groundwater in Bangladesh. In BGS Technical Report
387
WC/00/19; Kinniburgh, D. G., Smedley, P. L., Eds.; British Geological Survey: Keyworth, UK,
388
2001.
389 390
(11) Zheng, Y.; van Geen, A.; Stute, M.; Dhar, R.; Mo, Z.; Cheng, Z.; Horneman, A.; Gavrieli, I.;
391
Simpson, H. J.; Versteeg, R.; Steckler; Grazioli-Venier, A.; Goodbred, S.; Shahnewaz, M.;
392
Shamsudduha, M.; Hoque, M. A.; Ahmed, K. M. Geochemical and hydrogeological contrasts
393
between shallow and deeper aquifers in two villages of Araihazar, Bangladesh: Implications for
394
deeper aquifers as drinking water sources. Geochim. Cosmochim. Acta 2005, 69, 5203−5218.
395 396
(12) Norrman, J.; Sparrenbom, C. J.; Berg, M.; Dang, D. N.; Pham, Q. N.; Roskvist, H.; Jacks, G.;
397
Sigvardsson, E.; Baric, D.; Moreskog, J.; Harms-Ringdahl, P.; Nguyen, V. H. Arsenic mobilization
398
in a new well field for drinking water production along the Red River, Nam Du, Hanoi. Appl.
399
Geochem. 2008, 23, 3127–3142.
400 401
(13) Mukherjee, A.; Fryar, A. E.; Scanlon, B. R.; Bhattacharya, P.; Bhattacharya, A. Elevated
402
arsenic in deeper groundwater of the western Bengal basin, India: Extent and controls from regional
403
to local scale. Appl. Geochem. 2011, 26, 600–613.
404
18 ACS Paragon Plus Environment
Page 19 of 30
Environmental Science & Technology
405
(14) Knappett, P. S. K.; Mailloux, B. J.; Choudhury, I.; Khan, M. R.; Michael, H. A.; Barua, S.;
406
Mondal, D. R.; Steckler, M.S.; Akhter, S. H.; Ahmed, K. M.; Bostick, B.; Harvey, C. F.;
407
Shamsudduha, M.; Shuai, P.; Mihajlov, I.; Mozumder, R.; van Geen, A. Vulnerability of low-
408
arsenic aquifers to municipal pumping in Bangladesh. J. Hydrol. 2016,
409
doi.org/10.1016/j.jhydrol.2016.05. 035 .
410 411
(15) Darling, W. G.; Burgess, W G.; Hasan, M. K. Isotopic evidence for induced river recharge to
412
the Dupi Tila aquifer in the Dhaka urban area, Bangladesh. In: The application of isotope
413
techniques to the assessment of aquifer systems in major urban areas. IAEA-TECDOC, 2002, 1298,
414
IAEA, Vienna, pp 95–107.
415 416
(16) Hoque, M. A.; McArthur, J. M.; Sikdar, P. K.; Ball, J. D.; Molla, T. N. Tracing recharge to
417
aquifers beneath an Asian megacity with Cl/Br and stable isotopes: the example of Dhaka,
418
Bangladesh. Hydrogeol. J. 2014, 22, 1549–1560.
419 420
(17) Eckert, P.; Irmscher, R. Over 130 years of experience with Riverbank filtration in Düsseldorf,
421
Germany. J. Water Supply Res. T. 2006, 55, 283–291.
422 423
(18) Massmann, G.; Nogeitzig, A.; Taute, T.; Pekdeger, A. Seasonal and spatial distribution of
424
redox zones during lake bank filtration in Berlin, Germany. Environ. Geol. 2008, 54, 53–65.
425 426
(19) Bourg, A. C. M.; Bertin, C. Biogeochemical processes during the infiltration of river water into
427
an alluvial aquifer. Environ. Sci. Technol. 1993, 27, 661–666.
428
19 ACS Paragon Plus Environment
Environmental Science & Technology
429
(20) von Gunten, H. R.; Karametaxas, G.; Krahenbuhl, U.; Kuslys, M.; Giovanoli, R.; Hoehn, E.;
430
Keil, R. Seasonal biogeochemical cycles in riverborne groundwater. Geochim. Cosmochim. Acta
431
1991, 55, 3597–3609.
Page 20 of 30
432 433
(21) Farnsworth, C.E.; Hering, J. G. Inorganic geochemistry and redox dynamics in bank
434
infiltration settings. Environ. Sci. Technol. 2011, 45, 5079-5087.
435 436
(22) Hamann, E.; Stuyfzand, P. J.; Greskowiak, J.; Timmer, H.; Massmann, G. The fate of organic
437
micropollutants during long-term/long-distance river bank filtration. Sci. Total Environ. 2016, 545–
438
546 , 629–640.
439 440
(23) Trafford, J. M.; Lawrence, A. R.; Macdonald, D. M. J.; Nguyen Van Dan; Tran, D. N.; Ha, N.
441
T. The effect of urbanisation on the groundwater quality beneath the city of Hanoi. Technical
442
Report WC/96/22, 1996, British Geological Survey, Keyworth, Nottingham, UK.
443 444
(24) Norrman, J.; Sparrenbom, C. J.; Berg, M.; Dang, D. N.; Jacks, G.; Harms-Ringdahl, P.; Pham,
445
Q. N.; Rosqvist, H. Tracing sources of ammonium in reducing groundwater in a well field in Hanoi
446
(Vietnam) by means of stable nitrogen isotope (d15N) values. Appl. Geochem. 2015, 61, 248–258.
447 448
(25) Postma, D.; Larsen, F.; Nguyen, T. M. H.; Mai, T. D.; Pham, H. V.; Pham, Q. N.; Jessen, S.
449
Arsenic in groundwater of the Red River floodplain, Vietnam: Controlling geochemical processes
450
and reactive transport modeling. Geochim. Cosmochim. Acta 2007, 71, 5054–5071.
451
20 ACS Paragon Plus Environment
Page 21 of 30
Environmental Science & Technology
452
(26) Stahl, M. O.: Harvey, C. F.; van Geen, A.; Sun, J.; Pham, T. K. T.; Vi, M. L.; Thao, M. P.;
453
Pham, H. V.; Bostick, B. C. River bank geomorphology controls groundwater arsenic
454
concentrations in aquifers adjacent to the Red River, Hanoi Vietnam. Water Resour. Res. 2016, 52,
455
doi:10.1002/2016WR018891.
456 457
(27) van Geen, A.; Bostick, B. C.; Pham, T. K. T.; Vi, M. L.; Nguyen, N. M.; Phu, D. M.; Pham, H.
458
V.; Radloff, K.; Aziz, Z.; Mey, J. L.; Stahl, M.O.; Harvey, C.F., Oates, P., Weinman, B., Stengel,
459
C.; Frei, F.; Kipfer, R.; Berg, M. Retardation of arsenic transport through a Pleistocene aquifer.
460
Nature 2013, 501, 204-207.
461 462
(28) Mathers, S.; Zalasiewicz, J. Holocene sedimentary architecture of the Red River Delta,
463
Vietnam. J. Coastal Res.1999, 15, 314-325.
464 465
(29) Meng, X.; Wang, W. Speciation of arsenic by disposable cartridges. In: Proc. Third Int.Conf.
466
Arsenic Exposure and Health Effects. 1998, Soc. Environ. Geochem. and Health, Univ. Colorado at
467
Denver.
468 469 470
(30) Appelo, C. A. J.; Postma D. Geochemistry Groundwater and Pollution, 2nd ed. 2005, Balkema
471
Publ. 649 pp.
472 473
(31) Jacobs, L. A.; von Gunten, H. R.; Keil, R.; Kuslys, M. Geochemical changes along a river-
474
groundwater infiltration flow path—Glattfelden, Switzerland. Geochim. Cosmochim. Acta 1988, 52,
475
2693–2706. 21 ACS Paragon Plus Environment
Environmental Science & Technology
Page 22 of 30
476 477
(32) Postma, D.; Larsen, F.; Nguyen, T. T.; Pham, T. K. T.; Jakobsen, R.; Pham, Q. N.; Tran, V. L.;
478
Pham, H. V.; Murray, A. S. Groundwater arsenic concentration in Vietnam controlled by sediment
479
age. Nature Geoscience 2012, 5, 656-661.
480 481
(33) Postma, D; Pham, T. K. T.; Sø, H. U.; Hoang, V. H.; Vi, M. L.; Nguyen, T. T.; Larsen, F.;
482
Pham, H. V.; Jakobsen, R. A model for the evolution in water chemistry of an arsenic contaminated
483
aquifer over the last 6000 years, Red River floodplain, Vietnam. Geochim. Cosmochim. Acta 2016,
484
195, 277–292.
485 486
(34) Parkhurst, D. L.; Appelo, C. A. J. Description of input and examples for PHREEQC version
487
3—A computer program for speciation, batch-reaction, one-dimensional transport, and inverse
488
geochemical calculations: U.S. Geological Survey Techniques and Methods, 2013. book 6, chap.
489
A43, 497 p., available only at http://pubs.usgs.gov/tm/06/a43.
490 491
(35) Stute, M.; Zheng, Y.; Schlosser, P.; Horneman, A.; Dhar, R. K.; Datta, S.; Hoque, M. A.;
492
Seddique, A. A.; Shamsudduha, M.; Ahmed, K. M.; van Geen, A. Hydrological control of As
493
concentrations in Bangladesh Groundwater. Water Res. Res. 2007, 43, W09417,
494
doi:10.1029/2005WR004499.
495 496
(36) Nguyen, T. H. M.; Postma, D.; Pham, T. K. T.; Jessen, S.; Pham, H. V.; Larsen, F.; Adsorption
497
and desorption of arsenic to aquifer sediment on the Red River floodplain at Nam Du, Vietnam.
498
Geochim. Cosmochim. Acta 2014,142, 587-600.
499
22 ACS Paragon Plus Environment
Page 23 of 30
Environmental Science & Technology
500
(37) Jessen, S.; Postma, D.; Larsen, F.; Pham, Q. N.; Le, Q. H.; Pham, T. K. T.; Tran V. L.; Pham,
501
H. V., Jakobsen, R. Surface complexation modeling of groundwater arsenic mobility: Results of a
502
forced gradient experiment in a Red River flood plain aquifer, Vietnam. Geochim. Cosmochim.
503
Acta 2012, 98, 186–201.
504 505
(38) Radloff, K. A.; Zheng, Y.; Stute, M.; Weinman, B.; Bostick, B.; Mihajlov, I.; Bounds, M.;
506
Rahman, M. M.; Huq, M. R.; Ahmed, K. M.; Schlosser, P.; van Geen, A. Reversible adsorption and
507
flushing of arsenic in a shallow, Holocene aquifer of Bangladesh. Appl. Geochem. 2015, in press.
508 509
(39) Larsen, F.; Pham, Q. N.; Dang, D. N.; Postma, D.; Jessen, S.; Pham, H. V.; Nguyen, B. T.;
510
Duc, T. H.; Luu, T. T.; Nguyen, H.; Chambon, J.; Nguyen, V. H.; Hoang, D. H.; Nguyen, T. M. H.;
511
Mai, T. D.; Refsgaard, J. C. Controlling geological and hydrogeological processes in an arsenic
512
contaminated aquifer on the Red River flood plain, Vietnam. Appl. Geochem. 2008, 23, 3099–3115.
513 514
515 516
517
518 519
23 ACS Paragon Plus Environment
Environmental Science & Technology
520
Page 24 of 30
FIGURE CAPTIONS
521
522
Figure 1. Hydrogeological setting at the Nam Du field site. The upper layers of clay and fine-
523
medium grained sand belong to the Holocene, the coarse sand-gravel to the Pleistocene and the silt-
524
claystone to the Neogene period. Water levels in the different boreholes are indicated by ▼. The
525
low measured electrical conductivity (EC) values (200-270 µS/cm) indicate the presence of river
526
water in the whole system. Arsenic is low to absent in the river water and upper 17 m of the
527
groundwater but is mobilized further down.
528
529
Figure 2. Red River water infiltration into the aquifer: Tritium-Helium groundwater ages, plotted in
530
duplicate at all depths, show the presence of young water down to 40 m depth. The electrical
531
conductivity (EC) in the groundwater is close to the value of 205 µS/cm in the river water. The
532
partial pressure of CO2 is 10-3.3 in river water and 10-2.3 to 10-2.8 in the groundwater.
533
534
Figure 3. The effect of redox processes on the groundwater chemistry in the Holocene aquifer at
535
Nam Du. Symbols designate field data and lines reactive transport model predictions with
536
adsorption of As(III) disabled. In the model the upper cell consumes 0.27 mM DOC which reduces
537
O2 and NO3 in the infiltrating water. Further down the oxidation of organic carbon and reduction of
538
As-containing Fe-oxides are controlled by rate equations in the model of Postma et al. (2016).
539
Model lines are given for 70 years but are near identical for other times < 100 yr.
540
24 ACS Paragon Plus Environment
Page 25 of 30
Environmental Science & Technology
541
Figure 4. Leaching and retardation of As(III) from the Holocene aquifer. The symbols are field
542
data. Lines are leaching profiles calculated using a K’d of 12.6 L/kg and a flow velocity of 19 m/yr
543
derived from tritium-helium dating. The calculations were carried out for an initial groundwater
544
concentration of 1 µM As(III).
25 ACS Paragon Plus Environment
Environmental Science & Technology
Figure 1. ACS Paragon Plus Environment
Page 26 of 30
Page 27 of 30
Environmental Science & Technology
EC µS/cm 0
200
400
0
10
20
30
40
Figure 2.
ACS Paragon Plus Environment HeT_GEUSVietnam2015_26C (3)_dip_nam du.xls
Nam_Du_1_0510_Total_dp_newCH4_mod_v12.xls
Environmental Science & Technology
O2 mmol/L
DOC mmol/L
Depth, m
0
0.5
1
0
0
10
0
0.1
0.2
NO3 µmol/L 0.3
0
20
Page 28 of 30
Fe(II) µmol/L 40
0
50
As(III) µmol/L 100
0
0
0
0
10
10
10
10
20
20
20
20
20
30
30
30
30
30
40
40
40
40
40
1
2
Figure 3.
ACS Paragon Plus Environment Nam_Du_1_0510_Total_dp_newCH4_mod_v14.xls
Page 29 of 30
Environmental Science & Technology
As(III) µmol/L 0
1
2
0
60 yr
Depth, m
10
70 yr 80 yr
20
30
40
Figure 4.
ACS Paragon Plus Environment Nam_Du_1_0510_Total_dp_newCH4_mod_v13.xls
Environmental Science & Technology
47x26mm (300 x 300 DPI)
ACS Paragon Plus Environment
Page 30 of 30