Characterizing spatial diversity of passive sampling sites for

Aug 14, 2018 - More than 200 time series (length 5–11 years) were used to ... 14 sites were assessed using data from the MONET passive sampling netw...
0 downloads 0 Views 2MB Size
Subscriber access provided by UPSTATE Medical University Health Sciences Library

Environmental Measurements Methods

Characterizing spatial diversity of passive sampling sites for measuring levels and trends of semivolatile organic chemicals Ji#í Kalina, Martin Scheringer, Jana Boruvkova, Petr Kukucka, Petra P#ibylová, Ond#ej Sá#ka, Lisa Melymuk, Milan Vana, and Jana Klánová Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b03414 • Publication Date (Web): 14 Aug 2018 Downloaded from http://pubs.acs.org on August 19, 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.

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 29

Environmental Science & Technology

1

Characterizing spatial diversity of passive

2

sampling sites for measuring levels and trends of

3

semivolatile organic chemicals

4

Jiří Kalina1, Martin Scheringer*1, 2, Jana Borůvková1, Petr Kukučka1, Petra Přibylová1,

5

Ondřej Sáňka1, Lisa Melymuk1, Milan Váňa3, Jana Klánová1

6

1 Research Centre for Toxic Compounds in the Environment RECETOX,

7

Kamenice 5, 625 00 Brno, Czech Republic

8

2 Institute for Chemical and Bioengineering, ETH Zürich,

9

8093 Zürich, Switzerland, [email protected]

10

3 Czech Hydrometeorological Institute – Košetice Observatory,

11

394 22 Košetice, Czech Republic

12

KEYWORDS: Temporal trends, polychlorinated biphenyl, polycyclic aromatic hydrocarbon,

13

organochlorine pesticide, passive air sampling, MONET.

14

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 29

15

ABSTRACT: Passive air sampling of semivolatile organic compounds (SVOCs) is a

16

relatively inexpensive method that facilitates extensive campaigns with numerous sampling

17

sites. An important question in the design of passive-sampling networks concerns the number

18

and location of samplers. We investigate this question with the example of 17 SVOCs

19

sampled at 14 background sites across the Czech Republic. More than 200 time series (length

20

5–11 years) were used to characterize SVOC levels and trends in air between 2003 and 2015.

21

Six polychlorinated biphenyls (PCBs), 6 polyaromatic hydrocarbons (PAHs) and 5

22

organochlorine pesticides (OCPs) at 14 sites were assessed using data from the MONET

23

passive sampling network. Significant decreases were found for most PCBs and OCPs

24

whereas hexachlorobenzene (HCB) and most PAHs showed (mostly insignificantly)

25

increases. Spatial variability was rather low for PCBs and OCPs except for

26

dichlorodiphenyltrichloroethane (DDT) and rather high for PAHs. The variability of the

27

SVOC levels and trends depends on characteristics of the sites including their remoteness,

28

landscape, population and pollution sources. The sites can be grouped in distinct clusters,

29

which helps to identify similar and, thereby, potentially redundant sites. This information is

30

useful when monitoring networks need to be optimized regarding the location and number of

31

sites.

ACS Paragon Plus Environment

2

Page 3 of 29

Environmental Science & Technology

32

TOC Art

33 34

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 29

35

Introduction

36

Passive air sampling is a relatively inexpensive and efficient method for long-term

37

monitoring of semivolatile organic chemicals (SVOCs) such as persistent organic pollutants

38

(POPs) in air. Application of this technique enabled the development of the extensive passive

39

air monitoring networks on the regional and even global scales in the last decade. This effort

40

was largely driven by the Global Monitoring Plan of the Stockholm Convention. An

41

important question related to the set-up of passive air sampling networks is how the

42

minimum number of sampling sites can be found that are needed to represent the levels and

43

trends of an SVOC in the area of interest. The representativeness of a network is defined in

44

the Guidance on the Global Monitoring Plan for Persistent Organic Pollutants1 as a

45

“sufficient number of sampling sites to make general conclusions about POPs trends” in the

46

area.

47

Many studies have demonstrated that for relatively short sampling durations (weeks or

48

months) the variability between measured concentrations is rather high even in relatively

49

small regions and between sites classified as background sites.2–9 Long-term monitoring

50

makes it possible to reduce this variability because short-term fluctuations are averaged out.

51

Here, we address the question of the sufficient number of sampling sites in a network

52

using a large set of long-term passive air-sampling data. The data consist of mostly

53

unpublished results (only 17 of the 238 time series used here were previously published)

54

from the Czech Republic. Time series from 14 sites officially classified as background sites

55

according

56

Hydrometeorological Institute10,11,12 were analyzed for polychlorinated biphenyls (PCBs),

57

organochlorine pesticides (OCPs), and polycyclic aromatic hydrocarbons (PAHs). The length

58

of the time series collected in the period 2003-2015 was between 7 and 13 years. Annually

59

aggregated data for each chemical at each site were used to determine the average SVOC

to

the

Exchange-of-Information

(EoI)

classification

ACS Paragon Plus Environment

of

the

Czech

4

Page 5 of 29

Environmental Science & Technology

60

concentrations and the slopes of the time trends. Using the levels and long-term trends of the

61

17 SVOCs at the 14 sites, we demonstrate that it is possible to define clusters of sites with

62

similar SVOC pollution patterns. This makes it possible to define a minimum number of sites

63

and select a representative set of sites that is sufficient to draw conclusions about SVOC

64

levels and trends in the area. Defining a set of representative sites from a larger network

65

according to the approach presented here will be instrumental in securing the sustainability of

66

the monitoring networks built in the last decade in support of the Global Monitoring Plan of

67

the Stockholm Convention.

68 69

Materials and methods

70

Sampling design. In the 2002–2015 period, SVOC passive air sampling was conducted at

71

least once at 232 sites within the Monitoring Network (MONET) program.13 Out of these

72

sites, 40 sites operated by the Czech Hydrometeorological Institute are considered as

73

background sites according to the EoI classification based on the EU Council decision

74

97/101/EC and Commission decision 2007/752/EC,10,11,12 which means they are neither

75

industrial nor affected by traffic according to the detailed specification.10 A sub-set of 14

76

background sites provided time series sufficiently long to be included in this analysis.

77

Passive air samplers were deployed at all sites consisting of two stainless steel bowls to

78

form a protective chamber in which a polyurethane foam (PUF) disc was mounted.13 The

79

samplers were placed at 1.5 m height above the ground to allow the air to blow freely through

80

the gap between the bowls.

81 82

Experimental section. PUF discs of 15 cm diameter and 1.5 cm thickness were used for

83

sampling. After the sampling period (typically 28 days, see below) the discs were transported

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 29

84

to the accredited trace analytical laboratory at RECETOX, Brno, for extraction, clean up and

85

analysis. Analytical steps are described in detail by Kalina et al.14

86

From the total number of 30 monitored SVOCs, 17 were selected for the analysis because

87

they are not affected by sampling difficulties such as concentrations close to the detection

88

limit or breakthrough.14 The selection was made in order to represent three groups of SVOCs

89

with a similar number of substances: 5 OCPs, 6 indicator PCBs and 6 PAHs ranging from 3-

90

ring up to 6-ring. A list of the analytes is given in Table 1. Results of the analysis were given

91

in units of mass of chemical per PUF disc and the period of the exposure. Because neither the

92

temporal or spatial variability assessment requires working with air concentrations in units of

93

mass per air volume,14 the primary values in units of mass per PUF and 28 days were used,

94

which are not affected by known shortcomings of the conversion into air concentrations15

95

(see discussion). Measurement results converted into air concentrations are provided in the

96

SI file Kalina_monthly_data.xlsx.

97 98

Table 1: List of analyzed SVOCs, their abbreviations, and percentages of their values below

99

the limit of quantification. OCPs

PCBs

PAHs

α-HCH

A_HCH

3% PCB 28

PCB28

2% fluorene

FLU

0%

γ-HCH

G_HCH

2% PCB 52

PCB52

3% phenanthrene

PHE

0%

p,p‘-DDE

PP_DDE

0% PCB 101

PCB101

9% pyrene

PYR

0%

pentachloro benzene

PECB

5% PCB 137

PCB137

14% chrysene

CHRY

6%

hexachloro benzene

HCB

0% PCB 153

PCB153

4% benzo(a) pyrene

BAP

48%

PCB 180

PCB180

26% benzo(ghi) perylene

BGP

56%

100

ACS Paragon Plus Environment

6

Page 7 of 29

Environmental Science & Technology

101

Data treatment. In order to achieve robust results with the smallest possible influence of

102

random effects, only the background sites with at least 5-year time series (i.e. continuous

103

sampling took place in 5 consecutive calendar years) were taken for further processing. The

104

five-years span encompasses also incomplete initial/final years, when the sampling started or

105

ended within the year. Because these initial and final years of the time series are usually not

106

fully covered by the samples, the values from incomplete initial and final years were

107

excluded. This provided a set of 14 relevant sites with sampling durations from 7 to 11 years

108

(2004–2014) that were used for the analysis; a list of these 14 sites is provided in Table 2.

109 110

Table 2: List of background sites, their EoI classification12, spatial coordinates, and sampling

111

periods. EoI classification: classification according to criteria defined under the Exchange-of-

112

Information decision of the European Council.10,11 B: background, N: natural, R: countryside,

113

A: agricultural, S: suburban, AN: agricultural/natural, REG: regional. Sampling period

Name of the site

Longitude Latitude EoI Letter classification

Bily Kriz, Beskydy

X

2006– B/R/N-REG 49.50261 18.53856 2015

Decinsky Sneznik, Labske piskovce

D

2006– B/R/N-REG 50.78951 14.08684 2015

Churanov, Sumava, EMEP

C

2006– B/R/N-REG 49.06844 13.61488 2015

Jesenik, Jeseniky

J

B/R/N-NCI

2006– 50.24225 17.19022 2015

Kosetice, EMEP

E

B/R/ANREG

2003– 49.57345 15.08041 2015

L

B/R/N

2005– 50.72940 14.98790 2015

M

2006– B/R/A-REG 48.79175 16.72450 2015

Liberec, Jested Mikulov, Sedlec

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 29

2004– 50.00731 14.44620 2015

Praha, Libus

P

B/S/R

Primda, Sumava

I

2006– B/R/N-REG 49.66959 12.67785 2015

Rudolice v Horach, Krusne hory

R

2006– B/R/N-REG 50.57979 13.41922 2015

Rychory, Krkonose

Y

2006– B/R/N-REG 50.66046 15.85006 2015

Serlich, Orlicke hory

S

2006– B/R/N-REG 50.32804 16.38353 2013

Stitna nad Vlari-Popov, Planavy

T

2004– B/R/N-REG 49.04776 18.00781 2015

Svratouch

H

B/R/ANREG

2006– 49.73507 16.03413 2015

114 115

With PUF disc passive sampling, both the gas and particle phases are analyzed together,

116

which yields one value for each compound within the sample. Certain fractions of the results

117

were below the limit of quantification, LoQ, which leads to left-censored values. The LoQ is

118

calculated as 3.3 times the limit of detection, LoD; this methodology is described in detail by

119

Kalina et al.14 Percentages of values below the LoQ for individual compounds are listed in

120

Table 1.

121

Values below the LoQ were replaced by ½ of the LoQ as described in Kalina et al.14 The

122

most frequent values below the LoQ were found for heavier PAHs, which fluctuate strongly

123

between the summer and the winter. In such cases, the winter concentrations, which are

124

higher by several orders of magnitude, determine the results of the annual aggregation of the

125

concentrations.

126

Although the typical sampling period was 28 days for almost all samples, a part of the

127

samples differed in the length of the sampling interval (85% were sampled for 28 days, 3%

128

for 27 days, 3% for 29 days, 1% for 26 days, 1% for 30 days, and the remaining 7% were

129

sampled for different durations). It may be assumed that all sampling periods from minimum

ACS Paragon Plus Environment

8

Page 9 of 29

Environmental Science & Technology

130

of 14 days up to maximal period of 98 days fall into the linear phase of the uptake curve of

131

passive air sampling.15 Under these assumptions it was possible to standardize the amounts of

132

the analytes sequestered by the samplers to a sampling duration of 28 days, which provides

133

approximately 13 samples per year.

134

Annual aggregation. Considering the strong seasonal fluctuations of the SVOC air

135

concentrations,14 statistical evaluations that are sensitive to periodic changes cannot be

136

performed on the primary data. Therefore, we aggregated the values annually to average out

137

both the seasonal within-year fluctuations and possible random fluctuations. A standardized

138

method of the annual aggregation is taking the arithmetic mean.1 As shown in our previous

139

work,14 annual aggregation also avoids inconsistency of primary data caused by seasonal

140

differences in sampling rates. Under the assumption of long-term stable conditions, see

141

multi-year temperature data for site “E” (Kosetice) in Table S1 and Figure S1, and thus

142

similar within-year changes of the sampling rates, the annually aggregated values can be used

143

to establish the long-term time trend.

144

Temporal trends. It may generally be assumed that SVOCs in ambient air undergo

145

degradation processes that can be approximated by first-order kinetics.16–18 These reactions

146

are described by compound-specific half-lives. For background sites not affected by primary

147

emissions, this leads to an exponential decrease of the compounds’ levels in air.14

148

The exponential trend fit was performed by using the Theil-Sen estimator, a robust tool for

149

estimating trends used for regression of air pollution data.19–22 The trend was fitted as a linear

150

regression on log-transformed annually aggregated data on each time series obtained for one

151

SVOC at one sampling site. The assessment of significance was made for each of the 238

152

time series given by all combinations of sampling sites and chemicals; each time series was

153

characterized by the median level of the SVOC concentration (intercept), the half-life and

154

direction of the trend (slope) and its statistical significance.

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 29

155

Additionally, a median trend was computed for each compound from all individual trends

156

of the 14 sites with different levels of pollution, lengths, slopes and statistical significance. A

157

method of median slope and median intercept was used to determine these median trends.

158

Based on the principles of the Theil-Sen estimator on log-transformed data, this method

159

computes the median of the slopes and the median of the intercepts over the 14 sites and uses

160

these statistics to form a new trend. This makes it possible to compute a median half-life and

161

its significance for each compound. These half-lives are listed in Table 3.

162

Spatial variability. The selected 14 sites are spread over the whole area of the Czech

163

Republic; they are located mainly in the countryside, at higher elevations, or in the

164

surroundings of bigger municipalities. For each of the 17 compounds, these sites provide 14

165

time series of different length, starting in 2007 or earlier and ending in 2010 or later. Since all

166

of the time series include the year 2009, which is a middle date of the overall period (2004–

167

2014), this year was chosen as a reference to determine the levels of the compounds at the

168

sites. The levels were determined as the values of the Theil-Sen trend in the year 2009 in

169

units of mass per PUF and 28 days.

170

The interquartile ratio (IQR) was used to compare the spread of the levels of individual

171

compounds among the 14 sites. The value of an interquartile ratio is computed as a difference

172

between the 3rd and the 1st quartile (the interquartile range), divided by the median.

173

Cluster analysis. Each of the 17 SVOCs provides two characteristics of each site: the SVOC

174

level and its temporal trend. This allows us to construct a 34-dimensional (34-d) space in

175

which each site is represented as a point characterized by the 17 levels and the 17 trends of

176

the SVOCs (normalized values were used to handle the different magnitudes of SVOC levels

177

and trends). In this space, it is possible to characterize the sites according to their proximity

178

(measured by their Euclidean distance in the 34-d space) and to group them into several

ACS Paragon Plus Environment

10

Page 11 of 29

Environmental Science & Technology

179

clusters; a cluster is a set of sites that are closer to one another than to sites in a different

180

cluster.

181

The K-means clustering method was used to define three distinct clusters of sites based on

182

their mutual distances in the 34-d space.23,24 These clusters represent three typical types of

183

sampling sites within the Czech Republic. A principal-component analysis was then

184

performed to project the 34-d space to a plane spanned by two dimensions to visualize the

185

clusters.

186

Cluster analysis over data from passive air sampling was used in earlier studies25–27 to

187

separate monitoring sites influenced by different pollution sources and characterize local

188

fingerprints of groups of the sites. One of these studies25 also used primary data on pollution

189

levels to construct a multidimensional space for the clustering.

190 191

Results and discussion

192

The positions of the 14 background sites within the Czech Republic, along with their SVOC

193

levels, trend direction and statistical significance, are represented by colored dots in Figure 1.

194

Six characteristic compounds are shown to demonstrate the situation for lighter (PCB 28) and

195

heavier (PCB 180) PCBs, two main types of OCPs (γ-HCH and p,p’-DDE), and a lighter

196

(fluorene) and heavier (chrysene) PAHs. Maps for all 17 compounds are provided in

197

Figure S2 and individual SVOC levels, half-lives and rate constants of time trends are shown

198

in a colored heatmap as Table S5.

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 29

199

Figure 1: Sampling sites with different trends of SVOCs in ambient air. 14 background sites

200

are shown as colored dots of different size with the diameter depending logarithmically on

201

the pollution level at the site. Statistically significant decrease is represented by green,

202

increase by red, and insignificant change by gray color. Maps created with data from

203

http://www.diva-gis.org/Data and the ArcGIS ® software, v10.5, by ESRI, Environmental

204

System Research Institute (ESRI). ArcGIS® and ArcMap™ are the intellectual property of

205

Esri and are used herein under license.

ACS Paragon Plus Environment

12

Page 13 of 29

Environmental Science & Technology

206

The maps in Figure 1 show the spatial coverage and differences in SVOC levels and trends,

207

but do not facilitate a direct comparison of their magnitude. These are shown in time-

208

concentration plots for individual sites in Figure 2 (six of the total number of 238 plots are

209

shown) and all sites in one plot in Figure 3. The same six compounds as in Figure 1 were

210

chosen for Figures 2 and 3.

211

Figure 2: Individual time series of selected combinations of sites and chemicals (PCBs,

212

OCPs and PAHs). Insignificant exponential trends are drawn as gray lines, significantly

ACS Paragon Plus Environment

13

Environmental Science & Technology

Page 14 of 29

213

decreasing trends as green lines, and significantly increasing trends as red lines. The shaded

214

area represents the 95% confidence interval of the trend.

215

Figure 3: Exponential regression of time series of selected PCBs, OCPs and PAHs.

216

Insignificant trends are drawn as gray lines, significantly decreasing trends as green lines, and

217

significantly increasing trends as red lines (the only significant increase was observed for

218

chrysene). Letters indicate sites as in Table 2. The median trend computed by the method of

ACS Paragon Plus Environment

14

Page 15 of 29

Environmental Science & Technology

219

median slope and median intercept is depicted as a thick darker line. Six representative

220

compounds were selected, for others see Figure S3 in the Supporting Information.

221

Exponential trends, median levels of concentrations and half-lives of the median trends

222

computed for the 17 compounds by the method of median slope and median intercept are

223

listed in Table 3. Interquartile ratios are shown to indicate the variance of the levels. More

224

detailed values of SVOC levels and trends are provided in Table S5.

225 226

Table 3: Median trends, half-lives, and median levels with their interquartile ratio.

227

Statistically significantly decreasing trends are highlighted in green, increasing in red. compound

fluorene

overall exponential trend (∆T = years counted from 2004) [ng/PUF/28 days]

median annual increase level (2009) [%] [ng/ PUF/28 days]

interquartile ratio [–]

overall half-life [year]

9,58 × exp(-0.02 × ∆T)

-2.3

233.9

0.74

30.5

13.75 × exp(-0.04 × ∆T)

-4.0

651.1

0.78

17.2

6.78 × exp(-0.01 × ∆T)

-1.2

93.2

1.59

59.2

chrysene

3.33 × exp(0.06 × ∆T)

+6.4

7.86

2.05

-10.8

benzo(a)pyrene

1.05 × exp(0.03 × ∆T)

+3.4

0.77

0.93

-20.5

benzo(ghi)perylene

0.90 × exp(0.02 × ∆T)

+1.9

0.63

0.82

-37.3

α-HCH

0.77 × exp(-0.12 × ∆T)

-11.4

1.92

0.35

6.09

γ-HCH

1.54 × exp(-0.11 × ∆T)

-21.0

3.36

0.41

3.30

p,p‘-DDE

1.53 × exp(-0.02 × ∆T)

-1.9

3.27

1.20

36.3

pentachlorobenzene

0.88 × exp(-0.07 × ∆T)

-6.9

1.58

0.22

10.1

hexachlorobenzene

2.84 × exp(0.03 × ∆T)

+2.5

8.38

0.17

-27.5

PCB 28

0.65 × exp(-0.09 × ∆T)

-9.4

1.04

0.41

7.41

PCB 52

0.45 × exp(-0.15 × ∆T)

-15.2

0.85

0.45

4.56

PCB 101

0.42 × exp(-0.10 × ∆T)

-9.9

0.40

0.98

7.03

phenanthrene pyrene

ACS Paragon Plus Environment

15

Environmental Science & Technology

Page 16 of 29

PCB 138

0.28 × exp(-0.15 × ∆T)

-14.6

0.26

0.69

4.74

PCB 153

0.40 × exp(-0.13 × ∆T)

-13.0

0.49

0.70

5.32

PCB 180

0.28 × exp(-0.09 × ∆T)

-9.0

0.15

0.93

7.69

228 229

PAHs. The range of PAH levels among all sites is rather high compared to the other

230

compounds, namely up to two orders of magnitude (IQR 0.73–2.05). The lowest level was

231

found for BGP at site C (0.383 ng/PUF/28 days), the highest for PHE at site P (1.251

232

ng/PUF/28 days).

233

The sites D, T, H and P showed the highest levels for most of the PAHs. While site D is

234

presumably affected by sources in Germany,28 the site P located on southern edge of the

235

capital city of Prague could possibly be influenced by car traffic (total length of roads within

236

5 km diameter around the site is 999 km compared to median of 170 km over the whole set of

237

sites, for details see Table S3 (note that this site is still officially classified as a background

238

site and is not directly affected by local traffic emissions12). An explanation for increased

239

levels of PAHs for all the sites D, T, H and P is provided by an observed linear relationship

240

between the logarithm of population and the logarithm of lignite consumption which is the

241

most common solid fuel used for domestic heating in the Czech Republic: data from 145

242

Czech municipalities provided an expected annual consumption of approx. 2,000 Mg/year of

243

lignite for the sites T and H (2,700 Mg/year for D and 7,300 Mg/year for P), see Table S4 and

244

text on p. S7 in the SI. This is a significantly higher amount than the median estimated

245

consumption of the remaining sites, which is approx. 1,000 Mg. For details see Table S4.

246

In contrast, the lowest levels of PAHs were sampled at site C in the Sumava mountains in

247

southern Bohemia and at sites J and S in the mountains in the northern part of the country.

248

These mountains are far from cities and industrial centers with large areas covered by

249

coniferous forests.

ACS Paragon Plus Environment

16

Page 17 of 29

Environmental Science & Technology

250

The trend of PAHs pollution is in most cases insignificant, usually with negative estimated

251

half-lives, which means that both a concentration increase or decrease are possible. The only

252

significant increase was detected for chrysene at site E with a 10-years time series available.

253

The median trend calculated from all 14 sites reveals a statistically significant increase for

254

chrysene and decreases for fluorene and phenanthrene. For several PAHs the initial phase of

255

the time series reveals decreasing tendencies and changes to insignificant trends in the later

256

phase.

257

OCPs. In the cases of HCHs, PECB and HCB, the range of the levels is rather low with IQRs

258

between 0.17 and 0.41, which is less than one order of magnitude. For p,p’-DDE the spread

259

of the levels is intermediate with an IQR of 1.19, which is approximately one order of

260

magnitude. The highest levels of all 7 OCPs were detected in the same group of sites (D, H,

261

P, T) as for PAHs. High levels of p,p’-DDE were also found at site M. The lowest HCH

262

levels were found at mountainous sites such as C, L, S and X in the Sumava and northern

263

mountains.

264

At several sites, the half-lives of HCB and p,p’-DDE were detected as negative with

265

insignificantly increasing trends or as positive but with half-lives longer than 10 years. The

266

other OCPs were decreasing with different significance at all 14 sites, for γ-HCH significant

267

decreases were found at all 14 sites. Half-lives were between 2.3 and 11.5 years with the

268

exception of sites E, M and X with longer half-lives for PECB. This could be due to an

269

agricultural influence on sites E and M (63% (site E) and 67% (site M) of the area within 5

270

km diameter is agricultural land, compared to the median of 5% over all 14 sites, see Table

271

S3). The median trend showed a statistically significant decrease for both HCHs and PECB,

272

an insignificant trend for p,p’-DDE, and a significant increase for HCB.

273

PCBs. The range of PCBs is relatively low with IQRs between 0.41 and 0.97, which

274

corresponds to a variance of less than one order of magnitude. The highest levels of all 6

ACS Paragon Plus Environment

17

Environmental Science & Technology

Page 18 of 29

275

PCBs were found at site D, which could be affected by long-range transport from industrial

276

areas of northern Germany28 (2.2 ng/PUF/28 days for PCB28) and P (2.4 ng/PUF/28 days for

277

PCB28) in the suburb of the capital, Prague. The lowest levels of all PCBs were detected at

278

the mountainous sites in the northern part of the country (J, L, S, X).

279

With the only exception of PCB101 at site M, all PCBs were decreasing and their median

280

half-lives were about 6 years. The highest median half-life was detected for PCB 180 (7.7

281

years), the lowest for PCB 52 (4.6 years). These findings are in good accordance with many

282

previously published results.14,16–18 All 6 PCBs decreased significantly at site T, most of them

283

decreased significantly also at site P.

284

Clustering

285

The K-means clustering method was used to identify clusters of similar sites based on the

286

levels and trends of all 17 chemicals at all 14 sites. Three characteristic clusters were

287

identified, representing three types of sites that are typical by their pollution levels and trends

288

of the SVOCs. The results of the clustering with a characterization of the clusters and their

289

differences in SVOC levels and trends are provided in Table 4.

290 291

Table 4: Clusters of similar sites identified by the K-mean clustering method.

Number Character of the of the cluster cluster

Sites within the cluster

OCPs

PCBs

PAHs

1 (red)

Rural sites with C, E, I, M, increasing levels R, Y (rural/ of PAHs agricultural)

Mostly moderate levels and slow, mainly insignificant decrease.

Moderate levels and slow, mainly insignificant decrease.

Low levels, (sometimes insignificantly) increasing.

2 (green)

Mountain sites with low levels of all SVOC types and decreasing

Low levels and insignificant decrease except of HCB.

Low levels and partially significant decrease.

Low and (mostly insignificantly) decreasing levels.

J, L, S, X (remote/ mountainous)

ACS Paragon Plus Environment

18

Page 19 of 29

Environmental Science & Technology

levels of PAHs 3 (blue)

Sites with increased levels of all SVOC types

D, H, P, T (suburban/ rural)

Increased levels with different trends.

High and partially significantly decreasing levels.

High and rather stable levels.

292 293

In order to visualize the results of the K-mean clustering, a principal-component analysis

294

was conducted, which rotates the 34-d space and projects it into the plane of the 1st and 2nd

295

dimension, see Figure 4. In our case, these two dimensions explain in total 61% of the overall

296

variability of the dataset.

297

298 299

Figure 4: Clusters found in the 34-d space after rotation by PCA and projection onto the

300

plane spanned by the first two principal components. The three different clusters are depicted

301

by different colors. Similar sites are close to each other in the 34-d space, but not necessarily

ACS Paragon Plus Environment

19

Environmental Science & Technology

Page 20 of 29

302

in reality (see map in Figure 5). The first principal component is depicted on the x-axis,

303

correlating positively with the levels of all three substance groups, PAHs, PCBs and OCPs,

304

and explaining 42% of the variability; the second principal component is depicted on the y-

305

axis, correlating negatively with PAHs decrease rate constants. It explains 19% of the

306

variability. An overall centroid site is depicted as the yellow dot in the middle of 34-d space.

307 308

Figure 5: Location of sites from the three different clusters within the Czech Republic. Sites

309

from the cluster with increased SVOC pollution levels are shown in blue (cluster 3),

310

northern-mountain sites with decreasing PAHs levels in green (cluster 2), and rural sites with

311

moderate levels of SVOCs and even increasing levels of PAHs in red (cluster 1). Map from

312

https://www.czso.cz/csu/czso/13-2108-06-v_letech_2000_az_2005-h000.

313

The first principal component explains 42% of the total variability and correlates positively

314

with most SVOC concentrations. In fact, the cluster with the highest levels of PAHs (cluster

315

3) reveals also the highest levels of OCPs and PCBs and vice versa (cluster 1). The second

316

principal component, explaining another 19% of the total variability between the sites, is

317

positively correlated with the decrease in (mainly the lighter) PAHs (higher positive values of

ACS Paragon Plus Environment

20

Page 21 of 29

Environmental Science & Technology

318

PC2 means faster decrease, lower negative values of PC2 means faster increase). This means

319

PAHs are the most variable compounds between the sites, contributing mostly to the

320

differentiation of the clusters.

321

Sites D, H, P and T form the cluster with the highest pollution levels in OCPs, PCBs as

322

well as PAHs (cluster 3, blue). The remaining 10 sites are less polluted, with their main

323

difference in the PAHs half-lives. The cluster of sites J, L, S and X (cluster 2, green)

324

represents the least polluted sites with a clear and relatively fast decrease in PAHs levels

325

(half-lives around 10 years), compared to the cluster of sites C, E, I, M, R and Y (cluster 1,

326

red) with a slow increase in PAHs levels (median doubling time around 40 years).

327

The same coloring of the clusters is used in Figure 5, showing the distribution of the

328

clusters over the Czech Republic. The lowest levels of most SVOCs and decreasing PAH

329

trends are typical for the green cluster of sites in mountainous areas in the northern part of the

330

country (cluster 2). (The site Y (red, cluster 1) also located in that area has specific features

331

with lighter PAHs decreasing (FLU, PHE, PYR) and heavier ones increasing (CHRY, BAP,

332

BGP).) Considering a northwest-to-southeast prevailing wind direction, these sites are not

333

affected by PAH sources from the Czech Republic, but may be influenced by sources in

334

Germany and Poland.28 OCP and PCB levels at the sites of the green cluster are also the

335

lowest ones, which could be related to the high elevation of the sites within this cluster

336

(average elevation 840 m).

337

The blue cluster of sites with the highest PAH, OCP and PCB levels (cluster 3) is probably

338

affected by emissions from nearby settlements. These sites have the lowest elevation (550 m

339

on average) and lie close to densely populated areas (note that they are nevertheless officially

340

classified as background sites12). The levels of all SVOCs are relatively high and they change

341

only slowly with rather stable levels mainly for heavier PAHs.

ACS Paragon Plus Environment

21

Environmental Science & Technology

Page 22 of 29

342

The red cluster (cluster 1) contains moderately polluted sites in mostly rural areas differing

343

from cluster 2 (green) mainly due to increasing levels of PAHs and slower decrease of OCPs.

344

Most of the sites are located in southern parts of the Czech Republic with possible influences

345

of agricultural areas.

346

If a representative set of monitoring sites is to be set up, it is necessary to have at least one

347

site from each of the three clusters. This means that three is also the lowest possible number

348

of sampling sites when one intends to optimize the number of samplers while capturing most

349

of the variability of the SVOC levels and temporal trends found for the Czech Republic.

350

Using the primary values in units of mass per PUF and 28 days encompasses several

351

sources of uncertainty but also avoids possible shortcomings of the conversion into air

352

concentrations14,15. The main sources of uncertainty are different environmental conditions

353

(mostly temperature) changing between both seasons and sites. The seasonal changes of the

354

sampling rates are mostly averaged out by the annual aggregation of the data and the trends

355

from the passive samplers agree very well with trends determined with active samplers.14

356

Regarding the spatial differences, we used the recalculation model by Harner29 in

357

combination with site-specific temperatures to estimate sampling rates for all 14 sites. Long-

358

term monthly temperature averages30 were used to calculate the sampling rates; this showed

359

that differences between the sites were for most of the compounds less than 12% (with the

360

exceptions of PCB 138, PCB 180 and BAP, for which the differences are still below 30%),

361

see the SI file Kalina_monthly_data.xlsx. These uncertainties are small compared to the

362

differences of measured SVOC levels, see interquartile ratios above.

363 364

Environmental significance

365

Based on differences in pollution levels and trends, a set of 14 passive air monitoring sites in

366

the Czech Republic was classified into three clusters. There is a cluster of sites with high and

ACS Paragon Plus Environment

22

Page 23 of 29

Environmental Science & Technology

367

stable levels of pollution of all three groups of compounds and two clusters with lower levels

368

of SVOCs, differing in whether PAH levels are increasing or decreasing. For the cluster with

369

lowest SVOC levels a further decrease in PAH levels is typical whereas for the rest of the

370

sites PAHs are increasing.

371

Our results show that passive air sampling makes it possible to define clusters of

372

background sites with similar pollution profiles, which can help to optimize present

373

monitoring networks. In the case of the Czech Republic, three clusters of different type were

374

found (remote/mountain, rural/agricultural and sub-urban/rural), but the approach presented

375

here can be applied similarly to other passive sampling networks with long-term data series.

376

These data are currently being generated under the umbrella of the Global Monitoring Plan

377

(GMP) of the Stockholm Convention31. The networks included in the GMP cover areas on

378

larger (regional or continental) scales where other types of background sites (marine, polar,

379

remote) can also be included. A cluster analysis as presented here will reveal whether there

380

are several sites that are sufficiently similar so that they form a cluster and some of them may

381

be removed from the network or whether all sites present in the network are needed. In the

382

current discussion about the sustainability of sampling networks under the GMP (driven by

383

the costs of their long-term maintenance) it is highly important to develop strategies for

384

minimizing the sampling effort while the information gain from the networks is not

385

substantially reduced. The work presented in this study is a contribution to this process.

386 387

Acknowledgement

388

This work was supported by the Czech Ministry of Education, Youth and Sports (RECETOX

389

RI: LM2015051, ACTRIS_CZ RI: LM2015037) and the European Structural and Investment

390

Funds

(RECETOX

RI:

CZ.02.1.01/0.0/0.0/16_013/0001761,

ACS Paragon Plus Environment

ACTRIS_CZ

-

RI:

23

Environmental Science & Technology

Page 24 of 29

391

CZ.02.1.01/0.0/0.0/16_013/0001315, and CETOCOEN PLUS: CZ.02.1.01/0.0/0.0/15_003/

392

0000469).

393

Supporting Information

394

Average daily temperatures in Kosetice (Table S1); average daily temperatures in Kosetice

395

(Figure S1); list of sites suitable for trend estimation, their spatial coordinates, and sampling

396

periods (Table S2); characteristics of the sites including elevation, density of traffic networks,

397

expected PAHs emissions (Table S3); expected lignite consumption in the vicinity of the sites

398

(Table S4); sampling sites of different trends (Figure S2); exponential regression of time

399

series of selected PCBs, OCPs and PAHs (Figure S3); median SVOC levels, half-lives and

400

rate constants of SVOC time trends at individual sites (Table S5) (PDF).

401

Monthly averages of primary data from passive monitoring, estimated concentrations in air

402

and average sampling rates (XLS).

403 404

References

405

(1)

Conference of the Parties to the Stockholm Convention on Persistent Organic

406

Pollutants. Guidance on the Global Monitoring Plan for Persistent Organic Pollutants.

407

United Nations Environment Program 2015.

408

(2)

Melymuk, L.; Robson, M.; Helm, P. A.; Diamond, M. L. PCBs, PBDEs, and PAHs in

409

Toronto Air: Spatial and Seasonal Trends and Implications for Contaminant Transport.

410

Sci. Total Environ. 2012, 429, 272–280.

411

(3)

Khairy, M.; Muir, D.; Teixeira, C.; Lohmann, R. Spatial Trends, Sources, and Air–

412

Water Exchange of Organochlorine Pesticides in the Great Lakes Basin Using Low

413

Density Polyethylene Passive Samplers. Environ. Sci. Technol. 2014, 48 (16), 9315–

ACS Paragon Plus Environment

24

Page 25 of 29

Environmental Science & Technology

414 415

9324. (4)

Tombesi, N.; Pozo, K.; Harner, T. Persistent Organic Pollutants (POPs) in the

416

Atmosphere of Agricultural and Urban Areas in the Province of Buenos Aires in

417

Argentina Using PUF Disk Passive Air Samplers. Atmos. Pollut. Res. 2014, 5 (2),

418

170–178.

419

(5)

Roots, O.; Lukki, T.; Přibylová, P.; Borůvková, J.; Kukučka, P.; Audy, O.; Kalina, J.;

420

Klánová, J.; Holoubek, I.; Sweetman, A.; et al. Measurements of Persistent Organic

421

Pollutants in Estonian Ambient Air (1990–2013). Proc. Est. Acad. Sci. 2015, 64 (2),

422

184–199.

423

(6)

Pozo, K.; Estellano, V. H.; Harner, T.; Diaz-Robles, L.; Cereceda-Balic, F.; Etcharren,

424

P.; Pozo, K.; Vidal, V.; Guerrero, F.; Vergara-Fernández, A. Assessing Polycyclic

425

Aromatic Hydrocarbons (PAHs) Using Passive Air Sampling in the Atmosphere of

426

One of the Most Wood-Smoke-Polluted Cities in Chile: The Case Study of Temuco.

427

Chemosphere 2015, 134, 475–481.

428

(7)

Pozo, K.; Palmeri, M.; Palmeri, V.; Estellano, V. H.; Mulder, M. D.; Efstathiou, C. I.;

429

Sará, G. L.; Romeo, T.; Lammel, G.; Focardi, S. Assessing Persistent Organic

430

Pollutants (POPs) in the Sicily Island Atmosphere, Mediterranean, Using PUF Disk

431

Passive Air Samplers. Environ. Sci. Pollut. Res. 2016, 23 (20), 20796–20804.

432

(8)

Muñoz-Arnanz, J.; Roscales, J. L.; Ros, M.; Vicente, A.; Jiménez, B. Towards the

433

Implementation of the Stockholm Convention in Spain: Five-Year Monitoring (2008–

434

2013) of POPs in Air Based on Passive Sampling. Environ. Pollut. 2016, 217, 107–

435

113.

436

(9)

Estellano, V. H.; Pozo, K.; Přibylová, P.; Klánová, J.; Audy, O.; Focardi, S.

ACS Paragon Plus Environment

25

Environmental Science & Technology

Page 26 of 29

437

Assessment of Seasonal Variations in Persistent Organic Pollutants across the Region

438

of Tuscany Using Passive Air Samplers. Environ. Pollut. 2017, 222, 609–616.

439

(10)

97/101/EC: Council Decision of 27 January 1997 Establishing a Reciprocal Exchange

440

of Information and Data from Networks and Individual Stations Measuring Ambient

441

Air Pollution within the Member States; OPOCE, 1997.

442

(11)

2001/752/EC: Commission Decision of 17 October 2001 Amending the Annexes to

443

Council Decision 97/101/EC Establishing a Reciprocal Exchange of Information and

444

Data from Networks and Individual Stations Measuring Ambient Air Pollution within

445

the Member State; 2001.

446

(12)

Paličková, L.; Kotlík, B.; Sládeček, J.; Vlasáková, L.; Machálek, P.; Modlík, M.;

447

Škáchová, H.; Vlček, O.; Bäumelt, V.; Holubová Šmejkalová, A.; et al. Commentary

448

on

449

http://portal.chmi.cz/files/portal/docs/uoco/isko/tab_roc/2016_enh/pdf/03kom.pdf.

450

(13)

the

summary

annual

tabular

survey

Klánová, J.; Kohoutek, J.; Hamplová, L.; Urbanová, P.; Holoubek, I. Passive Air

451

Sampler as a Tool for Long-Term Air Pollution Monitoring: Part 1. Performance

452

Assessment for Seasonal and Spatial Variations. Environ. Pollut. 2006, 144 (2), 393–

453

405.

454

(14)

Kalina, J.; Scheringer, M.; Borůvková, J.; Kukučka, P.; Přibylová, P.; Bohlin-Nizzetto,

455

P.; Klánová, J. Passive Air Samplers As a Tool for Assessing Long-Term Trends in

456

Atmospheric Concentrations of Semivolatile Organic Compounds. Environ. Sci.

457

Technol. 2017, 51 (12), 7047–7054.

458 459

(15)

Holt, E.; Bohlin-Nizzetto, P.; Borůvková, J.; Harner, T.; Kalina, J.; Melymuk, L.; Klánová, J. Using Long-Term Air Monitoring of Semi-Volatile Organic Compounds

ACS Paragon Plus Environment

26

Page 27 of 29

Environmental Science & Technology

460

to Evaluate the Uncertainty in Polyurethane-Disk Passive Sampler-Derived Air

461

Concentrations. Environ. Pollut. 2017, 220 (Part B), 1100–1111.

462

(16)

Precipitation and Air at Chicago. Environ. Sci. Technol. 2006, 40 (4), 1178–1183.

463 464

Ping, S.; Basu, I.; Hites, R. A. Temporal Trends of Polychlorinated Biphenyls in

(17)

Schuster, J. K.; Gioia, R.; Sweetman, A. J.; Jones, K. C. Temporal Trends and

465

Controlling Factors for Polychlorinated Biphenyls in the UK Atmosphere

466

(1991−2008). Environ. Sci. Technol. 2010, 44 (21), 8068–8074.

467

(18)

Bruckmann, P.; Hiester, E.; Klees, M.; Zetzsch, C. Trends of PCDD/F and PCB

468

Concentrations and Depositions in Ambient Air in Northwestern Germany.

469

Chemosphere 2013, 93 (8), 1471–1478.

470

(19)

Bari, M. A.; Kindzierski, W. B.; Spink, D. Twelve-Year Trends in Ambient

471

Concentrations of Volatile Organic Compounds in a Community of the Alberta Oil

472

Sands Region, Canada. Environ. Int. 2016, 91, 40–50.

473

(20)

Concentration in the Po Valley. Atmos. Chem. Phys. 2014, 14 (10), 4895–4907.

474 475

Bigi, A.; Ghermandi, G. Long-Term Trend and Variability of Atmospheric PM10

(21)

Masiol, M.; Agostinelli, C.; Formenton, G.; Tarabotti, E.; Pavoni, B. Thirteen Years of

476

Air Pollution Hourly Monitoring in a Large City: Potential Sources, Trends, Cycles

477

and Effects of Car-Free Days. Sci. Total Environ. 2014, 494–495, 84–96.

478

(22)

Sanka, O.; Kalina, J.; Lin, Y.; Deutscher, J.; Futter, M.; Butterfield, D.; Melymuk, L.;

479

Brabec, K.; Nizzetto, L. Estimation of p,p’-DDT Degradation in Soil by Modeling and

480

Constraining Hydrological and Biogeochemical Controls. Environ. Pollut. 2018, 239,

481

179–188.

ACS Paragon Plus Environment

27

Environmental Science & Technology

Page 28 of 29

482

(23)

Legendre, P.; Legendre, L. Numerical Ecology; Elsevier, 2012.

483

(24)

Hogarh, J. N.; Seike, N.; Kobara, Y.; Habib, A.; Nam, J.-J.; Lee, J.-S.; Li, Q.; Liu, X.;

484

Li, J.; Zhang, G.; et al. Passive Air Monitoring of PCBs and PCNs across East Asia: A

485

Comprehensive Congener Evaluation for Source Characterization. Chemosphere 2012,

486

86 (7), 718–726.

487

(25)

Liu, X.; Wania, F. Cluster Analysis of Passive Air Sampling Data Based on the

488

Relative Composition of Persistent Organic Pollutants. Environ. Sci. Process. Impacts

489

2014, 16 (3), 453–463.

490

(26)

Wang, X.; Ren, J.; Gong, P.; Wang, C.; Xue, Y.; Yao, T.; Lohmann, R. Spatial

491

Distribution of the Persistent Organic Pollutants across the Tibetan Plateau and Its

492

Linkage with the Climate Systems: A 5-Year Air Monitoring Study. Atmos. Chem.

493

Phys. 2016, 16 (11), 6901–6911.

494

(27)

Pokhrel, B.; Gong, P.; Wang, X.; Gao, S.; Wang, C.; Yao, T. Sources and

495

Environmental Processes of Polycyclic Aromatic Hydrocarbons and Mercury along a

496

Southern Slope of the Central Himalayas, Nepal. Environ. Sci. Pollut. Res. 2016, 23

497

(14), 13843–13852.

498

(28)

Dvorská, A.; Lammel, G.; Holoubek, I. Recent Trends of Persistent Organic Pollutants

499

in Air in Central Europe - Air Monitoring in Combination with Air Mass Trajectory

500

Statistics as a Tool to Study the Effectivity of Regional Chemical Policy. Atmos.

501

Environ. 2009, 43 (6), 1280–1287.

502

(29)

Harner, T. 2017_v1_5_Template for Calculating Effective Air Sample Volumes for

503

PUF

and

SIP

Disk

Samplers_Sept_15.

504

https://www.researchgate.net/publication/319764519_2017_v1_5_Template_for_calcu

ACS Paragon Plus Environment

2017.

28

Page 29 of 29

Environmental Science & Technology

505

lating_Effective_Air_Sample_Volumes_for_PUF_and_SIP_Disk_Samplers_Sept_15

506

(30)

Future, F. the. Worldclim http://worldclim.org/version2.

507

(31)

Global Monitoring Plan of the Stockholm Convention on Persistent Organic

508

Pollutants: visualization and on-line analysis of data from the monitoring reports:

509

Contact http://www.pops-gmp.org/index.php?pg=contact (accessed Jun 14, 2018).

510

ACS Paragon Plus Environment

29