An Integrated Model for Input and Migration of Mercury in Chinese

Jan 28, 2019 - Coastal sediments are a major sink of the global mercury (Hg) biogeochemical cycle, bridging terrestrial Hg migration to the open ocean...
0 downloads 0 Views 2MB Size
Subscriber access provided by AUT Library

Environmental Processes

An Integrated Model for Input and Migration of Mercury in Chinese Coastal Sediments Mei Meng, Ruoyu Sun, Hongwei Liu, Ben Yu, Yongguang Yin, Ligang Hu, Jianbo Shi, and Guibin Jiang Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b06329 • Publication Date (Web): 28 Jan 2019 Downloaded from http://pubs.acs.org on January 29, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

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 36

Environmental Science & Technology

1

An Integrated Model for Input and Migration of Mercury in

2

Chinese Coastal Sediments

3 4 5

Mei Meng†, ‡, Ruo-yu Sun†, Hong-wei Liu‡, Ben Yu‡, Yong-guang Yin‡, Li-gang Hu‡,

6

Jian-bo Shi*, ‡, §, Gui-bin Jiang‡

7 8 9 10

† Institute

11



12

Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China

13

§ Institute

of Surface-Earth System Science, Tianjin University, Tianjin 300072, China

State Key Laboratory of Environmental Chemistry and Ecotoxicology, Research Center for

of Environment and Health, Jianghan University, Wuhan 430056, China

14 15 16 17

* Corresponding author (J.B. Shi)

18

Tel/fax: +86-10-62849129

19

E-mail: [email protected]

1

ACS Paragon Plus Environment

Environmental Science & Technology

20

TOC

21

22 23

2

ACS Paragon Plus Environment

Page 2 of 36

Page 3 of 36

Environmental Science & Technology

ABSTRACT

24 25

Coastal sediments are a major sink of global mercury (Hg) biogeochemical cycle,

26

bridging terrestrial Hg migration to the open ocean. It is thus of substantial interest to quantify

27

the Hg contributors to coastal sediments and the extents to which the Hg sequestered into

28

coastal sediments affects the ocean. Here, we measured concentrations and isotope

29

compositions of Hg in Chinese coastal sediments, and found that estuary sediments had

30

distinctly higher δ202Hg and lower Δ199Hg values than marine sediments. Hg isotope

31

compositions of marine sediments followed a latitudinal trend where δ202Hg decreases and

32

Δ199Hg increases from north to south. An integrated model was developed based on Hg

33

isotope mixing model and urban distance factor (UDF), which revealed a significant

34

difference in Hg source contributions among the estuary and marine sediments, and a gradual

35

change of dominant Hg sources from terrestrial inputs (riverine and industrial wastewater

36

discharges) to atmospheric deposition with the decrease of urban impact. An UDF value of

37

306 ± 217 was established as the critical point where dominant Hg sources started to change

38

from terrestrial inputs to atmospheric deposition. Our study help understand the input,

39

migration of Hg in Chinese marginal seas, and provide critical insights for targeted

40

environmental management.

3

ACS Paragon Plus Environment

Environmental Science & Technology

41

■ INTRODUCTION

42

Ocean serves as a central role in mediating the global cycling of mercury (Hg). It receives

43

atmospheric Hg dry/wet deposition (5000-6500 Mg yr-1) and riverine Hg discharges

44

(1000-5500 Mg yr-1), and meanwhile emits dissolved elemental Hg0 (4000-5000 Mg yr-1) into

45

the atmosphere.1-4 As estimated by Amos et al. (2014),1 about 70% of riverine Hg is buried in

46

the worldwide coastal sediments, with the remainder exported to the open ocean. Therefore,

47

burial of Hg in coastal sediments represents a major sink of global Hg biogeochemical cycle.

48

A small fraction of accumulated Hg in sediments and oceanic water columns (typically
500 km and an average water depth of 60 m.65

124

The Yangtze River, the largest inflowing river, discharges 5.0×108 tons of sediments

125

annually.65 The formed Yangtze River Estuary plays a vital role in delivering various

126

contaminants to the adjacent sea.66-68 Dominated by subtropical climate, the ECS has high

127

temperature (20-26 C) and salinity (> 34%). The SCS is characterized by a tropical marine

128

monsoon climate, with all-year-round high temperature (25-28 C), abundant rainfall 7

ACS Paragon Plus Environment

Environmental Science & Technology

129

(1100-2500 mm yr-1), and significant influence from typhoon. The Pearl River is the largest

130

inflowing river, bringing ~9.0×107 tons of sediments annually into the Pearl River Estuary, a

131

seriously contaminated estuary in China.58, 69, 70

132

Sample Collection and Pretreatment. The sampling campaigns were conducted during

133

two cruises in June 2011 and May 2012. A total of 141 surface sediments (0-3 cm) were

134

collected using a box corer from BS (n = 29, from 12.5 m to 38 m), YS (n = 55, from 20.7 m

135

to 79 m), and ECS (n = 37, from 39.5 m to 99.5 m) and its Yangtze River Estuary (n = 20)

136

(Figure 1). The sedimentation rates in the studied sea area varied from 0.23 to 2.86 cm

137

yr-1,71-76 and thus the collected sediments were expected to integrate deposition from 1998 to

138

2012, reflecting surrounding Hg inputs from recent 1-2 decades. The surface sediments from

139

SCS (n = 6, from 2918 m to 3739 m) and its Pearl River Estuary (n = 37) reported by Yin et al.

140

(2015)58 are also shown in Figure 1. These sampling sites cover most of Chinese marginal sea

141

areas. We distinguish surface sediments as estuary sediments and marine sediments for

142

samples collected from estuaries and seas, respectively. Besides, seven surface soil samples

143

(0-2 cm) were collected from the industrial areas of Dalian City in the northeast coast of BS

144

(Figure 1). Upon collection, all the samples were kept in a -20 C refrigerator immediately

145

and transported to the laboratory, where they were freeze-dried at -45 C, ground and

146

homogenized to less than 80 meshes before analysis.

147

Measurement of Total Hg and Hg Isotope Ratio. A Hydra-C Hg analyzer (Teledyne

148

Lemman Labs, USA) was first used for direct analysis of total Hg (THg) concentration in

149

sediments and soils, following the USEPA method 7473.77 Then, all the sediments and soils

150

were digested using aqua regia in a water bath of 95 C for 4 h. THg concentrations in 8

ACS Paragon Plus Environment

Page 8 of 36

Page 9 of 36

Environmental Science & Technology

151

digested sample solutions were determined with an inductively coupled plasma mass

152

spectrometry (ICP-MS, Thermo iCAP-Q, Thermo Fisher Scientific, USA), which agreed well

153

with THg concentrations measured by Hydra-C Hg analyzer. Before Hg isotope ratio

154

measurement by a multi-collector (MC) ICP-MS (Nu-Plasma Ⅱ, Nu Instruments, UK), the

155

digested sample solutions were diluted to Hg concentrations of 1-3 ng mL-1. Mass bias of

156

MC-ICP-MS was corrected by internal NIST-997 Tl standard using the exponential law and

157

by NIST-3133 Hg standard using the standard-sample bracketing method. Details on the

158

analytical procedures for THg concentrations and Hg isotope ratios in the sediments and soils

159

are provided in the Supporting Information (SI). We deem the larger 2SD uncertainties of Hg

160

isotope ratios in UM-Almadén and GBW07310 (sediment reference) as the typical analytic

161

uncertainties of Hg isotope ratios in sediments and soils, which are 0.14‰ for δ202Hg, 0.04‰

162

for Δ199Hg, 0.03‰ for Δ200Hg, and 0.04‰ for Δ201Hg. The 2SE uncertainties of Hg isotope

163

ratios in samples with multiple analysis are applied only when they are larger than the typical

164

analytic uncertainties.

165

Construction of An Integrated SC-UDF Model. To unravel the contributed percentage

166

of possible Hg sources in the sediments of the different sea areas and to depict how the urban

167

activities impact Hg migration from nearshore to offshore region, an integrated SC-UDF

168

model was constructed through combining the Hg isotope source contribution (SC) model and

169

urban distance factor (UDF) model.

170

The source contributions were modeled by a Monte Carlo simulation approach used by

171

Jiskra et al. (2015, 2017),44, 78 through the pseudorandom number generation function (i = 1 :

172

1,000,000) of the MatLab software (R2016a, MathWorks) based on the following equations: 9

ACS Paragon Plus Environment

Environmental Science & Technology

n

173

 f * IR = IR i

i

f

1

Page 10 of 36

(1)

sample

i 1 n

174

i

(2)

i 1

175

where i is the number of Hg source end-member, fi is the fraction of each end-member, IRi is

176

the isotope ratios of each end-member, and IRsample is the measured isotope ratios of samples.

177

The isotope ratios are δ202Hg and ΔxxxHg (xxx = 199, 200, 201, 204).

178

The degree of urban impact was modeled by an urban distance factor (UDF) based on the

179

city population (P) and the distances (D) between sampling sites and potential source cities.

180

Seven forms of UDF were calculated for each sampling site (see SI), and the UDF1 (= √P / D)

181

was selected here, considering the large population base in Chinese cities. We established a

182

database of UDF1 for the Chinese marginal sea areas, using the updated population data (i.e.,

183

the 2010 China Census) for all the 155 potential source cities (Table S2) along the whole

184

eastern and southern coastal regions in China (Figure S1). Details on the calculation process

185

(Table S3) and the final UDF1 database (Table S4) are provided in the SI.

186 187

■ RESULTS AND DISCUSSION

188

Variation of THg Concentrations in Surface Sediments. THg concentrations in all the

189

surface sediments ranged from 7.0 to 159.6 μg kg-1, with a mean value of 31.5 ± 19.1 μg kg-1

190

(n = 141) (Table 1). One sample (B7, Figure 1) collected at the near-lowest water depth (12.5

191

m) of BS was found to be highly enriched in Hg as compared to other samples (159.6 μg kg-1

192

vs. 7.0 to 84.3 μg kg-1). On average, the Yangtze River Estuary sediments had the highest

193

THg concentrations of 48.8 ± 18.7 μg kg-1 (n = 20), followed by BS of 38.7 ± 26.3 μg kg-1 (n 10

ACS Paragon Plus Environment

Page 11 of 36

Environmental Science & Technology

194

= 29), ECS of 27.3 ± 15.0 μg kg-1 (n = 37), and YS of 24.2 ± 9.8 μg kg-1 (n = 55). In

195

comparison, the SCS and its Pearl River Estuary sediments reported by Yin et al. (2015)58 had

196

average THg concentrations of 75.5 ± 6.3 μg kg-1 (n = 6) and 98.4 ± 39.4 μg kg-1 (n = 37),

197

respectively (Table 1; Figure S2-S3). This indicated that Hg levels in sediments from the

198

south coast were more elevated than those from the east coast of China. The contour map of

199

THg (Figure 2) clearly shows that the sediments near the estuaries (Yellow River Estuary and

200

Yangtze River Estuary) and along the ECS coast had obviously higher THg concentrations

201

than those of offshore areas, suggesting that terrestrial inputs by river runoffs and wastewater

202

discharges were the main sources of Hg in coastal areas. As estimated by Liu et al. (2016),12

203

the annual terrestrial Hg exported from mainland China into adjacent seas were 24 Mg for BS,

204

25 Mg for YS, and 125 Mg for ECS in 2012. By normalizing the annual terrestrial Hg input

205

by water area12 in each sea basin, we found that BS had the highest terrestrial Hg input rate of

206

0.3 kg m-2 yr-1, followed by ECS of 0.16 kg m-2 yr-1 and YS of 0.07 kg m-2 yr-1. The terrestrial

207

Hg input rates are positively correlated (R2 = 0.96, P = 0.12, n = 3) with mean THg

208

concentrations of marine sediments in the studied three marginal seas, highlighting the

209

importance of terrestrial inputs in contributing Hg enrichment in sediments.

210

Among the terrestrial Hg inputs, riverine runoffs have been shown to be the dominant Hg

211

sources.1, 12 The delivered Hg might combine with organic matter of river suspended loads, as

212

seen from the positive correlation between THg concentrations and total organic carbon (TOC)

213

contents in all sediments (R2 = 0.11, P < 0.01, n = 141). In addition, THg concentrations in

214

surface sediments have a negative correlation with co-located pH values (R2 = 0.08, P < 0.01,

215

n = 141), with high pH values generally corresponding to low sediment THg concentrations. 11

ACS Paragon Plus Environment

Environmental Science & Technology

216

The above observations corroborate previous studies demonstrating that clay sediments with

217

high TOC enriched in humic acid and fulvic acid are the main Hg sink of Chinese marginal

218

sea.16 However, caution should be taken to the obtained correlations, because TOC and pH

219

each only explains ~10% of THg variance in sediments. Despite less affected by terrestrial

220

inputs, the sediments in central YS (Y3, Y5, Y6, Y8, Y9, Y10, Y12, Y13, Y17, Y51, Y52)

221

had higher (t-test, P < 0.05) THg concentrations (37.1 ± 5.7 μg kg-1) than sediments near the

222

coast (21.0 ± 7.8 μg kg-1) (Figure 2; Figure S2). We attribute this contrast to the following

223

reasons. First, the low terrestrial Hg inputs (absence of large rivers)12 and large surface area of

224

YS could have diminished the contribution of terrestrial inputs to the Hg levels in the coastal

225

areas. Second, atmospheric Hg deposition accounts for a large fraction of Hg inputs in YS,

226

which has been observed to be enhanced in central YS.75 Third, an offshore mud area

227

containing significant high TOC (Figure 2; Figure S4) is frequently generated in central YS

228

under the influence of current circulation, serving as a scavenger of seawater Hg.

229

Variation of Hg Isotope Ratios in Surface Sediments. Surface sediments exhibited a

230

large variation in both δ202Hg (-2.43‰ to -0.09‰) and Δ199Hg (-0.08‰ to 0.31‰) (Table 1).

231

Δ200Hg values in surface sediments were very small from -0.04‰ to 0.09‰, and most were

232

within the measurement uncertainty. It is noted that the sample (B7) with the highest THg

233

concentration (159.6 μg kg-1) had the highest δ202Hg value of -0.09‰. Different from the

234

estuary sediments exhibiting small and negative Δ199Hg values, the marine sediments were

235

typically characterized by positive Δ199Hg values.

236

On average, δ202Hg values in surface sediments were high in BS (-0.93 ± 0.33‰, -1.48‰

237

to -0.09‰) and Yangtze River Estuary (-0.98 ± 0.15‰, -1.54‰ to -0.74‰), and low in ECS 12

ACS Paragon Plus Environment

Page 12 of 36

Page 13 of 36

Environmental Science & Technology

238

(-1.63 ± 0.30‰, -2.19‰ to -1.13‰) and YS (-1.82 ± 0.45‰, -2.43‰ to -0.47‰) (Table 1). In

239

contrast, Δ199Hg values in surface sediments were highest in YS (0.19 ± 0.05‰, 0.07‰ to

240

0.31‰) followed by ECS (0.12 ± 0.05‰, 0.03‰ to 0.22‰), BS (0.09 ± 0.05‰, -0.01‰ to

241

0.19‰) and Yangtze River Estuary (-0.02 ± 0.02‰, -0.08‰ to -0.01‰) (Table 1). Thus, it

242

appears that the surface sediments with elevated Hg were more enriched in heavy Hg isotopes

243

and had small or no MIF anomalies. Similar phenomenon has also been observed in SCS and

244

other marginal seas.54, 57, 58 For example, Yin et al. (2015)58 showed that the sediments from

245

SCS had low δ202Hg of -2.44‰ and high Δ199Hg of 0.35‰, whereas sediments from its Pearl

246

River Estuary had high δ202Hg of -1.60‰ and near-zero Δ199Hg of -0.01‰.

247

The linear regressions of δ202Hg versus 1/THg and Δ199Hg versus 1/THg in each estuary

248

and sea were performed and shown in Figure S5. Only THg concentrations in sediments from

249

Pearl River Estuary studied by Yin et al. (2015)58 showed significant correlations with both

250

δ202Hg (R2 = 0.65, P < 0.01, n = 37) and Δ199Hg (R2 = 0.33, P < 0.01, n = 37), suggesting a

251

mixing of binary Hg sources. THg concentrations in sediments from BS significantly

252

correlated with Δ199Hg (R2 = 0.14, P < 0.05, n = 29) rather than δ202Hg. THg concentrations in

253

sediments from other areas showed no correlations with either δ202Hg or Δ199Hg. Thus, we

254

suggest that the variations in Hg isotope compositions of our studied sediments can not be

255

convincingly explained by two sources.

256

Potential Hg Sources and Influencing Factors. The Yangtze River Estuary in ECS

257

primarily receives terrestrial inputs brought by the Yangtze River, the largest river in China.

258

The formed “Yangtze River Delta” is surrounded by the highly-developed industrial and

259

urbanized cities in Jiangsu and Zhejiang Provinces and Shanghai city. Discharges of industrial 13

ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 36

260

wastewater and domestic sewage into aquatic environments and soils would be carried away

261

by terrestrial runoffs and accumulated into the Yangtze River Estuary. Besides, the Yangtze

262

River Estuary is characterized by high total suspended particulates (0.7-51.5 g L-1),79 high

263

sedimentation rate (1.36 to 4.11 cm yr-1),80 and low water clarity (< 100 FTU).81 Thus,

264

Yangtze River Estuary sediments are expected to conserve the Hg isotope signatures of

265

terrestrial inputs containing large amounts of anthropogenic Hg. The similarity in δ202Hg

266

between Yangtze River Estuary sediments and BS marine sediments and their high THg

267

concentrations suggest that BS was seriously affected by terrestrial Hg inputs as well. Unlike

268

other Chinese marginal seas, BS is a nearly full-enclosed sea surrounded by a large number of

269

industrialized and densely populated cities in Liaoning, Hebei and Shandong Provinces, and

270

receives terrestrial discharges accumulated in the Yellow River Estuary. As shown in Figure

271

S4, the suspended loads along the coastlines and from the Yellow River Estuary could be well

272

dispersed into offshore areas of BS under the influence of Bohai Coastal Current (BCC).

273

Different from the Yangtze River Estuary sediments, the marine sediments of BS had

274

small, positive Δ199Hg of 0.09 ± 0.05‰ (-0.01‰ to 0.19‰). Positive but slightly higher

275

Δ199Hg values were also observed in the marine sediments of YS (0.19‰, 0.07‰ to 0.31‰)

276

and ECS (0.12‰, 0.03‰ to 0.22‰) of more open water areas in this study, and deep marine

277

sediments in SCS (0.35‰, 0.21‰ to 0.45‰) reported by Yin et al. (2015).58 In addition, the

278

marine sediments in other marginal seas were found with positive Δ199Hg as well.53,

279

Strock et al. (2015)39 showed that seawater samples from Canadian Arctic Archipelago were

280

typically characterized with positive Δ199Hg (0.1‰ to 0.4‰). Thus, Δ199Hg of marine

281

sediments was largely imparted by overlying seawater Δ199Hg. 14

ACS Paragon Plus Environment

56, 57

Page 15 of 36

Environmental Science & Technology

282

As shown in Figure S6, Δ199Hg values in marine sediments of Chinese seas had the

283

highest correlation (R2 = 0.26, P < 0.01, n = 141) with co-located water depths, suggesting

284

that the offshore and semi-pelagic seawater had higher Δ199Hg than coastal seawater. Hg in

285

offshore and semi-pelagic seawater largely derives from atmospheric wet Hg deposition with

286

large positive Δ199Hg, whereas Hg in coastal seawater is seriously affected by terrestrial Hg

287

with negative to near-zero Δ199Hg.57 We also found that Δ199Hg values in marine sediments

288

were positively correlated with TOC and negatively correlated with pH (Figure S6),

289

suggesting Hg in offshore seawater was probably sequestrated by organic matter in marine

290

sediments. Notably, the sediments in central YS with the highest TOC and lowest pH values

291

had very positive Δ199Hg values (Figure 2).

292

The linear regressions of Δ199Hg versus Δ201Hg (Figure S7) in sediments from BS (slope =

293

0.95; R2 = 0.54, P < 0.01, n = 29), YS (slope = 1.01; R2 = 0.32, P < 0.01, n = 55) and ECS

294

(slope = 1.07; R2 = 0.75, P < 0.01, n = 37) in this study, and from SCS/Pearl River Estuary

295

(slope = 1.05; R2 = 0.97, P < 0.01, n = 43) in Yin et al. (2015)58 all show a Δ199Hg/Δ201Hg

296

slope near the unity, demonstrating that Hg had undergone photo-reduction processes before

297

incorporation into the sedimentary particles in these sea areas. Photo-reduction of

298

post-depositional Hg in sediments is unlikely significant, because the sedimentary Hg is

299

thought to be less reactive than dissolved Hg, and has a strong complexion with organic

300

matter and sulfides in sediments. Additionally, the low clarity of coastal seawater would

301

greatly attenuate the intensity of light penetrated into the seawater, decreasing the extents of

302

photo-reduction of sedimentary Hg.

303

Integrated SC-UDF Model. Here, an integrated SC-UDF model was developed for 15

ACS Paragon Plus Environment

Environmental Science & Technology

304

Chinese marginal sea area using (i) Hg source contribution (SC) and (ii) urban distance factor

305

(UDF) database, as follows:

306

(i) Hg Source Contribution. Binary or ternary mixing model has been established to

307

identify Hg sources and quantify their contributions in sediments if different Hg sources have

308

characteristic and distinguishable isotope compositions.43, 58, 82 Since geochemical processes

309

could alter the Hg isotope compositions of original sources, the appointment of Hg

310

end-members in mixing models should reasonably correct Hg isotope fractionation attributed

311

to processes occurring before and after source Hg deposition.83 For example, Lepak et al.

312

(2015)43 assumed that adsorption of dissolved Hg onto particles could cause an enrichment of

313

lighter Hg isotopes in particles, and applied a negative MDF offset (-0.40‰ in δ202Hg) to

314

correct precipitation Hg deposited into the Great Lakes sediments in a ternary mixing model.

315

This MDF offset is based on a laboratory-controlled experiment on dissolved Hg adsorption

316

onto goethite reported by Jiskra et al. (2012).26 However, Washburn et al. (2017)84 have

317

recently observed a consistent and positive MDF offset (0.28‰ in δ202Hg) from dissolved Hg

318

to suspended particulate bound Hg in a contaminated river system of Virginia. This indicates

319

that the correction of source Hg isotope compositions might involve with considerable

320

uncertainty. For this reason, the use of the samples that integrate isotope signatures of both

321

potential Hg sources and the subsequent transport and transformation processes as the

322

end-members is more plausible than the use of sources itself or sources only corrected by a

323

process offset.85 Here, we constructed an improved ternary mixing model for surface

324

sediments of Chinese marginal seas (Figure 3), with the three end-members determined based

325

on soil samples in this study (representing direct discharge of industrial Hg into seawater), 16

ACS Paragon Plus Environment

Page 16 of 36

Page 17 of 36

Environmental Science & Technology

326

riverine sediments taken from a contaminated river49 (representing riverine Hg inputs), and

327

deep-sea sediments58 (representing atmospheric Hg deposition) (Table S6).

328

(1) Direct Discharge of Industrial Hg. THg concentrations in studied surface soils near

329

the coastal city varied over a wide range of 23.4-912.7 μg kg-1, with a mean value of 295.5 ±

330

386.9 μg kg-1 (Table 1). Two (So3 and So4, in bold, Table S5) of these soil samples with

331

extremely high THg levels (897.5 μg kg-1; 912.7 μg kg-1) were sampled near large wastewater

332

outflows of an industrial complex, which were probably severely contaminated by Hg in the

333

discharged wastewater. These surface coastal soils had a narrow variation range between

334

-0.67‰ and -0.37‰ in δ202Hg, and negligible MIF (-0.02‰ to 0.04‰ for Δ199Hg; -0.03‰ to

335

0.02‰ for Δ201Hg) (Table 1). The mean Hg isotope compositions (δ202Hg = -0.50 ± 0.10‰;

336

Δ199Hg = 0.01 ± 0.02‰) of these soils are comparable to those of commercial liquid Hg

337

(δ202Hg = -0.40‰; Δ199Hg = -0.02‰), Hg ores (δ202Hg = -0.59‰; Δ199Hg = 0.03‰), and

338

other non-ferrous metal ores (δ202Hg = -0.48‰, Δ199Hg = 0.01‰) as summarized in Sun et al.

339

(2016).86 We thus used Hg isotope compositions of these measured coastal soils to represent

340

those of discharged wastewater from industries.

341

(2) Riverine Hg Inputs. Hg isotope compositions in dissolved and particulate-bound Hg

342

in Chinese rivers are currently not available. Because particulate-bound Hg is the primary

343

form of Hg in rivers and most of particulate-bound Hg is buried in the coastal sediments,1 we

344

used Hg isotope compositions (δ202Hg = -1.65 ± 0.54‰; Δ199Hg = -0.11 ± 0.09‰) of surface

345

sediments from a contaminated river reported by Liu et al. (2011)49 to represent riverine Hg

346

inputs. The chosen river flows through urban, industrial, and background regions and is the

347

main inflowing river to Pearl River Estuary. Its sediments can be deemed as the best 17

ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 36

348

integrators of isotope signatures of mixed anthropogenic Hg sources. We expect that rivers

349

flowing into other Chinese estuaries were contaminated by similar anthropogenic Hg sources.

350

Our used values for riverine Hg inputs are comparable to Hg isotope compositions of

351

watershed soils (δ202Hg: -1.82 ± 0.39‰; Δ199Hg: -0.29 ± 0.12‰) that were compiled globally

352

and used in a ternary mixing model for Hg tracing in Chinese coastal sediment cores.82

353

(3) Atmospheric Hg Deposition. Atmospheric Hg deposition is an important Hg source

354

to oceans. However, the atmospheric Hg incorporated into sediments should have quite

355

different Hg isotope compositions from the original atmospheric Hg, because atmospheric

356

dry/wet deposition processes, the following Hg transformation (e.g., photo-reduction and

357

adsorption) in seawater columns and post-depositional processes are expected to cause

358

significant Hg isotope fractionation. Thus, we used Hg isotope compositions (δ202Hg = -2.44

359

± 0.25‰; Δ199Hg = 0.35 ± 0.08‰) in the deep-sea offshore sediments (2918-3739 m) of

360

SCS58 to represent those of atmospheric Hg deposition. The sediments in SCS were collected

361

from deep-sea offshore regions, and were mainly affected by Hg from atmospheric deposition

362

with limited terrestrial Hg contribution.

363 364

According to the isotope mixing model presented in equations 1-2, the fractions of the above three end-members (sources) were simulated according to the following equations:

365

f1 * δ202Hg1 + f2 * δ202Hg2 + f3 * δ202Hg3 = δ202Hgsample

(3)

366

f1 * Δ199Hg1 + f2 * Δ199Hg2 + f3 * Δ199Hg3 = Δ199Hgsample

(4)

367

f1 * Δ201Hg1 + f2 * Δ201Hg2 + f3 * Δ201Hg3 = Δ201Hgsample

(5)

368

f1 + f2 + f3 = 1

369

The simulated results (mean ± 1SD) for fractions of direct industrial Hg discharge (f1),

(6)

18

ACS Paragon Plus Environment

Page 19 of 36

Environmental Science & Technology

370

riverine Hg (f2), and atmospheric Hg (f3) were summarized in Table S7 and shown in Figure

371

S8. Terrestrial inputs were the main Hg sources in BS and Yangtze River Estuary and Pearl

372

River Estuary (Figure S8), with direct industrial wastewater Hg dominated in BS and riverine

373

Hg dominated in these two estuaries. SCS sediments were nearly exclusively contributed by

374

atmospheric Hg deposition (Figure S8), because we used them as the atmospheric Hg

375

end-member. Atmospheric Hg deposition was also the largest Hg source in most of the

376

sediments from YS, although Hg contribution from riverine inputs was significant as well

377

(Figure S8). Notably, some sediments in YS near Shandong/Liaoning Provinces had the

378

largest Hg contribution from industrial discharge. For ECS, the largest contributors changed

379

from terrestrial inputs in estuary sediments (on average, 92.7% in Yangtze River Estuary) to

380

atmospheric deposition in offshore sediments (on average, 50.9% in sites: E8, E9, E13, E14,

381

E18, E19, E26, E27, Figure 1; in bold, Table S5) (Figure S8). A transition of dominant Hg

382

sources from terrestrial inputs to atmospheric deposition in ECS was estimated to be 367 ± 24

383

km (Figure S8).

384

(ii) UDF Database. The Pearl River Estuary sampling sites showed the highest UDF1

385

values (738.1, 640.4 to 795.5), suggesting the greatest urban impact. The BS sampling sites

386

had the next high UDF1 values (473.3, 434.3 to 499.9), followed by sampling sites from

387

Yangtze River Estuary (425.2, 390.1 to 476.3), YS (407.0, 342.8 to 517.8), ECS (357.4, 276.2

388

to 456.9), and SCS (261.0, 228.0 to 283.7).

389

Based on the calculation results of SC and UDF1, an integrated SC-UDF model was

390

developed and shown in Figure 4. The left side of Figure 4 displayed both the Hg source

391

contribution percentage (pie-map) and the urban impact degree (UDF1, contour-map). In 19

ACS Paragon Plus Environment

Environmental Science & Technology

392

general, a decreasing degree of urban impact corresponds to a change in dominant Hg sources

393

from terrestrial inputs (i.e., the sum of industrial wastewater discharge and riverine inputs) to

394

atmospheric deposition. The right side of Figure 4 showed the relationships of UDF1 with

395

fractions of terrestrial Hg and atmospheric Hg. With the decrease of UDF1, the contribution

396

of terrestrial Hg showed a significantly decreasing trend (y = 0.163x + 0.10, R2 = 0.26, P