Elucidating the Behavior of Cyclic Volatile Methylsiloxanes in a

NILU − Norwegian Institute for Air Research, P.O. Box 100, 2027 Kjeller, Norway. ○ Department of Chemistry, University of Oslo, P.O. Box 1033, Bli...
1 downloads 9 Views 419KB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Elucidating the Behavior of Cyclic Volatile Methylsiloxanes in a Subarctic Freshwater Food Web: A Modeled and Measured Approach Ingjerd Sunde Krogseth, Emma Maria Undeman, Anita Evenset, Guttorm N. Christensen, Mick J. Whelan, Knut Breivik, and Nicholas A. Warner Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b03083 • Publication Date (Web): 05 Oct 2017 Downloaded from http://pubs.acs.org on October 9, 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 21

Environmental Science & Technology

3

Elucidating the Behavior of Cyclic Volatile Methylsiloxanes in a Subarctic Freshwater Food Web: A Modeled and Measured Approach

4

Dr. Ingjerd S. Krogseth,1,* Dr. Emma Undeman,2,3 Dr. Anita Evenset,4,5 Cand. Scient. Guttorm N.

5

Christensen,4 Dr. Mick J. Whelan,6 Dr. Knut Breivik,7,8 Dr. Nicholas A. Warner1

1 2

6 7 8 9 10 11 12 13

1 2

NILU – Norwegian Institute for Air Research, The Fram Centre, P.O. Box 6606 Langnes, 9296 Tromsø, Norway; Baltic Sea Centre, Stockholm University, SE-106 91 Stockholm, Sweden; 3Department of Environmental Science and Analytical Chemistry, Stockholm University, Svante Arrhenius väg 8, SE-106 91, Stockholm, Sweden; 4 5 Akvaplan-niva AS, The Fram Centre, P.O. Box 6606 Langnes, 9296 Tromsø, Norway; Faculty of Biosciences, Fisheries and Economics, UiT the Arctic University of Norway, Hansine Hansens veg 18, 9019 Tromsø, Norway; 6 School of Geography, Geology and the Environment, University of Leicester, Leicester LE1 7RH, United 7 8 Kingdom; NILU – Norwegian Institute for Air Research, P.O. Box 100, 2027 Kjeller, Norway; Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway.

14 15

*Corresponding author: NILU – Norwegian Institute for Air Research, The Fram Centre, P.O. Box 6606 Langnes, 9296 Tromsø, Norway, tel: +47 63 89 80 00, [email protected]

16 17

Abstract

18

environments through wastewater effluents, and their bioaccumulation potential is debated. Here, a

19

new bentho-pelagic version of the ACC-HUMAN model was evaluated for polychlorinated biphenyls

20

(PCBs) and applied to cVMS in combination with measurements to explore their bioaccumulation

21

behavior in a subarctic lake. Predictions agreed better with measured PCB concentrations in Arctic

22

char (Salvelinus alpinus) and brown trout (Salmo trutta) when the benthic link was included than in

23

the pelagic-only model. Measured concentrations of decamethylcyclopentasiloxane (D5) were 60 ±

24

1.2 (Chironomidae larvae), 107 ± 4.5 (pea clams Pisidium sp.), 131 ± 105 (three-spined sticklebacks:

25

Gasterosteus aculeatus), 41 ± 38 (char), and 9.9 ± 5.9 (trout) ng g-1 wet weight. Concentrations were

26

lower for octamethylcyclotetrasiloxane (D4) and dodecamethylcyclohexasiloxane (D6), and none of

27

the cVMS displayed trophic magnification. Predicted cVMS concentrations were lower than

28

measured in benthos, but agreed well with measurements in fish. cVMS removal through ventilation

29

was an important predicted loss mechanism for the benthic-feeding fish. Predictions were highly

30

sensitive to the partition coefficient between organic carbon and water (KOC) and its temperature

31

dependence, as this controlled bioavailability for benthos (the main source of cVMS for fish).

Cyclic volatile methylsiloxanes (cVMS) are used in personal care products and emitted to aquatic

1 ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 21

1. Introduction

32 33

Cyclic

34

decamethylcyclopentasiloxane (D5), and dodecamethylcyclohexasiloxane (D6), are used as

35

ingredients in personal care products and are emitted to aquatic environments through wastewater

36

effluents.1-3 D4 has been classified as a persistent, bioaccumulative, and toxic (PBT) substance and D5

37

as a very persistent and very bioaccumulative (vPvB) substance in a European Chemical Agency

38

Member State Committee opinion.4 A restriction on use in wash-off products is currently underway

39

in Europe.5 However, their hazardous properties, particularly their bioaccumulation potential, are

40

uncertain.6-8 Field studies of cVMS in aquatic food webs have demonstrated both trophic dilution9-12

41

and trophic magnification,13,14 and the explanation(s) for this variability remains unclear. Possible

42

influencing factors include behavioral and physiological characteristics of the investigated taxa, food

43

chain composition, and the nature of the physical environment (including possible confounding

44

effects such as spatial gradients in concentration).7,15

45

As well as being very volatile, cVMS are very hydrophobic, and are expected to partition significantly

46

to organic phases. Although it is now well established that the organic carbon (OC) to water partition

47

coefficients (KOC) for cVMS are lower than those expected from their octanol to water partitioning

48

(KOW),16-19 there is uncertainty in the temperature dependence of KOC which currently constrains our

49

understanding of the behavior of cVMS in cold aquatic environments.20 Although a positive

50

relationship between KOW and temperature has been reported (suggesting lower hydrophobicity at

51

low temperatures),19,21,22 a negative relationship between KOC and water was recently reported

52

(suggesting higher hydrophobicity at low temperatures).23 This will influence the behavior of these

53

compounds in systems with low temperatures. In aquatic systems without continuous wastewater

54

emissions, the exposure of pelagic organisms to these compounds is likely to be primarily via benthic

55

links in the food web, as cVMS may be rapidly removed from water between emission events.20

56

The combination of fugacity-based models and measurements has proved useful for developing an

57

understanding of the environmental behavior of cVMS and the identification of key knowledge

volatile

methyl

siloxanes

(cVMS),

such

as

octamethylcyclotetrasiloxane

2 ACS Paragon Plus Environment

(D4),

Page 3 of 21

Environmental Science & Technology

58

gaps.24-27 However, this approach has rarely been used for cVMS behavior in food webs. Exceptions

59

include applications of the AquaWeb6 and ACC-HUMAN models.24,28 ACC-HUMAN predicted trophic

60

dilution of cVMS in zooplankton, herring, and cod from the Inner Oslofjord (which receives high

61

wastewater emissions), in accordance with measurements.24 However, given the importance of

62

sediments as a reservoir of cVMS, there is a need for a benthic link in the model. This is particularly

63

important for systems such as Storvannet where fugacities in sediment often exceed those in the

64

water column due to intermittent emission events and rapid advection of the water column.20 Here,

65

we explore the bioaccumulation behavior of cVMS using the ACC-HUMAN model, which was

66

expanded to include a benthic link. Model performance was first evaluated for well-known PCBs, and

67

the model was subsequently applied to explain measured cVMS concentrations in biota sampled

68

from Lake Storvannet in Northern Norway. Lake Storvannet was chosen as (i) the lake and its food

69

web is relatively well-studied,29,30 (ii) measurements of PCBs from lake water, sediments, and biota

70

are available for model evaluation purposes,31 and (iii) a study of cVMS behavior in the physical

71

environment of the lake has already been carried out.20

72

2. Methods

73

2.1. Study site

74

Lake Storvannet is a small (0.2 km2, maximum depth 17 m), oligotrophic lake, located in a residential

75

area of the town of Hammerfest in Northern Norway (70°N 23°E; 7000 inhabitants).20 The lake has a

76

long history of local contamination with organic contaminants including PCBs.31 Untreated

77

wastewater from Hammerfest is now emitted into the sea. However, sewers run along both sides of

78

Storvannet, and the lake unintentionally receives wastewater through leaking pipes and combined

79

sewer overflow events (CSOs), particularly during the snowmelt season and after heavy rainfalls.20 As

80

a result, cVMS are present in sediments of the lake with concentrations comparable to much more

81

populated areas (207 ± 30, 3775 ± 973, and 848 ± 211 ng g-1 OC for D4, D5, and D6, respectively).20

82

However, cVMS were not detected in water samples from the lake (LOQ 19-40 ng L-1, 12-54 ng L-1,

83

and 42-74 ng L-1 for D4, D5, and D6, respectively).20 The lake is ice covered for approximately six

3 ACS Paragon Plus Environment

Environmental Science & Technology

84

months of the year and the summer season and the plankton bloom is short. Benthic Chironomidae

85

larvae (hereafter; chironomids) and pea clams (Pisidium sp.) are, consequently, important dietary

86

items for resident fish.29,32,33 Populations of both stationary and anadromous Arctic char (Salvelinus

87

alpinus) (hereafter; char) and brown trout (Salmo trutta) (hereafter; trout), as well as three-spined

88

sticklebacks (Gasterosteus aculeatus) (hereafter; sticklebacks) inhabit the lake.29,30 While char mostly

89

feed on benthos, the trout has a more pelagic diet including aquatic insects.29,34,35 Both the trout and

90

char may be cannibalistic, but the trout tends to be more piscivorous than the char and also feed on

91

sticklebacks (Figure 1, Table 1, Supporting Information (SI) 2.2). The sticklebacks feed on plankton

92

and benthos in the littoral zone of the lake. Only stationary fish were included in this study.

93 94

Char (n = 11) and trout (n = 13) were caught opportunistically in March, May, and June 2014 (Table

95

S1). In March and May, the fish were caught with lines and hooks. In June, gill nets were used. In

96

addition, traps, baited with fish roe (wrapped in fine mesh so the fish could not eat it), were

97

deployed in shallow areas of the lake for collection of sticklebacks (n = 5). Benthic fauna

98

(Chironomidae and Pisidium sp.) were collected by hauling a 500 µm trawl along the bottom. All

99

samples were wrapped in aluminum foil, placed in zip-lock bags, and transported in chilled boxes to

100

the lab where they were stored at -20 °C until sample preparation. All use of personal care products

101

was avoided by field and lab personnel. Muscle and liver tissues from Atlantic cod (Gadus morhua)

102

caught off the Norwegian coast were utilized as blank material. Prior to sample collection, the cod

103

tissues were split into three separate portions to serve as field blanks, dissection blanks, and

104

reference blanks, in order to assess if any contamination occurred during field work or sample

105

dissection.

106

Strict quality control procedures were enforced to limit cVMS contamination of samples (SI 1.1). The

107

bodyweight and total length of each fish were recorded, gender and maturity were determined

108

through visual inspection of the gonads, and stomach content was visually examined (Table S1, S2).

109

Otoliths were retrieved to determine the age of the char and trout (Table S1). Liver and dorsal

2.2. Sample collection, preparation and analysis

4 ACS Paragon Plus Environment

Page 4 of 21

Page 5 of 21

Environmental Science & Technology

110

muscle tissue was taken out for analysis of cVMS, extractable organic matter (assumed

111

representative for lipids),36 and stable isotope analysis (muscle only; δ15N, δ13C, δ14S) (SI 1.2, Table 2).

112

Prior to analyses, both liver and muscle samples were homogenized with a metal scalpel on a clean

113

glass surface. Due to the small size of the sticklebacks, the whole fish was homogenized and the

114

homogenate was analyzed for cVMS and stable isotopes. The collected chironomids were pooled and

115

homogenized in a mortar, giving two replicates for cVMS analysis and one for stable isotopes. The

116

same procedure was followed with whole Pisidium samples (including the carbonaceous shell).

117

All samples were extracted for cVMS using a biphasic cold solvent extraction, based on a previously

118

published method for sediments,20 but adapted to biotic tissues (SI 1.3). Analysis of cVMS was carried

119

out on an Agilent 7890A GC connected to an Agilent 5975C MS detector and a Gerstel MPS3

120

autosampler. Method details have previously been published.20 Calculation of limits of detection

121

(LOD) and quantification (LOQ) from blank samples is described in SI 1.1.

122

2.3. Bioaccumulation modeling

123 124

2.3.1. Development of the bentho-pelagic ACC-HUMAN The aquatic part of the original ACC-HUMAN model represents a pelagic Baltic Sea food web

125

consisting of zooplankton, a planktivorous fish (herring) and a piscivorous fish (cod).28 Here, a benthic

126

link was included in the model by adding a sediment compartment, a benthic filter feeder, a benthic

127

predator/deposit feeder, and up to three prey items for the latter. The model does not simulate

128

chemical transfer between lake water and sediment. Rather, time-resolved concentrations or

129

fugacities in sediment and lake water is provided as input to the model. Sediment pore water and

130

sediment solids are assumed to be in equilibrium.

131

The prey items are defined by their biochemical composition (fractions of water, lipid, and non-lipid

132

organic matter), the fraction of time they spend buried in sediments, and assumed thermodynamic

133

equilibrium with their surroundings. Bioaccumulation in the predator/deposit feeder and filter feeder

134

are predicted using a dynamic version of the steady-state model of Arnot and Gobas.37 The fugacity

135

capacity (Z-value) of the organism is calculated based on its biochemical composition. Chemical 5 ACS Paragon Plus Environment

Environmental Science & Technology

136

uptake occurs through feeding and ventilation (of sediment pore water and/or overlying water).

137

Elimination of chemicals can occur through ventilation, fecal egestion and/or biotransformation. The

138

predator/deposit feeder may feed on the prey items and/or sediment particles, while the filter

139

feeder may feed on suspended or settled sediment particles and/or zooplankton. The fish in the

140

model may prey on the benthic filter feeder and the deposit feeder, in addition to zooplankton and

141

fish. The model assumes annual constant growth rate, lipid content, and dietary composition.

142

However, a seasonal feeding pattern for the fish was included in the model, in which the fish are

143

assumed not to eat during wintertime (November – February) as is the case in Storvannet. All

144

equations, parameters, and assumptions for the new model are described in SI 2.1 and 2.2.

145

To adapt the bentho-pelagic ACC-HUMAN to Storvannet, the model organisms were parameterized

146

as key species of the Storvannet food web. Benthos were parameterized as chironomid larvae

147

(deposit feeder) and Pisidium clams (filter feeder). The pelagic fish was re-parametrized based on

148

sticklebacks, while the piscivorous fish was re-parametrized to represent either char or trout (SI 2.2,

149

Table 1). Dietary composition of the organisms and overall structure of the food web was based on

150

existing knowledge about the system aided by analysis of stable isotopes and gut content.

151

Measurements of water temperatures (Table S11),38 fraction of solids in sediments (18 % v/v) and

152

fraction of OC in sediment solids (5 % w/w)20 from Storvannet were used as model input.

153 154

2.3.2. Model evaluation and application Many model aspects are independent of the compounds studied and are best constrained on the

155

basis of data-rich chemicals like PCBs. To evaluate the performance of the new bentho-pelagic model

156

version, it was compared with the pelagic-only version of ACC-HUMAN parameterized for the taxa

157

present in Storvannet. Predicted concentrations in char and trout were compared to measurements

158

from Storvannet for five PCB congeners (SI 2.3).

159

Empirical data for exposure-relevant physical media are restricted to sediment for cVMS. Hence, only

160

the bentho-pelagic model version was applied to cVMS, as the pelagic-only model needs

6 ACS Paragon Plus Environment

Page 6 of 21

Page 7 of 21

Environmental Science & Technology

161

concentrations in water as input. Average measured concentrations of cVMS in surface sediments (0-

162

1 cm, OC-normalized) of Storvannet in 2014 were used as input to the model (Table S13).20 cVMS

163

concentrations in Storvannet lake water sampled at the same time were below limits of

164

quantification (LOQ).20 Given the lack of measured concentrations, a range of values for the fugacity

165

in water was evaluated, and expressed relative to the fugacity in sediment derived by the model

166

from the concentration in sediment (SI 2.4).

167

Measured partition coefficients (KOW and KAW (air-water partitioning)) and their temperature

168

dependencies (ΔUOW and ΔUAW) were used as input (Table S13).19,20,39 The KOC values for cVMS are

169

known to be significantly overestimated by empirical relationships with KOW.19,39 This deviation may

170

also apply for ΔUOC versus ΔUOW. The most recently measured values of KOC and ΔUOC for cVMS were,

171

therefore, used as direct inputs, rather than relying on relationships based on affinity to octanol

172

(Table S13).18,23 That said, partitioning to lipids was still assumed to be represented by partitioning to

173

octanol,40,41 although there is limited evidence to justify this assumption. cVMS may undergo

174

biotransformation in aquatic biota, but few measurements are available. Here, a fish

175

biotransformation rate constant (kM) of 4 · 10-4 h-1 (corresponding to a half-life of 72 days) was

176

applied for all cVMS (Table S13), following Whelan and Breivik.24,42 In the absence of evidence, cVMS

177

were assumed to not biotransform in benthos or zooplankton.

178

The model was run for 10 years with a 1 hour simulation time step. Predicted median lipid-

179

normalized concentrations in year 10 across all age classes at the time of year when the samples

180

were collected, as well as the annual range (min-max), were compared to measured lipid-normalized

181

concentrations in whole-body homogenates (chironomids, Pisidium, and sticklebacks) and in muscle

182

tissue (char and trout). As the model does not provide gender- and/or maturity-specific results for

183

fish, measurements were grouped irrespective of gender and maturity for comparison to model

184

results. Additional simulations were carried out to test the sensitivity of model results to changes in

7 ACS Paragon Plus Environment

Environmental Science & Technology

185

cVMS physical-chemical properties and biotransformation rates (properties were changed

186

individually by ± 10 % and ± 20 %).

187

3. Results and discussion 3.1. Measured concentrations of cVMS

188 189

All cVMS concentrations in benthos were above LOQ (Table S3). The D4, D5, and D6 whole body

190

concentrations (pooled duplicates) in chironomids were 9.9 ± 0.3, 60 ± 1.2, and 9.3 ± 0.1 ng g-1 wet

191

weight (ww), and in Pisidium 4.7 ± 0.4, 107 ± 4.5, and 11.8 ± 1.2 ng g-1 ww, respectively (Table 2). The

192

measured concentrations in benthos are within the range of measured concentrations in benthos

193

reported elsewhere, including chironomids and burrowing mayfly in Lake Pepin, ragworm in the

194

Humber Estuary, and polychaets from the Oslofjord (Table S8).9,10,43,44 Lipid content for benthos was

195

not measured due to insufficient sample material. Instead this was derived from literature (Table 2).

196

Hence, the lipid-normalized concentrations have added uncertainty, particularly for Pisidium which

197

has a very low lipid content. This does not matter for model comparison purposes as the lipid

198

content assumed for normalization was used as model input and hence errors in lipid-content will

199

cancel out in model-measurement comparisons. The benthic samples may also have contained traces

200

of sediment particles, as they were not depurated after sampling. This means that measured

201

concentrations of (i) D4 in chironomids and D5 in Pisidium could have been diluted by presence of

202

sediment particles (concbenthos > concsed (ng g-1 ww)), while (ii) D5, D6 in chironomids and D4, D6 in

203

Pisidium could have been contaminated by presence of sediment particles (concbenthos ≤ concsed (ng g-1

204

ww)).20 Hence, cVMS concentrations in benthos should be considered approximate and are only used

205

for model-measurement comparison.

206

All cVMS concentrations in whole-body homogenates of sticklebacks were above LOQ (Table S3),

207

with concentrations of 12.6 ± 9, 131 ± 105, and 6.3 ± 3.6 ng g-1 ww for D4, D5, and D6, respectively

208

(309 ± 226, 3203 ± 2584, and 155 ± 89 ng g-1 lipid weight (lw)) (n = 5) (Table 2). In muscle of char and

209

trout, D4, D5, and D6 were above LOQ in 37 %, 97 %, and 14 % of the samples, respectively (Table

210

S3). D4 and D6 were below LOQ in trout, while D5 ranged between 3.9 – 22.0 ng g-1 ww (383 – 2434 8 ACS Paragon Plus Environment

Page 8 of 21

Page 9 of 21

Environmental Science & Technology

211

ng g-1 lw) (n = 13). In char, D4, D5, and D6 ranged between