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