Climate Influence on Legacy Organochlorine Pollutants in Arctic

28 Jan 2019 - ... Influence on Legacy Organochlorine Pollutants in Arctic Seabirds. Karen L Foster , Birgit M. Braune , Anthony J. Gaston , and Mark L...
1 downloads 0 Views 569KB Size
Subscriber access provided by UNIV OF LOUISIANA

Environmental Processes

Climate Influence on Legacy Organochlorine Pollutants in Arctic Seabirds Karen L Foster, Birgit M. Braune, Anthony J. Gaston, and Mark L. Mallory Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b07106 • Publication Date (Web): 28 Jan 2019 Downloaded from http://pubs.acs.org on January 29, 2019

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

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

Page 1 of 34

Environmental Science & Technology

Hexachlorobenzene 0.08

ww) -1 (Conc. mg kg

Partial Residuals

0.03

0.06

0.02

0.04

0.01

0.02

0.00

0.00

-0.01

-0.02

-0.02

-0.04

-0.03 -2.0 -1.5 -1.0 -0.5

-0.06 0.0

NAO

0.5

1.0

1.5

20

40

60

80

100

Rain (mm)

ACS Paragon Plus Environment

120

140

Environmental Science & Technology

1

Climate Influence on Legacy Organochlorine

2

Pollutants in Arctic Seabirds

3

Karen L. Fostera,b*, Birgit M. Braunec, Anthony J. Gastonc, Mark L. Malloryd

4 5 6 7 8 9

aKaren

Page 2 of 34

Foster Environmental Research, Peterborough, ON, K9J 8L2, Canada of Modelling & Quantitative Methods (AMOD), Trent University, Peterborough, ON, K9L 0G2, Canada cEnvironment and Climate Change Canada, National Wildlife Research Centre, Carleton University, Ottawa, ON, K1A 0H3, Canada dBiology Department, Acadia University, Wolfville, NS, B4P 2R6, Canada bApplications

10 11

Keywords: climate change, Arctic, seabirds, PCBs, organochlorine pesticides, DDT,

12

chlorobenzenes, mirex, chlordanes

13

ABSTRACT

14

Changing climate can influence the transport of chemical pollutants into Arctic regions and their

15

fate once there. However, the influence of weather or climate variables on organochlorine

16

accumulation in Arctic wildlife, including seabirds, and associated timescale are poorly

17

understood.

18

variables for time lags of 0 to 10 yr and organochlorine pollutant concentrations spanning 1975-

19

2014 in eggs of two seabird species (northern fulmar Fulmarus glacialis, thick-billed murre Uria

20

lomvia) that breed in the Canadian High Arctic. The majority of variability in the data was

21

associated with declining organochlorine emissions (up to 70.2 % for murres and 77.4 % for

We assessed the interannual relationships between a suite of weather/climate

ACS Paragon Plus Environment

1

Page 3 of 34

Environmental Science & Technology

22

fulmars). Controlling for emissions using principal component ordination and general linear

23

modelling, correlations with the North Atlantic Oscillation (NAO) were found for fulmars, and

24

with rainfall for murres, after a time lag of 4-9 yr between weather/climate conditions and egg

25

collection. Our results suggest that with increasingly NAO+ conditions and increasing rainfall,

26

associated with climate change, concentrations of certain organochlorines such as

27

hexachlorobenzene and p,pʹ-DDE increased; dependent on seabird species and ecology, as well

28

as partitioning characteristics of the chemical. Analysis of a truncated version of the datasets

29

(2005-2014), consistent with typical time series lengths for environmental pollutants in Arctic

30

wildlife, found correlations with precipitation for murres, but not with NAO for fulmars,

31

suggesting that longer time series better elucidate relationships with broad-scale climate indices.

32

Organochlorine pollutant datasets spanning 40 years, which is rare for Arctic wildlife, for two

33

species of seabird were assessed and results highlight the association between weather/climate

34

and pollutant accumulation in Arctic food webs and the critical role of ongoing monitoring to

35

effectively elucidate these relationships.

36

ACS Paragon Plus Environment

2

Environmental Science & Technology

37

Page 4 of 34

INTRODUCTION

38

Climate change has been identified as a likely factor governing the fate of pollutants in

39

the environment. This has implications for Arctic ecosystems, which serve as sinks for persistent

40

organic pollutants (POPs), such as legacy organochlorines, predominantly produced and used at

41

lower latitudes but restricted or banned in most circumpolar countries in the 1970s and 1980s (1-

42

3). Widely agreed upon impacts of climate change on pollutants in Arctic ecosystems include

43

the increased volatilization from primary sources (4), increased re-mobilization from

44

environmental reservoirs or secondary sources such as soils, glaciers, sea ice and permafrost

45

(5,6), and changes to air-water exchange owing to reductions in sea ice cover, increased

46

precipitation and increased organic carbon inputs from terrestrial sources and altered carbon

47

cycling in the Arctic Ocean (1,7-9).

48

Ecological implications of a warming climate for Arctic wildlife have been documented

49

(10). For Arctic seabirds, climate-related ecological shifts include effects on population size and

50

predation (11,12), prey availability (13-15), diet (16-18), as well as reproduction and foraging

51

behaviour (19-23). In some cases, ecological shifts have led to changes in accumulated POPs

52

(24,25). A study of organochlorines in thick-billed murre (Uria lomvia) eggs from two colonies

53

found that in the northern Hudson Bay colony, a dietary shift from cold water Arctic cod

54

(Boreogadus saida) to sub-Arctic capelin (Mallotus villosus) resulted in a lower trophic position

55

and a decreased rate of decline of accumulated organochlorine pollutants (26). In the High

56

Arctic colony, an increase in stable nitrogen isotope (δ15N) levels, possibly owing to an increase

57

in fish consumption, led to more dramatic rates of decline of organochlorines but the magnitude

58

of declines were small compared to the effect of reduced emissions.

ACS Paragon Plus Environment

3

Page 5 of 34

Environmental Science & Technology

59

Large-scale climate patterns such as the Arctic Oscillation (AO) and North Atlantic

60

Oscillation (NAO) as well as local and regional weather, which are not always correlated (27),

61

are expected to influence pollutant fate globally, including transfer into and within Arctic regions

62

and accumulation in Arctic wildlife (1-3,7,28-33). AO and NAO are indices of atmospheric

63

circulation with positive (+) cyclonic and negative (–) anticyclonic phases, and they have been

64

found to correlate with POPs concentrations in various Arctic wildlife (e.g., 24,34,35) and in a

65

sub-Arctic terrestrial raptor (36). In seabirds, circulating blood levels of POPs in glaucous gulls

66

(Larus hyperboreus) from the Norwegian Arctic were found to positively correlate with AO

67

from the previous year and negatively correlate with AO from the current year, suggesting

68

increased background levels of POPs under AO+ circulation patterns (37).

69

common eiders (Somateria mollissima) from two colonies in the Norwegian Arctic (38),

70

researchers found higher concentrations of circulating POPs (hexachlorobenzene, PCB-153, and

71

p,pʹ-DDE) in blood when temperatures were low, suggesting a re-mobilization of hydrophobic

72

chemicals stored in fat reserves metabolized to meet the increased energetic demand of the birds

73

in colder conditions.

In a study of

74

Modelling assessments suggest that the maximum response of POPs concentrations in an

75

ecosystem to changes in climate (39) or emissions (40-43) can lag the change itself, with the

76

precise lag time dependent on the physical-chemical properties of the chemical. PCB-153 in an

77

Arctic marine ecosystem, for example, is predicted to have a response time of years between a

78

climate change-induced shift in exposure and the resultant change in a seabird food web (39).

79

Time lags of 4 to 7 yr between maximum PCB loadings to Lake Ontario and maximum PCB

80

accumulation in sediments have been observed (40).

81

declines in PCBs in Lake Ontario food web components have also been observed, including

Similarly, time lagged responses to

ACS Paragon Plus Environment

4

Environmental Science & Technology

Page 6 of 34

82

herring gull (Larus argentatus) eggs, with the estimated response time dependent on diet

83

composition (40). Time lags of years to decades between changes in POPs emissions and

84

resultant changes in accumulation in abiotic environmental compartments have also been

85

predicted for some chemicals in remote Arctic receptor regions (41), oceans (42), and a generic

86

evaluative environment (43). Yet, only zero or one yr time lags have been assessed for Arctic

87

wildlife (34,35,37,38), which may be optimistic given the slow response of the marine ecosystem

88

and food web to shifts in POPs emissions (39,41,42).

89

Prince Leopold Island Migratory Bird Sanctuary (74°02′N, 90°05′W) supports one of the

90

largest seabird colonies in the Canadian High Arctic, including approximately 100,000 breeding

91

pairs of thick-billed murres and 16,000 breeding pairs of northern fulmars (Fulmarus glacialis;

92

44). The ecology and POPs of murres and fulmars at Prince Leopold Island have been studied

93

since 1975. During this time, the diet of murres has remained consistent: 95 % fish (90 % of

94

which is Arctic cod) and 5 % marine invertebrates (16). An increase in the proportion of fish

95

was noted, but not significant. Similarly, the diet of fulmars has also remained fairly consistent:

96

fish (Arctic cod), cephalopods, polychaetes and crustaceans, though recent years suggest a higher

97

proportion of fish (45). Using previously reported POPs concentration data measured over 40

98

years in fulmar and murre eggs (46), and weather and climate data, we assessed potential

99

interannual climate effects on accumulated POPs and the timescale on which effects might occur

100

using principal component (PC) ordination and general linear modelling (GLM).

ACS Paragon Plus Environment

5

Page 7 of 34

101

Environmental Science & Technology

MATERIALS AND METHODS

102

Weather and Climate Data. We compiled an extensive suite of local and regional scale

103

weather variables pertaining to atmospheric and oceanic environmental compartments, and air-

104

to-ocean transfer (e.g., rain/snow), as well as large scale climate indices (e.g., AO, NAO), Arctic

105

sea ice extent, and temperature anomalies (Table 1). Local and regional data were compiled for

106

the Barrow Strait/Lancaster Sound region, or the nearest sampling station to Prince Leopold

107

Island for 1965 through 2014, where available, to enable time lag analyses of zero to 10 yr.

108

However, field- and satellite-based measurements were not always available for every year

109

(Supporting information, Table S1). Depending on the data available, seasonal (winter and

110

summer) and annual averages were compiled, or else data for ~June 20th; a date that reflects the

111

start of egg laying for the majority of murres, and the end of egg laying for fulmars.

112

POP

Concentrations

in

Seabird

Eggs.

Concentrations

of

organochlorine

113

pesticides/metabolites and polychlorinated biphenyls (PCBs) in seabird eggs were measured

114

between 1975 and 2014 and methods of collection and analysis are provided elsewhere (26), as

115

are summarized mean annual concentrations (46). Data for 18 yr were available for murres and

116

17 yr for fulmars. Two to five composite samples comprising equal aliquots of three eggs each

117

were analyzed per species per year, resulting in sample sizes of n = 80 for fulmars and n = 81 for

118

murres. However, all p-values were conservatively calculated for n = 18 for murres and n = 17

119

for fulmars to reflect the number of sampling years (rather than the total number of samples).

120

A subset of organochlorine chemicals was selected for each seabird species and consisted

121

of PCB-31/28, -74, -99, -118, -153, -105, -138, -187, -183, -128, -156 (fulmar only), -180, -170,

122

-199 (murre only), -194 (fulmar only), 1,2,3,4-tetrachlorobenzene (1,2,3,4-TCB; murre only),

ACS Paragon Plus Environment

6

Environmental Science & Technology

Page 8 of 34

123

pentachlorobenzene (PnCB), hexachlorobenzene (HCB), octachlorostyrene (OCS; murre only),

124

heptachlor epoxide (HE; murre only), oxychlordane, trans-nonachlor (fulmar only), cis-

125

nonachlor, p,p’-DDE, dieldrin, p,p’-DDD (fulmar only), p,p’-DDT (fulmar only), photomirex

126

(fulmar only), and mirex (fulmar only); a total of 22 chemicals for murres and 25 for fulmars.

127

Chemical selection was based on detection in > 90 % of samples for the species and

128

concentrations > PQL (practical quantification limit; five times the method detection limit,

129

MDL). Selected chemicals also had declining concentration trends with time. For statistical

130

analyses, non-detectable concentrations (i.e., no signal above baseline) and concentrations
0.8, |r| ≥ 0.6) and the time lags at which correlations were strongest.

147

Pearson correlations were also assessed for truncated versions of the fulmar and murre

148

datasets. Data were truncated to 2005-2014, a period during which concentration variability was

149

reduced, and correlations with weather/climate variables, time lagged as above, were assessed.

150

We deduced a priori that maximum POPs accumulation in seabird eggs might lag climate

151

by a minimum of three yr based on trophic exchange alone. The prey of fulmars and murres

152

breeding on Prince Leopold Island range from zooplankton to Arctic cod (16,45,49), thus POPs

153

likely pass through one to three trophic levels, including zooplankton food webs, which have

154

also been shown to biomagnify POPs (50). Additionally, murres prey on Arctic cod that, based

155

on cod size, are estimated to be more than two years old (49,51). Taken with predicted time lags

156

of years to decades between climate/emission changes and maximum POPs concentration

157

response reported in the literature as previously discussed (39-43), we conservatively opted to

158

assess time lags of zero to 10 yr.

159

Time lagged weather/climate variables identified in the preliminary PC correlation

160

assessments were included in GLM analyses and the most parsimonious models selected using

161

the difference in Akaike's Information Criteria (∆AICC) corrected for smaller sample size

162

following a practical information-theoretic approach (52). Two models were considered for each

163

chemical for each seabird species. The ‘base’ model consisted of Year + Constant. A second

164

model was then assessed: base + NAO (for fulmars) or base + Rain (for murres). Substantial

165

empirical support for both models was indicated when ∆AICC = 0-2, ∆AICC = 4-7 indicated

166

considerably less support for the second best model, and >10 essentially no support for the

ACS Paragon Plus Environment

8

Environmental Science & Technology

Page 10 of 34

167

second best model (52). A second set of exploratory GLMs were run on the most abundant

168

PCBs (i.e., -118, -138, -153, and -180) for both fulmars and murres, substituting global annual

169

emissions (“Emissions”, 53) for Year (i.e., GLMs were: Emissions + NAO + Constant for

170

fulmars, and Emissions + Rain + Constant for murres). Concentrations from this model were

171

compared to those from the Year-based models.

172

To assist in the interpretation of results, three environmentally relevant partition

173

coefficients: octanol-water (Log KOW), air-water (Log KAW), and octanol-air (Log KOA) were

174

compiled using EPISuite™ v.4.11 (54; Supporting information, Table S2). All p-values were

175

two-tailed for α = 0.05. Some weather variables had insufficient non-zero values to include in

176

statistical analyses; these included winter rain accumulation and multiyear sea ice concentrations.

177

Correlations between % lipid and stable nitrogen isotopes (δ15N) in eggs and weather/climate

178

variables at time lags of zero to 10 yr were included in preliminary assessments, but aside from a

179

moderate correlation between δ15N and summer NAO at a lag of 5 yr (r = -0.60) for fulmars no

180

other correlations of acceptable power were found and % lipid and δ15N were not considered

181

further. All statistics were performed in SYSTAT Version 13.

182

ACS Paragon Plus Environment

9

Page 11 of 34

183

Environmental Science & Technology

RESULTS

184

A decline in concentrations of legacy POPs since 1975 in the eggs of murres and fulmars

185

from Prince Leopold Island was evident (r = -0.2 to -0.91, also Figure 1), which corresponded

186

with periods of reduced POPs production. Exploratory GLMs that substituted Emissions for

187

Year yielded PCB-118, -138, -153 and -180 concentrations for both fulmars and murres that

188

were highly correlated with those using Year (r = 0.94 – 0.97, all p 0.8). PC1 was highly correlated

195

with sample collection year (r = -0.87, p < 0.001 for fulmars and r = -0.85, p < 0.001 for murres;

196

Supporting information, Figure S1) and captured the confounding influence of declining

197

emissions with time, thus PC1 was excluded from further analyses. The second extracted

198

principal component (PC2) accounted for 10.8 % and 11.2 % of the variability in the fulmar and

199

murre data, respectively. The third principal component (PC3) captured 5.7 % of the variability

200

in the murre data. For fulmars, PC3 accounted for only 2.4 % of the data variability and, as

201

might be expected, permitted little discrimination between chemical loadings. Thus, PC3 for

202

fulmars was not considered further. Neither PC2 nor PC3 (murres only) were strongly correlated

203

with Year for either species (r = -0.17 and -0.09 for murres, respectively, and r = -0.16 for

204

fulmars), making them suitable variables for further analyses.

ACS Paragon Plus Environment

10

Environmental Science & Technology

Page 12 of 34

205

PC2 loadings for fulmars grouped the chlorobenzenes (PnCB and HCB), chlordanes

206

(oxychlordane, trans-nonachlor, cis-nonachlor), dieldrin, photomirex and mirex in the left

207

quadrants, while PCBs and DDTs (DDE, DDD, DDT) grouped together in the right quadrants

208

(Figure 2a). Preliminary correlation analyses between fulmar PC2 and weather/climate variables

209

identified average annual NAO and local average annual atmospheric pressure as the only

210

comparisons that yielded an acceptable power (Supporting information, Table S3). Correlation

211

between fulmar PC2 and NAO was low (r < 0.06) at time lags of 0-3 yr (Supporting information,

212

Figure S2a), but during lags 4-9 yr they were moderately negatively correlated; the lag of 7 yr

213

provided the highest correlation (r = -0.61, p = 0.0093; Figure 2b). Similarly, correlations of

214

PC2 with atmospheric pressure were low (r < 0.25) at time lags of 0-3 yr, but moderately

215

positively correlated during 4-9 yr lags, with the highest correlation also at a time lag of 7 yr (r =

216

0.63, p = 0.0067; Supporting information, Table S3). Atmospheric pressure and NAO were

217

highly correlated from 1965 to 2014 (r = -0.83, p < 0.001), thus NAO alone was selected for

218

GLM analyses. Winter atmospheric pressure at a lag of 6 yr (r = 0.62, p = 0.0079; Supporting

219

information, Table S3) was also correlated with PC2.

220

PC2 and PC3 loadings for murres distinguished between chlorobenzenes (1,2,3,4-TCB,

221

PnCB, HCB) and OCS in the bottom left quadrant, PCBs and DDE in the right quadrants close to

222

the origin, and chlordanes (HE, oxychlordane, cis-nonachlor) and dieldrin in the upper left

223

quadrant (Figure 3a).

224

variables found that PC3 versus annual average rainfall was the only correlation that yielded an

225

acceptable power (Supporting information, Tables S4 and S5). Correlations between murre PC3

226

and rainfall were poor at time lags of 0-4 yr (all |r| < 0.32; Supporting information, Figure S2b),

Correlation analyses between murre PC2/PC3 and weather/climate

ACS Paragon Plus Environment

11

Page 13 of 34

Environmental Science & Technology

227

but moderately negative at lags of 5-8 yr; the strongest correlation being at a time lag of 6 yr (r =

228

-0.60, p = 0.0085; Figure 3b).

229

Groupings of organochlorine chemicals based on PC loadings for both fulmars and

230

murres (Gaussian bivariate ellipses, p = 0.68; Figure 2a and 3a) reflected intrinsic environmental

231

partitioning characteristics of the chemicals (Supporting information, Table S2). Fulmar PC2

232

loadings were moderately but significantly correlated with Log KOW values (Figure 4a, r = 0.40,

233

p = 0.049), whereas, murre PC3 loadings were highly negatively correlated with Log KAW,

234

stretching from chlordanes and dieldrin (Log KAW = -3.07 to -5.45), to PCBs and DDE (Log KAW

235

= -1.93 to -3.43), to chlorobenzenes and OCS (Log KAW = -1.16 to -2.03); representing a range

236

of four orders of magnitude of KAW (Figure 4b, r = -0.77, p r < 0.26;

240

Supporting information, Figure S3). Thus, average annual NAO at a time lag of 7 yr (“NAO”)

241

and annual rain accumulation at a time lag of 6 yr (“Rain”) were assessed as predictors of

242

individual chemical concentrations in fulmars and murres, respectively, using GLM.

243

General linear modelling (GLM). In fulmars, years of NAO+ conditions were followed

244

by increased concentrations of chlorobenzenes (PnCB, HCB), trans- and cis- nonachlor, dieldrin,

245

photomirex, and mirex in their eggs, as shown by the positive β-coefficients for NAO

246

(Supporting information, Table S6). For these chemicals, the NAO + Year + Constant model

247

was more parsimonious than the base model (∆AICC = 3.4 – 32.4). Partial residual plots

248

confirmed that when controlled for Year, concentrations of chlorobenzenes, trans-/cis-nonachlor,

ACS Paragon Plus Environment

12

Environmental Science & Technology

Page 14 of 34

249

dieldrin, photomirex and mirex concentrations increased with NAO+ (Figure 5, Supporting

250

information, Figure S4 and S5). The GLM for oxychlordane yielded a negative β-coefficient for

251

NAO, however, this was not a robust result – an outlier sample (studentized residuals = 3.33 to

252

8.63) for most chemicals was removed for the purpose of plotting partial residuals and the slope

253

for oxychlordane changed (from negative to positive). The base model was more parsimonious

254

for DDE and most PCBs in fulmar eggs (∆AICC = 1.1 - 2.2), however, partial residual plots of

255

DDTs (DDE, DDD, DDT) and PCBs exhibited a cyclical trend with NAO (Figure 5, Supporting

256

information, Figure S6 and S7).

257

In murres, years of increased rainfall were followed by decreased concentrations of

258

chlordanes (HE, oxychlordane, cis-nonachlor) as well as PCB-170 and -180 in eggs, as indicated

259

by negative β-coefficients for Rain (Supporting information, Table S7). However, only the

260

partial residual plot for oxychlordane yielded a convincingly negative trend (Figure 5).

261

Chlorobenzene (1,2,3,4-TCB, PnCB, HCB), OCS, DDE, dieldrin and most PCB concentrations

262

increased with rainfall and had positive β-coefficients for Rain. Partial residual plots illustrate

263

that when controlled for Year, concentrations of chlorobenzenes, OCS, DDE, dieldrin and most

264

PCBs increased in murre eggs after increased rainfall, while concentrations of oxychlordane

265

decreased (Figure 5, Supporting information, Figure S8 to S10).

266

Across the vast majority of organochlorine chemicals investigated, inclusion of NAO (for

267

fulmars) and Rain (for murres) greatly improved model fit. For fulmars, model fit for 13

268

chemicals indicated that the NAO model was most parsimonious (∆AICC = 0.66 – 32.44) and of

269

these, considerably less support for the base model was found for eight chemicals with ∆AICC ≥

270

4. For the 12 chemicals for which the base model was the better fit, low ∆AICC suggested

271

substantial empirical support for including NAO (∆AICC = 1.13 – 2.20). For murres, the model

ACS Paragon Plus Environment

13

Page 15 of 34

Environmental Science & Technology

272

that included Rain was more parsimonious for 12 organochlorine pollutants (∆AICC = 0.12 –

273

17.47) and of these, considerably less support was seen for the base model for eight chemicals

274

with ∆AICC ≥ 4. For 10 chemicals the base model was the better fit, however, low ∆AICC still

275

suggested substantial empirical support for the Rain model (0.67 ≥ ∆AICC ≥ 2.20).

276

Truncation (2005-2014). Truncated concentrations were, for the most part, not well

277

correlated with Year (-0.57 < r < 0.36), though PCB-170 in murres was moderately correlated (r

278

= -0.67) and not considered further. In murres, PnCB concentrations were positively correlated

279

with annual precipitation (i.e., the sum of rain and snow accumulation) after a 6 yr time lag (r =

280

0.76, p = 0.011), chlorobenzenes (1,2,3,4-TCB, PnCB, HCB) were moderately correlated (r =

281

0.67-0.76) with annual rain, while chlordanes (HE, oxychlordane, cis-nonachlor) and dieldrin

282

were less well correlated (r = 0.12-0.32) with annual rain, and DDE and PCBs were intermediate

283

(r = 0.20-0.68). The power of these correlations was low, nonetheless, these trends generally

284

followed those identified for murres using the full dataset. In fulmars, however, correlations

285

were not found with NAO, nor with other climate indices, with one exception. PCB-156 was

286

both negatively and positively correlated with a number of climate indices at various time lags,

287

but these correlations likely reflect the unusually high percent standard deviation (70 %) of PCB-

288

156.

ACS Paragon Plus Environment

14

Environmental Science & Technology

289

Page 16 of 34

DISCUSSION

290

We have found patterns of correlation between weather and climate and the

291

concentrations of some organochlorine legacy POPs in Arctic seabird eggs using unique, 40 yr

292

datasets for two species. Preliminary correlations between extracted PCs and all weather/climate

293

variables identified NAO (and sea level atmospheric pressure) for fulmars at time lags of 4-9 yr

294

and rainfall for murres at time lags of 5-8 yr as the variables of most relevance to the interannual

295

variability of organochlorine concentrations in eggs. Time lags of 7 and 6 yr, respectively, for

296

fulmars and murres yielded correlations of sufficiently high statistical power.

297

relationships were further examined for each individual chemical using GLMs. Since locally

298

measured sea level atmospheric pressure was well correlated with the broader scale regional

299

NAO for the time period studied (1965-2014), GLM investigations for fulmars focused on only

300

NAO as a main effect.

These

301

Across all organochlorine chemicals investigated, inclusion of NAO (for fulmars) and

302

rainfall (for murres) substantially improved GLM model fit and suggested a strong relationship

303

with organochlorine concentrations in seabird eggs, consistent with the PC correlation analyses.

304

Years of NAO+ conditions were followed by years of higher concentrations of chlorobenzenes,

305

trans-/cis-nonachlor, dieldrin, photomirex and mirex in fulmar eggs, and NAO- was associated

306

with the opposite trends. Concentrations of PCBs and DDTs were found to cycle with maxima

307

at approximately NAO = 0 and minima at NAO = -1 and +1; which is discussed below. Murre

308

GLM results suggested that years of higher rainfall were followed by years when concentrations

309

of chlorobenzenes, OCS, DDE, dieldrin and most PCBs were higher, while concentrations of

310

HE, oxychlordane, cis-nonachlor, PCBs-170 and -180 were lower.

ACS Paragon Plus Environment

15

Page 17 of 34

Environmental Science & Technology

311

Declining concentrations of organochlorine POPs in murre and fulmar eggs have been

312

shown previously, and are attributed to declines in historic production of these chemicals (e.g.,

313

26,46). We included Year in GLM models to control for this effect, in the absence of emissions

314

data for the majority of studied chemicals. However, using emissions data for PCBs that were

315

available (53), our exploratory emissions-based GLM for four PCBs yielded modelled

316

concentrations highly correlated with those from year-based GLMs for both fulmars and murres.

317

We conclude that historic emissions data that become available for other legacy organochlorine

318

chemicals should be included in predictive models of accumulation in Arctic wildlife, with

319

consideration for the time lag necessary for long-range transport to Arctic regions (41), but if

320

such data are not yet available year is a suitable proxy.

321

Time lags of zero up to 10 yr between weather/climate occurrences and seabird egg

322

collections were investigated to assess possible staggered effects of weather and climate on POPs

323

in seabird eggs. Lags of up to 10 yr we viewed as conservative, given that measurement- and

324

modelling-based assessments have suggested that time lags of years to decades may be relevant

325

for some chemicals, particularly with respect to marine ecosystems (39,41,42). Our results are in

326

agreement with these studies and support the consideration of a broader timeframe than currently

327

considered in assessments of climate change effects on POPs in Arctic wildlife (34,35,37,38),

328

particularly given the unknowns associated with climate driven shifts in background

329

concentrations of POPs via primary or secondary sources, and the complex processes that relay

330

these shifts into environmental compartments, into the base of the food web and up to top

331

predators.

332

Clues as to the mechanism for the associations between accumulated organochlorine

333

concentrations in fulmar and murre eggs and NAO and rainfall, respectively, may be gleaned

ACS Paragon Plus Environment

16

Environmental Science & Technology

Page 18 of 34

334

from partitioning properties and inter-chemical comparisons.

For murres, the majority of

335

organochlorine chemicals studied increased in concentration in eggs following years of heavier

336

rainfall, which supports findings that POPs increase in Arctic sea ice melt-ponds following

337

precipitation events (56). PC3 distinguished between chlorobenzenes and OCS, PCBs and DDE,

338

and chlordanes and dieldrin; chemical groupings that correspond with declining Log KAW values,

339

which was shown by the significant decline in PC3 loadings with Log KAW. Since PC3 was

340

found to decline with rainfall, our results suggest that years of heavier rainfall were followed by

341

higher concentrations of higher Log KAW chemicals, such as HCB and PnCB, which have an

342

increased likelihood of being found in air (55) and therefore available for atmospheric

343

scavenging. Oxychlordane, with the lowest Log KAW of the organochlorines investigated here,

344

supported this trend with a concentration decline in years following heavier rainfall, however, it

345

must be noted that oxychlordane is also a major metabolite of technical chlordane in murres (57)

346

and a source apportionment assessment of metabolism, background, and climate associated

347

sources is merited.

348

For fulmars, the concentrations of most organochlorines studied increased in eggs under

349

NAO+ conditions. NAO+ circulation is associated with the transport of POPs originating in

350

industrialized source regions in Eurasia and eastern North America into the Arctic, and is

351

particularly effective in winter (7,30), although, the colder temperatures of winter may hinder the

352

release of POPs from secondary sources (2). Associations between NAO+ and magnified POPs

353

accumulation in Arctic wildlife have been found (34,35,37). For fulmars, the yearly averaged

354

NAO explained more variability in the dataset than winter NAO, possibly owing to the increased

355

contribution of secondary sources during the warmer seasons - higher annual temperatures in

ACS Paragon Plus Environment

17

Page 19 of 34

Environmental Science & Technology

356

Arctic regions have been associated with increased re-mobilization of POPs from secondary

357

sources (5).

358

In our analyses, PC2 loadings for fulmar egg organochlorines distinguished between

359

chlorobenzenes, chlordanes, dieldrin, photomirex, mirex and PCB-187, and the DDTs and

360

remaining PCBs; a grouping that corresponds with Log KOW, and was confirmed by a moderate,

361

but significant, positive correlation between fulmar PC2 loadings and Log KOW. The negative

362

correlation between PC2 and NAO, therefore, indicated that NAO+ conditions favored chemicals

363

such as chlorobenzenes with lower Log KOW. Chlorobenzenes, based on Log KOW, Log KAW and

364

Log KOA, are considered “multiple hoppers;” or chemicals with increased long-range

365

atmospheric transport mobility to Arctic regions (41), and it follows that these would be the

366

chemicals most impacted by increased atmospheric advection into Arctic regions, which occurs

367

under NAO+ conditions (7,30). The negative correlation between PC2 and NAO also indicated

368

that NAO- conditions favored increased concentrations of higher Log KOW chemicals such as the

369

majority of PCBs in fulmar eggs. These higher Log KOW chemicals are mostly “single hoppers”

370

(41) and less readily mobile than the multiple hoppers, which supports the muted influence of

371

NAO+ as shown by GLM partial residual plots.

372

concentrations of higher Log KOW chemicals follow years of neutral NAO. It is probable that

373

given the multiplicity of factors associated with NAO phases, that competing factors are at play.

374

For example, depending on the source region of emissions, NAO+ conditions may also decrease

375

the transport of POPs into the Arctic (2). Future studies utilizing both POPs concentration time

376

series for Arctic wildlife and models capturing competing chemical fate processes under various

377

climate change scenarios are recommended to further elucidate the underlying mechanisms.

However, it is unclear why maximum

ACS Paragon Plus Environment

18

Environmental Science & Technology

Page 20 of 34

378

The association between organochlorine concentrations in the eggs of fulmars and murres

379

and weather/climate variables was particularly interesting considering individual seabird life

380

histories. Murres feed closer to the colony (generally 500 km) (58,59). The

382

importance of the more locally relevant variable of rainfall for murres and of the broader North

383

Atlantic Ocean conditions for fulmars aligns with distinctions in feeding range.

384

recorded locally may be a better indicator of conditions where the murres forage than broad-scale

385

climate indices (e.g., NAO, AO). Also, we recognize that fulmars overwinter in the North

386

Atlantic from the Labrador Sea to the northeast Atlantic (60), and murres in Davis Strait –

387

northern Labrador Sea and south to waters off Newfoundland (61) and there is some evidence

388

that the body burden of POPs in female migrating Arctic seabirds can be influenced by

389

overwintering location, and that influence can be passed on to chicks (62,63). However, we

390

contend that contributions of contaminants accrued overwinter to eggs are probably small

391

compared to contaminants in local food chains for these income breeders (e.g., 64).

Rainfall

392

Interannual assessment of weather/climate variables and pollutants in Arctic regions, for

393

example, contrasting chemical partitioning and fate under + as opposed to – modes of NAO or

394

AO in the Arctic, is thought to provide insight into likely climate-driven shifts over a longer

395

temporal scale (e.g., 7). The results of this study would, therefore, suggest that under a climate

396

change scenario where NAO is increasingly positive, concentrations of chlorobenzenes, trans-

397

/cis-nonachlor, dieldrin, photomirex and mirex are likely to increase in fulmar eggs, at least

398

initially, whereas concentrations of PCBs and DDTs cycle with NAO.

399

associated with climate change (65) is likely to lead to increased concentrations of

Increased rainfall

ACS Paragon Plus Environment

19

Page 21 of 34

Environmental Science & Technology

400

chlorobenzenes, OCS, DDE, dieldrin and PCBs in the eggs of murres, at least initially, and

401

decreased concentrations of oxychlordane.

402

To identify greater than decadal variability in climate cycles, time series data should be

403

longer than 30 yr (66). Climate cycling in the Arctic is thought to occur on a timeline of six

404

years to decades (67). Yet, the typical length of environmental pollutant monitoring time series

405

in Arctic biota is around a decade (1-3), and very few monitoring datasets for POPs in Arctic

406

wildlife exceed the span of 40 yr covered by the datasets used here. Assessments of a truncated

407

version of the murre dataset (2005-2014) found positive correlations between annual

408

precipitation and chlorobenzenes (most notably PnCB) at a time lag of 6 yr, as well as low power

409

correlations with the other organochlorine chemicals, which, in part, reflects the trends identified

410

with the full dataset. For fulmars, however, no correlations with NAO, or indeed with any

411

climate index, were found. Thus, we conclude that a time series of shorter, more typical length

412

for POPs would be too short to effectively detect the trends with weather and climate identified

413

using the full dataset.

414

The single largest barrier to studies of weather and climate influence on environmental

415

pollutants in the Arctic is limited availability of consistent, long-term monitoring data; an issue

416

that has been noted previously (7). Long-term pollutant monitoring datasets for Arctic wildlife

417

are valuable and comparatively rare because of the greater logistical challenges and operational

418

costs for field work than at more southerly latitudes. However, the importance of these data is

419

clear and the optimization of methods to use them effectively to assess the impact of changing

420

climate on pollutants within Arctic ecosystems can only increase their value in the future.

ACS Paragon Plus Environment

20

Environmental Science & Technology

Page 22 of 34

murres fulmars 0.8 0.20 0.16 0.12 0.08 0.04 0.00

PCB

0.6 0.4 0.2 0.0 0.15

Oxychlordane

700 600 500 400 300 200 100 0

Global Emissions (Tonnes)

421

0.12 0.09

Concentration (mg kg-1 ww)

0.06 0.03 0.00 0.15

HCB

0.12 0.09 0.06 0.03 0.00 Dieldrin

0.03 0.02 0.01 0.00

DDE

0.4 0.3 0.2 0.1 0.0 1970

422 423 424 425

1980

1990

2000

2010

Year Figure 1 Concentrations of selected contaminants in fulmar (○) and murre (●) eggs from Prince Leopold Island (mean ± SE; fulmar n = 2-5/yr, murre n = 3-5/yr), with predominant periods of POPs production shaded grey (68-71). Red line shows global Σ4PCB emissions (53). Σ4PCBs is the sum of PCBs-118, -138, -153 and -180.

426

ACS Paragon Plus Environment

21

Page 23 of 34

Environmental Science & Technology

1.0

a)

0.5

0.0

Oxychlordane

b)

r = -0.61, p = 0.0093

r = 0.63, p = 0.0067

4 2

DDD nonachlors (cis/trans) 194 DDE 156 138 105 118 99 180 170 183 31/28 74 128 DDT PMirex Mirex 187 153 Dieldrin HCB

PC2

PC3 Loadings

6

0 -2

-0.5

PnCB

-4 -1.0 -1.0

427 428 429 430

-6 -0.5

0.0

0.5

PC2 Loadings

1.0

-2

-1

0

NAO

1

1012

1014

Pressure (

Figure 2 Northern fulmar PCA results: a) Principal component (PC2 and PC3) loadings for organochlorine chemicals with Gaussian bivariate ellipses for chemical groups shown (p = 0.68), b) PC2 scores with average annual NAO lagged to 7 yr. Chemical groups identified by loadings correspond with Log KOW (Figure 4a).

431

ACS Paragon Plus Environment

22

Environmental Science & Technology

0.5

433 434 435

4

-0.5

b)

2

Dieldrin

0.0

r = -0.60, p = 0.0085

3

cis-nonOxychlordane HE

-1.0 -1.0

432

a)

180 170 138183 153 199 187 105 118 74DDE 99 31/28 OCS 128 HCB 1,2,3,4-TCB PnCB

PC3

PC3 Loadings

1.0

Page 24 of 34

1 0 -1 -2

-0.5

0.0

0.5

PC2 Loadings

1.0

40

60

80

100

Rain (mm)

120

Figure 3 Thick-billed murre PCA results: a) Principal component (PC2 and PC3) loadings for organochlorine chemicals with Gaussian bivariate ellipses for chemical groups shown (p = 0.68), b) PC3 scores with annual rain accumulation lagged to 6 yr. Chemical groups identified by loadings correspond with Log KAW (Figure 4b).

ACS Paragon Plus Environment

23

Page 25 of 34

Environmental Science & Technology

0.4

0.6

Loadings PC3

0.2

Loadings PC2

0.8

a)

r = 0.40, p = 0.049

0.0 -0.2 -0.4 -0.6

0.4 0.2 0.0 -0.2

-0.8

-0.4

-1.0

-0.6 4

5

6

7

Log KOW

8

9

b)

r = -0.77, p < 0.001

-6

-5

-4

-3

-2

-1

0

Log KAW

437 438 439

Figure 4 Principal component loadings versus relevant equilibrium partition coefficients for a) fulmar PC2 loadings versus Log KOW, and b) murre PC3 loadings versus Log KAW.

440

ACS Paragon Plus Environment

24

Oxychlordane

Partial Residuals (Concentration, mg kg-1 ww) PCB 153 HCB

DDE

Environmental Science & Technology

0.2

0.15

0.1

0.10

0.0

0.05

-0.1

0.00

-0.2

-0.05

-0.3 0.03

-0.10 0.08

0.02 0.04

0.01 0.00

0.00

-0.01 -0.02

-0.04

-0.03 0.12

0.03

0.08

0.02

0.04

0.01

0.00 0.00

-0.04 -0.08

-0.01

-0.12 0.09

-0.02 0.015

0.06

0.010

0.03

0.005

0.00

0.000

-0.03

-0.005

-0.06 -2.0 -1.5 -1.0 -0.5

441 442 443 444 445

Page 26 of 34

-0.010 0.0

0.5

1.0

1.5

20

40

60

NAO

80

100

120

140

Rain (mm)

Figure 5 Partial residual plots of the concentrations (mg kg-1 ww) of four chemicals in the eggs of fulmars (○) and murres (●) when year is controlled for versus average annual NAO (lagged 7 yr) and annual accumulation of rain (lagged 6 yr) for fulmars and murres, respectively. One outlier was removed for fulmars for all chemicals (studentized residuals = 3.33 to 8.63).

ACS Paragon Plus Environment

25

Page 27 of 34

Environmental Science & Technology

Table 1 Climate and weather time series data compiled for this study include atmospheric, oceanic, and sea ice variables, as well as air-ocean transfer (i.e., atmospheric scavenging), and teleconnection indices. Data compiled for this study were for 1965-2014 as available. Annual, summer (June, July, August), and winter (December, January, February, March) averages compiled except where noted. Category Air

Variable Surface Air Temperature Wind Speed Temperature Anomalies

Air-Ocean Transfer

Sea Level Pressure Precipitation

Unit °C km h-1 °C

hPa mm

Ocean

Sea Surface Temperature (SST)

°C

Sea Ice

Sea Ice Conc.

%

Sea Ice Conc.

%

Teleconnections

Arctic Sea Ice Extent North Atlantic Oscillation (NAO) Arctic Oscillation (AO)

Pacific North America index (PNA) Multivariate El Niño-Southern Oscillation (MEI)

106 km2

Data Type and Treatment Seasonal means computed from monthly means Seasonal means computed from monthly means Temperature departures from long-term average of reference period (1901-2000) for June Seasonal means computed from monthly means Total precipitation, snow and rain accumulation computed as the sum of monthly totals

Region Resolute

Source Environment and Climate Change Canada, http://ec.gc.ca/dccha-ahccd/

Resolute

Environment and Climate Change Canada, http://ec.gc.ca/dccha-ahccd/

Northern Hemisphere

https://www.ncdc.noaa.gov/, (72)

Resolute

Environment and Climate Change Canada, http://ec.gc.ca/dccha-ahccd/

Resolute

Environment and Climate Change Canada, http://ec.gc.ca/dccha-ahccd/

Computed from satellite imagery; yearly 15 day composite July 20th-Aug 3rd (ice cover in some years prevented SST measurements June 20th) Total, multiyear, and first year ice concentrations computed from ice charts; used data for the week of June 18th Satellite brightness measurements; average June values for six regions

Lancaster Sound

St. Lawrence Global Observatory, http://slgo.ca/en/remotesensing/data.html

Eastern Parry: Lancaster to Eastern Barrow region Six marine regions: all data within 25 km of colonies Northern Hemisphere North Atlantic

Canadian Ice Service Digital http://iceweb1.cis.ec.gc.ca/IceGraph/

Summer sea ice extent in June Principal component based index Wintertime sea level pressure (SLP) based northern annular mode (NAM) index; seasonal means computed Historical index standardized to 1981 to 2010 climatology; means computed Bimonthly values; seasonal and yearly means computed

Arctic

Western Canada associated with Asian jet stream Tropical Pacific

Archive,

IceGraph

2.5,

National Snow & Ice Data Center (NSIDC), Polaris, http://nsidc.org/data/polaris/, (73)

http://nsidc.org/data/seaice_index/archives.html; (74) Climate Analysis Section, National Center for Atmospheric Research (NCAR), https://climatedataguide.ucar.edu/climate-data, (75) Climate Prediction Center, National Oceanic and Atmospheric Administration (NOAA), http://www.cpc.ncep.noaa.gov/products/precip/CWlink/daily_ao_index/ao.shtml, (76) Climate Prediction Center, National Oceanic and Atmospheric Administration (NOAA), http://www.cpc.ncep.noaa.gov/data/teledoc/pna.shtml https://www.esrl.noaa.gov/psd/enso/mei/, (77)

ACS Paragon Plus Environment

26

Environmental Science & Technology

Supporting Information.

Page 28 of 34

POPs concentration data are entered into the Environment and

Climate Change Canada Data Catalogue and will be available through the Government of Canada’s Open Government site (open.canada.ca/en) in the near future. Additionally, data are available from the authors upon request. Tables: data availability for POPs and weather/climate variables, POPs partition coefficients, preliminary PC correlations, GLM results. Figures: PC1 versus sample collection year, r values at time lags 0 to 10 yr, climate/weather variable values (1965-2014), GLM partial residual plots, comparison of Year- and Emission-based GLM models. This material is available free of charge via the Internet at http://pubs.acs.org. Corresponding Author *[email protected]; (705) 768-2081 ACKNOWLEDGMENT The authors thank D. Mackay (Trent University), S. Findlay (University of Ottawa) and R. Macdonald (Department of Fisheries and Oceans, DFO, Emeritus) for insightful suggestions early in the project and anonymous reviewers for constructive comments and suggestions. We also thank those individuals maintaining the climate variable databases and public access websites, and the field assistants over the years who endured often challenging conditions to acquire eggs for analyses. Financial support provided by the Northern Contaminants Program (NCP) of Indigenous and Northern Affairs Canada, Environment and Climate Change Canada, and Acadia University.

ACS Paragon Plus Environment

27

Page 29 of 34

Environmental Science & Technology

REFERENCES 1) AMAP. Influence of Climate Change on Transport, Levels, and Effects of Contaminants in Northern Areas – Part 2. By: P. Carlsson, J.H. Christensen, K. Borgå, R. Kallenborn, K. Aspmo Pfaffhuber, J.Ø. Odland, L.-O. Reiersen, and J.F. Pawlak. Arctic Monitoring and Assessment Programme (AMAP), Oslo, 2016, 52 pp. 2) NCP. Canadian Arctic Contaminants Assessment Report On Persistent Organic Pollutants. Muir D, Kurt-Karakus P, Stow J (Eds). Northern Contaminants Program, Aboriginal Affairs and Northern Development, Ottawa, 2013, pp. 487 3) UNEP/AMAP Climate change and POPs: Predicting the impacts. Report of the UNEP/AMAP expert group. Secretariat of the Stockholm Convention, Geneva, 2011, pp. 62. 4) Zhao, S.; Breivik, K.; Liu, G.; Zheng, M.; Jones, K.C.; Sweetman, A.J. Long-term temporal trends of polychlorinated biphenyls and their controlling sources in China. Environ Sci Technol 2017, 51, 2838-2845. 5) Ma, J.; Hung, H.; Tian, C.; Kallenborn, R. Revolatilization of persistent organic pollutants in the Arctic induced by climate change. Nat Clim Chang 2011, 1, 256–260. 6) Nizzetto, L.; Macleod, M.; Borgå, K.; Cabrerizo, A.; Dachs, J.; Di Guardo, A.; Ghirardello, D.; Hansen, K.M.; Jarvis, A.; Lindroth, A.; Ludwig, B.; Monteith, D.; Perlinger, J.A.; Scheringer, M.; Schwendenmann, L.; Semple, K.T.; Wick, L.Y.; Zhang, G.; Jones, K.C. Past present and future controls on levels of persistent organic pollutants in the global environment. Environ Sci Technol 2010, 44, 6526-6531. 7) Macdonald, R.W.; Harner, T.; Fyfe, J. Recent climate change in the Arctic and its impact on contaminant pathways and the interpretation of temporal trend data. Sci Total Environ 2005, 342, 5-86. 8) Armitage, J.M.; Wania, F. Exploring the potential influence of climate change and particulate organic carbon scenarios on the fate of neutral organic contaminants in the Arctic environment. Environ Sci, Processes Impacts 2013, 15, 2263-2272. 9) Macdonald, R.W.; Kuzyk, Z.A.; Johannessen, S.C. It is not just about the ice: a geochemical perspective on the changing Arctic Ocean. J Environ Stud Sci 2015, 5, 288-301. 10) Descamps, S.; Aars, J; Fuglei, E; Kovacs, K.M.; Lydersen, C.; Pavlova, O.; Pedersen, Å.Ø.; Ravolainen, V.; Strøm, H. Climate change impacts on wildlife in a High Arctic archipelago – Svalbard, Norway. Glob Chang Biol 2017, 23, 490-502. 11) Irons, D.B.; Anker-Nilssen, T.; Gaston, A.J.; Byrd, G.V.; Falk, K.; Gilchrist, G.; Hario, M.; Hjernquist, M.; Krasnov, Y.V.; Mosbech, A.; Olsen, B.; Petersen, A.; Reid, J.B.; Robertson, G.J.; Strøm, H.; Wohl, K.D. Fluctuations in circumpolar seabird populations linked to climate oscillations. Glob Chang Biol 2008, 14, 1455-1463. 12) Prop, J.; Aars, J.; Bårdsen, B.-J.; Hanssen, S.A.; Bech, C.; Bourgeon, S.; de Fouw, J.; Gabrielsen, G.W.; Lang, J.; Noreen, E.; Oudman, T.; Sittler, B.; Stempniewicz, L.; Tombre, I.; Wolters, E.; Moe, B. Climate change and the increasing impact of polar bears on bird populations. Front Ecol Evol 2015, 3, 33. 13) Kortsch, S.; Primicerio, R.; Beuchel, F.; Renaud, P.E.; Rodrigues, J.; Lønne, O.J.; Gulliksena, B. Climate-driven regime shifts in Arctic marine benthos. Proc Nat Acad Sci 2012, 109, 14052-14057. 14) Kortsch, S.; Primicerio, R.; Fossheim, M.; Dolgov, A.V.; Aschan, M. Climate change alters the structure of arctic marine food webs due to poleward shifts of boreal generalists. Proc R Soc B 2015, 282, 31-39.

ACS Paragon Plus Environment

28

Environmental Science & Technology

Page 30 of 34

15) Steiner, N.; Azetsu-Scott, K.; Hamilton, J.; Hedges, K.; Hu, X.; Janjua, M.Y.; Lavoie, D.; Loder, J.; Melling, H.; Merzouk, A.; Perrie, W.; Peterson, I.; Scarratt, M.; Sou, T.; Tallmann, R. Observed trends and climate projections affecting marine ecosystems in the Canadian Arctic. Environ Rev 2015, 23, 191-239. 16) Provencher, J.F.; Gaston, A.J.; O’Hara, P.D.; Gilchrist, H.G. Seabird diet indicates changing Arctic marine communities in eastern Canada. Mar Ecol Prog Ser 2012, 454, 171-182. 17) Regular, P.M.; Shuhood, F.; Power, T.; Montevecchi, W.A.; Robertson, G.J.; Ballam, D.; Piatt, J.F.; Nakashima, B. Murres capelin and ocean climate: inter-annual associations across a decadal shift. Environ Monit Assess 2009, 156, 293-302. 18) Smith, P.A.; Gaston A.J. Environmental variation and the demography and diet of thickbilled murres. Mar Ecol Prog Ser 2012, 454, 237-249. 19) Gaston, A.J.; Gilchrist, H.G.; Mallory, M.L. Variation in ice conditions has strong effects on the breeding of marine birds at Prince Leopold Island Nunavut. Ecography 2005, 28, 331344. 20) Grémillet, D.; Welcker, J.; Karnovsky, N.J.; Walkusz, W.; Hall, M.E.; Fort, J.; Brown, Z.W.; Speakman, J.R.; Harding, A.M.A. Little auks buffer the impact of current Arctic climate change. Mar Ecol Prog Ser 2012, 454, 197-206. 21) Karnovsky, N.; Harding, A.; Walkusz, W.; Kwasniewski, S.; Goszczko, I.; Wiktor, Jr. J.; Routti, H; Bailey, A.; McFadden, L.; Brown, Z.; Beaugrand, G.; Grémillet, D. Foraging distributions of little auks Alle alle across the Greenland Sea, implications of present and future Arctic climate change. Mar Ecol Prog Ser 2010, 415, 283-293. 22) Moe, B.; Stempniewicz, L.; Jakubas, D.; Angelier, F.; Chastel, O.; Dinessen, F.; Gabrielsen, G.W.; Hanssen, F.; Karnovsky, N.J.; Rønning, B.; Welcker, J.; Wojczulanis-Jakubas, K.; Bech, C. Climate change and phenological responses of two seabird species breeding in the high-Arctic. Mar Ecol Prog Ser 2009, 393, 235-246. 23) Shultz, M.T.; Piatt, J.F.; Harding, A.M.A.; Kettle, A.B.; van Pelt, T. Timing of breeding and reproductive performance in murres and kittiwakes reflect mismatched seasonal prey dynamics. Mar Ecol Prog Ser 2009, 393, 247-258. 24) McKinney, M.A.; Pedro, S.; Dietz, R.; Sonne, C.; Fisk, A.T.; Roy, D. Jenssen, B.M.; Letcher, R.J. A review of ecological impacts of global climate change on persistent organic pollutant and mercury pathways and exposures in arctic marine ecosystems. Curr Zool 2015, 61, 617-628. 25) McKinney, M.A.; McMeans, B.C.; Tomy, G.T.; Rosenberg, B.; Ferguson, S.H.; Morris, A.; Muir, D.C.G.; Fisk, A.T. Trophic transfer of contaminants in a changing Arctic marine food web: Cumberland Sound, Nunavut, Canada. Environ Sci Technol 2012, 46, 9914-9922. 26) Braune, B.M.; Gaston, A.J.; Hobson, K.A.; Gilchrist, H.G.; Mallory, M.L. Changes in trophic position affect rates of contaminant decline at two seabird colonies in the Canadian Arctic. Ecotox Environ Safe 2015, 115, 7-13. 27) Stenseth, N.C.; Ottersen, G.; Hurrell, J.W.; Mysterud, A.; Lima, M.; Chan, K.S.; Yoccoz, N.G.; Adlandsvik, B. Studying climate effects on ecology through the use of climate indices, the North Atlantic Oscillation El Nino Southern Oscillation and beyond. Proc R Soc Lond B 2003, 270, 2087-96. 28) Armitage, J.M.; Quinn, C.L.; Wania, F. Global climate change and contaminants-an overview of opportunities and priorities for modelling the potential implications for long-term human exposure to organic compounds in the Arctic. J Environ Monit 2011, 13, 1532-46.

ACS Paragon Plus Environment

29

Page 31 of 34

Environmental Science & Technology

29) Kallenborn, R.; Halsall, C.; Dellong, M.; Carlsson, P. The influence of climate change on the global distribution and fate processes of anthropogenic persistent organic pollutants. J Environ Monit 2012, 4, 2854–2869. 30) Ma, J.; Hung, H.; Macdonald, R.W. The influence of global climate change on the environmental fate of persistent organic pollutants, A review with emphasis on the Northern Hemisphere and the Arctic as a receptor. Glob Planet Change 2016, 146, 89-108. 31) Nadal, M.; Marquès, M.; Mari, M.; Domingo, J. Climate change and environmental concentrations of POPs, A review. Environ Res 143,177-185. 32) Wang X, Sun D, Yao T (2016) Climate change and global cycling of persistent organic pollutants, A critical review. Sci China Earth Sci 2015, 59, 1899-1911. 33) Alava, J.J.; Cisneros-Montemayor, A.M.; Sumaila, U.R.; Cheung, W.W.L. Projected amplification of food web bioaccumulation of MeHg and PCBs under climate change in the Northeastern Pacific. Sci Rep 2018, 8, 13460-13472. 34) Cabrerizo, A.; Muir, D.C.G.; Köck, G.; Iqaluk, D.; Wang, X. Climatic influence on temporal trends of polychlorinated biphenyls and organochlorine pesticides in landlocked char from lakes in the Canadian High Arctic. Environ Sci Technol 2018, 52, 10380-10390. 35) Rigét, F.; Vorkamp, K.; Hobson, K.A.; Muir, D.C.; Dietz, R. Temporal trends of selected POPs and the potential influence of climate variability in a Greenland ringed seal population. Environ Sci Process Impacts 2013, 15, 1706-1716. 36) Bustnes, J.O.; Yoccoz, N.G.; Bangjord, G.; Herzke, D.; Ahrens, L.; Skaare, J.U. Impacts of climate and feeding conditions on the annual accumulation (1986−2009) of persistent organic pollutants in a terrestrial raptor. Environ. Sci. Technol. 2011, 45, 7542−7547. 37) Bustnes, J.O.; Gabrielsen, G.W.; Verreault, J. Climate variability and temporal trends of persistent organic pollutants in the Arctic, A study of glaucous gulls. Environ Sci Technol 2010, 44, 3155–3161. 38) Bustnes, J.O.; Moe, B.; Hanssen, S.A.; Herzke, D.; Fenstad, A.A.; Nordstad, T.; Borgå, K.; Gabrielsen, G.W. Temporal dynamics of circulating persistent organic pollutants in a fasting seabird under different environmental conditions. Environ Sci Technol 2012, 46, 10287–10294. 39) Borgå, K.; Saloranta, T.M.; Ruus, A. Simulating climate change-induced alterations in bioaccumulation of organic contaminants in an Arctic marine food web. Environ Toxicol Chem 2010, 29, 1349-1357. 40) Gobas, F.A.P.C.; Z'graggen, M.N.; Zhang, X. Time response of the Lake Ontario ecosystem to virtual elimination of PCBs. Environ Sci Technol 1995, 29, 2038-2046. 41) Gouin, T.; Wania, F. Time trends of Arctic contamination in relation to emission history and chemical persistence and partitioning properties. Environ Sci Technol 2007, 41, 5986-5992. 42) Stemmler, I.; Lammel, G. Evidence of the return of past pollution in the ocean, A model study. Geophys Res Lett 2013, 40, 1373–1378. 43) Choi, S.D.; Wania, F. On the Reversibility of Environmental contamination with persistent organic pollutants. Environ Sci Technol 2011, 45, 8834-8841. 44) Gaston, A.J.; Mallory, M.L.; Gilchrist, H.G. Populations and trends of Canadian Arctic seabirds. Polar Biol 2012, 35, 1221-1232. 45) Mallory, M.L.; Karnovsky, N.J.; Gaston, A.J.; Hobson, K.A.; Provencher, J.F.; Forbes, M.R.; Hunt, G.L.; Byers, T.; Dick, T.A. Temporal and spatial patterns in the diet of northern fulmars Fulmarus glacialis in the Canadian High Arctic. Aquat Biol 2010, 10, 181-191.

ACS Paragon Plus Environment

30

Environmental Science & Technology

Page 32 of 34

46) Braune, B.M.; Gaston, A.J.; Mallory, M.L. Temporal trends of legacy organochlorines in eggs of Canadian Arctic seabirds monitored over four decades. Sci Total Environ 2019, 646, 551-563. 47) U.S. Environmental Protection Agency. Guidance for data quality assessment: Practical methods for data analysis. 2000. EPA/600/R-96/084. Office of Environmental Information,Washington, DC. 48) Helsel, D. Much ado about next to nothing: Incorporating nondetects in science. Ann Occup Hyg 2010, 54, 257–262. 49) Elliott, K.H.; Gaston, A.J. Mass-length relationships and energy content of fishes and invertebrates delivered to nestling thick-billed murres Uria lomvia in the Canadian Arctic 19812007. Mar Ornithol 2008, 36, 25-34. 50) Hallanger, I.G.; Ruus, A.; Warner, N.A.; Herzke, D.; Evenset, A.; Schøyen, M.; Gabrielsen, G.W.; Borgå, K. Differences between Arctic and Atlantic fjord systems on bioaccumulation of persistent organic pollutants in zooplankton from Svalbard. Sci Total Environ 2011, 409, 2783-2795. 51) Lønne, O.J.; Gulliksen, B. Size, age and diet of polar cod, Boreogadus saida (Lepechin 1773), in ice covered waters. Polar Biol 1989, 9, 187-191. 52) Burnham, K.P.; Anderson, D.R. Model selection and multimodel inference, A practical information-theoretic approach, 2nd ed.; Springer Science+Business Media LLC, New York, NY, 2002. 53) Breivik, K.; Armitage, J.M.; Wania, F.; Sweetman, A.J.; Jones, K.C. Tracking the global distribution of persistent organic pollutants accounting for e-waste exports to developing regions. Environ Sci Technol 2016, 50, 798-805. 54) US EPA Estimation Programs Interface Suite™ for Microsoft® Windows, v 4.11. United States Environmental Protection Agency, Washington, DC, 2017. 55) Mackay, D. Multimedia Environmental Models, The Fugacity Approach, 2nd ed.; CRC Press LLC, Boca Raton, FL, USA, 2001. 56) Pućko, M.; Stern, G.A.; Burt, A.E.; Jantunen, L.M.; Bidleman, T.F.; Macdonald, R.W.; Barber, D.G.; Geilfus, N.X.; Rysgaard, S. Current use pesticide and legacy organochlorine pesticide dynamics at the ocean-sea ice-atmosphere interface in resolute passage, Canadian Arctic, during winter-summer transition. Sci Total Environ 2017, 580, 1460-1469. 57) Fisk, A.T.; Moisey, J.; Hobson, K.A.; Karnovosky, N.J.; Norstrom, N.J. Chlordane components and metabolites in seven species of Arctic seabirds from the Northwater Polynya: relationships with stable isotopes of nitrogen and enantiomeric fractions of chiral compounds. Environ Pollut 2001, 113, 225-238. 58) Gaston, A.J.; Elliott, K.H.; Ropert-Coudert, Y.; Kato, A.; Macdonald, C.A.; Mallory, M.L.; Gilchrist, H.G. Modeling foraging range for breeding colonies of thick-billed murres Uria lomvia in the Eastern Canadian Arctic and potential overlap with industrial development. Biol Conserv 2013, 168, 134-143. 59) Mallory, M.L.; Braune, B.M. Tracking contaminants in seabirds of Arctic Canada, temporal and spatial insights. Mar Pollut Bull 2012, 64, 1475-1484. 60) Mallory, M.L.; Akearok, J.A.; Edwards, D.B.; O'Donovan, K.; Gilbert, C.D. Autumn migration and wintering of northern fulmars (Fulmarus glacialis) from the Canadian high Arctic. Polar Biol 2008, 31, 745–750. 61) McFarlane Tranquilla, L.; Montevecchi, W.A.; Hedd, A.; Regular, P.M.; Robertson, G.J.; Fifield, D.A.; Devillers, R. Ecological segregation among thick-billed murres (Uria lomvia) and

ACS Paragon Plus Environment

31

Page 33 of 34

Environmental Science & Technology

common murres (Uria aalge) in the Northwest Atlantic persists through the nonbreeding season. Can J Zool 2015, 93, 447-460. 62) Baert, J.M.; Janssen, C.R.; Borgå, K.; De Laender, F. Migration and opportunistic feeding increase PCB accumulation in Arctic seabirds. Environ Sci Technol 2013, 47, 1179311801. 63) Bourgeon, S.; Leat, E.K.H.; Furness, R.W.; Borgå, K.; Hanssen, S.A.; Bustnes, J.O. Dietary versus maternal sources of organochlorines in top predator seabird chicks, An experimental approach. Environ Sci Technol 2013, 47, 5963-5970. 64) Mallory, M.L.; Forbes, M.R.; Ankney, C.D.; Alisauskas, R.T. Nutrient dynamics and constraints on the pre-laying exodus of high Arctic northern fulmars. Aquatic Biology 2008, 211223. 65) Bintanja, R.; Andry, O. Towards a rain-dominated Arctic. Nat Clim Change 2017, 7, 263267. 66) Le Treut, H.R.; Somerville, U.; Cubasch, Y.; Ding, C.; Mauritzen, A.; Mokssit, T.; Peterson; Prather, M. Historical Overview of Climate Change In, Climate Change 2007, The Physical Science Basis Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change [Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL (eds)] Cambridge University Press Cambridge United Kingdom and New York NY, 2007. 67) Venegas, S.A.; Mysak, L.A. Is there a dominant timescale of natural climate variability in the Arctic? J Climate 2000, 13, 3412-3434. 68) AMAP AMAP assessment 2002, Persistent organic pollutants in the Arctic (Arctic Monitoring and Assessment Program (AMAP), Oslo, 2004. Xvi+310 pp. 69) Barber, J.L.; Sweetman, A.J.; van Wijk, D.; Jones, K.C. Hexachlorobenzene in the global environment, emissions levels distribution trends and processes. Sci Total Environ 2005, 349, 144. 70) Barrie, L.A.; Gregor, D.; Hargrave, B.; Lake, R.; Muir, D.; Shearer, R.; Tracey, B.; Bidleman, T. Arctic contaminants, sources occurrence and pathways. Sci Total Environ 1992, 122, 1-74. 71) Li, Y.F.; Macdonald, R.W. Sources and pathways of selected organochlorine pesticides to the Arctic and the effect of pathway divergence on HCH trends in biota, a review. Sci Total Environ 2005, 342, 87-106. 72) Lawrimore, J.H.; Menne, M.J.; Gleason, B.E.; Williams, C.N.; Wuertz, D.B.; Vose, R.S.; Rennie, J. Global Historical Climatology Network - Monthly (GHCN-M) Version 330 NOAA National Centers for Environmental Information, 2011. 73) Cavalieri, D.J.; Parkinson, C.L.; Gloersen, P.; Zwally, H. Sea Ice Concentrations from Nimbus-7 SMMR and DMSP SSM/I-SSMIS Passive Microwave Data Boulder Colorado USA, NASA National Snow and Ice Data Center Distributed Active Archive Center, 1996. http://dxdoiorg/105067/8GQ8LZQVL0VL. 74) Fetterer, F.; Knowles, K.; Meier, W.; Savoie, M. Sea Ice Index Version 2 Boulder Colorado USA NSIDC, National Snow and Ice Data Center doi, 2016. http://dxdoiorg/107265/N5736NV7 Accessed 20/02/ 2017. 75) Hurrell, J. NAO Index Data provided by the Climate Analysis Section National Center for Atmospheric Research (NCAR) Boulder USA, 2003 https://climatedataguideucaredu/climatedata. Accessed 03/09/2015.

ACS Paragon Plus Environment

32

Environmental Science & Technology

Page 34 of 34

76) Higgins, R.W. AO Index Data provided by the Climate Prediction Center National Oceanic and Atmospheric Administration (NOAA) Camp Springs USA, 2002. http://wwwcpcncepnoaagov/products/precip/CWlink/daily_ao_index/aoshtml. Accessed 03/09/2015. 77) Wolter, K.; Timlin, M.S. Monitoring ENSO in COADS with a seasonally adjusted principal component index. Proc of the 17th Climate Diagnostics Workshop Norman OK NOAA/NMC/CAC NSSL Oklahoma Clim Survey CIMMS and the School of Meteor Univ of Oklahoma, 1993, pp. 52-57.

ACS Paragon Plus Environment

33