Subscriber access provided by READING UNIV
Article
Aggregation Behavior of Dissolved Black Carbon: Implications for Vertical Mass Flux and Fractionation in Aquatic Systems Fanchao Xu, Chenhui Wei, Qingqing Zeng, Xuening Li, Pedro J. J. Alvarez, Qilin Li, Xiaolei Qu, and Dongqiang Zhu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b04232 • Publication Date (Web): 13 Nov 2017 Downloaded from http://pubs.acs.org on November 15, 2017
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 27
Environmental Science & Technology
1
Aggregation Behavior of Dissolved Black Carbon:
2
Implications for Vertical Mass Flux and Fractionation in
3
Aquatic Systems
4
Fanchao Xu†, Chenhui Wei†, Qingqing Zeng†, Xuening Li¶, Pedro J.J. Alvarez§, Qilin Li§, Xiaolei Qu†*,
5
Dongqiang Zhu‡
6
†
State Key Laboratory of Pollution Control and Resource Reuse, School of the Environment, Nanjing
7
University, Jiangsu 210023, China ¶
8
State Key Laboratory of Petroleum Pollution Control, CNPC Research Institute of Safety and
9 10
Environmental Technology, Beijing 102206, China §
Department of Civil and Environmental Engineering, Rice University, Houston TX 77005, United
11 12 13
States ‡
School of Urban and Environmental Sciences, Peking University, Beijing 100871, China
* Corresponding author Xiaolei Qu, phone: +86-025-8968-0256; email:
[email protected] 14 15
ACS Paragon Plus Environment
1
Environmental Science & Technology
16
Page 2 of 27
Abstract
17
The fluvial export of dissolved black carbon (DBC) is a major land-ocean flux in the global black
18
carbon cycle, affecting the size of refractory carbon pool in the oceans. The aggregation behavior of
19
DBC is a significant determinant of its transport and vertical mass flux. In this study, the aggregation
20
kinetics and interaction energy of DBC leached from biochar were investigated. DBC was mainly
21
stabilized by hydration force and underwent structural compacting in divalent cation solutions. Na+ and
22
Mg2+ had limited impact on the colloidal stability of DBC due to the strong hydration of these cations.
23
Ca2+ and Ba2+ readily destabilized DBC by forming inner-sphere complexes, reducing its hydrophilicity.
24
Consistently, charge reversal of DBC was observed with high concentrations of Ca2+ and Ba2+.
25
Simulated sunlight exposure led to photo-oxidation of DBC, increasing its colloidal stability. DBC
26
behaved non-conservatively in laboratory mixing experiments using estuary water samples due to
27
aggregation/sedimentation; while model aquatic humic acid behaved conservatively. Our results infer
28
that there is a vertical mass flux of DBC and possible fractionation from the dissolved organic matter
29
pool in the fluvial and estuarine systems, which have been overlooked in efforts to determine global
30
carbon budgets and associated climate change implications.
31
ACS Paragon Plus Environment
2
Page 3 of 27
Environmental Science & Technology
32 33
Introduction
34
Black carbon is the carbonaceous residue formed during incomplete combustion of biomass or fossil
35
fuel.1 It can release the water soluble fraction, usually referred to as dissolved black carbon (DBC),
36
when in contact with water during infiltration or surface runoff of storm water, which comprises around
37
10% of the dissolved organic carbon (DOC) in the rivers and 2% in the oceans.2-6 Global fluvial
38
discharge of DBC is estimated to be about 26.5 million tons per year, being the major flux of black
39
carbon to the oceans.5 Owing to the recalcitrant nature of DBC, this flux is an important contributor to
40
the oceanic pool of refractory carbon, which plays a significant role in oceanic carbon cycling,
41
atmospheric CO2 levels, and climate change.2,
42
domains and oxygen-containing functional groups,3, 4 making it a good sorbent for organic compounds
43
and metals. Thus, it can serve as a carrier for priority pollutants,9 affecting their fate and transport.
5-8
Furthermore, DBC contains abundant aromatic
44
The rate and extent of aggregation of DBC and subsequent formation of fast-settling particles in
45
fluvial systems can significantly affect its flux to the oceans, influencing the global budget of refractory
46
carbon, its impact on climate change, and the environmental fate of priority contaminants. However,
47
little is known about the aggregation behavior and vertical mass flux of DBC in aquatic systems. Two
48
previous studies investigated the deposition behavior of DBC in the saturated granular media,10, 11 and
49
described the deposition behavior adequately using the extended Derjaguin-Landau-Verwey-Overbeek
50
(XDLVO) theory to consider electrostatic and Lewis acid-base interactions between DBC and quartz
51
sand.10 The presence of humic acid was reported to enhance the transport of DBC in sand columns,
52
while iron oxyhydroxide coatings on sand grains decreased DBC transport.11 Nevertheless, the
53
aggregation behavior of DBC has not been addressed in the literature. Furthermore, DBC was reported
54
to undergo structural changes during solar irradiation,3, 12-14 and the impact of solar exposure on DBC
55
aggregation behavior is still unknown.
56
The removal of DOM from aquatic systems is determined by its degradation or aggregation into large
57
particulates.15 Several studies have reported that the riverine DOM pool is transported conservatively ACS Paragon Plus Environment
3
Environmental Science & Technology
Page 4 of 27
58
(i.e., without addition or loss) to the oceans.16-18 On the other hand, several studies pointed out that
59
DOM behaves non-conservatively in several estuaries due to multiple biogeochemical processes
60
including aggregation/sedimentation, degradation, autochthonous production, and solute-particle
61
interactions.19-21 The vertical mass flux of DBC was largely unaccounted for in the current estimation of
62
its land-ocean flux.2, 5, 8, 22 The non-conservative behavior of DBC due to aggregation/sedimentation
63
may potentially result in the fractionation of DBC from the DOM pool and its accumulation in terrestrial
64
sediments as important intermediate pools, which has not been addressed yet.
65
In the present study, DBC leached from biomass-derived black carbon was used as a proxy for DBC
66
in natural aquatic systems.3, 12 It is a continuum of macromolecules with high aromaticity, which is
67
different from black carbon nanoparticles generated by ball-milling in a previous study.4, 12, 23 The initial
68
aggregation kinetics and electrokinetic properties of DBC were characterized in the presence of
69
naturally abundant mono- and di-valent cations, including Na+, Mg2+, Ca2+, and Ba2+. The interaction
70
energy of London-van der Waals force, electrostatic force, and Lewis-acid-base interactions was
71
calculated based on the XDLVO model. The impact of solar exposure on its structural properties and
72
aggregation behavior was also examined. The behavior of DBC during mixing of natural water samples
73
collected from two end-members in the Yangtze River Estuary was compared with a theoretical
74
conservative mixing scenario to discern the significance of non-conservative processes (i.e.,
75
aggregation/sedimentation in this case). Our main objectives were to 1) determine the effect of ionic
76
strength, ionic composition, and sunlight exposure on the aggregation behavior of DBC; 2) discern the
77
aggregation mechanisms; and 3) examine possible vertical mass flux and fractionation of DBC along the
78
salinity gradient due to aggregation/sedimentation processes, which are important knowledge gaps to
79
inform global black carbon budgets and associated climate change implications.
80 81
Materials and Methods
82
Materials. CaCl2 (> 99%) was obtained from Alfa Aesar, UK. MgCl2 (> 99.5%), NaCl (> 99.5%),
83
and BaCl2 (> 99%) were purchased from Sigma-Aldrich, USA. Deuteroxide (D2O, 99.8 atom % D) was ACS Paragon Plus Environment
4
Page 5 of 27
Environmental Science & Technology
84
purchased from Tokyo Chemical Industry, Japan. Sodium-3-trimethyl silyl propionate (TMSP-2,2,3,3-
85
D4) was purchased from Cambridge Isotope Laboratories, USA. Suwannee River humic acid (SRHA)
86
was provided by the International Humic Substances Society (IHSS, St. Paul, USA). Deionized water
87
(18.2 MΩ·cm resistivity at 25 °C) used in the experiments was produced by an ELGA Labwater system
88
(PURELAB Ultra, ELGA LabWater Global Operations, UK).
89
Preparation of DBC. DBC was released from bulk black carbon by gentle sonication in water as
90
described in a previous study.3 Briefly, bamboo shavings collected from Lishui, Zhejiang Province,
91
China was pulverized into powder (high-speed pulverizer FW 100, Tianjin Taisite Instrument, China)
92
and pyrolyzed in a muffle furnace under oxygen-limited condition. The pyrolysis temperature was set to
93
increase from 20 to 400 °C in 2 h and kept at 400 °C for 3 h. The resulting black carbon was pulverized
94
and passed through a 100-mesh sieve. The sieved black carbon powder (30 g) was mixed with 500 mL
95
deionized water in a 1000-mL glass beaker; sonicated in a sonication bath (KH-800TDB, Kunshan
96
Hechuang Ultrasonic Instrument, China) at 50 W for 30 min; and filtered through a 0.45-µm membrane
97
(Pall, USA). The residue retained on the membrane was collected and subjected to another round of
98
sonication and filtration. The sonication-filtration cycle was repeated three times. The filtrate from each
99
cycle was collected and freeze-dried to yield DBC powder, which was stored in a desiccator at room
100
temperature until use. The DBC stock solution was made by dissolving the DBC powder in deionized
101
water. The pH of the stock solution was adjusted to 6.8 using HCl.
102
Irradiation experiments. DBC solution was irradiated with simulated sunlight to examine the impact
103
of sunlight exposure on the colloidal stability of DBC. A xenon lamp (CEL-HXF300, AULTT, China)
104
was used to simulate the solar irradiation. The lamp spectrum was similar to natural sunlight as
105
measured by a spectrometer (USB2000+, Ocean Optics, FL, USA) (Figure S1). The DBC solution (20
106
mL, 100 mg/L) was stirred at 200 rpm in a 50-mL glass vial which was immersed in a water-circulating
107
bath at 20 ± 0.1 °C. The DBC solution was irradiated by the xenon lamp from the top at a distance of 0.2
108
m for 24 and 48 h respectively. Water loss due to evaporation during the experiment was compensated
109
using deionized water. ACS Paragon Plus Environment
5
Environmental Science & Technology
Page 6 of 27
110
Aggregation kinetics of DBC. Measurement of Particle Aggregation Rate. The aggregation kinetics
111
of the original and irradiated DBC were studied in various solution chemistry by time-resolved dynamic
112
light scattering (DLS) measurements using a ZEN 3500 Zetasizer Nano ZS (Malvern, Worcestershire,
113
UK) equipped with a 532-nm laser. The refractive index of DBC was set to be 2.20.24 The DBC stock
114
solution was sonicated in a sonication bath at 50 W for 5 min before the aggregation experiments. In the
115
measurements of aggregation kinetics, 0.2 mL of 100 mg/L DBC solution and a pre-determined amount
116
of electrolyte stock solution and deionized water were introduced into a disposable polystyrene cuvette
117
(Sarstedt, Germany) to yield a sample with a total volume of 1 mL. The sample was briefly mixed and
118
immediately analyzed at 25 °C. The time-resolved DLS measurements monitor the average
119
hydrodynamic diameter of DBC with the autocorrelation function being accumulated over a period of 7
120
s. If the particle size increased or decreased more than 50% in 7 s, the measurement was removed from
121
the dataset.
122
Determination of Aggregation Kinetics. The initial aggregation rate of DBC was determined by
123
applying linear least square regression to the dh(t) data from the initial value dh(0) to 1.5 dh(0). The
124
attachment efficiency α is a key matric to quantify the tendency of particle aggregation in a given water
125
chemistry. It was determined by normalizing the initial aggregation rate in the solution of interest by
126
that in the diffusion-limited regime.25, 26 Stability curves of DBC were drawn by plotting the α value as
127
a function of the electrolyte concentration. The initial aggregation rate in the diffusion-limited regime
128
was the average of the initial aggregation rates at the plateau of the stability curve. The critical
129
coagulation concentration (CCC) was approximated by the electrolyte concentration where the
130
diffusion-limited regime was achieved.
131 132
Characterization of DBC. Transmission electron microscope (TEM) images of DBC were taken using a JEM-200CX, JEOL, Japan.
ACS Paragon Plus Environment
6
Page 7 of 27
Environmental Science & Technology
133
Electrokinetic properties of DBC. The electrophoretic mobility (EPM) and ζ-potential of DBC in
134
different solution chemistry were measured by phase analysis light scattering (PALS) using the ZEN
135
3500 Zetasizer Nano ZS (Malvern, Worcestershire, UK). Each sample was measured five times at 25 °C
136
in a folded capillary cell (Malvern, Worcestershire, UK).
137
Contact angle measurements. The DBC solution with different water chemistry was filtered through
138
an ultrafiltration membrane (3 kDa, Millipore Corporation, USA) using an ultrafiltration cell (Model
139
8010, Millipore Corporation, USA) to form a DBC layer on the membrane. The DBC layer on top of the
140
ultrafiltration membrane was dried by nitrogen gas under room temperature and then subjected to static
141
contact angle measurements. Static contact angles of probe liquids including water, glycerol, and n-
142
decane on the surface of DBC layer were measured at least three times at different locations by an
143
optical contact angle measuring device (OCA30, Dataphysics Instruments Gmb, Germany).
144
Proton nuclear magnetic resonance (1H-NMR) analysis. Ten milliliters of 1000 mg/L DBC in D2O
145
was irradiated using the same setting as irradiation experiments. Samples were withdrawn from the
146
reactor at pre-determined irradiation time and analyzed by NMR spectroscopy. The 1H-NMR spectra
147
were acquired on a Bruker AVANCE DRX-600 NMR spectrometer (Bruker, Massachusetts, USA). The
148
chemical shifts (δ, ppm) were internally referenced to a standard, TMSP-2,2,3,3-D4.
149
Water sample collection and mixing experiments. The freshwater sample was collected from
150
Yangtze River near Shanghai (N 31°33.143', E 121°44.645'). The seawater sample was collected from
151
the beach at Gouqi Island in the East China Sea (N 30°43.176', E 122°47.573'). Both water samples
152
were filtered through 0.45-µm membranes (Pall, USA) after collected. To investigate the estuarine
153
mixing behavior of DBC, the laboratory mixing experiments were carried out using these two water
154
samples as the end-members of Yangtze River Estuary. Two milligrams of DBC or SRHA was
155
dissolved in 200 mL freshwater sample to simulate a high-DBC/SRHA-containing river water. It was
156
mixed with the seawater sample in different proportions to form a salinity gradient. The resulting
157
mixtures were kept quiescent for 24 h in the dark. The absorption spectra of the mixtures after 24 h ACS Paragon Plus Environment
7
Environmental Science & Technology
Page 8 of 27
158
settling were measured by UV−vis spectroscopy in a 1-cm quartz cuvette. The absorption coefficient at
159
254 nm, a254, was used to quantify the abundance of DBC/SRHA in water samples.22
160
Result and Discussion
161
Characterization of DBC.
162
DBC forms a stable solution of brown color in deionized water. TEM images suggest that DBC
163
assembles into irregular-shaped particles with diameters of hundreds of nanometers (Figure S2). It was
164
negatively charged with a ζ-potential of -25.8 ± 1.9 mV (i.e., electrophoretic mobility of -1.8 ± 0.1
165
×10-8 m2 V-1 s-1) in 1 mM NaCl solution at pH 6.85. This is lower than the ζ-potential of DBC released
166
from black carbon made from the slow pyrolysis of pine needle and wheat straw, -28.9 ~ -36.6 mV.10, 11
167
The elemental and structural properties of DBC were thoroughly characterized in our previous studies.3,
168
4
169
international humic substances society, 0.62.4,
170
carbon, 32.4 % carboxyl/ester/quinone carbon, 5.5 % aromatic C-O carbon, and 18.5 % aliphatic carbon
171
as determined by quantitative solid-state
172
carbons, but less aliphatic carbon than SRHA.3, 25 The elemental and structural analyses suggest the
173
abundance of polar moieties, especially carboxyl groups, in DBC. Note that DBC in this study refers to
174
all carbon leached from biomass-derived black carbon. It contains a wider range of carbon materials
175
than that in previous spatiotemporal distribution studies which considered DBC as only highly
176
condensed aromatic structures in the DOM pool.2, 28, 29
177
Aggregation behavior of DBC in monovalent cation solutions.
The (O+N)/C ratio of DBC was 0.66, slightly higher than SRHA, a reference humic acid from the 27
DBC used in this work contains 43.7 % aromatic
13
C NMR.3 It has more carboxyl/ester/quinone and aromatic
ACS Paragon Plus Environment
8
Page 9 of 27
Environmental Science & Technology
1500 1200
(b)
20
10 mM NaCl 800 mM NaCl
10 ζ− potential (mV)
Particle Diameter (d.nm)
(a)
900 600 300
0 -10 -20 -30
0
-40 1000
2000
3000
1
10
Time (s)
178
100
1000
NaCl Concentration (mM)
179
3000
3000
2000
(d)
AB Interaction EDL Interaction LW Interaction TOT Interaction
Interaction Energy (kB T)
Interaction Energy (kB T)
(c)
1000
0
-1000
2000
AB Interaction EDL Interaction LW Interaction TOT Interaction
1000
0
-1000
0
10
20
Separation Distance h (nm)
0
10
20
Separation Distance h (nm)
180 181
Figure 1. (a) Aggregation profile of 20 mg/L DBC in 10 and 800 mM NaCl solutions at pH 6.5 ± 0.4
182
measured by time-resolved dynamic light scattering; (b) ζ-potential of DBC as a function of NaCl
183
concentration. Error bar represents one standard deviation of five measurements. The interaction energy
184
of London-van der Waals force (LW), electrostatic force (EDL), Lewis-acid-base interactions (AB), and
185
the total interaction energy (TOT) of DBC in (c) 10 mM and (d) 800 mM NaCl solutions.
186
The aggregation kinetics of DBC in 10 mM and 800 mM NaCl solutions are shown in Figure 1a.
187
There was no significant change of hydrodynamic size of DBC within 1 h with NaCl concentrations up
ACS Paragon Plus Environment
9
Environmental Science & Technology
Page 10 of 27
188
to 800 mM, which was higher than most ionic strength found in natural aquatic systems. Classical
189
DLVO theory points out that the thermodynamics of particle attachment is controlled by the interplay of
190
London-van der Waals force and electrostatic force;30, 31 and the aggregation behavior is mostly decided
191
by the surface charge as reflected by its electrokinetic properties.32-34 The ζ-potential of DBC was
192
examined over a wide range of NaCl concentrations as shown in Figure 1b. DBC became less
193
negatively charged from -25.8 ± 1.9 mV to -4.5 ± 2.6 mV as the NaCl concentration increased from 1
194
mM to 800 mM owing to charge screening. In the 800 mM NaCl solution, the ζ-potential of DBC was -
195
4.5 ± 2.6 mV, close to neutral. Thus, electrostatic repulsion alone cannot account for the strong colloidal
196
stability of DBC (i.e., its ability to resist aggregation). We further examined the aggregation behavior of
197
DBC in the framework of XDLVO theory. The interaction energy of London-van der Waals force,
198
electrostatic force, Lewis-acid-base interactions, and the total interaction energy of DBC particles in 10
199
mM and 800 mM NaCl solutions were summarized in Figure 1c and 1d (see detailed calculation35, 36 and
200
dataset in Supporting Information). The free energy of Lewis-acid-base interactions between DBC
201
particles (GAB) was determined using contact angle measurements and Young-Dupré equation.10, 35 The
202
calculated ΔGAB were 37.3 mJ m−2 and 35.7 mJ m−2 in 10 mM and 800 mM NaCl solutions respectively,
203
indicating the strong hydrophilic nature of DBC particle surfaces.35 This is consistent with the high
204
polarity index (i.e., (O+N)/C ratio) and the abundant oxygen-containing functional groups of DBC. The
205
positive nature of ΔGAB also suggests repulsive forces between DBC particles due to surface cohesion
206
water, i.e., hydration force.35,
207
stabilization of DBC at both low and high NaCl concentrations (Figure 1c and 1d). The aggregation
208
behavior of DBC is different from black carbon nanoparticles which can be destabilized by Na+ with a
209
CCC of 250 mM,23 consistent with its macromolecular nature.
210
Aggregation behavior of DBC in divalent cation solutions and charge reversal
37
Clearly, the hydration force played a predominant role in the
ACS Paragon Plus Environment
10
Page 11 of 27
Environmental Science & Technology
1500
20
(b)
10 mM MgCl2
1200
800 mM MgCl2
ζ− potential (mV)
Particle Diameter (d.nm)
(a)
900 600 300 0
Mg2+ Ba 2+ Ca 2+
10
0
-10
-20 1000
2000
3000
1
10
Time (s)
211
1000
Concentration (mM)
3000
3000
2000
(d)
AB Interaction EDL Interaction LW Interaction TOT Interaction
Interaction Energy (kB T)
(c) Interaction Energy (kB T)
100
1000
0
-1000
2000
AB Interaction EDL Interaction LW Interaction TOT Interaction
1000
0
-1000 0
10
20
Separation Distance h (nm)
0
10
20
Separation Distance h (nm)
212 213
Figure 2. (a) Aggregation profile of 20 mg/L DBC in 10 and 800 mM MgCl2 solutions at pH 6.2 ± 0.5
214
measured by time-resolved dynamic light scattering; (b) ζ-potential of DBC as a function of MgCl2,
215
CaCl2, and BaCl2 concentration. Error bar represents one standard deviation of five measurements. The
216
interaction energy of London-van der Waals force (LW), electrostatic force (EDL), Lewis-acid-base
217
interactions (AB), and the total interaction energy (TOT) of DBC in (c) 10 mM and (d) 800 mM MgCl2
218
solutions.
ACS Paragon Plus Environment
11
Environmental Science & Technology
1 mM BaCl2 10 mM BaCl 2 1000
200 mM CaCl2
(a)
Particle Diameter (d.nm)
Particle Diameter (d.nm)
1000
1 mM CaCl 2 10 mM CaCl2 100 mM CaCl2
800 600 400 200
100 mM BaCl2
(b)
200 mM BaCl2
800 600 400 200 0
0 200
400
600
800
200
1000 1200
400
600
800
1000 1200
Time (s)
Time (s)
(d) Attachment Efficiency, α
(c) Attachment Efficiency, α
Page 12 of 27
1
.1
1
.1 10
100
CaCl2 Concentration (mM)
10
100
BaCl2 Concentration (mM)
ACS Paragon Plus Environment
12
Page 13 of 27
Environmental Science & Technology
50
(e)
(f) Interaction Energy (kB T)
Interaction Energy (kB T)
50 0 -50 -100
AB Interaction EDL Interaction LW Interaction TOT Interaction
-150 -200
0 -50 -100
AB Interaction EDL Interaction LW Interaction TOT Interaction
-150 -200
0
10 Separation Distance h (nm)
20
0
10
20
Separation Distance h (nm)
219 220
Figure 3. Aggregation profile of 20 mg/L DBC in (a) CaCl2 and (b) BaCl2 solutions at pH 6.0 ± 0.5
221
measured by time-resolved dynamic light scattering; and attachment efficiencies of DBC as a function
222
of (c) CaCl2 concentration and (d) BaCl2 concentration. The interaction energy of London-van der
223
Waals force (LW), electrostatic force (EDL), Lewis-acid-base interactions (AB), and the total
224
interaction energy (TOT) of DBC in 100 mM (e) CaCl2 and (f) BaCl2 solutions.
225
Colloidal stability depends strongly on the valence of counter ions, as described by the Schulze-Hardy
226
Rule.38 Mg2+, Ca2+, and Ba2+ are divalent cations found in natural waters, which are expected to impact
227
the colloidal stability of DBC more significantly than monovalent cations. The aggregation kinetics of
228
DBC in 10 mM and 800 mM MgCl2 solutions are shown in Figure 2a. The particle size of DBC
229
decreased in the initial stage after mixing with MgCl2 solutions and then stabilized at 422.7 ± 77.3 nm
230
and 229.2 ± 22.4 nm in 10 mM and 800 mM MgCl2 solutions, respectively. DBC contains abundant
231
carboxyl and phenolic groups, providing electrostatic repulsion between functional moieties and
232
consequently leads to a stretched configuration. The presence of Mg2+ reduced the intramolecular
233
electrostatic repulsion, resulting in a more compact configuration of DBC. Thus, the initial decrease in
234
particle size can be attributed to structural compacting, which was previously observed for DOM and
235
polyelectrolytes.39-41 No significant aggregation was observed in 800 mM MgCl2 solution in 1 h. The ACS Paragon Plus Environment
13
Environmental Science & Technology
Page 14 of 27
236
XDLVO model calculation of the interaction energy of DBC in 10 mM and 800 mM MgCl2 solutions is
237
present in Figure 2c and 2d. The high colloidal stability of DBC in MgCl2 solutions is consistent with
238
the high energy barrier provided by hydration force.
239
In contrast to that observed in NaCl and MgCl2 solutions, Ca2+ and Ba2+ can readily induce the
240
aggregation of DBC. Representative aggregation profiles of DBC in CaCl2 and BaCl2 solutions are
241
present in Figure 3a and 3b. The attachment efficiency α, ranging from 0 to 1, represents the probability
242
of forming an aggregate during a collision between two DBC particles in the given water chemistry. It
243
was determined by normalizing the initial aggregation rate in the solution of interest by that in the
244
diffusion-limited regime.25, 26 Higher attachment efficiency suggests a less stable colloidal system. The
245
attachment efficiencies of DBC in CaCl2 and BaCl2 solutions are summarized in Figure 3c and 3d. The
246
stability plot has two distinct regimes. At low salt concentrations, the attachment efficiency increased
247
with salt concentration. Then the attachment efficiency reached a plateau after the salt concentration
248
exceeded the CCC. The CCC values for Ca2+ and Ba2+ are 75 mM and 60 mM, respectively. Overall, the
249
destabilization capacity of the cations tested follows the order of Ba2+ > Ca2+ >> Mg2+, Na+.
250
The electrokinetic properties of DBC in MgCl2, CaCl2, and BaCl2 solutions are shown in Figure 2b.
251
The ζ-potential of DBC became less negatively charged as salt concentration increased at the initial
252
stage. Among the three divalent cations tested, the efficiency in reducing the surface potential of DBC
253
follows Ca2+ > Ba2+ > Mg2+ (Figure 2b), and they were all much more efficient in reducing the surface
254
charge of DBC than monovalent ions. As the salt concentration further increased, charge reversal was
255
observed in CaCl2 and BaCl2 solutions, which was in contrast to that in NaCl and MgCl2 solutions and
256
cannot be predicted by electric double layer models. The charge reversal concentration, cR, (i.e., the
257
concentration of electrolyte at the isoelectric point) was 150 mM for Ca2+ and 190 mM for Ba2+. Charge
258
reversal of particles and macromolecules was previously attributed to the formation of adsorbed layer of
259
cationic polyelectrolytes, charge correlation effect, or specific sorption.42-47 We hypothesize that Ca2+
260
and Ba2+ can specifically bond to DBC, which leads to the overcharging of the surface.
ACS Paragon Plus Environment
14
Page 15 of 27
Environmental Science & Technology
261
The drastically different aggregation and electrokinetic properties of DBC in Na+, Mg2+, Ca2+, and
262
Ba2+ solutions can be attributed to the different interaction mechanisms between cations and DBC. Na+
263
is surrounded by a compact water shell in the aqueous phase and behaves as indifferent counter ions that
264
only affect the electrical double layer thickness. Thus, Na+ has little impact on the hydrophilicity of
265
DBC particles and consequently its aggregation behavior. Consistently, the Δ GAB only slightly
266
decreased from 37.3 mJ m−2 to 35.7 mJ m−2 as NaCl concentration increased from 10 mM to 800 mM.
267
Mg2+ was reported to either form no complex with carboxyls in natural organic matter48 or weak
268
complexes.49, 50 It is also reported that the stability constants for complexes between Mg2+ and fulvic
269
acids were much lower than that for Ca2+.51 The ΔGAB decreased from 30.1 mJ m−2 to 9.8 mJ m−2 as
270
MgCl2 concentration increased from 10 mM to 800 mM. This decrease is more significant than that in
271
NaCl solutions. Nevertheless, the energy barrier for aggregation was still high at 800 mM MgCl2,
272
1.49×10-18 J (i.e., 368 kBT, Figure 2d). Thus, Mg2+ also had limited impact on the aggregation behavior
273
of DBC.
274
On the other hand, Ca2+ and Ba2+ were previously reported to form inner-sphere complexes with the
275
carboxylic groups of humic substances.48-50, 52, 53 The strong complexation between Ca2+/Ba2+ and the
276
oxygen-containing functional groups in DBC was consistent with the charge reversal observed (Figure
277
2b) and the XDLVO model calculation (Figure 3e and 3f). The ΔGAB decreased to -4.6 mJ m−2 and -5.3
278
mJ m−2 in 100 mM CaCl2 and BaCl2 solutions receptively, indicating that the Lewis-acid-base
279
interactions shifted from repulsive (i.e., hydration force) to attractive (i.e., hydrophobic effect) at
280
elevated Ca2+/Ba2+ concentrations. Previous studies reported that highly charged biosurfaces changed
281
from hydrophilic to hydrophobic as their charge was partially or completely neutralized by Ca2+.35, 54, 55
282
The complexation with Ca2+/Ba2+ brings polar groups together, leading to a compact configuration of
283
DBC. In this configuration, less polar groups but more hydrocarbon structures of DBC were exposed to
284
the water phase, resulting in increased hydrophobicity.54 This is consistent with the significantly
285
reduced electron-donor surface tension component (γ-) and the increased electron-acceptor surface ACS Paragon Plus Environment
15
Environmental Science & Technology
Page 16 of 27
286
tension component (γ+) of DBC after the addition of Ca2+/Ba2+ (Table S1). In addition, Ca2+ and Ba2+
287
can form complexes with carboxyl groups in different DBC particles, potentially leading to enhanced
288
aggregation through bridging effect.56
289
Sunlight exposure increases the colloidal stability of DBC
0
0h 24h 48h .1 10
290
(b)
-10
1
ζ− potential (mV)
Attachment Efficiency, α
(a)
-20 -30 -40 -50
100
1000
0
Concentration (mM)
10
20
30
40
50
60
Irradiated time (h)
291 292
Figure 4. (a) Attachment efficiencies of original and irradiated DBC as a function of CaCl2
293
concentration at pH 6.0 ± 0.4; (b) ζ-potential of original and irradiated DBC in deionized water. Error
294
bar represents one standard deviation of five measurements.
295 296
Table 1. Composition of H atoms in different moieties in DBC obtained by 1H-NMR.
Moieties Location (ppm) )
Original DBC 24h-irradiated DBC 48h-irradiated DBC
0.0-1.2
1.2-1.9
1.9-3.1
3.1-4.7
5.9-8.6
Methyl
Aliphatic
Functionalized aliphatic
Heteroatom substituted
Aromatic
8.7 7.7 9.0
7.4 8.7 8.5
29.9 30.0 29.3
11.2 12.0 13.2
42.8 41.6 39.9
297 298
The colloidal stability of DBC is closely related to its structural properties, which would determine its
299
hydrophilicity and surface charge. Interactions with natural elements such as sunlight are expected to ACS Paragon Plus Environment
16
Page 17 of 27
Environmental Science & Technology
300
change the structure of DBC and consequently influence its colloidal stability. The aggregation profiles
301
and electrokinetic properties of DBC before and after simulated sunlight exposure are compared in
302
Figure 4. The colloidal stability of DBC in CaCl2 solutions increased with irradiation time (Figure 4a).
303
The CCC of the 48-h irradiated DBC was 400 mM, more than 5 times higher than that of the original
304
DBC, 75 mM.
305
We investigated the structural evolution and hydrophilicity of DBC during irradiation using 1H NMR.
306
The 1H NMR spectra were integrated over five regions: 0-1.2 ppm, methyl groups and unfunctionalized
307
aliphatics; 1.2-1.9 ppm, aliphatic protons β and γ connected to keto or aromatic groups; 1.9-3.1 ppm,
308
functionalized aliphatics, aliphatic protons α connected to keto or aromatic groups; 3.1-4.7 ppm,
309
aliphatic protons connected to oxygen and nitrogen atoms; and 5.9-8.6 ppm, aromatic protons.57, 58 The
310
percentages of protons in different moieties in DBC before and after irradiation from NMR analysis
311
were summarized in Table 1.
312
The aromatic protons decreased with increasing irradiation time from 42.8 % to 39.9 % after 48 h
313
irradiation (Table 1), which was attributed to the substitution on aromatic rings and ring cleavage. On
314
the other hand, the percentage of protons connected to oxygen and nitrogen atoms (i.e., heteroatom
315
substituted H) increased from 11.2 % to 13.2 % after 48 h irradiation. As the formation of nitrogen-
316
containing functional groups was expected to be minimal, the increase in the presence of these protons
317
was attributed mostly to the oxidation processes that led to the formation of oxygen-containing
318
functional groups. The percentage of aliphatic protons also increased, owing to the selective
319
preservation and the generation of aliphatics during aromatic ring cleavage. Thus, DBC lost aromatic
320
moieties but gained oxygen-containing and aliphatic moieties during sunlight exposure. This argument
321
is in general consistent with a previous report using solid-state
322
(including methyl, aliphatics, and functionalized aliphatics), which favor the hydrophobicity of DBC,
323
grew 0.8 % after 48 h. However, the loss of hydrophobic aromatic structures was more significant at
324
2.9 %. DBC became more hydrophilic due to the loss of aromatic units and the addition of oxygen-
325
containing moieties during irradiation, contributing to stronger hydration force. Meanwhile, the increase
13
C NMR technique.3 Total aliphatics
ACS Paragon Plus Environment
17
Environmental Science & Technology
Page 18 of 27
326
in oxygen-containing functional groups also leads to the increase in surface charge. The ζ-potential of
327
DBC became more negative with irradiation time (Figure 4b). As a result, DBC samples with longer
328
irradiation time underwent higher electrostatic repulsion and hydration force which facilitate their
329
colloidal stability. The role of solar exposure on the aggregation behavior of DBC will be affected by
330
the turbidity of the water which influences the number of photons received by DBC.
331
DBC in the estuarine mixing
332
The aggregation/sedimentation processes of DBC during Yangtze River Estuarine mixing were
333
investigated by laboratory mixing experiments. The mixing curves for a254 of SRHA and DBC between
334
two end-members, Yangtze River freshwater and East China Sea seawater, are presented in Figure 5a
335
and 5b. The theoretical conservative mixing line represents the theoretical dilution scenario in which the
336
DBC/SRHA concentration decreases linearly from 10 mg/L to ~ 0 mg/L with the mixing of water. It is a
337
straight line connecting the two end-members. The mixing curve of SRHA coincided with the
338
theoretical conservative mixing line as shown in Figure 5a, implying its conservative behavior during
339
the mixing. The mixing curve of DBC was significantly below the theoretical dilution line, indicating
340
the removal of DBC during the mixing (Figure 5b). The non-conservative behavior of DBC was
341
observed over the entire salinity gradient including freshwater samples. The presence of 269 mg/L
342
suspended particulate matter in the freshwater sample had limited impact on the sedimentation pattern
343
of SRHA/DBC (Figure S3).
344
The non-conservative behavior of DBC in the laboratory mixing experiments can be attributed to its
345
aggregation and sedimentation in the presence of Ca2+ and Ba2+. The water chemistry of the natural
346
water samples is summarized in Table S2. The aggregation profiles of SRHA and DBC in water
347
samples collected from Yangtze River and the East China Sea, and their mixture are summarized in
348
Figure 5c and 5d. The freshwater end-member contained 1.04 mM Ca2+ and no Ba2+, and therefore was
349
only able to induce slow aggregation of DBC (Figure 5d). The seawater end-member contained 12.82
350
mM Ca2+ and 0.09 mM Ba2+, which led to fast aggregation of DBC (Figure 5d). Ca2+ is more important
ACS Paragon Plus Environment
18
Page 19 of 27
Environmental Science & Technology
351
than Ba2+ due to its natural abundance. EDTA can effective chelate with Ca2+, competing for their
352
bonding with DBC. The addition of 13.68 mM EDTA to the seawater drastically inhibited the
353
aggregation behavior of DBC (Figure 5d), further confirming the predominant role of Ca2+. The
354
aggregation kinetics increased with increasing seawater proportion which introduced more Ca2+ cations
355
into the system. The aggregation behavior of DBC was also tested by sedimentation experiments in
356
freshwater, seawater, and their mixture. The sedimentation rate increased with increasing seawater
357
proportion (Figure S4), agrees with the aggregation kinetics. On the contrary, the size of SRHA
358
remained unchanged in all the water samples within 50 min (Figure 5c). These results are consistent
359
with their behavior in laboratory mixing experiments. Theoretical Conservative Mixing SRHA 40
Theoretical Conservative Mixing DBC
(a)
(b)
30 25
30
a 254
a 254
20
20
15 10
10 5
0
360
5
10
15
20
Salinity( ‰)
25
30
0
10
20
30
Salinity( ‰)
361 362
ACS Paragon Plus Environment
19
Environmental Science & Technology
2000 Particle Diameter (d.nm)
1500
Seawater SRHA Sea : Fresh=1: 1 SRHA Freshwater SRHA
1000
500
0
Particle Diameter (d.nm)
1200
(c)
Page 20 of 27
Seawater DBC Sea : fresh=1: 1 DBC Freshwater DBC Seawater+EDTA DBC
(d)
1000 800 600 400 200 0
0
500 1000 1500 2000 2500 3000
500
1000
Time (s)
363
1500
2000
2500
Time (s)
364 365
Figure 5. The a254 as a function of the salinity for the mixture of seawater endmember and 10 mg/L (a)
366
SRHA- or (b) DBC-containing freshwater endmember with different ratios after 24h settling in dark
367
condition at pH 7.7 ± 0.2 and 7.9 ± 0.2. The dashed line represents the theoretical conservative mixing
368
line. The aggregation profile of (c) SRHA and (d) DBC in natural freshwater, seawater, their 1:1
369
mixture, and seawater with 13.68 mM EDTA at pH 7.0 ± 0.2 measured by time-resolved dynamic light
370
scattering.
371
Environmental Implications.
372
Our results suggest that the aggregation behavior of DBC is mainly controlled by the concentration of
373
cations in aquatic systems that can complex with the polar functional groups in DBC, predominantly
374
Ca2+ owing to its natural abundance. The electrokinetic properties commonly used for the prediction of
375
colloidal stability, such as ζ-potential and electrophoretic mobility, don’t apply for DBC as it is
376
stabilized mainly through hydration force. Solar irradiation enhanced the colloidal stability of DBC in
377
aquatic systems. Thus the photo-exposed DBC is expected to be more resistant to aggregation than that
378
initially released from soils. Meanwhile, photodegradation processes were suggested to be one of the
379
mechanisms for the loss of DBC in estuary waters.59 Nevertheless, these effect is influenced by the light
380
transmittance of the water column and will be hampered in high-turbidity waters. Results in this work ACS Paragon Plus Environment
20
Page 21 of 27
Environmental Science & Technology
381
imply that a portion of DBC flux will be transformed into the particulate matter by aggregation, leading
382
to vertical mass flux in the estuaries. It is worth noting that the vertical mass flux of DBC is also
383
expected in the freshwater systems based on our data. Thus a portion of DBC released from the soil
384
system, along with the contaminants associated, might accumulate in the riverine and estuarine
385
sediments and will, therefore, be lost before it enters the oceans. The source material and combustion
386
conditions will affect the properties and colloidal behavior of DBC. DBC with more hydrophilic
387
structures is expected to have higher colloidal stability due to the predominant role of Lewis-acid-base
388
interactions. The structure-activity relationships for the colloidal stability of DBC should be addressed
389
in future research. The different behavior of DBC and other DOM fractions in the fluvial and estuarine
390
systems also may lead to the fractionation of DBC from the DOM pool, affecting the DBC/DOC ratio
391
which is often assumed to be a stable characteristic for major rivers around the world and is used to
392
estimate the global DBC flux.5, 22 Our results imply that this approach might lead to an over-estimation
393
of the land-ocean flux of DBC through fluvial transport as the DBC/DOC ratio will decrease along the
394
salinity gradient due to the fractionation process. The significance of the fractionation process as well as
395
riverine/estuarine sediments as a sink of DBC is unclear so far and needs further investigation at
396
multiple scales.
397 398
ASSOCIATED CONTENT
399
Supporting Information
400
The experimental setup for irradiation experiments, the spectra of the Xenon lamp and sunlight, TEM
401
images of DBC particles, the lab mixing experiment with 269 mg/L suspended particulate matter, the
402
calculation of the Lewis-acid-base interaction energy (Δ GAB), London-van der Waals force (Δ GLW),
403
electrostatic interaction energy (ΔGEDL) of DBC in water, the contact angle and the surface tension data,
404
sedimentation profiles of DBC in water samples, and water chemistry of the natural water samples can
ACS Paragon Plus Environment
21
Environmental Science & Technology
Page 22 of 27
405
be found in the Supporting Information. This material is available free of charge via the Internet at
406
http://pubs.acs.org/.
407 408
AUTHOR INFORMATION
409
Corresponding Author
410
*Xiaolei Qu; Phone: +86-025-8968-0256; email:
[email protected] 411 412
ACKNOWLEDGMENTS
413
We thank Yuxuan Yao for her assistance in measuring DBC aggregation kinetics. This work was
414
supported by the National Key Basic Research Program of China (Grant 2014CB441103), the National
415
Natural Science Foundation of China (Grant 21407073 and 21622703), the NSF ERC on
416
Nanotechnology-Enabled Water Treatment (EEC-1449500), and the State Key Laboratory of
417
Environmental Chemistry and Ecotoxicology, Research Center for Eco-Environmental Sciences,
418
Chinese Academy of Sciences (KF2014-11).
419 420
REFERENCES
421
(1)
422 423 424
(2) Dittmar, T.; de Rezende, C. E.; Manecki, M.; Niggemann, J.; Coelho Ovalle, A. R.; Stubbins, A.; Bernardes, M. C., Continuous flux of dissolved black carbon from a vanished tropical forest biome. Nat. Geosci. 2012, 5, 618-622.
425 426 427
(3) Fu, H.; Liu, H.; Mao, J.; Chu, W.; Li, Q.; Alvarez, P. J. J.; Qu, X.; Zhu, D., Photochemistry of dissolved black carbon released from biochar: Reactive oxygen species generation and phototransformation. Environ. Sci. Technol. 2016, 50, 1218-1226.
428 429
(4) Qu, X.; Fu, H.; Mao, J.; Ran, Y.; Zhang, D.; Zhu, D., Chemical and structural properties of dissolved black carbon released from biochars. Carbon 2016, 96, 759-767.
430 431 432
(5) Jaffe, R.; Ding, Y.; Niggemann, J.; Vahatalo, A. V.; Stubbins, A.; Spencer, R. G. M.; Campbell, J.; Dittmar, T., Global charcoal mobilization from soils via dissolution and riverine transport to the oceans. Science 2013, 340, 345-347.
433 434
(6) Dittmar, T.; Paeng, J., A heat-induced molecular signature in marine dissolved organic matter. Nat. Geosci. 2009, 2, 175-179.
Goldberg, E., Black carbon in the environment. Wiley and Sons, New York, NY: 1985.
ACS Paragon Plus Environment
22
Page 23 of 27
Environmental Science & Technology
435 436
(7) Hopkinson, C. S.; Vallino, J. J., Efficient export of carbon to the deep ocean through dissolved organic matter. Nature 2005, 433, 142-145.
437 438 439
(8) Santin, C.; Doerr, S. H.; Kane, E. S.; Masiello, C. A.; Ohlson, M.; Maria de la Rosa, J.; Preston, C. M.; Dittmar, T., Towards a global assessment of pyrogenic carbon from vegetation fires. Global Change Biol. 2016, 22, 76-91.
440 441
(9) McCarthy, J.; Zachara, J., Es&t features: Subsurface transport of contaminants. Environ. Sci. Technol. 1989, 23, 496-502.
442 443 444
(10) Wang, D. J.; Zhang, W.; Hao, X. Z.; Zhou, D. M., Transport of biochar particles in saturated granular media: Effects of pyrolysis temperature and particle size. Environ. Sci. Technol. 2013, 47, 821828.
445 446 447
(11) Wang, D. J.; Zhang, W.; Zhou, D. M., Antagonistic effects of humic acid and iron oxyhydroxide grain-coating on biochar nanoparticle transport in saturated sand. Environ. Sci. Technol. 2013, 47, 51545161.
448 449
(12) Ward, C. P.; Sleighter, R. L.; Hatcher, P. G.; Cory, R. M., Insights into the complete and partial photooxidation of black carbon in surface waters. Environ. Sci.-Process Impacts 2014, 16, 721-731.
450 451
(13) Wagner, S.; Jaffe, R., Effect of photodegradation on molecular size distribution and quality of dissolved black carbon. Org. Geochem. 2015, 86, 1-4.
452 453
(14) Stubbins, A.; Niggemann, J.; Dittmar, T., Photo-lability of deep ocean dissolved black carbon. Biogeosciences 2012, 9, 1661-1670.
454
(15)
455 456
(16) Stedmon, C. A.; Markager, S., Behaviour of the optical properties of coloured dissolved organic matter under conservative mixing. Estuar. Coast. Shelf Sci. 2003, 57, 973-979.
457 458 459
(17) Guo, W.; Stedmon, C. A.; Han, Y.; Wu, F.; Yu, X.; Hu, M., The conservative and nonconservative behavior of chromophoric dissolved organic matter in chinese estuarine waters. Mar. Chem. 2007, 107, 357-366.
460 461 462
(18) Mantoura, R. F. C.; Woodward, E. M. S., Conservative behaviour of riverine dissolved organic carbon in the severn estuary: Chemical and geochemical implications. Geochim. Cosmochim. Ac. 1983, 47, 1293-1309.
463 464 465
(19) Yang, L. Y.; Guo, W. D.; Hong, H. S.; Wang, G. Z., Non-conservative behaviors of chromophoric dissolved organic matter in a turbid estuary: Roles of multiple biogeochemical processes. Estuar. Coast. Shelf Sci. 2013, 133, 285-292.
466 467
(20) Benner, R.; Opsahl, S., Molecular indicators of the sources and transformations of dissolved organic matter in the mississippi river plume. Org. Geochem. 2001, 32, 597-611.
468 469 470
(21) Yamashita, Y.; Jaffe, R.; Maie, N.; Tanoue, E., Assessing the dynamics of dissolved organic matter (dom) in coastal environments by excitation emission matrix fluorescence and parallel factor analysis (eem-parafac). Limnol. Oceanogr. 2008, 53, 1900-1908.
471 472 473
(22) Stubbins, A.; Spencer, R.; Mann, P. J.; Holmes, R. M.; McClelland, J.; Niggemann, J.; Dittmar, T., Utilizing colored dissolved organic matter to derive dissolved black carbon export by arctic rivers. Front. Earth Sci. 2015, 3, 63.
Wells, M. L., Marine colloids: A neglected dimension. Nature 1998, 391, 530-531.
ACS Paragon Plus Environment
23
Environmental Science & Technology
Page 24 of 27
474 475
(23) Yi, P.; Pignatello, J. J.; Uchimiya, M.; White, J. C., Heteroaggregation of cerium oxide nanoparticles and nanoparticles of pyrolyzed biomass. Environ. Sci. Technol. 2015, 49, 13294-13303.
476 477
(24) Brant, J. A.; Labille, J.; Bottero, J. Y.; Wiesner, M. R., Characterizing the impact of preparation method on fullerene cluster structure and chemistry. Langmuir 2006, 22, 3878-3885.
478 479
(25) Chen, K. L.; Elimelech, M., Aggregation and deposition kinetics of fullerene (c-60) nanoparticles. Langmuir 2006, 22, 10994-11001.
480 481 482
(26) Qu, X.; Hwang, Y. S.; Alvarez, P. J. J.; Bouchard, D.; Li, Q., Uv irradiation and humic acid mediate aggregation of aqueous fullerene (nc60) nanoparticles. Environ. Sci. Technol. 2010, 44, 78217826.
483
(27)
484 485 486
(28) Kim, S. W.; Kaplan, L. A.; Benner, R.; Hatcher, P. G., Hydrogen-deficient molecules in natural riverine water samples - evidence for the existence of black carbon in dom. Mar. Chem. 2004, 92, 225234.
487 488
(29) Dittmar, T., The molecular level determination of black carbon in marine dissolved organic matter. Org. Geochem. 2008, 39, 396-407.
489 490 491
(30) Derjaguin, B.; Landau, L., Theory of stability of highly charged liophobic sols and adhesion of highly charged particles in solutions of electrolytes. Zhurnal Eksperimentalnoi Teor. Fiz. 1945, 15, 663682.
492 493
(31) 49.
494 495
(32) Hotze, E. M.; Phenrat, T.; Lowry, G. V., Nanoparticle aggregation: Challenges to understanding transport and reactivity in the environment. J. Environ. Qual. 2010, 39, 1909-1924.
496 497 498
(33) Petosa, A. R.; Jaisi, D. P.; Quevedo, I. R.; Elimelech, M.; Tufenkji, N., Aggregation and deposition of engineered nanomaterials in aquatic environments: Role of physicochemical interactions. Environ. Sci. Technol. 2010, 44, 6532-6549.
499 500
(34) Chen, K. L.; Elimelech, M., Relating colloidal stability of fullerene (c60) nanoparticles to nanoparticle charge and electrokinetic properties. Environ. Sci. Technol. 2009, 43, 7270-7276.
501 502
(35) van Oss, C. J., Hydrophobicity of biosurfaces — origin, quantitative determination and interaction energies. Colloids Surf., B 1995, 5, 91-110.
503 504
(36) Zhang, Y.; Chen, Y.; Westerhoff, P.; Hristovski, K.; Crittenden, J. C., Stability of commercial metal oxide nanoparticles in water. Water Res. 2008, 42, 2204-2212.
505 506
(37) van Oss, C. J., Acid—base interfacial interactions in aqueous media. Colloids Surf., A 1993, 78, 1-49.
507 508
(38) Hardy, W. B., "A preliminary investigation of the conditions which determine the stability of irreversible hydrosols ". Proc. R. Soc. London 1900, 66, 110-125.
509 510
(39) Chin, W.-C.; Orellana, M. V.; Verdugo, P., Spontaneous assembly of marine dissolved organic matter into polymer gels. Nature 1998, 391, 568-572.
511 512
(40) Engebretson, R. R.; Amos, T.; vonWandruszka, R., Quantitative approach to humic acid associations. Environ. Sci. Technol. 1996, 30, 990-997.
IHSS International humic substance society website. www.humicsubstances.org/
Verwey, E. J. W., Theory of the stability of lyophobic colloids. Philips Res. Repts. 1945, 1, 33-
ACS Paragon Plus Environment
24
Page 25 of 27
Environmental Science & Technology
513 514
(41) Hsiao, P.-Y.; Luijten, E., Salt-induced collapse and reexpansion of highly charged flexible polyelectrolytes. Phys. Rev. Lett. 2006, 97, 148301.
515 516
(42) Besteman, K.; Van Eijk, K.; Lemay, S. G., Charge inversion accompanies DNA condensation by multivalent ions. Nat. Phys. 2007, 3, 641-644.
517 518 519
(43) Hierrezuelo, J.; Vaccaro, A.; Borkovec, M., Stability of negatively charged latex particles in the presence of a strong cationic polyelectrolyte at elevated ionic strengths. J. Colloid Interface Sci. 2010, 347, 202-208.
520 521 522
(44) Szilagyi, I.; Polomska, A.; Citherlet, D.; Sadeghpour, A.; Borkovec, M., Charging and aggregation of negatively charged colloidal latex particles in the presence of multivalent oligoamine cations. J. Colloid Interface Sci. 2013, 392, 34-41.
523 524 525
(45) Schneider, C.; Hanisch, M.; Wedel, B.; Jusufi, A.; Ballauff, M., Experimental study of electrostatically stabilized colloidal particles: Colloidal stability and charge reversal. J. Colloid Interface Sci. 2011, 358, 62-67.
526 527
(46) Lyklema, J., Overcharging, charge reversal: Chemistry or physics? Colloids Surf., A 2006, 291, 3-12.
528 529 530
(47) Borkovec, M.; Szilagyi, I.; Popa, I.; Finessi, M.; Sinha, P.; Maroni, P.; Papastavrou, G., Investigating forces between charged particles in the presence of oppositely charged polyelectrolytes with the multi-particle colloidal probe technique. Adv. Colloid Interface Sci. 2012, 179–182, 85-98.
531 532
(48) Kalinichev, A. G.; Kirkpatrick, R. J., Molecular dynamics simulation of cationic complexation with natural organic matter. Eur. J. Soil Sci. 2007, 58, 909-917.
533 534 535
(49) Ahn, W.-Y.; Kalinichev, A. G.; Clark, M. M., Effects of background cations on the fouling of polyethersulfone membranes by natural organic matter: Experimental and molecular modeling study. J. Membr. Sci. 2008, 309, 128-140.
536 537 538
(50) Iskrenova-Tchoukova, E.; Kalinichev, A. G.; Kirkpatrick, R. J., Metal cation complexation with natural organic matter in aqueous solutions: Molecular dynamics simulations and potentials of mean force. Langmuir 2010, 26, 15909-15919.
539 540
(51) Schnitzer, M.; Hansen, E., Organo-metallic interactions in soils: 8. An evaluation of methods for the determination of stability constants of metal-fulvic acid complexes. Soil Sci. 1970, 109, 333-340.
541 542
(52) Hering, J. G.; Morel, F. M. M., Humic-acid complexation of calcium and copper. Environ. Sci. Technol. 1988, 22, 1234-1237.
543 544
(53) Celebi, O.; Kilikli, A.; Erten, H. N., Sorption of radioactive cesium and barium ions onto solid humic acid. J. Hazard. Mater. 2009, 168, 695-703.
545 546
(54) Ohki, S., A mechanism of divalent ion-induced phosphatidylserine membrane fusion. Biochimica et Biophysica Acta (BBA) - Biomembranes 1982, 689, 1-11.
547 548
(55) Van Oss, C. J., Hydrophobic, hydrophilic and other interactions in epitope-paratope binding. Mol. Immunol. 1995, 32, 199-211.
549 550
(56) Kloster, N.; Brigante, M.; Zanini, G.; Avena, M., Aggregation kinetics of humic acids in the presence of calcium ions. Colloids Surf., A 2013, 427, 76-82.
ACS Paragon Plus Environment
25
Environmental Science & Technology
Page 26 of 27
551 552 553
(57) Schmitt-Kopplin, P.; Hertkorn, N.; Schulten, H. R.; Kettrup, A., Structural changes in a dissolved soil humic acid during photochemical degradation processes under O2 and N2 atmosphere. Environ. Sci. Technol. 1998, 32, 2531-2541.
554 555
(58) L. Malcolm, R., The uniqueness of humic substances in each of soil, stream and marine environments. Anal. Chim. Acta 1990, 232, 19-30.
556 557 558 559
(59) Stubbins, A.; Spencer, R. G. M.; Chen, H. M.; Hatcher, P. G.; Mopper, K.; Hernes, P. J.; Mwamba, V. L.; Mangangu, A. M.; Wabakanghanzi, J. N.; Six, J., Illuminated darkness: Molecular signatures of congo river dissolved organic matter and its photochemical alteration as revealed by ultrahigh precision mass spectrometry. Limnol. Oceanogr. 2010, 55, 1467-1477.
560 561 562
ACS Paragon Plus Environment
26
Page 27 of 27
563
Environmental Science & Technology
TOC
564
ACS Paragon Plus Environment
27