Heightened Biological Uptake of Polybrominated

May 20, 2014 - patterns and help identify potential source and sink regions.6. Uptake of POPs ... Vancouver, and Seattle in the adjacent Puget Sound),...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/est

Heightened Biological Uptake of Polybrominated Diphenyl Ethers Relative to Polychlorinated Biphenyls Near-Source Revealed by Sediment and Plankton Profiles along a Coastal Transect in British Columbia Jean-Pierre W. Desforges,† Neil Dangerfield,‡ Patrick D. Shaw,§ and Peter S. Ross*,† †

Ocean Pollution Research Program, Vancouver Aquarium, P.O. Box 3232, Vancouver BC V6B 3X8, Canada Institute of Ocean Sciences, Fisheries and Oceans Canada, P.O. Box 6000, Sidney BC V8L 4B2, Canada § Environment Canada, 201-401 Burrard Street, Vancouver, BC V6C 3S5, Canada ‡

S Supporting Information *

ABSTRACT: Polychlorinated biphenyl (PCB) and polybrominated diphenyl ether (PBDE) concentrations and profiles in paired sediment− plankton samples were determined along a 500 km transect in coastal British Columbia, Canada. PCB and PBDE levels in sediment were both greater in the industrialized Strait of Georgia than in remote northern sites and exhibited parallel spatial trends. In plankton, recent-use PBDE levels were higher near-source, while levels of legacy PCBs were uniform across sites. Principal component analysis of 95 PCB congeners illustrated the influence of proximity to source (i.e., latitude) on congener patterns for both matrices (sediment, r2 = 0.52, p = 0.012; plankton, r2 = 0.59, p = 0.016). The PCB pattern in plankton grew lighter with latitude, but the opposite pattern in sediments suggested that temperature-related fractionation, sediment processes, and basin-wide oceanography had divergent effects on each matrix. Biota-sediment accumulation factors (BSAFs) were greater for PBDEs than PCBs, but spatial profiles were similar; PCBs and PBDEs were near equilibrium in remote atmospherically driven sites (BSAF = 1.7 and 1.3) but accumulated preferentially in sediments at source-driven sites (BSAF = 0.2 and 0.4). The influences of particle-binding and hydrophobicity on the aquatic fate of PCBs and PBDEs was evident by the strong influence of log KOW on congener-specific BSAFs (PCBs, r2 = 0.18 p < 0.001; PBDEs, r2 = 0.61 p < 0.001). While biotic uptake of PCBs has become spatially uniform in coastal BC because of dilution over time, biomagnification of PBDEs remains higher in industrialized waters.



INTRODUCTION

reflects thermodynamic equilibrium between the organic carbon fraction in sediment and the lipid fraction in biota.5 Given that sediment and biota may not be tightly coupled and that equilibrium conditions are rarely observed in the field, the BSAF approach may be more useful as a tool to evaluate deviation from equilibrium (i.e., BSAF ∼1 to 2) in order to garner insight into chemical behavior and fate. In this manner, the “equilibrium status” may be used to shed light on spatial patterns and help identify potential source and sink regions.6 Uptake of POPs by phytoplankton is governed by equilibrium partitioning between tissue and water.6 Uptake in phytoplankton is particularly important as their settling to sediment drives regional and global POP fate and because plankton represent the entry point of POPs to marine food

Chemicals released into the environment cycle through various media, with pathways and ultimate fate dependent upon the interaction of chemical properties and environmental conditions.1 The global distribution of persistent organic pollutants (POPs) has been largely attributed to long-range atmospheric transport (LRAT), a temperature-driven process of volatilization and deposition which influences chemical fractionation.2 Despite the importance of LRAT, POP fate in the marine environment may be ultimately driven by biogeochemical processes.3 Air−water exchange, uptake by phytoplankton, and sedimentation of particles are among the most important processes affecting the fate of POPs in aquatic environments. As POPs sorb preferentially to organic matter, sequestration by sinking particles and burial in deep ocean sediment is a major sink for these lipophilic contaminants.3,4 The calculation of a biota-sediment accumulation factor (BSAF) provides insight into the propensity of a contaminant to biomagnify in food webs. A BSAF between 1 and 2 typically © 2014 American Chemical Society

Received: Revised: Accepted: Published: 6981

January 14, 2014 April 16, 2014 May 20, 2014 May 20, 2014 dx.doi.org/10.1021/es500218b | Environ. Sci. Technol. 2014, 48, 6981−6988

Environmental Science & Technology

Article

Figure 1. Total concentration of PCBs and PBDEs in sediment (ng/g TOC) and plankton (ng/g lw) samples along a coastal transect in the northeast Pacific Ocean from Hecate Strait to the Strait of Georgia, British Columbia, Canada. Sediment PBDEs include BDE-209 (displayed in stacked gray bar) while plankton do not (see Methods). Sediment samples from sites 3 and 6 were not collected (nc) due to the rocky nature of the seabottom.

transect of coastal British Columbia from Hecate Strait to the southern Strait of Georgia (SoG). The SoG is urbanized (∼4 million people) and receives industrial and wastewater discharges from multiple sources (e.g., urban Victoria and Vancouver, and Seattle in the adjacent Puget Sound),10 while the northern part of the transect is remote and has no major point-sources of contaminants. The SoG is semienclosed and subject to high sedimentation rates and elevated inputs of organic carbon through net primary production, the Fraser River and anthropogenic sources (i.e., wastewater discharge and pulp mills).10 In contrast, Hecate Strait/Queen Charlotte Sound is characterized by open ocean conditions with moderate to high exposure and low sedimentation rates and organic carbon.23 Our aim was to examine the geographical pattern and bioaccumulation of PCBs and PBDEs along a marine transect in order to garner insight into the uptake of these two priority pollutants in marine food webs.

webs. The accumulation of POPs in plankton with subsequent food web biomagnification leads to high concentrations in higher trophic levels, such as the heavily PCB contaminated killer whales (Orcinus orca) in the NE Pacific Ocean.7 Consequently, sediment and plankton POP concentrations can be used to predict uptake by higher trophic level organisms such as killer whales and provide guidance on different sediment management or ocean disposal scenarios.8 Polychlorinated biphenyls and polybrominated diphenyl ethers (PBDEs) are two contaminants of particular concern since both are persistent, bioaccumulative, and toxic. PCBs are legacy pollutants (banned in 1970s in Canada and the USA), whereas PBDEs have only recently faced regulatory scrutiny.9 The contrasting history of these pollutants has implications for their current environmental distribution: PCBs largely reflect cycling of legacy reservoirs while PBDEs continue to reflect point-source releases.10,11 In addition, PBDEs have a wider range of molecular size and KOW and are more influenced by environmental transformation than PCBs, resulting in differences in their behavior and fate.12−14 While studies have been conducted along transects for contaminants in air,6,15,16 water,16,17 soil,6,18 vegetation,19,20 and biota,21,22 few such studies compare different contaminants in biota and sediment in the marine environment. Here were present the first large spatial transect study of PCB and PBDE profiles in paired sediment−plankton samples. We investigated congener-specific concentrations and patterns along a 500 km



METHODS Sample Collection. Sediment samples were taken from nine sites and plankton samples from 11 equidistant sites along a transect from Hecate Strait, Queen Charlotte Sound (QCS), Johnstone Strait through to the southern Strait of Georgia (BC, Canada) in September 2006 aboard the CCGS Vector (Figure 1). Surface sediments were collected using a Smith-McIntyre grab sampler (∼15 cm depth), and the top few centimeters were transferred to solvent-rinsed amber glass jars after 6982

dx.doi.org/10.1021/es500218b | Environ. Sci. Technol. 2014, 48, 6981−6988

Environmental Science & Technology

Article

Table 1. Concentrations of total PCBs and total PBDEs in Sediments and Plankton (ng/g wet weight)a sediment I.D.

site

plankton

depth (m)

TOC (%)

lipid (%)

PCB

PBDE

PCB

PBDE

1 Queen’s Sound 2 W. Calvert Island 3 S. Calvert Island 4 N. Malcolm Island 5 Havannah Channel 6 S. Hardwicke Island North (x)̅

167 126 151 175 367 200 198 (±87)

6.2 0.5

1.4 0.2

0.5 0.1

1.2 0.3

0.4 0.2

3.8 (±2.4)

1.9 1.8 3.8 1.6 3.9 0.4 2.2 (±1.4)

0.8 (±0.6)

0.3 (±0.2)

0.5 0.7 0.6 0.4 0.8 0.2 0.6 (±0.2)

0.2 0.2 0.2 0.2 0.4 0.3 0.2 (±0.1)

7 E. of Oyster River 8 N. Hornby Island 9 E. Gabriola Island 10 Active Pass 11 Patricia Bay South (x)̅

288 180 380 154 185 237 (±95)

1.5 3.5 1.8 0.5 2.9 2.0 (±1.2)

2.6 1.3 2.6 1.6 3.2 2.3 (±0.8)

2.9 5.1 2.6 1.2 14.9 5.3 (±5.5)

0.7 0.7 1.6 0.5 1.5 1.0 (±0.5)

1.2 0.6 0.8 0.5 0.6 0.7 (±0.3)

0.7 0.4 0.6 0.5 0.2 0.5 (±0.2)

0.49

0.19

0.98

0.021

0.035

0.27

0.015

p-value (North vs South)

3.8 4.7

a

Samples were collected along a transect in the northeast Pacific Ocean from Hecate Strait to the Strait of Georgia, British Columbia, Canada. Sites are divided geographically into north and south groupings via delineation by the oceanographically constrained Discovery Passage, with mean ± standard deviation for each region presented in italics.

removing visible biota. Seafloor depths ranged from 126 to 380 m along the transect (Table 1). One sample from each site was placed in a cooler during transport to shore and stored at −20 °C until analyzed for organic carbon and contaminants. Plankton samples were collected as vertical tows from 20 m above the seafloor to the surface using a SCOR plankton net (0.5 m) with 335 μm cod end. The contents of the cod end were collected into solvent-rinsed 500 mL glass jars and stored frozen. A separate tow was conducted, and contents were collected in solvent-rinsed jars and stored with 10% buffered formalin for species identification and plankton counts. Organic Carbon Determination. Sediment organic carbon was determined using a Leemans 440 Elemental Analyzer after removal of carbonates by acid fuming with HCl (1 N) for 24 h in a closed container. Results of replicates of acetanilide standard (71.09% carbon) were 70.61 ± 0.41%, and the pooled standard deviation for sample replicates of total organic carbon (TOC) (%) was 0.03 (n = 3 pairs). Contaminant Analysis. Sediment samples were analyzed by the Laboratory of Expertise in Aquatic Chemical Analysis (LEACA) at the Institute of Ocean Sciences (Sidney, BC), and the plankton samples were analyzed by AXYS Analytical Services (Sidney, BC). Both laboratories analyzed PCBs and PBDEs using high-resolution gas chromatography/high-resolution mass spectrometry according the U.S. Environmental Protection Agency methods 1668 (PCBs) and 1614 (PBDEs). Detailed methodology for sample extraction, lipid determination, cleanup, and quantitation is reported elsewhere.24,25 Sample analysis batches included all samples, two procedural blanks, one standard reference material, and one duplicate sample. For sediment, the average BDE-209 blank concentration (666.1 ± 55.5 pg/g dry weight) fell within the lower range of sample concentrations (421−3608 pg/g dw) and the average surrogate recovery for BDE-209 was 64%. For plankton, BDE-209 detection limits and blank concentrations were greater than sample concentrations of BDE-47 and all other congeners. Plankton BDE-209 concentrations ranged from detection limit (325 pg/g) to 978 pg/g dw (mean: 414.3 ± 249.3 pg/g dw), with four samples below detection, six near

detection, and one above detection. BDE-209 in plankton was highly variable between samples (no spatial pattern) and represented 30 to 80% of sumPBDEs in plankton (mean: 51.1 ± 14.7%). All figures and tables herein include BDE-209 in the sum of PBDEs for sediment samples but exclude BDE-209 in plankton samples due to low detection and poor analytical results in this matrix. PCB and PBDE concentrations in sediment and plankton were all blank-corrected. Sediment concentrations were further corrected with TOC and plankton samples with lipid content for better comparisons between samples. To reduce the influence of nondetected congeners on total concentrations, detection limit substitutions were made for congeners when 70% or more of samples had detectable concentrations. Where the 70% threshold was not reached, those congeners were not considered for further analysis. Data Analysis. Principal component analysis (PCA) was performed using combined PCB data from sediment and plankton using Pirouette software (Infometrix, Inc., Washington, USA). Each congener was evaluated for the number of detection limit substitutions before inclusion in the final PCA data set; congeners were excluded if over 25% of samples were below detection. The final data set included 95 PCB congeners. The centered log-ratio procedure was applied to the data set to avoid the effects of negative bias and closure;22 data were normalized to total PCB concentration followed by geometric mean and then log-transformed and autoscaled. Hierarchal cluster analysis (HCA) was performed on TOCand lipid-corrected sediment and plankton PCB and PBDE concentrations to examine the similarity/differences between sites. Differences between northern and southern regions were determined using independent-samples t test. The BSAF was calculated by dividing the lipid weight concentration of PCB or PBDE in plankton by the TOC-corrected contaminant concentration in sediment.



RESULTS AND DISCUSSION PCB and PBDE Concentrations: Influence of Source History. PCBs and PBDEs were detected in sediment from all

6983

dx.doi.org/10.1021/es500218b | Environ. Sci. Technol. 2014, 48, 6981−6988

Environmental Science & Technology

Article

Figure 2. Principal component analysis (PCA) of 95 PCB congeners in sediment and plankton samples collected along a coastal transect in the northeast Pacific Ocean from Hecate Strait to the Strait of Georgia, British Columbia, Canada. Lower plots illustrate those variables that drive the principal components. Sites are divided geographically into north and south via delineation by the oceanographically constrained Discovery Passage (North = 1−6; South = 7−11).

plankton are relatively uniform over larger spatial scales due to biogeochemical cycling and environmental dilution over time, while spatial PBDE trends reflect those of a more recent-use contaminant. PCB and PBDE Patterns: Influence of Sources, Physicochemical Properties, and Receiving Environments. Dominant PCB congeners in sediment consisted of PCB congeners 138 > 153 > 118 > 101 > 110 > 149. The PCB concentration-weighted log KOW (herein referred as weighted KOW)an average KOW value weighted as a function of congener composition13varied with latitude but did not differ between northern and southern sites (t test, p = 0.10; Supporting Information, Figure S2). Sediments in southern sites had higher contributions of tri- and tetrachlorinated homologue groups, while sediment from the north had higher contributions of hexa- and heptachlorinated homologues (Figure S3). The dominant PCB congeners in plankton consisted of PCB congeners 153 > 138 > 101 > 149 > 99 > 70/76. The geographical pattern in plankton diverged from that of sediment; northern sites were compositionally higher in diand trichlorinated congeners whereas the southern sites had greater hexa- and heptachlorinated congeners (Figure S3). This resulted in a decrease of the concentration-weighted log KOW with latitude (Figure S2). A PCA of 95 PCB congeners in sediments and plankton (Figure 2) clearly illustrates the influence of a number of key processes, including bioaccumulation of select congeners, latitudinal fractionation, and biogeochemical influences on environmental partitioning. The first principal component (p1) related to the BSAF (r2 = 0.68, p < 0.001) and distinguished plankton samples (midweight congeners) from sediment (light and heavy congeners). The second component (p2) contrasted light (2−5 Cl) and heavy congeners (6−8 Cl), which associated

nine stations, and concentrations ranged over several orders of magnitude across sites (Figure 1, Table 1). Cluster analysis (HCA) of PCBs and PBDEs separated sites in the north from the south on the basis of total concentrations, with the delineation occurring at the hydrographically constrained Discovery Passage (Figure 1). Sediment PCB concentrations were 7 and 11 times greater in the south compared to the north for raw and TOC-corrected data respectively, while PBDE concentrations were three and five times greater (Figure 1, Table 1). This was not unexpected, as these trends follow known emissions histories in the NE Pacific. The SoG and adjacent water bodies (Juan de Fuca and Puget Sound) are home to large urban centers, including Vancouver, Victoria, and Seattle, and receive industrial discharges from sewage treatment facilities and pulp and paper mills.10,11,13 In comparison, Hecate Strait and QCS are remote, relatively unpopulated, and have no major documented point sources of contaminants. Cluster analysis of PBDE concentrations in plankton separated northern from southern sites, excluding site 6 which fell separately as an outlier (low lipid content and % carbon). Plankton from southern sites had higher PBDE levels than those in the north (ww t test p = 0.015; lw t test, p = 0.019). In contrast, PCB concentrations did not differ between northern and southern regions (ww t test, p = 0.27; lw t test, p = 0.77). The increasing PBDE−PCB ratio with proximity to wastewater discharge sources in the SoG reflects the contrasting usage history of the two contaminants (Supporting Information, Figure S1). In the SoG, PBDE tissue levels in plankton were found to rapidly respond to ongoing PBDE inputs to the water column via wastewater effluent, whereas PCBs appeared to be less available due to high particle scavenging of already low background concentrations.12 Our results are consistent with this model in the SoG and further show that PCB levels in 6984

dx.doi.org/10.1021/es500218b | Environ. Sci. Technol. 2014, 48, 6981−6988

Environmental Science & Technology

Article

Biota-Sediment Accumulation Factors (BSAF). Given the persistence of PCBs and PBDEs in the environment, evaluating their biota-sediment equilibrium status (BSAF ≈ 1− 2) can provide valuable context on accumulation history and enable the assessment of chemical fate in the ocean.6 For both PCBs and PBDEs, observed BSAFs were at, or above, the equilibrium reference value for northern sites (PCB, 1.7 ± 1.3; PBDE, 1.3 ± 1.0) and below equilibrium in southern sites (PCB, 0.2 ± 0.1; PBDE, 0.4 ± 0.4) (Figure 3). The difference

with latitudinal groupings of sediment and plankton samples (north vs south). Although p2 correlated with latitude for both sediment (r2 = 0.52, p = 0.012) and plankton (r2 = 0.59, p = 0.016), the direction of the relationship was inverted: sediments in the southern and northern groupings plotted with light and heavy congeners, respectively, while the opposite occurred for plankton (Figure, 2). It is interesting to note that proximity to source (i.e., in this case, latitude) had little effect on PCB concentrations in plankton but did influence PCB patterns. The lighter PCB pattern with latitude in plankton is likely the result of the combined effect of fractionation during atmospheric transport and biogeochemical processes, while biogeochemical processes alone appear to drive the pattern in sediment. This divergent latitudinal pattern in sediment and plankton likely reflects the difference in response time as well as the contrasting influence of sediment processes on each matrix. Plankton are more strongly associated with patterns in the atmosphere which respond quickly to changes in emissions, while sediments respond much more slowly and can retain legacy patterns.6 We propose that the heavier PCB pattern in plankton of the SoG is a response to continuous terrestrial inputs (i.e., Fraser River and anthropogenic sources) which act to replenish PCBs in the water column, whereas the lack of replenishment in Hecate Strait/QCS contributes to the depletion of mid- and heavy-weight congeners in plankton. Source history and biogeochemical processes have the opposite effect on patterns in sediment; new terrestrial inputs of PCBs are mixing with legacy deposits in the SoG (e.g., bioturbation) and dilute the historically heavier profile, while this legacy signal is retained in the north as a consequence of minimal new inputs and continued erosion of seafloor sediments. Overall, the importance of sediment−contaminant property interactions in governing the PCB fate is highlighted by the relationship between log KOW (i.e., hydrophobicity) and congener patterns in both matrices (log KOW vs p2: r2 = 0.33, p < 0.001) (Figure 2). A total of 49 PBDE congeners were detected in sediments, with BDE-209 contributing 66 ± 18% of ∑PBDE on average. BDE-47 (10.5 ± 6.6%), BDE-99 (7.0 ± 7.7), and BDE-100 (2.5 ± 2.2) were the next most dominant congeners. Forty congeners were detected in plankton, and the dominant congeners were BDE-47 (36.4 ± 9.9%), BDE-99 (27.9 ± 7.1%), and BDE-100 (7.4 ± 1.7%). The dominance of these congeners in sediment and biota correspond well with patterns observed in other studies and reflect the partitioning of congeners over a log KOW range of 4.5−7.5.10,12,26 PBDE patterns in sediment did not follow a clear geographical trend (Supporting Information, Figures S2 & S3). This lack of a trend for PBDEs across coastal BC highlights a dichotomy with PCBs and underscores the influence of physicochemical properties on the fate of these two contaminants. PBDE fate is strongly linked to particle processes and organic carbon content,10,11,13,27 and indeed we found that patterns in sediment (weighted log KOW) were strongly related to sediment TOC (r2 = 0.53 p = 0.025). In plankton, the PBDE-weighted log KOW increased with latitude and was negatively correlated with water depth (r2 = 0.44 p = 0.036) (Supporting Information, Figure S2). The relatively lighter congener patterns in southern plankton may result from the sedimentation and burial of heavier congeners in the SoG and/ or the influence of particle composition and transport on the fate of hydrophobic PBDEs in the marine environment.

Figure 3. Biota−sediment accumulation factors (BSAFs) for PCBs and PBDEs at sites along a coastal transect in the northeast Pacific Ocean from Hecate Strait to the Strait of Georgia, British Columbia, Canada. The vertical dotted line separates northern from southern sites, which are delineated by the oceanographically constrained Discovery Passage. The horizontal dashed line represents theoretical equilibrium conditions between plankton and sediment (BSAF = 1). BSAFs were calculated using lipid- and TOC-corrected concentrations for plankton and sediment, respectively.

between the northern and the southern portion of the transect clearly illustrates the effect of source-driven vs atmospherically driven systems. Continued terrestrial and point-source release of pollutants in the SoG has resulted in a source-driven system where sediments are a major sink of PCBs and PBDEs. This is compounded by intense particle scavenging which can reduce the bioavailability of contaminants in the water column.28 The dominant input of pollutants for northern sites is presumably via atmospheric deposition, and combined with low sedimentation rates and strong mixing of the water column,23 it appears as though both PCBs and PBDEs have reached some level of equilibrium between sediment and biota. In other words, the sediment compartment does not act as a strong contaminant source to plankton, and vice versa. Deviation from equilibrium was not only a function of source-proximity but was also influenced by physicochemical properties (Figure 4). Congener-specific BSAFs were correlated with Log KOW, with a parabolic fit best describing the relationship (PCB: r2=0.18, p < 0.001; PBDE: r2=0.61, p < 0.001). BSAFs peaked for congeners with log KOW ≈ 6−7. Unlike PCBs, there was no dramatic decline in BSAFs for the heaviest PDBEs (Figure 4); however, BDE-209 was excluded in 6985

dx.doi.org/10.1021/es500218b | Environ. Sci. Technol. 2014, 48, 6981−6988

Environmental Science & Technology

Article

local source history and characteristics of the receiving environment is essential for accurate modeling of contaminant fate in marine environments. The definition of BSAF includes a biotic as well as a sediment term and resulting BSAFs are therefore sensitive to the features of both matrices. Plankton collected in this study were dominated by calanoid copepods and euphausiids, representing 70 and 8.2% of all organisms, respectively. While plankton composition in northern and southern sites were similar and overall lipid content did not differ between regions (t test; p = 0.98), species-specific differences in growth rates, lipid composition, feeding rate, feeding preferences, and uptake and elimination may have contributed to some of the variability in resulting contaminant patterns.12,35 Biotransformation of some PCB and PBDE congeners, through oxidation or debromination, has been observed in invertebrates and other bottom food web organisms.36−38 However, exposure from water and diet appear to more strongly drive chemical accumulation in plankton while biotransformation becomes more important in higher trophic level consumers.39 For sediment, although no general TOC pattern was evident along the transect (Table 1), site-specific differences in particle and organic carbon composition can influence BSAFs.27 For instance, high affinity carbonaceous geosorbents, such as black carbon, strongly bind to hydrophobic contaminants and considerably reduce their bioavailability.29 Models that include feeding parameters and sediment composition have been shown to accurately predict bioaccumulation of contaminants from sediment.29,30,35,40 In a north−south transect spanning much of coastal BC, we show that highest BSAFs do not correspond with areas receiving the highest contaminant inputs. Nonetheless, higher trophic organisms accumulate pollutants in proportion to absolute concentrations. Despite higher BSAFs in remote waters, higher PBDEs but not PCBs in plankton near-source highlight the dichotomy between current-use and legacy pollutants. This trend is extended up the food web to harbor seals which have greater blubber levels of PBDEs but not PCBs in the SoG compared to Queen Charlotte Strait; however, the much longer lived killer whales appear to retain the legacy PCB spatial profile.41 Evidence from plankton in this study has shown that while the bioavailability of PCBs has become spatially uniform in coastal BC, the risk of food web magnification of PBDEs is explicitly higher in the industrialized southern waters.

Figure 4. Biota−sediment accumulation factors (BSAFs) for PCB and PBDE congeners are greater for northern sites than southern sites and correlate with log KOW. The horizontal black line represents theoretical equilibrium conditions between plankton and sediment (BSAF = 1). BSAF values were calculated using lipid- and TOC-corrected concentrations for plankton and sediment, respectively.

plankton due to high detection limits (see Methods) but trends did reveal a decline in biotic uptake relative to octa- and nonaBDEs (not shown). High BSAFs for those congeners expected to not be bioavailable12,27 may be due to plankton stomach contents or cell adsorption rather than retained tissue burden. The parabolic log KOW curve has been widely reported in studies of biota−water accumulation,12 BSAFs,29,30 maternal transfer to offspring,31 and biomagnification in food webs,32 making evident the universal influence of physicochemical properties on the fate of POPs. All but two PBDE congeners were above the equilibrium reference value (i.e., BSAF > 1) in the north and approximately half were above in the south. As with total BSAF, the BSAF vs log KOW slopes were higher for PBDEs than PCBs for both regions (Figure 4). Annual atmospheric deposition in urban and remote regions of coastal BC delivers over twice as much PBDEs as PCBs.33 Therefore, despite shared physicochemical properties for PCBs and PBDEs, the greater input of PBDEs to the environment results in increased uptake at the base of the food web. The BSAF can be used predictively in models to forward calculate contaminant concentrations in biota from the concentration in sediment (i.e., CBiota = BSAF × CSed) or can be used to back calculate sediment concentrations protective of biological thresholds (i.e., CSed = CBiota/BSAF).8,34 However, the utility of predictive models using literature BSAF values will depend on the applicability of values from one site to another. Our findings using oceanic plankton indicate that total and congener-specific BSAF vary among sites as a function of source history, physicochemical properties, sedimentation rates, sediment composition, and water depth. Thus, incorporating



ASSOCIATED CONTENT

S Supporting Information *

Additional information on spatial patterns of contaminant ratios, concentration weighted log KOW and homologue groupings. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*Phone: 604-659-3563, fax: 604-659-3562, e-mail: peter.ross@ vanaqua.org. Notes

The authors declare no competing financial interest. 6986

dx.doi.org/10.1021/es500218b | Environ. Sci. Technol. 2014, 48, 6981−6988

Environmental Science & Technology



Article

(15) Meire, R. O.; Targino, A. C.; Torres, J. P. M. Bulk Atmospheric Deposition of Persistent Toxic Substances (PTS) along Environmental Gradients in Brazil. Environ. Sci. Pollut. Res. Int. 2013, 20, 2046−2053. (16) Gioia, R.; Nizzetto, L.; Lohmann, R.; Dachs, J.; Temme, C.; Jones, K. C. Polychlorinated Biphenyls (PCBs) in Air and Seawater of the Atlantic Ocean: Sources, Trends and Processes. Environ. Sci. Technol. 2008, 42, 1416−1422. (17) Gioia, R.; Lohmann, R.; Dachs, J.; Temme, C.; Lakaschus, S.; Schulz-Bull, D.; Hand, I.; Jones, K. C. Polychlorinated Biphenyls in Air and Water of the North Atlantic and Arctic Ocean. J. Geophys. Res. 2008, 113, D19302. (18) Meijer, S. N.; Steinnes, E.; Ockenden, W. A.; Jones, K. C. Influence of Environmental Variables on the Spatial Distribution of PCBs in Norwegian and U. K. Soils: Implications for Global Cycling. Environ. Sci. Technol. 2002, 36, 2146−2153. (19) Calamari, D.; Tremolada, P.; Di Guardo, A.; Vighi, M. Chlorinated Hydrocarbons in Pine Needles in Europe: Fingerprint for the Past and Recent Use. Environ. Sci. Technol. 1994, 28, 429−434. (20) Davidson, D. A.; Wilkinson, A. C.; Blais, J. M.; Kimpe, L. E.; McDonald, K. M.; Schindler, D. W. Orographic Cold-Trapping of Persistent Organic Pollutants by Vegetation in Mountains of Western Canada. Environ. Sci. Technol. 2003, 37, 209−215. (21) Muir, D. C.; Ford, C. A.; Grift, N. P.; Metner, D. A.; Lockhart, W. L. Geographic Variation of Chlorinated Hydrocarbons in Burbot (Lota Lota) from Remote Lakes and Rivers in Canada. Arch. Environ. Contam. Toxicol. 1990, 19, 530−542. (22) Ross, P. S.; Jeffries, S. J.; Yunker, M. B.; Addison, R. F.; Ikonomou, M. G.; Calambokidis, J. C. Harbor Seals (Phoca Vitulina) in British Columbia, Canada, and Washington State, USA, Reveal a Combination of Local and Global Polychlorinated Biphenyl, Dioxin, and Furan Signals. Environ. Toxicol. Chem. 2004, 23, 157−165. (23) Zacharias, M. A.; Howes, D. E.; Harper, J. R.; Wainwright, P. The British Columbia Marine Ecosystem Classification: Rationale, Development, and Verification. Coast. Manage. 1998, 26, 105−124. (24) Christensen, J. R.; MacDuffee, M.; Macdonald, R. W.; Whiticar, M.; Ross, P. S. Persistent Organic Pollutants in British Columbia Grizzly Bears: Consequence of Divergent Diets. Environ. Sci. Technol. 2005, 39, 6952−6960. (25) Ikonomou, M.; Fraser, T.; Crewe, N.; Fischer, M. B.; Rogers, I. H.; He, T.; Sather, P.; Lamb, R. A Comprehensive Multiresidue Ultratrace Analytical Method, Based on HRGC/HRMS, for the Determination of PCDDs, PCDFs, PCBs, PBDEs, PCDEs, and Organochlorine Pesticides in Six Different Environmental Matrices. Can. Technol. Reports Fish. Aquat. Sci. 2001, 2389, 1−95. (26) DeWit, C. A. An Overview of Brominated Flame Retardants in the Environment. Chemosphere 2002, 46, 583−624. (27) Dinn, P. M.; Johannessen, S. C.; Ross, P. S.; Macdonald, R. W.; Whiticar, M. J.; Lowe, C. J.; Roodselaar, A. Van. PBDE and PCB Accumulation in Benthos near Marine Wastewater Outfalls: The Role of Sediment Organic Carbon. Environ. Pollut. 2012, 171, 241−248. (28) Fisher, J. B.; Petty, R. L.; Lick, W. Release of Polychlorinated Biphenyls from Contaminated Lake Sediments: Flux and Apparent Diffusivities of Four Individual PCBs. Environ. Pollut. 1983, 5, 121− 132. (29) Moermond, C. T. A.; Zwolsman, J. J. G.; Koelmans, A. A. Black Carbon and Ecological Factors Affect in Situ Biota to Sediment Accumulation Factors for Hydrophobic Organic Compounds in Flood Plain Lakes. Environ. Sci. Technol. 2005, 39, 3101−3109. (30) Thomann, R. V.; Connolly, J. P.; Parkerton, T. F. An Equilibrium Model of Organic Chemical Accumulation in Aquatic Food Webs with Sediment Interactions. Environ. Toxicol. Chem. 1992, 11, 615−629. (31) Desforges, J.-P. W.; Ross, P. S.; Loseto, L. L. Transplacental Transfer of Polychlorinated Biphenyls and Polybrominated Diphenyl Ethers in Arctic Beluga Whales (Delphinapterus Leucas). Environ. Toxicol. Chem. 2012, 31, 296−300. (32) Fisk, A. T.; Hobson, K. A.; Norstrom, R. J. Influence of Chemical and Biological Factors on Trophic Transfer of Persistent

ACKNOWLEDGMENTS This study was supported by Fisheries and Oceans Canada (DFO) and the Georgia Basin Ecosystem Initiative of Environment Canada. Thanks to the crew of the CCGS Vector, the staff at the former Laboratory of Expertise for Aquatic Chemical Analysis (LEACA) in Sidney, BC (DFO) and to Axys Analytical Services in Sidney, BC. The authors gratefully acknowledge comments from Robie W. Macdonald, Sophia Johannessen, and Moira Galbraith.



REFERENCES

(1) Wania, F.; Mackay, D. The Evolution of Mass Balance Models of Persistent Organic Pollutant Fate in the Environment. Environ. Pollut. 1999, 100, 223−240. (2) Wania, F.; Mackay, D. A Global Distribution Model for Persistent Organic Chemicals. Sci. Total Environ. 1995, 160−161, 211−232. (3) Dachs, J.; Lohmann, R.; Ockenden, W. A.; Méjanelle, L.; Eisenreich, S. J.; Jones, K. C. Oceanic Biogeochemical Controls on Global Dynamics of Persistent Organic Pollutants. Environ. Sci. Technol. 2002, 36, 4229−4237. (4) Galbán-Malagón, C. J.; Berrojalbiz, N.; Gioia, R.; Dachs, J. The “Degradative” and “Biological” Pumps Controls on the Atmospheric Deposition and Sequestration of Hexachlorocyclohexanes and Hexachlorobenzene in the North Atlantic and Arctic Oceans. Environ. Sci. Technol. 2013, 47, 7195−7203. (5) Morrison, H. A.; Gobas, F. A.; Lazar, R.; Haffner, G. D. Development and Verification of a Bioaccumulation Model for Organic Contaminants in Benthic Invertebrates. Environ. Sci. Technol. 1996, 30, 3377−3384. (6) Li, Y.-F.; Harner, T.; Liu, L.; Zhang, Z.; Ren, N.-Q.; Jia, H.; Ma, J.; Sverko, E. Polychlorinated Biphenyls in Global Air and Surface Soil: Distributions, Air-Soil Exchange, and Fractionation Effect. Environ. Sci. Technol. 2010, 44, 2784−2790. (7) Ross, P.; Ellis, G. M.; Ikonomou, M. G.; Barrett-Lennard, L. G.; Addison, R. F. High PCB Concentrations in Free-Ranging Pacific Killer Whales, Orcinus Orca: Effects of Age, Sex and Dietary Preference. Mar. Pollut. Bull. 2000, 40, 504−515. (8) Alava, J. J.; Ross, P. S.; Lachmuth, C.; Ford, J. K. B.; Hickie, B. E.; Gobas, F. A. P. C. Habitat-Based PCB Environmental Quality Criteria for the Protection of Endangered Killer Whales (Orcinus Orca). Environ. Sci. Technol. 2012, 46, 12655−12663. (9) DeWit, C. A.; Herzke, D.; Vorkamp, K. Brominated Flame Retardants in the Arctic EnvironmentTrends and New Candidates. Sci. Total Environ. 2009, 408, 2885−2918. (10) Johannessen, S. C.; Macdonald, R. W.; Wright, C. A.; Burd, B.; Shaw, D. P.; van Roodselaar, A. Joined by Geochemistry, Divided by History: PCBs and PBDEs in Strait of Georgia Sediments. Mar. Environ. Res. 2008, 66 (Suppl), S112−20. (11) Dinn, P. M.; Johannessen, S. C.; Macdonald, R. W.; Lowe, C. J.; Whiticar, M. J. Effect of Receiving Environment on the Transport and Fate of Polybrominated Diphenyl Ethers near Two Submarine Municipal Outfalls. Environ. Toxicol. Chem. 2012, 31, 566−573. (12) Frouin, H.; Dangerfield, N.; Macdonald, R. W.; Galbraith, M.; Crewe, N.; Shaw, P.; Mackas, D.; Ross, P. S. Partitioning and Bioaccumulation of PCBs and PBDEs in Marine Plankton from the Strait of Georgia, British Columbia, Canada. Prog. Oceanogr. 2013, 115, 65−75. (13) Grant, P. B. C.; Johannessen, S. C.; Macdonald, R. W.; Yunker, M. B.; Sanborn, M.; Dangerfield, N.; Wright, C.; Ross, P. S. Environmental Fractionation of PCBs and PBDEs during Particle Transport as Recorded by Sediments in Coastal Waters. Environ. Toxicol. Chem. 2011, 30, 1522−1532. (14) Ross, P. S.; Couillard, C. M.; Ikonomou, M. G.; Johannessen, S. C.; Lebeuf, M.; Macdonald, R. W.; Tomy, G. T. Large and Growing Environmental Reservoirs of Deca-BDE Present an Emerging Health Risk for Fish and Marine Mammals. Mar. Pollut. Bull. 2009, 58, 7−10. 6987

dx.doi.org/10.1021/es500218b | Environ. Sci. Technol. 2014, 48, 6981−6988

Environmental Science & Technology

Article

Organic Pollutants in the Northwater Polynya Marine Food Web. Environ. Sci. Technol. 2001, 35, 732−738. (33) Noël, M.; Dangerfield, N.; Hourston, R. S.; Belzer, W.; Shaw, P.; Yunker, M. B.; Ross, P. S. Do Trans-Pacific Air Masses Deliver PBDEs to Coastal British Columbia, Canada? Environ. Pollut. 2009, 157, 3404−3412. (34) Gobas, F. A. P. C.; Arnot, J. A. Food Web Bioaccumulation Model for Polychlorinated Biphenyls in San Francisco Bay, California, USA. Environ. Toxicol. Chem. 2010, 29, 1385−1395. (35) Wong, C. S.; Capel, P. D.; Nowell, L. H. National-Scale Evaluation of the Biota−Sediment Accumulation Factor Model. Environ. Sci. Technol. 2001, 35, 1709−1715. (36) Tian, S.; Zhu, L.; Liu, M. Bioaccumulation and Distribution of Polybrominated Diphenyl Ethers in Marine Species from Bohai Bay, China. Environ. Toxicol. Chem. 2010, 29, 2278−2285. (37) Bartrons, M.; Grimalt, J. O.; de Mendoza, G.; Catalan, J. Pollutant Dehalogenation Capability May Depend on the Trophic Evolutionary History of the Organism: PBDEs in Freshwater Food Webs. PLoS One 2012, 7, e41829. (38) Hakk, H.; Letcher, R. J. Metabolism in the Toxicokinetics and Fate of Brominated Flame RetardantsA Review. Environ. Int. 2003, 29, 801−828. (39) Borgå, K.; Fisk, A. T.; Hoekstra, P. E.; Muir, D. C. G. Biological and Chemical Factors of Importance in the Bioaccumulation and Trophic Transfer of Persistent Organochlorine Contaminants in Arctic Marine Food Webs. Environ. Toxicol. Chem. 2004, 23, 2367−2385. (40) Mackenbach, E. M.; Jing, Y.; Mills, M. A.; Landrum, P. F.; Lydy, M. J. Application of a Tenax Model to Assess Bioavailability of PCBs in Field Sediments. Environ. Toxicol. Chem. 2012, 31, 2210−2216. (41) Ross, P. S. Fireproof Killer Whales (Orcinus Orca): FlameRetardant Chemicals and the Conservation Imperative in the Charismatic Icon of British Columbia, Canada. Can. J. Fish. Aquat. Sci. 2006, 63, 224−234.

6988

dx.doi.org/10.1021/es500218b | Environ. Sci. Technol. 2014, 48, 6981−6988