Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE
Article
Atmospheric arsenic deposition in the Pearl River Delta region, South China: influencing factors and speciation Minjuan HUANG, Haoran Sun, Hongtao Liu, Xuemei Wang, Baomin Wang, and Dan Zheng Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b04427 • Publication Date (Web): 15 Feb 2018 Downloaded from http://pubs.acs.org on February 16, 2018
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.
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 34
Environmental Science & Technology
1
Atmospheric arsenic deposition in the Pearl River Delta region, South
2
China: influencing factors and speciation
3
Minjuan Huang
4
1,2
5
1
6
China
7
2
8
Studies, Sun Yat-sen University, Guangzhou, 510275, P.R. China
9
3
1,2*
, Haoran Sun 3,5, Hongtao Liu 4, Xuemei Wang
1,6*
, Baomin Wang
, Dan Zheng 1,2
School of Atmospheric Sciences, Sun Yat-sen University, Guangzhou, 510275, P.R.
Guangdong Province Key Laboratory for Climate Change and Natural Disaster
School of Environmental Science and Engineering, Sun Yat-sen University,
10
Guangzhou 510275, P.R. China
11
4
12
5
13
Remediation Technology, Sun Yat-sen University, Guangzhou 510275, P.R. China
14
6
15
510632, P.R. China
16
*Corresponding to: School of Atmospheric Sciences, Sun Yat-sen University,
17
Guangzhou, 510275, PR China. Tel: 020-84112493. Email: Mnjuan Huang:
18
[email protected]; Xuemei Wang:
[email protected] Instrumental Analysis & Research Center, Guangzhou 510275, P.R. China Guangdong Province Key Laboratory of Environmental Pollution Control and
Institute for Environmental and Climate Research, Jinan University, Guangzhou,
1
ACS Paragon Plus Environment
Environmental Science & Technology
19
Abstract: This is a comprehensive study on mobilization/speciation and
20
temporal/spatial variation of atmospheric arsenic (As) deposition in the Pearl River
21
Delta (PRD) region. A set of experimental procedure was established for measuring
22
the deposition fluxes of individual As species. The deposition carrying inorganic AsIII %
23
was significantly higher than that contained in atmospheric particles. Compared with
24
dry deposition, wet deposition was much more harmful to the regional ecosystem, as
25
it contributed the majority of bulk deposition (>75%), and carried most of the
26
mobilized iAsIII compounds. Stepwise linear regression model was utilized to identify
27
the factors influencing total As deposition (wet: precipitation and PM2.5, dry: relative
28
humidity, wind speed and PM10, bulk: precipitation, PM2.5 and wind speed). By
29
examining the representativeness of the study sites and comparison with the literature
30
data, the statistic models were verified to explain the temporal/spatial variation of
31
total As deposition in the entire PRD region, where significant seasonal variation was
32
only found for wet deposition (wet season > dry season). The annual As load
33
contributed from regional atmospheric deposition was increasing from 2013 to 2015,
34
when the contributions of individual cities varied annually.
35 36
Keywords:
37
Atmospheric arsenic deposition, mobilization/speciation, temporal/spatial variation,
38
PM2.5/PM10, meteorological factors
2
ACS Paragon Plus Environment
Page 2 of 34
Page 3 of 34
Environmental Science & Technology
39
1 INTRODUCTION
40
Arsenic (As) is a widely distributed toxic element in the global environment. Its
41
toxicities vary with different species. In general, the inorganic As species are more toxic
42
than most of the organic ones (except for the trivalent organoarsenicals), and inorganic
43
trivalent arsenicals (iAsIII) are more toxic than inorganic pentavalent arsenicals (iAsV)
44
Arsenic and its inorganic compounds have been classified as Group 1 human carcinogens
45
by the International Agency for Research on Cancer (IARC) 4 and dangerous substances for
46
the environment by the European Union (EU) 5. Recently, the China Food and Drug
47
Administration (CFDA) also re-organized the data from IARC and issued the
48
carcinogenicity of As 6. Furthermore, arsenic could also cause many other non-carcinogenic
49
human health effects (e.g., cardiovascular disease, neurological disorders, gastrointestinal
50
disturbances, liver disease and renal disease, reproductive health effects, dermal changes…)
51
7
1-3
.
.
52
East Asia (anthropogenic sources: 15.5 Gg/year, natural sources: 0.3 Gg/year) is the
53
largest As emission contributor around the world, and the atmospheric As concentration is
54
extremely high in eastern China 8. Up to 56.3% of atmospheric concentration (5.9-53 ng/m3)
55
and 58.0% of total deposition in this section are attributed to the anthropogenic Asian
56
emissions 8. However, the modeled deposition flux results were not reported in the same
57
study due to lack of deposition measurement data for their evaluation. In China, some
58
studies have quantified the atmospheric As emissions from some major anthropogenic
59
sources in recent years, such as coal combustion (1564.5 tons/year) 3
ACS Paragon Plus Environment
9, 10
, cement plants
Environmental Science & Technology
Page 4 of 34
60
(166.4 tons/year)11, as well as iron and steel industry (130 tons/year) 12. The aforementioned
61
sources account for 12% of the whole emission in East Asia. And Guangdong province is
62
ranked as one of the top As emission contributors within China 9-12.
63
The Pearl River Delta (PRD) region, located in Guangdong province, south China, is
64
one of the most industrialized, urbanized and populated regions in China, and has been long
65
suffering from severe air pollution
66
(50 ng/m3) are much higher than the area outside the PRD region
67
quality standards (6 ng/m3) issued by the Ministry of Environmental Protection of P.R.
68
China (2012) and European Union (EU) Directive (2013).
13-15
. The atmospheric As concentrations in this region 16
, far exceeding the air
69
In atmosphere, more than 90% of As exists in particulate form 17. Human can expose to
70
atmospheric As through inhalation of fine particles as well as non-dietary ingestion of and
71
dermal contact with settled particles. Besides, the atmospheric As can also affect human
72
health by enhancing As level in food items and surface water through deposition. The
73
anomalous concentration in surface waters and its availability to living beings is the major
74
environmental concern of As
75
carrying speciation need to be accounted to fully understand and evaluate the environmental
76
and human health risks caused by atmospheric As contamination.
18
. Accordingly, the concentrations, deposition fluxes and its
77
There have been increasing concerns about the atmospheric As and its speciation
78
contained in particles 19-21, however, the studies on atmospheric As deposition measurement,
79
especially its carrying speciation and mobilization of individual species, are extremely
80
limited. Lack of deposition measurement data has hindered the evaluation and further 4
ACS Paragon Plus Environment
Page 5 of 34
Environmental Science & Technology
81
improvement of numerical As deposition models 8.
82
In this study, an automated wet−dry sampler (a water surface sampler for dry
83
deposition measurement), which is one of the method widely used to measure atmospheric
84
deposition fluxes of metal(loid)s
85
fluxes of atmospheric As at two representative sampling sites in the PRD region. The
86
sampling sites and measurement procedure will be introduced in detail in the Material and
87
Methodology Section.
22, 23
, is employed to measure both wet and dry deposition
88
This study firstly attempts to establish a set of experimental procedure to stabilize and
89
determine different As species of low level contained in the aqueous deposition samples.
90
Based on the trace determination results, we further (1) investigate the characteristics of both
91
wet and dry deposition of atmospheric As; (2) statistically identify the factors influencing the
92
total As deposition at the sampling sites; afterwards, based on the recognized influencing
93
factors, (3) attempt to imply the temporal/spatial variation around the PRD region; lastly, (4)
94
compare the relative contributions respectively from wet and dry deposition to the
95
environmental risks on the basis of As speciation and its mobilization.
96
2 MATERIALS AND METHODOLOGY
97
2.1 Study Area, Sampling Sites and Atmospheric Deposition Measurements
98
The PRD region is located in adjacent to the South China Sea, lying at both sides of the
99
Pearl River Estuary (Figure 1) and subject to a typical East Asian monsoon climate.
100
Southwesterly wind from the sea prevails and brings abundant precipitation in monsoon
101
season (April to September), while northeasterly wind from the mainland dominates and 5
ACS Paragon Plus Environment
Environmental Science & Technology
102
brings little precipitation in non-monsoon season (October to March). The mean annual
103
temperature from 1981 to 2010 recorded at Wushan meteorological station in Guangzhou city
104
was 22.48℃, the mean annual relative humidity was 75%. The annual precipitation ranged
105
from 1240 to 2679 mm with a mean value of 1801 mm 24.
106
Most of atmospheric particulate-As exists in the particles less than or equal to 3.5 µm 17,
107
and the PM2.5 emissions in the PRD region are located primarily in the central area, especially
108
in Guangzhou and Foshan, rather than in the east or west
109
urban site (GZ-SYSU, representing the typical urban surface area with abundant PM2.5
110
emissions) and the Mt. Dinghu Natural Reserve site (DH: representing the suburban surface
111
area with fewer emissions) were selected as two representative study sites in the present study
112
(Figure 1). The Guangzhou sampling site (GZ: 113°17’E, 23°06‘N, at the altitude of about 50
113
m) is located on the roof of the Di Huan Building (7 floors) at the Haizhu Campus of Sun
114
Yat-sen University in Guangzhou, the center of PRD region, while the Mt. Dinghu sampling
115
site (DH: 112°33’ E, 23°10’N, at the altitude of about 100 m) is located on the roof of the
116
office building of Dinghu Mountain Forest Ecosystem Research Station at the foothill of the
117
biosphere reserve, the northwest of PRD. During the study period, the annual mean
118
temperatures were 23.41℃ and 22.30℃, and the total precipitations were 2257.81 mm and
119
2070.66 mm at GZ site and DH site, respectively.
25
. Accordingly, the Guangzhou
120
A total of 119 deposition samples were collected respectively at GZ site (35 wet
121
deposition samples and 20 dry deposition samples) and DH site (40 wet deposition samples
122
and 24 dry deposition samples) continuously from March 2015 to February 2016, using the 6
ACS Paragon Plus Environment
Page 6 of 34
Page 7 of 34
Environmental Science & Technology
123
automated wet-dry samplers (Tianhong Instrument Factory, China, ASP-2), which is the
124
water surface sampler for dry deposition. The sampling procedure is based on those reported
125
in our previous studies 24, 26. In brief, the sampler was equipped with a movable polyethylene
126
cover that alternately covered the dry or wet deposition sampling dish and was regulated by a
127
rain sensor. The wet deposition samples were collected during the rainy time, whereas the dry
128
deposition samples were collected every 15-day. For wet deposition, the soluble As fraction
129
was analyzed for each collected samples, while its insoluble fraction was analyzed for the
130
aggregate 15-day sample. For dry deposition, both fractions were analyzed for each collected
131
sample. The dish for dry deposition sampling was filled with Milli-Q water to maintain a
132
water depth of approximately 2.5 cm manually. Ahead of sampling each time, the sampling
133
dishes were rinsed with 10% HCl solution and Milli-Q water for 3 times, respectively .
134
The deposition fluxes were calculated using Eqs. 1 and 2, where Fw and Fd are the wet
135
and dry deposition fluxes (µg/m2), respectively, R is the annual rainfall (mm), T is the annual
136
period without rainfall (h), ri is the recorded rainfall specific to each wet deposition sample
137
(mm), tj is the sampling period for each dry deposition sample (h), Hj is the water depth
138
specific to each dry deposition sample (cm), Ci/Cj is the total As and its species
139
concentrations in the wet/dry deposition samples (µg/l), n/m is the number of wet/dry
140
deposition samples. Information on precipitation time and amount of rainfall at both sites
141
during study period were provided respectively by the Dinghu Mountain Forest Ecosystem
142
Research Station, CAS and the Resource Platform of Atmosphere and Environmental Science,
143
Sun Yat-sen University. 7
ACS Paragon Plus Environment
Environmental Science & Technology
=
144
=
145
∑
∑
Page 8 of 34
× 10
Equation 1
× 10
Equation 2
∑ ∑
146 147 148
2.2 Preservation and Extraction of Soluble and Insoluble Compounds Based on the reported methods of preservation and stabilization of As species contained 27-30
149
in natural and polluted water
150
stabilize different As species of trace level contained in the aqueous deposition samples.
, a series of modified experiments were established to
151
For the investigation of As solubility, both wet and dry deposition samples (50-250 ml)
152
were immediately filtered with 0.45 µm membrane (Mixed Cellulose Esters, 47mm,
153
ADVANTEC)
154
ethylenediaminetetraacetic acid (EDTA) (Beyotime, GR, 0.5M) were immediately added
155
into the sample filtrate once after sampling and filtration to inhibit the As species
156
interconversion by decreasing free metal ions concentrations
157
soluble inorganic As species could be stabilized within 28 days with 1.25 mM EDTA spiked
158
(SI Figure S1). The experiments for the effective EDTA amount to stabilize As compounds
159
are described in detail in Supporting Information.
once
collected.
For
the
soluble
compounds,
aliquots
of
31, 32
. More than 90% of the
160
For the insoluble compounds, a modified orthophosphoric acid (H3PO4) (Fulka, HPLC
161
grade, 85%) microwave assistant method based on a previous study 21 was applied to extract
162
the As compounds contained on the filtered membranes. Specific concentration of 40mM
163
L-ascorbic
164
the potential transformation of iAsIII to iAsV in the extraction process. Avoid of the potential
acid (Fulka, HPLC grade, 99.9998%), as an anti-oxidant, was applied to prevent
8
ACS Paragon Plus Environment
Page 9 of 34
Environmental Science & Technology
165
photochemical reactions of iAsIII/iAsV, all of the experiments were carried out in dark
166
environment. By applying the extraction method, the extractability of total As varied from
167
80.6% to 82.4% for SRM2711a Montana reference soil (National Institute of Standards and
168
Technology, NIST). The extraction procedure and the experiments to obtain effective
169
L-ascorbic
170
Supporting Information (SI Figure S2).
acid amount spiked to stabilize As speciation are described in detail in
171
Both the preserved filtrates (with soluble As compounds) and filtered membranes (with
172
insoluble As compounds) were transported and stored in dark and cold condition (4oC), and
173
all of them have to be determined within 1 month.
174
2.3 Concentration and Determination of trace As Species
175
The As contents were determined with ICP-MS (ICAP Qc, Thermo Fisher). The
176
separation of different As species in the eluents was conducted using the HPLC system
177
(Agilent Technologies, Santa Clara, CA, USA, 1260 series) based on their retention times.
178
Injection volume was set as 10µl of eluent and the HPLC mobile phase flow rate was
179
maintained at 1.0 ml/min. The mobile phase for anion exchange chromatography included
180
12.5 mM (NH4)2CO3 and 60 mM (NH4)2CO3 (Sigma-Aldrich Chemical Co.). The outlet of
181
the HPLC column was connected directly to a concentric nebulizer of ICP-MS, allowing
182
continuous transportation of the determinants to the argon plasma of ICP-MS. Retention
183
time for the As species was determined using mixed standards of iAsIII, iAsV, cacodylic acid
184
(DMA), methylarsonic acid (MMA), arsenocholine (AsC), and arsenobetaine (AsB)
185
(National Institute of Metrology, China). Peaks of different As species were identified by 9
ACS Paragon Plus Environment
Environmental Science & Technology
186
comparison with the retention times of individual standard compounds (SI Figure S3). The
187
instrumental detection limits were 0.05, 0.05, 0.1, 0.2, 0.3 and 0.6 µg/l for AsC, AsB, iAsIII,
188
DMA, MMA and iAsV, respectively. The matrix detection limits for the soluble As
189
compounds were equal to their instrumental detection limits, while the matrix detection
190
limits were 0.02 and 0.12 µg/l respectively for the insoluble iAsIII and iAsV.
191
For those filtrates and membrane extracts under the detection limits, solid phase
192
extraction (SPE) method was applied to concentrate the As compounds 33. The experiments
193
and the QA/QC results are introduced in detail in Supporting Information. The SPE
194
concentration method was able to bring down the detection limits by ten times.
195
2.4 Data Sources
196
The meteorological data at the sampling sites, including precipitation time, rainfall
197
amount, ambient temperature, pressure, relative humidity, wind speed and irradiance, for the
198
calculation of deposition fluxes and linear regression analysis were provided by the Dinghu
199
Mountain Forest Ecosystem Research Station, CAS and the Resource Platform of
200
Atmosphere and Environmental Science, Sun Yat-sen University. Furthermore, the
201
PM2.5/PM10 levels for the linear regression analysis at sampling sites were provided by the
202
China Air Quality Online Monitoring and Analysis Platform.
203
For the temporal/spatial analysis of atmospheric deposition of total As, the monthly
204
concentrations of PM2.5/PM10 in the PRD region, Hong Kong and Macao were cited from
205
the monitoring results reports of Guangdong-Hong Kong(-Macao) Pearl River Delta
206
Regional Air Quality Monitoring Network
34-37
, where the monitoring locations were also
10
ACS Paragon Plus Environment
Page 10 of 34
Page 11 of 34
Environmental Science & Technology
207
given. . The meteorological data and their monitoring locations around the PRD region were
208
downloaded from the Web of China Met. Data Services 38, which in Hong Kong and Macao
209
were provided respectively on the websites of Hong Kong Observatory and Macao
210
Meteorological and Geophysical Bureau. All of the monitoring stations for air quality and
211
meteorological factors are marked in Figure 1, which was created using R (with the
212
packages of maptools and ggplot2). Likewise, the regional spatial variations of both wet and
213
dry deposition were also created in R based on their annual projected z scores, which would
214
be introduced and discussed in detail in Section of Implication of seasonal, annual and
215
spatial variation around the PRD region.
216
3 RESULTS AND DISCUSSION
217
3.1 Total As level in rain
218
The total As level is the sum of levels of all the detected As species. The research on
219
the levels of total As in rain and snow is very limited in the past decades. Different from the
220
limited studies around the world (SI Table S1), the concentrations of total As in rain varied
221
in a tremendous range (GZ: 0.14 -7.91µg/l, DH: 0.15 - 5.88µg/l) (SI Table S1) in the present
222
study. And the average concentrations (GZ: 0.90 µg/l, DH: 1.12 µg/l) were the highest
223
compared with Singapore (0.14 µg/l)
224
California, USA (0.46µg/l) 41, the Shale bedrock areas in Southeastern Nigeria (0.56µg/l) 42
225
and the urban area of central Poland (0.74µg/l)
226
found between two study sites, where the high As level in rain could be explained by the
227
high concentrations of atmospheric arsenic in PRD region 16.
39
, Mt. Mansfield in Vermont, USA (0.1µg/l)
40
,
43
. No significant difference (p>0.05) was
11
ACS Paragon Plus Environment
Environmental Science & Technology
228
Page 12 of 34
3.2 Atmospheric deposition of total As
229
The total As deposition flux is the sum of fluxes of all the detected As species. During
230
the study period, the atmospheric As monthly deposition fluxes varied from 23.49 to 329.22
231
µg.m-2.month-1, with the mean of 135.25 µg.m-2.month-1 at GZ site, and varied from 54.31
232
to 471.25 µg m-2.month-1, with the mean of 185.09 µg.m-2.month-1 at DH site. No significant
233
difference (p>0.05) was observed between two sites. In the comparison with other studies
234
around the world, the total deposition fluxes in the present study were extraordinarily higher
235
than those reported in US 44, Europe
236
Delta (YRD) region
237
largest coal basin in China 50, as well as the Beijing-Tianjin-Hebei region in China 51, where
238
is one of the regions with the highest atmospheric As concentration around the world
239
(SI Table S2).
47
, Taiwan
48
45
and Australia 46; and comparable to Yangtze River
and Japan
49
; whereas, lower than the Shanxi Basin, the
8, 52
240
The majority of the annual As deposition at both sites was contributed by wet
241
deposition (GZ: 78.85%, DH: 84.05%). The percentages were significantly higher than the
242
results studied in Chialy, South Taiwan 48, North China
243
Moreover, dry deposition is even found to be the more important mechanism for the annual
244
As deposition in North China
245
precipitation ranges from 430 to 1390 mm, much lower than the PRD region, and most is
246
concentrated in June – September
247
contribution of wet deposition and its relatively high contribution in summer 51. In Chialy,
248
South Taiwan, precipitation also plays an important role in the relative differences between
51
51
and most cities in Japan
and the cities across Japan
49, 53
.
49, 53
. In North China, the annual
54
, which could explain the generally less annual
12
ACS Paragon Plus Environment
Page 13 of 34
Environmental Science & Technology
249
wet and dry deposition, and the wet deposition contribution (6.9-47.1%) in relatively dry
250
months (January to April 2011, September to October 2011, December 2011 to March /2012,
251
with precipitation less than 50mm) were significantly lower than those (63.6-87.6%) in wet
252
months (May to August 2011, November 2011, April to June 2012, with precipitation more
253
than 100mm) 48. Similarly, the annual wet deposition contributions in Noshiro, Akita (1713
254
mm), Matsuura, Nagasaki (1960mm) and Kashima, Ishikawa (2157mm ) were much higher
255
than other cities with less precipitation throughout Japan
256
wet deposition played the dominating role in both wet and dry seasons throughout the study
257
period, its average contributed percentages were significantly lower in 03/2015, 04/2015,
258
11/2015 and 02/2016 (GZ: 49.73%, DH: 64.72%) , compared with other months with more
259
precipitation (GZ: 86.08%, DH: 81.31%). Hence, the relative significance of the two
260
deposition mechanisms could be attributed to the seasonal and local availability of
261
precipitation.
53
. In the present study, although
262
In addition to precipitation, it has been verified that wet deposition dominates the total
263
flux of trace elements existing in fine particles 49, 51, and the atmospheric As tends to exist in
264
fine airborne particles 17, 55. Accordingly, existence of atmospheric As in fine particles might
265
be the other reason to explain the higher partition of wet deposition. The influencing factors
266
for both wet and dry deposition would be described in more detail in the next section.
267
In general, atmospheric As deposition fluxes (wet and dry) and its concentration in rain
268
are being influenced by many complicated factors, such as emission, ambient air
269
concentration, precipitation, scavenging, deposition velocity and transport processes in 13
ACS Paragon Plus Environment
Environmental Science & Technology
270
atmosphere. And its transport and deposition processes can be indicated by airborne
271
particulate matters (PM)
272
PRD region has been reported contributed from northeastern and southwestern Guangdong as
273
well as the super- region
274
deposition flux and its concentration in rain are found between the two study sites, even
275
though the PM2.5 emission level at GZ site are relatively higher compared with DH site 25.
276
3.3 Factors influencing atmospheric deposition of total As
8, 56
. Except the local emission, most of the ambient PM2.5 in the
57
, so it is not surprising that no significant differences of As
277
When we attempted to investigate the factors influencing the wet and dry deposition of
278
total As, the levels of particulate matters (PM) were assumed to explain the ambient air As
279
contamination in the present study.
280
3.3.1 Wet deposition
281
At both sites, no significant correlation was found between the fluxes and total As
282
levels in rain (p>0.05). However, the monthly wet deposition fluxes varied consistent with
283
the precipitation (r=0.813, p