Emission Changes Dwarf the Influence of Feeding Habits on Temporal

Sep 18, 2017 - Tools & Sharing. Add to Favorites · Download Citation · Email a Colleague · Order Reprints · Rights & Permissions · Citation Alerts · A...
2 downloads 19 Views 974KB Size
Subscriber access provided by Binghamton University | Libraries

Article

Emission changes dwarf the influence of feeding habits on temporal trends of per- and polyfluoroalkyl substances in two Arctic top predators Heli Routti, Jon Aars, Eva Fuglei, Linda Hanssen, Karen Lone, Anuschka Polder, Åshild Ø. Pedersen, Sabrina Tartu, Jeffrey M. Welker, and Nigel G. Yoccoz Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b03585 • Publication Date (Web): 18 Sep 2017 Downloaded from http://pubs.acs.org on September 19, 2017

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

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

Page 1 of 43

Environmental Science & Technology

1

Emission changes dwarf the influence of feeding

2

habits on temporal trends of per- and polyfluoroalkyl

3

substances in two Arctic top predators

4

Heli Routti,†* Jon Aars,† Eva Fuglei,† Linda Hanssen,‡ Karen Lone,† Anuschka Polder,§

5

Åshild Ø. Pedersen,† Sabrina Tartu,† Jeffrey M. Welker,∥ and Nigel G. Yoccoz,⊥

6

† Norwegian Polar Institute, Fram Centre, Tromsø, Norway

7

‡ Norwegian Institute for Air Research, Fram Centre, Tromsø, Norway

8

§ Norwegian University of Life Science, Campus Adamstua, Oslo, Norway

9

∥ University of Alaska Anchorage, Department of Biological Sciences, Anchorage, AK, USA

10

⊥ UiT-The Arctic University of Norway, Department of Arctic and Marine Biology, Tromsø,

11

Norway

12

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 43

13

ABSTRACT

14

We monitored concentrations of per- and polyfluoroalkyl substances (PFASs) in relation to

15

climate-associated changes in feeding habits and food availability in polar bears (Ursus

16

maritimus) and arctic foxes (Vulpes lagopus) (192 plasma and 113 liver samples, respectively)

17

sampled from Svalbard, Norway, during 1997-2014. PFASs concentrations became greater with

18

increasing dietary trophic level, as bears and foxes consumed more marine as opposed to

19

terrestrial food, and as the availability of sea ice habitat increased. Long-chained perfluoroalkyl

20

carboxylates (PFCAs) in arctic foxes decreased with availability of reindeer carcasses. The ~9-

21

14% yearly decline of C6-8 perfluoroalkyl sulfonates (PFSAs) following the cease in C6-8 PFSA

22

precursor production in 2001 indicates that the peak exposure was mainly a result of atmospheric

23

transport of the volatile precursors. However, the stable PFSA concentrations since 2009-2010

24

suggest that Svalbard biota is still exposed to ocean-transported PFSAs. Long-chain ocean-

25

transported PFCAs increased 2-4% per year and the increase in C12-14 PFCAs in polar bears

26

tended to level off since ~2009. Emerging short-chain PFASs showed no temporal changes.

27

Climate-related changes in feeding habits and food availability moderately affected PFAS trends.

28

Our results indicate that PFAS concentrations in polar bears and arctic foxes are mainly affected

29

by emissions.

30

ACS Paragon Plus Environment

2

Page 3 of 43

Environmental Science & Technology

31

INTRODUCTION

32

Per- and polyfluoroalkyl substances (PFASs) have been produced at large quantities for a wide

33

variety of purposes.1,2 In particular perfluoroalkyl carboxylates (PFCAs) and perfluoroalkyl

34

sulfonates (PFSAs) (and their precursors) are recognized as global pollutants of emerging

35

concern because they are persistent, bioaccumulative and toxic.3-5 PFASs have complex emission

36

history.6,7 Production of perfluorooctane sulfonyl fluoride (POSF) based products, including

37

perfluorooctane sulfonate (PFOS), grew rapidly in Europe and North-America from the 1970s

38

until the major producer, 3M, voluntarily phased-out the production in 2001-2002.8,9 After the

39

production ceased in Western countries, China started to produce POSF-based products in the

40

early 2000s.6 However, the estimated yearly production volumes in China have been less than a

41

tenth compared to the yearly production volumes in the 1990s.6 China has exported large

42

quantities of POSF to Brazil where it is used to synthetize N-ethyl perfluorooctane sulfonamide

43

(N-EtFOSA), used as a pesticide.10 An amendment of the Stockholm Convention listing PFOS,

44

its salts and POSF was ratified by most countries in 2010 and, in 2014 by China (chm.pops.int).

45

They are included in Annex B of the Convention, which means that they are restricted, but not

46

completely banned.11

47

Large quantities of C9-13 PFCAs were emitted from Western countries and Japan until 2002,

48

but since 2003, the emissions have increased from continental Asia.7 Furthermore, the release of

49

fluorotelomer-based precursors, including fluorotelomer alcohols (FTOHs) and –sulfonates

50

(FTSAs), increasingly contributes to global PFCA emissions.7

51

(PFOA), its salts and PFOA-related compounds are proposed for listing under the Stockholm

52

Convention.12 C8-14 PFCAs including their salts and precursors are currently proposed or

53

intended for restrictions by European Union (www.echa.europa.eu).

Today, perfluorooctanoate

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 43

54

PFSAs and PFCAs, that are released directly into the aquatic system, may slowly reach the

55

Arctic via ocean currents,13-15 whereas, their volatile precursors, for example FTOHs and

56

perfluoroalkane sulfonamides (FASAs), undergo rapid transport via air currents. These

57

precursors can be degraded to PFSAs and PFCAs in the atmosphere.16-19 The transport of PFCAs

58

via ocean currents is expected to be major compared to atmospheric transport13,20 whereas the

59

contribution of oceanic vs atmospheric transport of PFSAs and PFCAs may vary locally.21 The

60

European Arctic, including Svalbard, is exposed to high input of both oceanic and atmospheric

61

contaminants as the area is largely affected by winds and ocean currents from Europe and North-

62

America.22,23

63

Arctic mammalian predators are exposed to high concentrations of PFAS because PFASs

64

biomagnify in high latitude food webs.3,24,25 PFAS exposure in Arctic biota, in particularly in

65

polar bears (Ursus maritimus), is alarming. Concentrations of PFOS in polar bear plasma are

66

similar to those reported in humans living at the proximity of a PFAS manufacturing plant in

67

China.26-28 Documented adverse health effects of PFASs in experimental animals and humans are

68

numerous.5,29,30 Modelling studies suggest that current PFAS concentrations in polar bears are

69

high enough to contribute to reproductive, immunotoxic and genotoxic effects31 and correlative

70

studies have linked PFAS concentrations in polar bears to changes in steroid and thyroid

71

hormone concentrations.32,33 Furthermore, PFAS concentrations in Arctic mammals may

72

increase during seasonal food deprivation periods.28,34

73

Temporal trends of PFASs are highly variable in Arctic predators,35-38 which is likely related to

74

the complex emission sources and variance in transport mechanisms of these compounds. It is

75

thus essential to quantify the temporal trends (year to year) of PFAS concentrations in Arctic

76

wildlife that may reflect source changes toward Asia and South America from traditional sources

ACS Paragon Plus Environment

4

Page 5 of 43

Environmental Science & Technology

77

of Western countries.7,10 However, emission history is not the only factor that may influence

78

temporal changes in contaminant exposure in Arctic predators. Climate-driven changes in ocean

79

circulation patterns may affect long-range transport of PFAS15and, furthermore food web

80

structures and predator-prey interactions are changing in the Arctic ecosystems.39 This likely

81

influences contaminant exposure in Arctic biota as documented recently for polar bears.40-42

82

Polar bears and arctic foxes (Vulpes lagopus) are among the highest contaminated Arctic apex

83

predators.43 Polar bears and arctic foxes from Svalbard are species that can provide information

84

about temporal trends of PFASs in the Arctic environment as well as how PFAS accumulation

85

may be related to diet choices. For instance, both species feed on variable food items from both

86

the marine and terrestrial ecosystem,44-46 which may affect their exposure to contaminants. Polar

87

bears feed preferably on ringed seals (Pusa hispida), but they also ingest other marine food items

88

such as whale carcasses, seabirds and terrestrial prey like geese and reindeer, particularly in the

89

absence of sea ice.44,47-49 The arctic fox in Svalbard is a predator and scavenger foraging on both

90

marine and terrestrial food items including reindeer, ptarmigans, geese, seabirds and seals.45,46,50

91

However, feeding habits of polar bears and arctic foxes in Svalbard have likely changed over

92

time because climate warming is especially pronounced in the Svalbard area compared to other

93

Arctic areas.51 For instance, the maximum extent of sea ice in the Barents Sea, including

94

Svalbard, has decreased by ~ 50% since the late 1990s,52 the length of the summer season (i.e.

95

time with reduced ice) has increased by over 20 weeks between 1979 and 201353 and winter sea

96

ice is retreating.54 Furthermore, duration of snow cover on land has decreased 0.5 day per year

97

since 1976 whereas rain-on-snow events are predicted to increase by 40 % at the end of this

98

century.51,55

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 43

99

The aim of the study was to investigate temporal trends of PFASs in relation to climate-

100

associated changes in feeding habits and food availability in two Arctic top predators from

101

Svalbard, namely the polar bear and the arctic fox during the period 1997-2014. Furthermore, we

102

discussed about the role of emission history and long-range transport mechanisms on PFAS

103

trends in polar bears and arctic foxes. We also analyzed several emerging PFASs.

104 105

MATERIALS AND METHODS

106

Field sampling

107

Adult female polar bears from the Barents Sea subpopulation were captured annually between

108

25th March and 5th May in the period 2000-2014. The 192 samples collected opportunistically

109

throughout Svalbard represented 137 individual females. Twenty-eight females were captured 2-

110

8 times. Details of immobilization, determination of age, body condition, breeding status and

111

blood sampling are given in supporting information. Breeding status, age and body condition for

112

each sampling year are given in Table S1.

113

Arctic foxes (n=113) were trapped by local trappers on western Spitsbergen, Svalbard, around

114

the Isfjorden area and Nordenskiöld Land during the annual harvest between 1st November and

115

15th March in the period 1997-2014. All trapped individuals were less than two years old and

116

none of the females had given birth. Among the sampled foxes, number of individuals, as well as

117

sex, age and body condition were balanced over the eleven trapping seasons (Table S2). Details

118

of sampling as well as determination of body condition (ranging from 1-barely measurable fat to

119

4-extensive fat) have been described previously.56,57 Samples of skeletal muscle and liver were

120

packed in aluminium foil and stored at -20°C until analysed.

121

Proxies for feeding habits

ACS Paragon Plus Environment

6

Page 7 of 43

Environmental Science & Technology

122

Stable isotope values of nitrogen and carbon (δ15N and δ13C, respectively) were used as

123

proxies for feeding habits of polar bears and arctic foxes. Due to the availability of different

124

body compartments for the two species, δ15N and δ13C were measured in polar bear red blood

125

cells and arctic fox muscle tissue. Half-lives for δ15N and δ13C in polar bear red blood cells are

126

over 3 months and ~1.5 months, respectively,58 whereas half-lives for δ15N and δ13C for muscle

127

tissue of a small mammal is approximately a month.59 As polar bears were captured in April and

128

arctic foxes collected in November-March, stable isotope values for both species thus reflect

129

mostly their winter diet. Analytical procedures are described in supporting information.

130

Proxies for food availability

131

Polar bear habitat preference likely reflects the occurrence and availability of seals to polar

132

bears.60 We thus used sea ice habitat quality available for individual polar bears as a proxy for

133

availability of seals to polar bears. First, we categorized surrounding areas of Svalbard according

134

to the abundance of preferred sea ice habitat. The yearly (2000-2014) distribution and duration of

135

preferred sea ice habitat was based on resource selection function (RSF) models as described in

136

detail in supporting information. According to the maps produced by RSF models, preferred

137

habitat was consistently available for more days of the year on the east side of Svalbard than in

138

the west side of Svalbard over the whole study period (Figure S1). Among the 137 individual

139

polar bears studied, we used telemetry data (n=83) or capture position (n=54) to assign them to

140

the better habitat on the eastern side or the poorer habitat on the western side as detailed in

141

supporting information.

142

We used reindeer mortality and sea ice cover as proxies for availability of Svalbard reindeer

143

carcasses and seals, respectively, to arctic foxes.57 Reindeer mortality was derived from the

144

ongoing annual long-term summer monitoring of structural composition and mortality of

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 43

145

Svalbard reindeer in Adventdalen, Nordenskiöld, and expressed as number of reindeer that died

146

during the past 12 months in Advendtdalen.57,61,62 As reindeer mainly die in late winter and

147

spring63 they are available as food items for foxes mainly in the late winter and spring.

148

Therefore, we used the number of reindeer carcasses observed during the summer preceding the

149

arctic fox trapping season as a proxy.

150 151

To determine availability of sea ice cover for individual foxes, we calculated an index for each foxes based on daily sea ice maps of Isfjorden as described in supporting information.

152

Analyses of PFAS

153

Due to the availability of different body compartments for the two species, we monitored

154

PFASs in plasma of polar bears and liver of arctic foxes. The 17 PFASs monitored in polar bear

155

plasma and arctic fox liver are listed in Table 1. Methods for clean-up, separation, quantification

156

and quality assurance are explained in supporting information.

157

Data analyses

158

We used the statistical program R version 3.3.1. for data analyses.64 We used additive mixed

159

models (GAMM) on R-package mgcv65 to analyse the effects of year, feeding habits and food

160

availability, in addition to biological variables, on PFAS concentrations in polar bears and arctic

161

foxes. To study PFAS concentrations in polar bears, we defined 26 candidate models with year,

162

δ15N and δ13C values in red blood cells, habitat quality (high vs. low), breeding status, age and

163

body condition as fixed predictor variables, and individual ID as random variable to account for

164

repeated measurements (Table S6). As perfluoroundecanoate (PFUnDA) concentration was

165

considerably higher in samples analysed in 2012 compared to those analysed in 2014-15 (median

166

49 and 20 ng/g, respectively) for unknown reasons, we included batch (2012 vs. 2014-15) as a

167

random variable in the polar bear PFUnDA model. A smaller set of candidate models (n=5) were

ACS Paragon Plus Environment

8

Page 9 of 43

Environmental Science & Technology

168

applied

for

perfluorobutane

sulfonate

(PFBS),

perfluoroheptanoate

169

perfluorotetradecanoate (PFTeDA) due to the limited number of polar bear samples (n=70)

170

analysed for these compounds. To study liver PFAS concentrations in arctic foxes, we defined 18

171

candidate models including the effects of year, muscle δ15N and δ13C, reindeer carcasses, sea ice

172

cover and body condition (Table S7). As we hypothesize non-linear PFAS trends due to their

173

complex emission history, we entered year as a smoothed term using penalized regression splines

174

and Gaussian distribution.65 The other predictor variables were applied as linear numerical or

175

categorical variables. The models were fitted using the restricted maximum likelihood (REML)

176

method. Highly correlated predictor variables (e.g. δ15N and δ13C) were not included in the same

177

model as recommended.66 To make inference from all candidate models and subsequently from

178

all predictor variables, we used model averaging based on Akaike’s Information Criterion

179

(AIC)67 using the R-package MuMIn.68 More specifically, we ranked the models according to

180

AIC, which was further used to calculate AIC weight (e(0.5(AICmin

181

divided by the sum of all relative likelihoods). Further, we calculated model averaged estimates

182

for all the predictor variables in the candidate model list weighted using AIC weights. We used

183

95% confidence intervals of the model-averaged estimates to determine whether the parameters

184

were significantly different from 0 at the 5% confidence level. We used diagnostic plots (Figure

185

S2-S3) to identify that the distribution of model residuals met the model assumptions.69 The

186

residual plots revealed an outlier in models explaining perfluorohexane sulfonate (PFHxS)

187

concentrations in arctic foxes. We excluded this outlier, which did not lead to significant changes

188

in parameter estimates.

– AICi))

(PFHpA)

and

); relative likelihood

189

We used plots from the highest ranked GAMMs to evaluate during which period PFAS

190

concentrations linearly changed in polar bear and arctic fox food web. Estimates for yearly

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 43

191

changes in PFAS concentrations were derived from linear mixed effects models (LMEM) (R

192

package lme470) and linear models (LM) for polar bears and arctic foxes, respectively. Polar bear

193

individual ID was included in LMEMs as random variable. To obtain yearly changes in polar

194

bear/arctic fox food web, we adjusted the changes for biological and environmental variables i.e.

195

we included the same covariates in the LMEMs and LMs as in the top GAMM for the given

196

compound in the given species. The non-adjusted trends were derived from models that excluded

197

the continuous predictors reflecting feeding habits and food availability from the top GAMMs

198

(δ13C, δ15N in both species, and sea ice cover and reindeer mortality in arctic foxes). The annual

199

changes (%) in the median concentrations were calculated from 100 * (eestimate

200

Temporal changes in feeding habits were investigated by GAMMs, and the yearly linear changes

201

were obtained using LMEMs or LMs.

for year

– 1).

202 203

RESULTS AND DISCUSSION

204

Concentrations of PFASs in polar bears and arctic foxes

205

PFOS was quantitatively the most abundant PFAS detected in both polar bears and arctic foxes

206

(Table 1). The contribution of PFOS to ΣPFAS decreased from over 80% to ~60% during the

207

study period. Linear and branched PFOS, investigated only in arctic fox liver samples, were

208

highly correlated (r=0.94, 95% CI [0.91, 0.96], pLOD Fluorotelomer sulfonates (FTSA) 6:2 FTSA C8 n.a. 8:2 FTSA C10 n.a. Perfluoroalkane sulfonamido substances FOSA C8 0.02