Aggregation Behavior of Dissolved Black Carbon: Implications for

The fluvial export of dissolved black carbon (DBC) is a major land-ocean flux in the global black. 17 carbon cycle, affecting the size of refractory c...
4 downloads 0 Views 662KB Size
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