Mercury Release to Aquatic Environments from Anthropogenic

Jul 5, 2016 - Based on an analysis of measured data and distribution factors, we developed the China Aquatic Mercury Release (CAMR) model, which we us...
0 downloads 9 Views 3MB Size
Subscriber access provided by La Trobe University Library

Article

Mercury Release to Aquatic Environments from Anthropogenic Sources in China from 2001 to 2012 Maodian Liu, Wei Zhang, Xuejun Wang, Long Chen, Huanhuan Wang, Yao Luo, Haoran Zhang, Huizhong Shen, Yindong Tong, Langbo Ou, Han Xie, Xuejie Ye, and Chunyan Deng Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b01386 • Publication Date (Web): 05 Jul 2016 Downloaded from http://pubs.acs.org on July 9, 2016

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

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32

Environmental Science & Technology

1

Mercury Release to Aquatic Environments from Anthropogenic Sources in China

2

from 2001 to 2012

3

Maodian Liu,† Wei Zhang,*,‡ Xuejun Wang,*,† Long Chen,† Huanhuan Wang,† Yao

4

Luo,† Haoran Zhang,† Huizhong Shen,† Yindong Tong,§ Langbo Ou,† Han Xie,†

5

Xuejie Ye† and Chunyan Deng†

6



7

Environmental Sciences, Peking University, Beijing 100871, China

8



9

100872, China

Ministry of Education Laboratory of Earth Surface Processes, College of Urban and

School of Environment and Natural Resources, Renmin University of China, Beijing

10

§

11

300072, China

12

Corresponding authors:

13

Wei Zhang. School of Environment and Natural Resources, Renmin University of

14

China, Beijing, 100872, China. Tel: +86-10-62514030. E-mail: [email protected]

15

Xuejun Wang. Ministry of Education Laboratory of Earth Surface Processes, College

16

of Urban and Environmental Sciences, Peking University, Beijing, 100871, China.

17

Tel: +86-10-62759190. E-mail: [email protected]

18

Word count for text:

4475

19

Word count for 6 figures:

6×300=1800

20

Total Word count:

6275

School of Environmental Science and Engineering, Tianjin University, Tianjin

1

ACS Paragon Plus Environment

Environmental Science & Technology

21

ABSTRACT: Based on an analysis of measured data and distribution factors, we

22

developed the China Aquatic Mercury Release (CAMR) model, which we used to

23

calculate an inventory of mercury (Hg) that was released to aquatic environments

24

from primary anthropogenic sources in China. We estimated a total release of 98 tons

25

of Hg in 2012, including coal-fired power plants (17%), nonferrous metal smelting

26

(33%), coal mining and washing (25%), domestic sewage (17%), and other sectors

27

(8.3%). The total primary anthropogenic Hg released to aquatic environments in

28

China decreased at an annual average rate of 1.7% between 2001 and 2012, even

29

though GDP grew at an annual average rate of 10% during this period. In addition to

30

the Hg that was released to aquatic environments in China’s provinces, we estimated

31

the Hg release amounts and intensities (in g/km2·yr) for China’s 58 secondary river

32

basins. The highest aquatic Hg release intensities in China were associated with

33

industrial wastewater on the North China Plain and domestic sewage in eastern China

34

and southern China. We found that the overall uncertainty of our inventory ranges

35

from -22% to 32%. We suggest that the inventory provided by this study can help

36

establish a more accurate map of regional and global Hg cycling; it also has

37

implications for water quality management in China.

2

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

Environmental Science & Technology

38 39

40

Abstract art

41

3

ACS Paragon Plus Environment

Environmental Science & Technology

42

INTRODUCTION

43

Mercury (Hg) can be released from both anthropogenic and natural sources1

44

Anthropogenic mercury can alter the natural global Hg cycle, and because mercury is

45

transported through the atmosphere and oceans, it can cause worldwide

46

contamination.2-4 Most previous studies that examined Hg inventories have focused

47

on atmospheric anthropogenic emissions,5 which amount to 2000-3200 tons/yr

48

globally.6, 7 Hg release to aquatic environments is also important because this Hg can

49

be converted into a bioavailable and more toxic form (MeHg), which increases the

50

accumulation and biomagnification of Hg in aquatic food webs.5, 8 However, little

51

research has examined Hg release to aquatic environments, including both

52

anthropogenic and natural sources.5 AMAP/UNEP first produced a preliminary

53

inventory of Hg releases from anthropogenic sources to aquatic environments; this

54

was estimated to be 1100 tons globally in 2013.5 Mason et al.3 found that the current

55

annual global Hg discharge from rivers to oceans was 5 times higher than during

56

preindustrial times, reaching 1000 to 5500 tons/yr.4,

57

recognized riverine Hg discharge to oceans as an important pathway for global Hg

58

cycling,9-11 but detailed information on Hg sources and discharge to aquatic

59

environments is not yet available. Therefore, a detailed estimate of Hg released from

60

primary anthropogenic sources to the aquatic environment is important for

61

understanding global Hg cycling.

9, 10

Other researchers have

62

Previous studies have indicated that China is the largest contributor of atmospheric

63

mercury in the world.5, 12 Streets et al.13 developed the first complete anthropogenic 4

ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32

Environmental Science & Technology

64

Hg emission inventory for China, estimating its atmospheric Hg emissions from

65

anthropogenic sources to be 536 tons in 1999. Similarly, Pacyna et al.12 found that

66

atmospheric Hg emissions from primary anthropogenic sources in China were 605

67

tons in 2000 and 825 tons in 2005. However, to our knowledge, no inventory of

68

aquatic Hg release from primary anthropogenic sources in China exists. The

69

AMAP/UNEP inventory indicated that Hg release to aquatic environments from

70

primary anthropogenic sources in East and Southeast Asia reached 520 tons/yr (49%

71

of the total global release) and that atmospheric Hg emissions from anthropogenic

72

sources in this region were 780 tons/yr (40% of the total global emissions).5 Another

73

study7 reported that the ratio of commercial liquid Hg (from the intentional use of Hg

74

in products and processes) released to water and air was approximately 1:2 in 2010.

75

Therefore, primary anthropogenic Hg released to aquatic environments could be an

76

important pathway of Hg release in China.

77

Non-ferrous metal production, vinyl chloride monomer production, artisanal and

78

small-scale gold mining (ASGM), and oil refining were the major sectors considered

79

in the AMAP/UNEP inventory.5 Unlike atmospheric Hg emissions, there are few

80

records of Hg release to aquatic environments from anthropogenic sources. The

81

AMAP/UNEP inventory ignored several important sectors that largely release Hg to

82

aquatic environments rather than air.5,

83

estimates did not consider the coal washing industry, but this industry releases

84

significant amounts of Hg to water.15 This is especially true in China, which is the

85

largest coal producer in the world, accounting for 45% of the world’s production in

14

For example, their aquatic Hg emission

5

ACS Paragon Plus Environment

Environmental Science & Technology

86

2011.16 Domestic sewage discharge from treatment plants can also be an important

87

source of primary anthropogenic Hg to aquatic environments;17-19 again, this sector

88

was not considered in the AMAP/UNEP inventory.5

89

The purpose of this study was to develop the first complete inventory of Hg release

90

to aquatic environments from primary anthropogenic sources in China, including the

91

temporal distribution from 2001 to 2012 and the spatial distribution. In this inventory,

92

primary anthropogenic Hg release was divided into industrial wastewater and

93

domestic sewage, and we included all major sectors that discharge Hg in wastewater.

94

Hg release to water was further divided into particulate and dissolved phases. The

95

uncertainty associated with different sectors was calculated. The inventory produced

96

in this study will help establish a more accurate map of regional and global Hg

97

cycling and improve water quality management in China in the future.

98

METHODS

99

Model Description. The model developed in this study to estimate the primary

100

anthropogenic aquatic Hg released in China is called the China Aquatic Mercury

101

Release (CAMR) model. Parameterization of the uncertainty function was

102

accomplished with Monte Carlo simulations, following the method of Zhang et al.20

103

Aquatic Hg released from various sectors was divided into two groups (Table S1,

104

Supporting Information). The first group (Group 1) was chosen based on measured

105

Hg concentration data and the volume of wastewater released from various sectors.5, 9

106

The volume of wastewater discharge for each sector in Group 1 was end-of-pipe

6

ACS Paragon Plus Environment

Page 6 of 32

Page 7 of 32

Environmental Science & Technology

107

discharge data, which we collected from the China Environmental Statistics

108

Yearbook.21 Details regarding the concentration data are shown in Table S2 of the

109

Supporting Information. In contrast to the method of estimating emission factors for

110

atmospheric Hg emissions, aquatic Hg release from various sectors is estimated based

111

on the volume of wastewater discharge multiplied by the Hg concentration. The Hg

112

concentration data used in this study are end-of-pipe data following the recycling and

113

treatment of wastewater. For those sectors that lack measured Hg concentrations

114

(Group 2), we used a second estimation method. For Group 2, Hg release was

115

estimated using the method developed by AMAP/UNEP,5 in which UNEP Toolkit

116

distribution factors (the ratio of Hg released to air, water or soil from a given source)

117

were used to estimate the Hg release to aquatic environments based on the air

118

emission methods developed in previous studies.7 The CAMR model was established

119

based on the probabilistic distribution and includes these two groups as follows:

120 121

where

is the probabilistic distribution of the aquatic Hg released;

122

the probabilistic distribution of the Hg concentration in industrial wastewater and

123

domestic sewage (Tables S2 and S3, Supporting Information);

124

industrial wastewater or domestic sewage discharged (Table S4, Supporting

125

Information);

126

of Hg emissions from a given industrial sector (Table S2, Supporting Information),

127

based on the inventory of Zhang et al.20;

is the unit conversion factor;

is

is the volume of

is the probabilistic distribution

and

are the distribution factors of

7

ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 32

128

Hg released to water and air, respectively, for a given industrial sector (Table S5,

129

Supporting Information);

is the province;

is the release sector of Group 1; and

130

is the release sector of Group 2. The ratio of dissolved to particulate aquatic Hg

131

released from anthropogenic sources was determined based on observational data

132

from previous studies; it was set to be 1:4.22-25

133

This inventory includes 20 aquatic Hg release sectors, which are categorized as

134

industrial wastewater discharge or domestic sewage discharge. The industrial

135

wastewater sectors include the energy industry, nonferrous metal smelting, iron and

136

steel production, intentional use, and other low-release industries (Table S1,

137

Supporting Information). We estimated the aquatic Hg released from the coal washing

138

industry (part of the energy industry) based on the method used for Group 1. Given

139

the lack of previous research and distribution factors for the aquatic Hg released from

140

coal washing,5 the Hg concentration of wastewater from coal washing was estimated

141

based on the Hg removal efficiency of coal washing and the coal washing rate, as

142

follows:13, 26 (2)

143

144

where

is the probabilistic distribution of the Hg concentration in coal washing

145

wastewater;

146

determined in previous study;20

147

the China Energy Statistical Yearbook;27

148

washing;

is the probabilistic distribution of the Hg content of raw coal, is the activity level of raw coal production from is the Hg removal efficiency of coal

is the coal washing rate in China;

is the unit conversion factor;

8

ACS Paragon Plus Environment

Page 9 of 32

Environmental Science & Technology

149

is the water consumption volume of coal washing from the China Environmental

150

Statistics Report;28 and

151

ranged from 10% to 80%,29-33 and most estimates were between 25% and 60%. Based

152

on previous research,13 30% was selected as the average and 25% to 60% as the range

153

of Hg removal efficiency in this study. The coal washing rate in China was 16.6% in

154

2002 and 50.9% in 2010.26, 34 Because the wastewater from coal washing is always

155

directly discharged,35 treatment rates were not considered. The reuse rate of

156

wastewater from the coal washing industry was considered based on the volumes of

157

water consumption and wastewater discharge for the industry.

is the province. The Hg removal efficiency of coal washing

158

We developed a new method of estimating the aquatic Hg release from domestic

159

sewage discharge. This sector has rarely been considered in previous inventories.5

160

Based on the development level of different regions, domestic sewage discharge was

161

divided into treated versus non-treated discharge. The treatment rates of domestic

162

sewage in different regions from 2001 to 2012 are shown in the Supporting

163

Information (Tables S6 and S7) and come from the China Environmental Statistics

164

Yearbook.21 For domestic sewage, the Hg removal efficiencies of sewage treatment

165

plants were estimated based on the measured Hg data for different regions (Table S3,

166

Supporting Information). Therefore, the Hg concentration of domestic sewage

167

discharge is calculated as follows: (3)

168

169

where

is the probabilistic distribution of the Hg concentration in domestic

9

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 32

170

sewage;

is the probabilistic distribution of the Hg concentration in domestic

171

sewage discharged after treatment (Table S3, Supporting Information);

172

probabilistic distribution of the Hg concentration in domestic sewage discharged

173

before treatment (Table S3, Supporting Information);

174

domestic sewage; and

is the

is the treatment rate of

is the province.

175

Basic River Basin Unit. To understand the contributions of Hg sources at the river

176

basin level and to assist in river basin water quality management, the aquatic Hg

177

release from anthropogenic sources in China was recalculated at the scale of

178

individual river basins.36, 37 Mainland China was divided into 58 basins (Figure S1,

179

Supporting Information) based on the Industry Standard of China; this number

180

includes all of the secondary rivers in China.38 The geographic information layers for

181

China’s administrative regions and river basin system were obtained from the

182

National Geomatics Center of China (http://sms.webmap.cn/). The map of China’s

183

river basins was created using ArcGIS 10.3 software.

184

Mass Flow of Aquatic Hg. The mass flow of aquatic Hg to the ocean from primary

185

anthropogenic sources in China was calculated. To accomplish this, we simplified the

186

processes of Hg re-release from Hg-contaminated places, including erosion of

187

Hg-contaminated land surfaces, river bank erosion, and resuspension of legacy

188

anthropogenic Hg sources.5 This was necessary because, compared to atmospheric Hg

189

emissions, the Hg released to aquatic environments has multiple sources and transport

190

routes, making it more complex and difficult to quantify.5 Water resource data were

191

obtained from the China Water Resources Bulletin39 and the China Environmental 10

ACS Paragon Plus Environment

Page 11 of 32

Environmental Science & Technology

192

Statistics Yearbook.21 The China Water Resources Bulletin shows that, in 2012, 65%

193

of China’s surface water resources flowed into seas adjacent to China through rivers,

194

19% flowed to neighboring countries through rivers, 13% was stored in lakes and

195

reservoirs, and 3.8% flowed into boundary rivers.

196

To determine the contribution from aquatic Hg released from primary anthropogenic

197

sources in China, we estimated the Hg discharged from rivers into seas adjacent to

198

China in 2012. The total riverine Hg discharge into the seas was obtained from our

199

previous study,40 and the estimation method is similar to that of Sunderland and

200

Mason,4 who estimated the Hg discharge into oceans from rivers in the particulate

201

phase (Hg(P)) and the dissolved phase (Hg(D)).

202

of the Hg(P) discharged from a river is deposited on the river mouth and continental

203

shelf, and only 28% reaches the open ocean.4,

204

contribution of anthropogenic sources, we also calculated the natural sources of

205

riverine Hg discharge. We estimated the Hg contributions from natural sources in

206

rivers based on published discharge-Hg concentration relationships for Arctic rivers as

207

previous studies,9, 44 due to the lack of observational data on the background of Hg

208

concentrations in rivers. Rivers in the Arctic region are rarely affected by human

209

activity, therefore, they should reflect natural sources. The contributions of legacy

210

anthropogenic sources were obtained by subtracting the anthropogenic sources and

211

natural sources from total riverine Hg discharge into Chinese adjacent seas.

According to previous studies, 72%

42, 43

To estimate the relative

212

Uncertainty Analysis. A Monte Carlo simulation (10,000 runs) was used to

213

characterize the uncertainty of the inventory. The uncertainty originated from the 11

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 32

214

probabilistic distributions of Hg concentrations in Group 1, Hg emissions in Group 2,

215

discharge volumes, and other parameters. The probabilistic distributions of Hg

216

concentrations in Group 1 were assumed to follow a log-normal distribution based on

217

measured data of previous studies.20,

218

emissions found in previous studies were used; these had skewed distributions (e.g.,

219

log-normal and Weibull distributions).20 For the discharge volumes and activity levels,

220

uniform distributions with a fixed coefficient of deviation (5%) were assumed5, 48

221

because the data were derived from official statistics.28,

222

presented by the UNEP were used to estimate the Hg release from wastewater in

223

Group 2,5 because Hg releases to aquatic environments from anthropogenic point

224

sources are very poorly documented. We note that using these distribution factors

225

might also involve the uncertainties associated with atmospheric emission estimates

226

and increase the estimation uncertainties in the inventories of AMAP/UNEP and this

227

study.5 Mean values and P20-80 confidence intervals of the statistical distribution of

228

aquatic Hg release were calculated to quantify the aquatic Hg release and characterize

229

the uncertainty.29

230

RESULTS AND DISCUSSION

45-47

In Group 2, the uncertainties of Hg

49

Distribution factors

231

Aquatic Hg Release from Primary Anthropogenic Sources and Trends. Figure

232

1 summarizes the total aquatic Hg released by sectors associated with industrial

233

wastewater and domestic sewage in China. The total aquatic Hg released from

234

anthropogenic sources in China was 98 tons in 2012, with 78 tons of Hg(P) and 20

235

tons of Hg(D). From 2001 to 2012, the maximum Hg release was 123 tons in 2001 12

ACS Paragon Plus Environment

Page 13 of 32

Environmental Science & Technology

236

and the minimum was 87 tons in 2010. The total aquatic Hg release decreased

237

discontinuously during this period but increased from 2002 to 2003 and from 2010 to

238

2011. We note that the primary anthropogenic Hg emissions to the atmosphere

239

increased over this time period;20 this is not inconsistent with the decrease in aquatic

240

Hg releases that was observed in this study because the Chinese government has

241

regulated wastewater discharge more stringently in recent years.21 Although air

242

pollution control measures have reduced atmospheric Hg emissions in China, primary

243

anthropogenic Hg emissions to the atmosphere still grew by an average of 4.2% per

244

year from 2000 to 2010.20

245

Because Hg releases to aquatic environments from primary anthropogenic sources

246

are poorly documented,5 the total national aquatic Hg release was compared with the

247

aquatic Hg release inventory developed by AMAP/UNEP and atmospheric Hg

248

emissions estimated by previous studies.5, 12, 13, 20 In this study, we found that the ratio

249

of anthropogenic Hg released to water and air in China is approximately 1:5,12, 13, 20

250

which is lower than the global ratio of 1:2 estimated by previous studies.5,

251

suggest that this result is inconsistent because (1) the specific industrial structure in

252

China gives rise to greater atmospheric Hg emissions, including significant

253

contributions from coal-fired power generation and nonferrous metal smelting;13, 20

254

and (2) the aquatic Hg release intensities of various sectors in China estimated in this

255

study differ from the estimates made by AMAP/UNEP. AMAP/UNEP estimated that

256

the aquatic Hg released from artisanal and small-scale gold mining (ASGM)

257

accounted for 83% of the total global aquatic Hg release and 87% of the release in 13

ACS Paragon Plus Environment

7

We

Environmental Science & Technology

258

East and Southeast Asia. However, following the implementation of strict regulations

259

in China, ASGM has decreased sharply in recent years to only 1−3% of the total gold

260

output in 2010.20, 50 As a primary pathway in the AMAP/UNEP inventory, aquatic Hg

261

released from ASGM in China was only 2.7 tons in 2012 (2.8% of the total Hg

262

released), which represents a decrease from 10 tons in 2001 (8.1%) at an average rate

263

of 11% per year.

14

ACS Paragon Plus Environment

Page 14 of 32

Page 15 of 32

Environmental Science & Technology

264 265

Figure 1. Release of mercury to the aquatic environment from various sectors in China from 2001 to 2012. 15

ACS Paragon Plus Environment

Environmental Science & Technology

266

From 2001 to 2011, there was a variable amount of aquatic Hg released from other

267

important sources, such as coal-fired power plants, nonferrous metal smelting, mining

268

and washing of coal, and domestic sewage (Figure 2 (A)). As the leading source of

269

anthropogenic Hg emissions in China (releasing 68 to 110 tons/yr),13, 20, 29 the release

270

from coal-fired power plants decreased discontinuously. Due to the gradual decrease

271

in wastewater discharged from coal-fired power plants,28 aquatic Hg release from this

272

sector accounted for 31% of the overall release in 2001 and 17% in 2012, with an

273

annual average decrease of 5.9%. Because these data were unavailable or incomplete,

274

AMAP/UNEP did not consider the aquatic Hg release from coal-fired power plants.5

275

However, there is now abundant measured data for aquatic Hg release from this sector

276

in China.45 Aquatic Hg release from domestic sewage represents another sector that

277

might have been ignored in the past; this also showed a gradual decrease (Figure 2

278

(B)), accounting for 22% of the overall release in 2001 and 17% in 2012, with an

279

annual average decrease of 3.0%. Because of the increase in the number of sewage

280

treatment plants and the improved treatment rates in China (19% in 2001 and 83% in

281

2012),28 the aquatic Hg release from treated domestic sewage increased gradually.

282

The aquatic Hg released from nonferrous metal smelting (including Hg, Al, Zn, Pb,

283

Cu, large-scale gold mining, and ASGM) also decreased gradually after 2003 at an

284

average rate of 0.6%. In 2012, the aquatic Hg released from nonferrous metal

285

smelting accounted for 33% of the total Hg release. Previous studies estimated that

286

Hg emissions to the atmosphere were 90 to 240 tons/yr from nonferrous metal

287

smelting in China,13, 20 and in this study, the release to aquatic environments was 16

ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32

Environmental Science & Technology

288

estimated to be 33.4 tons in 2012. Zn smelting was the largest nonferrous smelting

289

source of aquatic Hg release; this was estimated to be 12.3 tons in 2012 (accounting

290

for 13% of the total aquatic Hg release from anthropogenic sources in China), and it

291

decreased at an average rate of 2.5% over the study period. The Hg releases from Pb

292

smelting and Cu smelting were 6.9 tons (7.0%) and 0.7 tons (0.6%) in 2012, with

293

average growth rates of 5.5% and -0.65%, respectively.

294

Aquatic Hg release from coal mining and washing increased rapidly and accounted

295

for 1.9% of the total release in 2001 and 25% in 2012, with an annual average growth

296

rate of 25%. Aquatic Hg release from coal mining and washing is now the second

297

largest anthropogenic source in China. This can be attributed to the increase in coal

298

production18 and the wastewater discharged from this sector to the environment.29, 36

299

If there is lack of strict governmental supervision and wastewater treatment and reuse

300

rates are not improved, aquatic Hg release from coal mining and washing could soon

301

become the largest anthropogenic source of aquatic Hg in China.

17

ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 32

302 303

Figure 2. Variations in aquatic mercury release from the four major sectors in China from 2001 to 2012 from (A) coal-fired power plants, domestic sewage,

304

nonferrous metal smelting, and coal mining and washing; (B) aquatic Hg release from treated versus untreated discharge from domestic sewage.

18

ACS Paragon Plus Environment

Page 19 of 32

Environmental Science & Technology

305

In addition to the four major aquatic Hg release sectors mentioned above, the

306

manufacture of wearable textiles and apparel and the intentional use of Hg (including

307

vinyl chloride monomer production, reagent production, thermometer production,

308

fluorescent lamp production, and battery production) contributed 4.2% (4.1 tons) and

309

3.3% (3.2 tons), respectively, to the overall release in 2012. Aquatic Hg released from

310

the intentional use of Hg decreased at an annual average rate of 14% due to the

311

decreased use of Hg in China.20 However, aquatic Hg release from the manufacture of

312

wearable textiles and apparel increased at an annual average rate of 20% due to the

313

growth in wastewater discharge from this sector. Other small sources of aquatic Hg

314

release, including oil refining, iron and steel production, the printing industry, and

315

waste resource utilization, together accounted for 1.0 tons (1.0%) of Hg in 2012, and

316

had an annual average growth rate of 6.1% from 2001 to 2012.

317

Spatial Distribution of Aquatic Hg Release from Primary Anthropogenic

318

Sources in China. Figure 3 shows the distribution of aquatic Hg released from

319

various primary anthropogenic sources by province in China in 2012. The five

320

provinces with the highest aquatic Hg releases were Henan (8.8 tons in 2012 and

321

8.6% of the national total), Shandong (8.2 tons, 8.0%), Hunan (7.3 tons, 7.2%),

322

Shanxi (6.5 tons, 6.4%), and Hebei (5.5 tons, 5.4%). These values are similar to

323

previously reported provincial distributions of atmospheric Hg emissions.20 The high

324

aquatic Hg release from Henan could be attributed to abundant Pb smelting

325

(accounting for 29% of the total aquatic Hg released in the province) and coal mining

326

and washing (14%). Shandong’s numbers are likely due to coal mining and washing 19

ACS Paragon Plus Environment

Environmental Science & Technology

327

(accounting for 28%) and large-scale gold mining (21%). Zn smelting was the

328

primary Hg release source in Hunan (40%). In Shanxi and Hebei, coal mining and

329

washing contributed 73% and 67%, respectively.

330

We found that the major sources of aquatic Hg release varied by province (Figure

331

3). In Guangdong, domestic sewage (45%) was the major source of aquatic Hg

332

because Guangdong has the largest population and is one of the most developed

333

provinces in China (Figure S3). Hg mining (29%) was the major source in Guizhou.

334

The lowest aquatic Hg release occurred in Tibet (0.07% of the national total release),

335

where domestic sewage (79%) was the major source of Hg.

20

ACS Paragon Plus Environment

Page 20 of 32

Page 21 of 32

Environmental Science & Technology

336 337

Figure 3. Provincial distribution of mercury release to the aquatic environment in China in 2012.

21

ACS Paragon Plus Environment

Environmental Science & Technology

338

Figure 4 (A) and (C) show the distributions of aquatic Hg released from industrial

339

wastewater and domestic sewage into 58 Chinese secondary river basins in 2012.

340

There are significant differences in the aquatic Hg release from industrial wastewater

341

and domestic sewage among the basins. The Yellow River basin (ID: 18) in northern

342

China showed the largest Hg release from industrial wastewater, with a total of 9.0

343

tons, and Dongting Lake basin (ID: 31) in central China showed the largest release

344

from domestic sewage, with a total of 1.1 tons. In contrast, the Senggecangbu River

345

basin (ID: 50), located on the western Tibetan Plateau, released the smallest amount

346

of Hg from both industrial wastewater and domestic sewage, with 0.9 × 10-3 ton/yr

347

and 3.1 × 10-3 ton/yr, respectively.

348

Figure 4 (B) and (D) show the aquatic Hg release density (g/km2·yr) of industrial

349

wastewater and domestic sewage in 58 Chinese secondary basins in 2012. On the

350

whole, among all Chinese basins, the highest aquatic Hg release density of industrial

351

wastewater occurred on the North China Plain (IDs: 12, 13, 14, 15, 16, 19, 21, 22, 23,

352

24, and 34; Figure 4 (B)). Because aquatic Hg release from coal-fired power plants is

353

associated with economic level and population density, the release density in eastern

354

China (IDs: 16, 21, 23, 34, 35, and 36) was higher than in other regions (Figure S4

355

(F)). Aquatic Hg release from coal mining and washing and nonferrous metal smelting

356

is associated with the distribution of coal and non-ferrous metal resources. Therefore,

357

the highest release density from coal mining and washing occurred on the North

358

China Plain (IDs: 11, 12, 13, 14, 15, 16, 19, 21, and 23; Figure S4 (D)) and the

359

highest release densities from nonferrous metal smelting occurred in central China 22

ACS Paragon Plus Environment

Page 22 of 32

Page 23 of 32

Environmental Science & Technology

360

and the North China Plain (IDs: 16, 21, 22, 23, 30, and 31; Figure S4 (H)).

361

In general, aquatic Hg release from domestic sewage is affected by economic level

362

and population density. Thus, the aquatic Hg release densities of domestic sewage

363

were highest in eastern China (IDs: 23, 24, 34, 35, and 36) and southern China (IDs:

364

39, 41, 42, 43, and 45), the two most populous regions in China (Figure 4 (D)).

365

Mainland China can be separated into two parts (east and west) by the famous “Hu

366

Huanyong line” (Figure 4 (D)); this division highlights the impact of population

367

density on the distribution of aquatic Hg release. Aquatic Hg release from domestic

368

sewage in the eastern part of China is an order of magnitude higher than in the

369

western part. Similar results have been reported for other pollutants, such as

370

antibiotics, steroids, black carbon, and personal care products.36, 37, 51, 52

23

ACS Paragon Plus Environment

Environmental Science & Technology

371 372

Figure 4. Distribution of mercury release to aquatic environments from industrial wastewater

373

(A and B) and domestic sewage (C and D) in Chinese secondary river basins in 2012. (A) and

374

(C) show the total release in units of ton/yr; (B) and (D) show the release densities in units of

375

g/m2·yr. The black line in (D) is the “Hu Huanyong line” through China. Basin Ids are shown

376

in Table S1 in the Supporting Information.

377

Mass Flow of Aquatic Mercury to Rivers and the Ocean from Anthropogenic

378

Sources and Environmental Implications. Figure 5 shows the mass flow of aquatic

379

Hg from primary anthropogenic sources in China in 2012. In our previous study, we

380

estimated that the amount of Hg discharged from rivers into seas adjacent to China

381

was 157 tons in 2012.40 Of this total, 16% (27 tons) came from industrial wastewater 24

ACS Paragon Plus Environment

Page 24 of 32

Page 25 of 32

Environmental Science & Technology

382

and domestic sewage, which are direct anthropogenic sources. An additional 23% (36

383

tons) came from natural sources (such as leaching, runoff and erosion processes from

384

land), a total which is consistent with the results of previous studies that estimated a

385

natural Hg release of 20% to 46% globally.3-5, 53 Finally, 61% (94 tons) came from

386

legacy anthropogenic sources, including (1) erosion of Hg-contaminated land surfaces;

387

(2) resuspension of Hg-contaminated river bed sediments; and (3) land and water

388

management practices.5 Therefore, anthropogenic sources contributed up to 77% of

389

the riverine Hg discharge, which is similar to previously published values.3 Hg(D)

390

discharge from Chinese rivers into ocean margins accounted for 45% of the total Hg

391

discharge;40 this would result in more Hg transported to the open ocean.4 Because of

392

the limited information that is currently available, more studies involving extensive

393

sampling and monitoring are needed.

394

With its large population and rapid industrialization, China releases abundant Hg

395

from anthropogenic sources into the aquatic environment. Although aquatic Hg

396

releases from primary anthropogenic sources in China decreased from 2001 to 2012,

397

there was still a release of 98 tons in 2012. Aquatic Hg releases from anthropogenic

398

sources threatens not only local environments but also global marine ecosystems.

399

Therefore, it is essential to reduce the discharge of industrial wastewater in the future

400

by increasing the recycling rate and enhancing the use of environmentally-friendly

401

technologies in production and wastewater treatment processes.

25

ACS Paragon Plus Environment

Environmental Science & Technology

Page 26 of 32

402

403

404 405

Figure 5. Mass flow of aquatic mercury release from anthropogenic sources to seas adjacent to China in 2012. (A) Mass flow of aquatic primary

406

anthropogenic Hg; (B) Contributions of anthropogenic and natural sources to the riverine Hg discharge into the seas. Note: Hg deposition in the river channels

407

was estimated based on the distribution of Hg in the dissolved versus particulate phases near the river mouth areas in our previous study.40 In that study, we

26

ACS Paragon Plus Environment

Page 27 of 32

Environmental Science & Technology

408

found a 45%:55% ratio of Hg in the dissolved versus particulate phases in riverine Hg discharge into the seas, while previous studies have found a ratio of 1:4

409

in wastewater (see the Method section). We assumed that all the Hg(D) released from anthropogenic sources reached the adjacent seas, which were 12 tons

410

calculated using the ratio of 1:4. The ratio of Hg(P) released from anthropogenic sources and reached the adjacent seas was calculated by (1/4) × (55%/45%),

411

which was 31% (15 tons). Therefore, we assumed that 69% (33 tons) of Hg(P) released from anthropogenic sources deposited in river channels before

412

reaching the river mouths. 4.7 tons Hg (4.8% of wastewater) were discharged into the seas directly and 93 tons were released into aquatic environment of

413

mainland China, according to the China Marine Statistical Yearbook.54 The contribution of Hg remobilization in rivers is involved in estimates of the

414

contributions of legacy anthropogenic Hg release and natural Hg release.

27

ACS Paragon Plus Environment

Environmental Science & Technology

415

Uncertainties Associated with the Inventory. We estimated the uncertainties

416

associated with the aquatic Hg released from anthropogenic activities in China and in

417

each sector based on the probabilistic distribution of Hg concentrations and emission

418

factors (Figure 6). The overall uncertainty ranged from -22% to 32% (P20-80). The

419

smallest uncertainties were associated with coal-fired power plants (-35%, 52%) and

420

domestic sewage (-42%, 77%), while the largest uncertainty was associated with

421

ASGM (-80%, 400%). This can be attributed to the large uncertainty in production

422

data from this sector in China. In addition, when the statistical distribution of a

423

parameter in this sector was unknown, we assumed a uniform distribution. The

424

uncertainties associated with coal mining and washing (-58%, 260%) were second

425

only to ASGM; these uncertainties largely mainly from the probabilistic distribution

426

of Hg concentrations in the raw materials from each province.20

427 428

Figure 6. Uncertainties (%) associated with the inventory by sector. 28

ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32

Environmental Science & Technology

429

ASSOCIATED CONTENT

430

Supporting Information

431

The Supplemental Information includes the Hg concentrations of wastewater in Group

432

1, distribution factor in Group 2, activity levels of each sector, treatment rates of

433

domestic sewage in China, Hg concentrations collected at river mouths, and detailed

434

descriptions of some of the data used to estimate Hg release from anthropogenic

435

sources in China. The trend of detailed anthropogenic Hg release to the aquatic

436

environment in China from 2001 to 2012, provincial distribution of aquatic Hg release

437

from domestic sewage in China, Hg discharge from major Chinese rivers into seas

438

adjacent to China are also provided in the supporting information.

439

ACKONWLEDGEMENTS

440

This work was funded by National Natural Science Foundation of China (41571484,

441

41130535 and 41471403). The authors would like to thank the editor and reviewers

442

for their insight comments and suggestions on the manuscript.

443

REFERENCES

444 445 446 447 448 449 450 451 452 453 454

1.

Driscoll, C. T.; Mason, R. P.; Chan, H. M.; Jacob, D. J.; Pirrone, N., Mercury as a global pollutant:

sources, pathways, and effects. Environmental science & technology 2013, 47, (10), 4967-4983. 2.

Corbitt, E. S.; Jacob, D. J.; Holmes, C. D.; Streets, D. G.; Sunderland, E. M., Global source–receptor

relationships for mercury deposition under present-day and 2050 emissions scenarios. Environmental science & technology 2011, 45, (24), 10477-10484. 3.

Mason, R. P.; Sheu, G. R., Role of the ocean in the global mercury cycle. Global biogeochemical

cycles 2002, 16, (4), 40-1-40-14. 4.

Sunderland, E. M.; Mason, R. P., Human impacts on open ocean mercury concentrations. Global

Biogeochemical Cycles 2007, 21, (4), 177-180. 5.

AMAP/UNEP, Technical Background Report for the Global Mercury Assessment 2013. In Arctic

Monitoring and Assessment Programme, Oslo, Norway/UNEP Chemicals Branch Geneva, Switzerland: 29

ACS Paragon Plus Environment

Environmental Science & Technology

455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498

2013; pp vi-263. 6.

Streets, D. G.; Devane, M. K.; Lu, Z.; Bond, T. C.; Sunderland, E. M.; Jacob, D. J., All-Time Releases

of Mercury to the Atmosphere from Human Activities. Environmental Science & Technology 2011, 45, (24), 10485-10491. 7.

Horowitz, H. M.; Jacob, D. J.; Amos, H. M.; Streets, D. G.; Sunderland, E. M., Historical mercury

releases from commercial products: Global environmental implications. Environmental science & technology 2014, 48, (17), 10242-10250. 8.

Mahaffey, K. R.; Sunderland, E. M.; Chan, H. M.; Choi, A. L.; Grandjean, P.; Mariën, K.; Oken, E.;

Sakamoto, M.; Schoeny, R.; Weihe, P., Balancing the benefits of n-3 polyunsaturated fatty acids and the risks of methylmercury exposure from fish consumption. Nutrition reviews 2011, 69, (9), 493-508. 9.

Amos, H. M.; Jacob, D. J.; Kocman, D.; Horowitz, H. M.; Zhang, Y.; Dutkiewicz, S.; Horvat, M.;

Corbitt, E. S.; Krabbenhoft, D. P.; Sunderland, E. M., Global biogeochemical implications of mercury discharges from rivers and sediment burial. Environmental science & technology 2014, 48, (16), 9514-9522. 10. Cossa, D.; Coquery, M.; Gobeil, C.; Martin, J.-M., Mercury fluxes at the ocean margins. In Global and regional mercury cycles: sources, fluxes and mass balances, Springer: 1996; pp 229-247. 11. Fisher, J. A.; Jacob, D. J.; Soerensen, A. L.; Amos, H. M.; Steffen, A.; Sunderland, E. M., Riverine source of Arctic Ocean mercury inferred from atmospheric observations. Nature Geoscience 2012, 5, (7), 499-504. 12. Pacyna, E. G.; Pacyna, J.; Sundseth, K.; Munthe, J.; Kindbom, K.; Wilson, S.; Steenhuisen, F.; Maxson, P., Global emission of mercury to the atmosphere from anthropogenic sources in 2005 and projections to 2020. Atmospheric Environment 2010, 44, (20), 2487-2499. 13. Streets, D. G.; Hao, J.; Wu, Y.; Jiang, J.; Chan, M.; Tian, H.; Feng, X., Anthropogenic mercury emissions in China. Atmospheric Environment 2005, 39, (40), 7789-7806. 14. UNEP, Toolkit for Identification and Quantification of Mercury Releases. In United Nations Environment Programme (UNEP) Chemicals Geneva, Switzerland: 2005. 15. Tang, Y.; Chang, C.; Zhang, Y.; Li, W., Migration and distribution of fifteen toxic trace elements during the coal washing of the Kailuan Coalfield, Hebei Province, China. Energy, Exploration & Exploitation 2009, 27, (2), 143-152. 16. IEA, Coal Information 2012. International Energy Agency (IEA): Paris, 2012. 17. Gbondo-Tugbawa, S. S.; McAlear, J. A.; Driscoll, C. T.; Sharpe, C. W., Total and methyl mercury transformations and mass loadings within a wastewater treatment plant and the impact of the effluent discharge to an alkaline hypereutrophic lake. Water research 2010, 44, (9), 2863-2875. 18. Kang, D.; Luo, L.; Wang, X. C.; Du, J. In Investigation of Cu, Zn, Mn, Cr, Hg and As migration in domestic wastewater treatment plant, Mechanic Automation and Control Engineering (MACE), 2010 International Conference on, 2010; IEEE: 2010; pp 2086-2089. 19. Zhao, J. L.; Jiang, Y. X.; Yan, B.; Wei, C.; Zhang, L. J.; Ying, G. G., Multispecies acute toxicity evaluation of wastewaters from different treatment stages in a coking wastewater‐treatment plant. Environmental toxicology and chemistry 2014, 33, (9), 1967-1975. 20. Zhang, L.; Wang, S.; Wang, L.; Wu, Y.; Duan, L.; Wu, Q.; Wang, F.; Yang, M.; Yang, H.; Hao, J., Updated Emission Inventories for Speciated Atmospheric Mercury from Anthropogenic Sources in China. Environmental science & technology 2015, 49, (5), 3185-3194. 21. NSB, China Environmental Statistics Yearbook. National Statistical Bureau of China (NBS): Beijing, China, 2002-2013. 30

ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32

Environmental Science & Technology

499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542

22. Zhang, G.; Liu, C.-Q.; Wu, P.; Yang, Y., The geochemical characteristics of mine-waste calcines and runoff from the Wanshan mercury mine, Guizhou, China. Applied geochemistry 2004, 19, (11), 1735-1744. 23. Li, P.; Feng, X.; Qiu, G.; Shang, L.; Wang, S., Mercury pollution in Wuchuan mercury mining area, Guizhou, Southwestern China: The impacts from large scale and artisanal mercury mining. Environment international 2012, 42, 59-66. 24. Li, P.; Feng, X.; Chou, G., The Preliminary Study of Mercury Pollution in Wuchuan Mercury Area Guizhou Province. Environmental Chemistry 2008, 27, (1), 96-99 (in Chinese). 25. Li, H.; Mao, Y.; Li, Y.; Ma, Y.; Wang, X.; Song, D., Occurrence and Mass Balance of Mercury at a Sewage Treatment Plant. Part 1: Total Mercury. Environmental Chemistry 2014, 33, (7), 1059-1065 (in Chinese). 26. Jiang, J.; Ming, J.; Wu, Y.; Streets, D. G.; Duan, L.; Tian, H., Development of Mercury Emission Inventory from Coal Combustion in China. Environmental Sciences 2005, 26, (2), 34-39 (in Chinese). 27. NESA, China Energy Statistical Yearbook. National Energy Statistical Agency of China (NESA): Beijing, China, 2002-2013. 28. MEP China Environmental Statistics Report; China Ministry of Environmental protection (MEP): Beijing, China, 2001-2012. 29. Wu, Y.; Streets, D.; Wang, S.; Hao, J., Uncertainties in estimating mercury emissions from coal-fired power plants in China. Atmospheric Chemistry and Physics 2010, 10, (6), 2937-2946. 30. Wang, L.; Peng, S., Effectiveness of coal preparation in removing hazardous trace elements in air. Environmental Science and Technology 2003, 26, (1), 5-7 (in Chinese). 31. Feng, X.; Hong, Y.; Hong, B.; Ni, J.; Zhu, Y., Modes of occurrence of mercury in coal. Bulletin of Mineralogy, Petrology and Geochemistry 2001, 20, (2), 71-78 (in Chinese). 32. Luttrell, G. H.; Kohmuench, J. N.; Yoon, R.-H., An evaluation of coal preparation technologies for controlling trace element emissions. Fuel Processing Technology 2000, 65, 407-422. 33. Akers, D., Coal cleaning controls HAP emissions. Power Engineering (Barrington) 1996, 100, (6), 33-36. 34. Gu, Y.; Yu, W.; Ma, X., Comparative Study on Current Situation of Coal Preparation and Policy at Home and Abroad. Coal Preparation Technology 2012, 4, 110-116 (in Chinese). 35. Li, Y. Study on Treatment and Reuse of Technology of Concentrated Wastewater of Coal Washing. Northeastern University, 2006 (in Chinese). 36. Zhang, Q. Q.; Zhao, J.-L.; Ying, G.-G.; Liu, Y.-S.; Pan, C.-G., Emission estimation and multimedia fate modeling of seven steroids at the river basin scale in China. Environmental science & technology 2014, 48, (14), 7982-7992. 37. Zhang, Q.-Q.; Ying, G.-G.; Pan, C.-G.; Liu, Y.-S.; Zhao, J.-L., A comprehensive evaluation of antibiotics emission and fate in the river basins of China: Source analysis, multimedia modelling, and linkage to bacterial resistance. Environmental science & technology 2015, 49, (11), 6772-6728. 38. MWR, Code for China River Name. In Ministry of Water Resources of China (MWR): Beijing,China, 2000. 39. MWR, China Water Resources Bulletin. Ministry of Water Resources of China (WMR): Beijing: China, 2012. 40. Liu, M.; Chen, L.; Wang, X.; Zhang, W.; Tong, Y.; Ou, L.; Xie, H.; Shen, H.; Ye, X.; Deng, C.; Wang, H., Mercury Export from Mainland China to Adjacent Seas and Its Influence on the Marine Mercury Balance. Environmental Science & Technology 2016, 50, (12), 6224-6232. 31

ACS Paragon Plus Environment

Environmental Science & Technology

543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574

41. MWR, China river sediment Bulletin. Ministry of water resources of China (MWR): Beijing,China, 2012. 42. Chester, R.; Jickells, T., The transport of material to the oceans: Relative flux magnitudes. Marine Geochemistry 2003, 92-124. 43. Cossa, D.; Martin, J.-M.; Takayanagi, K.; Sanjuan, J., The distribution and cycling of mercury species in the western Mediterranean. Deep Sea Research Part II: Topical Studies in Oceanography 1997, 44, (3), 721-740. 44. Emmerton, C. A.; Graydon, J. A.; Gareis, J. A.; St. Louis, V. L.; Lesack, L. F.; Banack, J. K.; Hicks, F.; Nafziger, J., Mercury export to the Arctic Ocean from the Mackenzie River, Canada. Environmental science & technology 2013, 47, (14), 7644-7654. 45. Tang, Q.; Liu, G.; Yan, Z.; Sun, R., Distribution and fate of environmentally sensitive elements (arsenic, mercury, stibium and selenium) in coal-fired power plants at Huainan, Anhui, China. Fuel 2012, 95, 334-339. 46. Chen, C.; Wang, H.; Zhang, W.; Hu, D.; Chen, L.; Wang, X., High‐resolution inventory of mercury emissions from biomass burning in China for 2000–2010 and a projection for 2020. Journal of Geophysical Research: Atmospheres 2013, 118, (21), 12,248-12,256. 47. Wang, Y.; Sun, F.; Lu, J.; Yuan, H.; Liu, X.; Shi, Y.; Jin, X.; Zhang, H., Monitoring and Evaluation of Heavy Metals Pollution in the Sewage Pool of Karamay Petrochemical Plant. Arid Environmental Monitoring 2010, 24, (3), 133-137 (in Chinese). 48. Shen, H.; Huang, Y.; Wang, R.; Zhu, D.; Li, W.; Shen, G.; Wang, B.; Zhang, Y.; Chen, Y.; Lu, Y., Global atmospheric emissions of polycyclic aromatic hydrocarbons from 1960 to 2008 and future predictions. Environmental science & technology 2013, 47, (12), 6415-6424. 49. NBS, China Statistical Yearbook. National Bureau of Statistics (NBS): Beijing, China, 2013. 50. CGA, China Gold Yearbook. China Gold Association (CGA): Beijing, China, 2002-2012. 51. Wang, R.; Tao, S.; Wang, W.; Liu, J.; Shen, H.; Shen, G.; Wang, B.; Liu, X.; Li, W.; Huang, Y., Black carbon emissions in China from 1949 to 2050. Environmental science & technology 2012, 46, (14), 7595-7603. 52. Hodges, J.; Holmes, C.; Vamshi, R.; Mao, D.; Price, O., Estimating chemical emissions from home and personal care products in China. Environmental Pollution 2012, 165, 199-207. 53. Zhang, Y.; Jaeglé, L.; Thompson, L.; Streets, D. G., Six centuries of changing oceanic mercury. Global Biogeochemical Cycles 2014, 28, (11), 1251-1261. 54. NBS, China Marine Statistical Yearbook. National Bureau of Statistics (NBS): Beijing, China, 2012.

575

32

ACS Paragon Plus Environment

Page 32 of 32