Time Trends of Persistent Organic Pollutants in Benthic and Pelagic Indicator Fishes from Puget Sound, Washington, USA

We modeled temporal trends in polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), and dichlorodiphenyltrichloroethane and its metabolites (DDTs) in two indicator fish species representing benthic and pelagic habitats in Puget Sound, Washington, USA. English sole (Parophrys vetulus, benthic) index sites and larger-scale Pacific herring (Clupea pallasii, pelagic) foraging areas represented a wide range of possible contamination conditions, with sampling locations situated adjacent to watersheds exhibiting high, medium and low development. Consistency in analytical data throughout the study was maintained by either calculating method-bias-correction factors on paired samples as methods evolved or by analyzing older archived samples by current methods. PCBs declined moderately in two herring stocks from a low-development basin (2.3 and 4.0% annual rate of decline) and showed no change in the highly developed and moderately developed basins during a 16- to 21-year period. PCBs increased in English sole from four of ten sites (2.9–7.1%), and the remaining six exhibited no significant change. PBDEs and DDTs declined significantly in all herring stocks (4.2–8.1%), although analytical challenges warrant caution in interpreting DDT results. PBDEs declined in English sole from two high-development and one low-development site (3.7–7.2%) and remained unchanged in the remaining seven. DDTs increased in English sole from one high-development site (Tacoma City Waterway) and declined in two high-development and one low development site. As with herring, analytical challenges warrant caution in interpreting the English sole DDT results. It is likely that source controls and mitigation efforts have contributed to the declines in PBDEs and DDTs overall, whereas PCBs appear to have persisted, especially in the pelagic food web, despite bans in PCB production and use.

Abstract We modeled temporal trends in polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), and dichlorodiphenyltrichloroethane and its metabolites (DDTs) in two indicator fish species representing benthic and pelagic habitats in Puget Sound, Washington, USA. English sole (Parophrys vetulus, benthic) index sites and larger-scale Pacific herring (Clupea pallasii, pelagic) foraging areas represented a wide range of possible contamination conditions, with sampling locations situated adjacent to watersheds exhibiting high, medium and low development. Consistency in analytical data throughout the study was maintained by either calculating method-bias-correction factors on paired samples as methods evolved or by analyzing older archived samples by current methods. PCBs declined moderately in two herring stocks from a low-development basin (2.3 and 4.0% annual rate of decline) and showed no change in the highly developed and moderately developed basins during a 16-to 21-year period. PCBs increased in English sole from four of ten sites (2.9-7.1%), and the remaining six exhibited no significant change. PBDEs and DDTs declined significantly in all herring stocks (4.2-8.1%), although analytical challenges warrant caution in interpreting DDT results. PBDEs declined in English sole from two high-development and one low-development site (3.7-7.2%) and remained unchanged in the remaining seven. DDTs increased in English sole from one high-development site (Tacoma City Waterway) and declined in two high-development and one low development site. As with herring, analytical challenges warrant caution in interpreting the English sole DDT results. It is likely that source controls and mitigation efforts have contributed to the declines in PBDEs and DDTs overall, whereas PCBs appear to have persisted, especially in the pelagic food web, despite bans in PCB production and use.
Although the production and use of a number of persistent organic pollutants (POPs) has been curtailed to protect environmental and human health in many countries (Tierney et al. 2014), it can be difficult to evaluate the success of these efforts over the long term, especially in aquatic environments. POPs include a wide variety of toxic industrial chemicals such as polychlorinated biphenyls (PCBs) and flame-retardant polybrominated diphenyl ethers (PBDEs), as well as chlorinated pesticides such as dichlorodiphenyltrichloroethane and its metabolites (DDTs), chlordanes, and hexachlorobenzene. These POPs can accumulate in sediments and biota where they may persist for many years (Jones and de Voogt 1999;Tierney et al. 2014). Earliest efforts to measure the level of POP contamination in marine animals (DDTs and PCBs; Jensen et al. 1969) combined with more recent time trend analyses suggest that many aquatic systems exhibited a substantial drop in DDTs and PCBs through the 1970s and 1980s [e.g., the Canadian arctic (Braune et al. 2005), the Baltic Sea (Bignert et al. 1998), and the Great Lakes of North America (Stow et al. 1994)]. Variation in biological covariates, such as fish age, lipid levels, sex, and trophic level, can confound analysis of temporal and spatial trends in biota (West et al. 2008), as well as comparisons between studies (this report). Moreover, it is difficult to maintain consistency in measuring POPs as analytical methods evolve over time.
The production and use of DDTs and PCB technical mixtures in the United States were sharply curtailed from 1972 to 1979, resulting from federal legislation banning these chemicals. In addition, Washington State passed legislation in 2008 to restrict the use of polybrominated diphenylether (PBDE) flame retardants. Washington State monitors the concentration of these three chemical groups in two marine indicator fish species representing benthic (English sole, Parophrys vetulus) and pelagic (Pacific herring, Clupea pallasii) habitats as part of a directed effort to track the health the Puget Sound ecosystem. Tissue residue recovery targets using the concept of critical POP body residues (Meador et al. 2011) are used to evaluate the success of pollution remediation efforts, on a large, ecosystem scale. The two Puget Sound fish contaminant indicators (species) reported here are part of a larger portfolio of Vital Sign indicators designed to simplify science reporting and link it to regional policy, ultimately resulting in easily communicated policy statements or targets that define the desired condition or goals (Puget Sound Partnership 2016).
Puget Sound comprises a complex, fjord-like, interconnected series of inland marine and estuarine waters in Washington State (Burns 1985; Moore et al. 2008) that exchange with the Pacific Ocean primarily via a relatively narrow waterway, the Strait of Juan de Fuca (Fig. 1). Oceanographic stratification and geographically restricted connections between deep Puget Sound basins result in recirculation of basin waters (Thomson 1994), which increases retention of regionally produced pollutants (Harrison et al. 1994). These inland marine waters drain watersheds that range from highly developed or urbanized with large cities (central Puget Sound) to moderately developed (southern Puget Sound), to rural, undeveloped lands (northern Puget Sound), resulting in marine waters that range from highly polluted to nearly pristine.
Based on results from these studies, POPs in these species illustrate broad, basin-wide contaminant conditions in the pelagic food web. Ross et al. (2013) identified strongly declining trends in concentrations of PCBs from 1984 to 2010 and one other POP class (polychlorinated naphthalenes) from 1984 to 2003 in harbor seals from Southern Puget Sound. PBDEs in the same harbor seals showed a different temporal pattern, with a strong increase from 1984 to 2003 and a subsequent 40% decline in 2009. More recent POP trends in Puget Sound marine mammals remain unknown, although Hickie et al. (2007) predicted it would take between 14 and 57 years for PCB concentration to fall below an effects threshold of 17 mg total PCBs/kg blubber lipids for the populations of fish-eating killer whales that typically forage in Puget Sound.
Contamination of pelagic species with POPs has occurred at a large, basin scale in Puget Sound, whereas POP contamination of benthic and sessile species appears to be limited to urbanized areas, which is where the greatest sediment contamination occurs (Long et al. 2005). This is evident from POPs reported in blue mussels (Mytilus trossulus; Lanksbury et al. 2013) and English sole in Puget Sound (this report and Ylitalo et al. 1999), and consistent with other species from the California Coast (Dodder et al. 2013) and San Francisco Bay (Davis et al. 2007;Greenfield and Allen 2013).
The Washington Department of Fish and Wildlife (WDFW) has tracked POPs in indicator species representing key food-web pathways including English sole, a bottom-dwelling flatfish, as well as Pacific herring, a small, schooling mid-water planktivore, or forage fish. English sole consume infaunal invertebrates and so represent the sediment-to-biota contaminant link for contaminants that tend to accumulate in sediments. Adult Pacific herring primarily consume zooplankton and are themselves prey to virtually every large piscivorous species in Puget Sound, and so generally represent the pathway of contaminants to higher level predators including Pacific salmon, seabirds, harbor seals and killer whales.
A significant challenge for long-term toxics monitoring programs is collecting data that are comparable over long time periods either by maintaining consistent analytical methods through time, or by conducting comparison studies as methods change. This includes selection of extraction solvents and extraction methods, and subtle changes in sample handling. In addition, comparisons of contaminant levels among studies can be difficult because differing methods or summation algorithms are often used to estimate total PCBs, PBDEs, or other POP groups. Monitoring programs often rely on subsets of congeners or analytes that are readily measured using low-resolution methods, to avoid the high cost of measuring all 209 congeners using more comprehensive analytical methods such as high-resolution gas chromatography/mass spectrometry.
In the current study, we report concentrations of PCBs, PBDEs and DDTs in two primary indicator species, English sole and Pacific herring, from Puget Sound,  Washington over a 16-to 21-year time period. Our objectives were (1) to apply correction factors to POP data to mitigate biases related to changes in analytical methodology and quantitation methods as techniques evolved during the monitoring periods, (2) to evaluate whether contaminant conditions are improving or worsening in Puget Sound, (3) to compare current contaminant conditions for these species for ten index sites (English sole) and three oceanographic basins (Pacific herring) representing a wide range of contamination, (4) to compare PCB summations used in our monitoring program with P ICES 7 , a PCB congener summation method sanctioned by the International Council for the Exploration of the Sea, and used by many European monitoring programs, (5) to compare our POP results with published critical body residues (see, Meador et al. 2011) to infer the likelihood of toxicopathic health effects on fish, and (6) to identify cases where POP body residues in fish were great enough to generate consumption advisories to protect human health.

Study Area and Sample Collection
Study areas for both indicator species were chosen to satisfy a wide range of needs for long-term monitoring of contaminants in Puget Sound's pelagic and benthic food webs (Fig. 1). Three of the most abundant stocks of winterspawning Pacific herring were selected to represent the pelagic food web in three oceanographic basins characterized by a wide range of land-development in their watersheds. The Port Orchard stock was sampled from central Puget Sound, which is characterized by highly developed or urbanized watersheds (termed high-development). The Semiahmoo stock was sampled from the southern Strait of Georgia, which is characterized by rural watersheds with little development, and a relatively direct connection to oceanic waters (termed low-development). The Squaxin stock was sampled from south Puget Sound, which exhibited relatively low land development, however its connection to oceanic waters was through the highly developed central Puget Sound, so its level of development was characterized as moderate. A fourth herring stock, Cherry Point, also occurring in the low-development southern Strait of Georgia was tracked because its unique springtime spawn-timing separates it from the other stocks, and its long history of precipitous decline has made it a priority for protection and recovery. All herring were collected in the winter (January through March) as they aggregated for spawning in predictable locations, except for Cherry Point herring, which were taken in May or June. Fish were captured using gill-nets along shores where spawning activity was occurring, or using midwater trawls in nearby off-shore waters where pre-spawners aggregated. To minimize the effects of reproduction and fish age on POP concentration, only male herring were selected, from sizes that were likely to be 3-year-old fish. Fifty male herring were sampled from each location and year. Ten composite samples, each containing five male herring were created by homogenizing whole bodies together using a food-grade stainless steel and aluminum grinder.
Ten English sole index sites were selected to represent a wide range of potential contaminant inputs for this benthic species, across the full geographic range of Puget Sound (Fig. 1). Degree of impervious surface in watershed uplands was used as a proxy for land development and also to classify index sites as follows: high development (within 500 m of land with [50% impervious surface), moderate development (within 500 m of land with[25% impervious surface or within 2000 m of land with [50% impervious surface), or low development (more than 2000 m away from land with [25% impervious surface). High development sites were typically close to shore (\500 m), whereas some moderate-and low-development sites were more than 2000 m from the shoreline used to classify them. Watershed impervious surface was estimated using the % imperviousness data layer from the 2006 National Land Cover Database (Fry et al. 2011;Wickham et al. 2013) for closest catchment areas draining into Puget Sound.
Although English sole are known to make long-distance movements for spawning in winter months before homing back to their feeding ground in the spring and summer months (Moser et al. 2013), they are typically sedentary during spring and summer months (Day 1976), where they exhibit small, predictable foraging ranges (O'Neill et al. 2007). Because they are thought to feed little during the winter (Day 1976), their POP tissue burdens likely represent local conditions of their spring/summer habitats. All English sole were captured using a bottom trawl, and skinoff muscle tissue was later resected from either fresh/iced or frozen/thawed specimens. An equal mass of muscle (typically 10 g) was resected from each of 60-120 specimens per site per year and combined randomly to create three or six composite groups of 10-20 fish each. Composited tissues were homogenized using hand-held stainless steel electric mixers.

Chemical Analyses
Two analytical methods, two extraction methods, and two extraction solvents were used to measure POPs in English sole and Pacific herring samples over the duration of this study. The most recently collected samples were analyzed by gas chromatography/mass spectrometry (GC/MS) and accelerated solvent extraction (ASEx) with methylene chloride (MeCl) according to Sloan et al. (2014). Specifications for instruments, supplies and materials used in all GC/MS runs are detailed in Sloan et al. (2004Sloan et al. ( , 2014. Earlier collected samples were analyzed via high-performance liquid chromatography with photodiode array (HPLC/PDA) detection according to Krahn et al. (1994). The HPLC/PDA analytical method was conducted with either manual solvent extraction (MSEx) and equal volumes of pentane and hexane (PenHex), or ASEx using MeCl.
Tissue samples were analyzed using GC/MS with ASEx according to Sloan et al. (2014) from 2005 to 2015 for PCBs and DDTs and from 1994 to 2015 for PBDEs (using a combination of real-time and archived samples). In brief, this method comprised three steps: (1) ASEx of tissue using MeCl, (2) two-step cleanup of the MeCl extract by silica/aluminum columns and size-exclusion high-performance liquid chromatography (SEC HPLC) to remove lipids and other biogenic compounds, and (3) quantitation of POPs using GC/MS with selected-ion-monitoring (SIM). The ASEx provided an extract that was used for analyte recovery and gravimetric lipid evaluation. Accuracy of the instrument was improved by including chemical ionization filaments (used to increase source temperature), a cool oncolumn injection system in the GC, and a guard column before the analytical column. Point-to-point calibration was used improve data fit over the full range of GC/MS calibration standards (Sloan et al. 2004). Forty-seven PCBs, 11 PBDEs, and 5 DDTs were routinely detected in all samples.
Total PCBs analyzed by GC/MS were estimated using a simple algorithm based on a subset of 17 commonly detected PCB congeners representing homologs containing three to ten chlorine atoms (IUPAC numbers 18,28,44,52,66,101,105,118,128,138,153,170,180,187,195,206,and 209) wherein the sum of detected values for these 17 congeners was multiplied by 2 to estimate total PCBs. This algorithm was developed in the context of a US nationwide toxics monitoring program from empirical data (Lauenstein and Cantillo 1993), as a cost-effective way to estimate total PCBs from a small subset of congeners using relatively inexpensive low-resolution analytical techniques. Total PCBs were estimated by Lauenstein and Cantillo (1993) using 18 congeners (those above, plus PCB 8); however, PCB 8 was not quantitated in the current study, so the algorithm employed herein summed only 17 congeners, with the estimation of total PCBs hereafter referred to as 2 P 17 PCBs. Although PCB 8 was detected in 30 samples of fish tissue in another study of pelagic Puget Sound fish using high-resolution GC/MS methods (West et al. 2011b), concentrations ranged only from 0.0014 to 0.012 ng/g wet wt, suggesting that it would be undetectable by the methods used in this study. Lauenstein and Cantillo's (1993) 2 P 18 PCBs provided excellent agreement with total PCBs estimated from a total chlorination method in their study, and 2 P 17 PCBs in the current study agreed well with total PCBs estimated from the sum of 209 congeners measured by high resolution methods (linear regression of P 209 congeners by 2 P 17 PCBs, r 2 = 0.997, p = 0.003, slope coefficient =0.991, intercept not significant, p = 0.892) for a subset of five English sole and herring from this study.
To facilitate broader interpretation of our 2 P 17 PCBs results, we also compared them with P ICES 7 PCBs, a commonly calculated summed PCB estimate in European biota (International Council for the Exploration of the Sea 1987). In a linear regression of 1089 fish and invertebrate samples analyzed by WDFW from 1994 to 2014, the 2 P 17 PCBs were approximately 2.6 times higher than P ICES 7 PCBs calculated from the same samples (linear regression, r 2 = 0.998, p \ 0.0001, slope coefficient = 2.58, intercept = 0.61).
Summed PBDEs were calculated as the sum of detected values for the 11 identified PBDE congeners (IUPAC numbers 28,47,49,66,85,99,100,153,154,155,183), hereafter referred to as P 11 PBDEs, and summed DDTs were calculated as the sum of detected values for five identified DDT compounds (o,p'-DDD, p,p'-DDD, p,p'-DDE, o,p'-DDT, and p,p'-DDT), hereafter referred to as P 5 DDTs. In cases where no single PBDE or DDT analyte was detected in a sample for its group, a value of the greatest LOQ for all analytes in a group was substituted for the total. This value was typically less than 1 ng/g wet wt. At least one PCB congener was detected in every sample.
Although all PBDEs were analyzed by the most current method by combining archived tissues and real-time samples, before 2005 PCBs and DDTs were analyzed using HPLC/PDA detection according to Krahn et al. (1994) (see correction procedure in ''Method Performance'' section). In addition the method of extraction changed from MSEx using PenHex from 1997 to 2003 (Sloan et al. 2005) to ASEx using MeCl (2003. Briefly, for the HPLC/ PDA, sample extracts were reduced in volume to approximately 1 ml and the POPs were separated from interfering compounds (e.g., lipids and aromatic compounds) on a gravity flow cleanup column that contained neutral, basic, and acidic silica gels eluted with PenHex (1997PenHex ( -2003 or MeCl (2003MeCl ( -2005. PCB congeners (IUPAC numbers 77,105,101,110,118,126,128,138,153,156,157,169,170/194,180,and 189) were resolved from five chlorinated pesticides (o,p'-DDD, p,p'-DDD, p,p'-DDE, o,p'-DDT, p,p'-DDT) by HPLC on two Cosmosil PYE analytical columns, connected in series and cooled to 16°C. Compounds were measured by an ultraviolet (UV) PDA and were identified by comparing their UV spectra (200-310 nm) and retention times to those of reference standards in a library. Analyte purity was confirmed by comparing Arch Environ Contam Toxicol (2017) 73:207-229 211 spectra within a peak to the apex spectrum. The concentration of total PCBs analyzed by HPLC/PDA was estimated as the sum of the concentrations of the 16 PCBs listed above (based on individual response factors) plus the sum of the concentrations of other unidentified PCBs calculated by summing areas of peaks identified as PCBs and using an average PCB response factor. Summed DDTs were calculated as the sum of detected values of the five DDT isomers. If no DDT isomer was detected in a sample, the greatest LOQ for the five isomers (typically \0.5 ng/g, wet wt) was substituted as the summed DDT value for that sample.
Lipid concentration was determined gravimetrically in all samples from all years. For PCBs and DDTs from 1997 to 2003, lipids were extracted using MSEx with PenHex and homogenized with a tissue grinder (Sloan et al. 2004). Lipids for PCBs and DDTs from 2003 onward and for all PBDEs (1994PBDEs ( -2015 were extracted using ASEx with MeCl (Sloan et al. 2014).

Method Performance
Recoveries of analytes from the National Institute of Standards and Technology (NIST) Standard Reference Materials (SRMs) were considered the primary measure of analytical performance for all methods, especially over the long-term, as equipment and supplies can vary through time, even for standardized methods. Additional quality control measures, including sample replication and addition of surrogate standards, were used to identify problematic batches according to Sloan et al. (2004Sloan et al. ( , 2014; however, these metrics typically passed quality control limits and so are not further reported here. A sample of NIST SRM 1974a (blue mussel) was analyzed with each HPLC/PDA sample batch; NIST SRM 1947 (Lake Michigan fish tissue) samples were run with GC/MS batches before 2015 and NIST SRM 1946 (Lake Superior fish tissue) was used after 2015. All batches contained 12-14 field samples. A NIST SRM quality control limit was used to flag poorly analyzed batches for reanalysis; batches were rejected if C30% of individual analytes measured in the NIST SRM exceeded 30% of the 95% confidence interval range of the published NIST certified concentrations and reanalyzed until they complied with recovery controls.
SRM recovery target concentrations were further calculated for each of the three POP classes by summing the NIST certified or reference values for analytes that were analyzed in field samples. Measuring recovery of SRM analytes through time provided an estimate of variance for method performance, which we calculated as the 95% confidence interval (CI) for the year factor in a linear regression of the SRM analyte recovery by year (or the CI of the mean analyte % recovery if the year coefficient was not significant). Significant year coefficients from field sample regressions falling within the 95% confidence interval from SRM % recovery regressions were deemed indistinguishable from method variability and censored as not significant as described below.
For the 17 congeners used in 2 P 17 PCB measured by GC/MS, certified values for PCB 66 and 195 were unavailable from NIST, so % recovery was calculated using a sum of the other 15 congeners from the 2 P 17 PCB list. Ninety-six percent of the resulting P 15 PCB concentrations recovered from SRM 1947 and SRM 1946 were within 15% of the NIST certified 95% confidence interval, indicating a high degree of accuracy.
P 15 PCB recoveries in SRMs from the GC/MS showed no significant timerelated bias, based on a linear regression of SRM % recovery by year for English sole (p = 0.49) and herring (p = 0.68). Overall, the 95% confidence interval for the mean P 15 PCB SRM recovery was only ±2.1% for English sole and ±1.9% for herring; accounting for this variability in the % recovery of the SRM, significant time trend year coefficients (i.e., slopes) from field sample regressions of 2 P 17 PCBs by year within these ranges were considered indistinguishable from analytical variability and subsequently censored as not significant. Eighty-six percent of summed PBDE SRM recoveries were within 15% of the NIST-certified 95% confidence interval. PBDE SRM recoveries did not vary significantly with time for English sole and only weakly for herring batches (linear regression of PBDE recovery by time, p = 0.057, adjusted r 2 = 0.084 for English sole; p = 0.012, adjusted r 2 = 0.22 for herring) with slight annual rates of decline of 0.86 and 0.83%. PBDE time trends in field samples falling within the 95% confidence interval of the PBDE SRM % recovery regression slope (-1.7 to -0.024% for English sole and -1.4 to -0.24% for herring) were considered indistinguishable from analytical variability, and subsequently censored as not significant. Eighty-nine percent of summed DDT SRM recoveries were within 15% of the NIST-certified 95% confidence interval; however, SRM recoveries of DDTs decreased consistently through time (linear regression of summed DDT recovery by year, p \ 0.001, adjusted r 2 = 0.52 for English sole; p \ 0.001, adjusted r 2 = 0.67 for herring). DDT time trends reported herein for field collected samples were considered indistinguishable from this bias if they fell within the 95% confidence interval of the slope coefficient from the SRM time regressions (-2.6 to -1.3% for English sole; -2.9 to -1.6%) and subsequently censored as not significant.
SRM recoveries for PCBs analyzed by HPLC/PDA were lower than the GC/MS method; 16% of PCBs, calculated as a sum of 8 congeners in common between the field samples and SRM certified values, were within 15% of the NISTcertified 95% confidence interval. Recoveries of DDTs from SRMs analyzed by the HPLC/PDA were hampered by the presence of unidentified interfering compounds, and so SRM recoveries were unusable for that analyte group. PBDEs were not measured by the HPLC/PDA method. However, DDT and PCB results from HPLC/PDA analyzed field samples were corrected using a robust sampleto-sample comparison described below. This method corrected HPLC/PDA results to GC/MS standards, essentially precluding the need to evaluate or correct for SRM recoveries from the HPLC/PDA results. Moreover, the SRM recoveries for PCBs and DDTs appeared to be consistent through time for each species, based on a visual inspection of the SRM recoveries; low sample size prevented computing valid regression models for the HPLC/ PDA methods.

Method-Bias Corrections for Time Trends
To correct for bias in estimated total PCB estimates associated with analytical method changes, 76 samples that had been analyzed by the older methods were reanalyzed using the GC/MS with ASEx described above: 14 samples by the HPLC/PDA method with MSEx and 62 samples by the HPLC/PDA with ASEx. Correction factors were calculated by regressing 2 P 17 PCBs from the GC/MS with ASEx method against estimated total PCBs calculated from each of the HPLC/PDA groups in the same samples. In both cases, the HPLC/PDA methods underestimated 2 P 17 PCBs; by approximately 33% for the MSEx group (linear regression of 2 P 17 PCBs by total PCBs, r 2 = 0.90, p \ 0.001, slope coefficient = 1.331, intercept not significant, p = 0.97), and 27% for the ASEx group (r 2 = 0.94, p \ 0.001, slope coefficient = 1.267, intercept not significant, p = 0.83).
The same five DDT analytes were measured across both analytical and extraction methods, precluding a need for method-corrections based on summation of differing chemicals in the group. However, a comparison of P 5 DDTs across the methods revealed a difference in quantitation of individual DDTs related to the solvent extraction method. P 5 DDTs measured by the HPLC/PDA method with ASEx were similar to results from the GC/MS with ASEx and therefore required no adjustment (linear regression of paired P 5 DDTs by both methods; r 2 = 0.95, p \ 0.001, slope coefficient = 0.95, intercept not significant, p = 0.77); however, P 5 DDTs measured by HPLC/ PDA with MSEx were approximately 45% lower than the GC/MS method (linear regression; r 2 = 0.95, p \ 0.001, slope coefficient = 1.45, intercept not significant, p = 0.68). Hence, P 5 DDTs were adjusted upwards according to this latter regression model, for samples that were only analyzed by the HPLC/PDA and MSEx.
Beginning in 2005, PBDEs were always measured in both species of fish by the GC/MS method described above with ASEx, and all PBDE in prior years were determined from archived tissue samples. Thus, no adjustments were made for any PBDE data.

Data Analysis
For evaluation of time trends, we used multiple linear regressions of 2 P 17 PCBs, P 11 PBDEs, and P 5 DDTs over time, with year as the independent variable, to evaluate whether observed time trends in POPs were statistically significant. The contribution of biological covariates in explaining POP variability also was evaluated. Compositemean fish length (MCL) was used as a proxy for fish age, which can be considered an estimate of total potential exposure time to contamination. Because the POPs measured in this study are considered lipophilic, total nonvolatile extractibles (percent lipids) were included to help control for variability in that measure. The sex ratio (proportion of male fish) was included to help control for potential differences related to loss of POPs in reproducing female fish. Automatic stepwise regression models were run using SYSTAT (2007), for each log-transformed POP total in separate runs, with the covariates listed above and their interactions. Although some interaction terms were occasionally significant, all interaction terms explained little POP variation (usually less than 3%), and including them did not appreciably change results. Hence, for simplicity they were omitted from the final regression models. Geometric mean POP concentrations and confidence intervals were adjusted for significant covariate terms by predicting the POP concentration through time from the multiple linear regression models, using the grand mean value for MCL, lipids, or sex ratio, when they were included in the final model. Partial coefficients of determination (r 2 ) were calculated at each step of the forward stepwise multiple linear regression to estimate the contribution of each significant covariate in explaining POP concentration in the final model. A significant (p [ 0.05) year coefficient in the final model was accepted as evidence for a significant time trend; however, as an additional guard against overinterpreting weak trends, year coefficients with r 2 \ 0.10 (i.e., time explaining \10% of POP variability) were censored as trivial (and models presented as dashed lines in graphs).
To evaluate the most current status of POPs in English sole and herring, we used analysis of variance (ANOVA) to test the significance of observed log-transformed POP differences between the ten English sole index sites collected in 2015 and between the four herring stocks collected in 2014. Although lipids, fish size, and sex ratio were tested as ANOVA covariates, none was a significant contributor in explaining POP variation, so POPs were not adjusted by these factors for the current status comparison.

Results
Time series data for herring comprised 12 sampling years spanning a 16-year period from 1999 to 2014 for two stocks (Port Orchard and Squaxin), 9 years over the same period for Semiahmoo and 6 years of sampling over a 21-year period for Cherry Point (Table 1). The time series for English sole covered 9 to 12 sampling years over an 18or 19-year span for all sites except the Duwamish River, which was sampled in 1997 and then not again until 2007-2015 (Table 1). The most recent years for both species were analyzed below to represent current status.
Of the three POP classes reported in our study, 2 P 17 PCBs dominated in concentration for both fish species, followed by P 11 PBDE and P 5 DDTs. The 2 P 17 PCB concentrations from the most current sampling were 3-6 times greater than P 11 PBDE and 3-16 times greater than P 5 DDTs in 2014 Pacific herring (Fig. 2) and 10-63 times greater than P 11 PBDE and 10-67 times greater than P 5 DDTs in 2015 English sole (Fig. 3). Percent lipids in whole body herring ranged from 2.5 to 10%, with highest values determined in herring from Squaxin (4.4-10%) and lowest in those from Cherry Point (2.5-4.2%). Size of the male herring ranged relatively narrowly (by study design) from 143 to 181 mm SL (Table 1). Lipid values for English sole muscle were consistently low across all years and all sites, ranging from 0.16 to 0.88% (Table 1). Size of English sole (as mean composite length) ranged from 253 to 380 mm TL, and the proportion of male fish (sex ratio) in the composites ranged widely from 1 to 92%.

Pacific Herring Current Status
PCB concentrations in male Pacific herring in 2014 were highest in fish from the Squaxin stock (210 ng/g, moderately developed basin) and Port Orchard stock (190 ng/g, highly developed basin), followed by 46 ng/g in fish from Cherry Point and 18 ng/g from Semiahmoo stocks, both from a basin characterized by low-development (Table 1; Fig. 2). Levels of P 11 PBDEs and P 5 DDTs showed a similar pattern, with decreasing concentrations from Squaxin to Port Orchard to Cherry Point to Semiahmoo. Fish length and lipid covariates (or their interactions) were sometimes significant in the ANCOVA models comparing locations in the current status analysis; however, they were dropped from the final model, because they consistently accounted for a minor amount (\3%) to the variability in concentration in any of the POPs. Geometric means and confidence intervals (shown in Fig. 2) were not adjusted for any covariates.

Pacific Herring Time Trends
Concentrations of 2 P 17 PCBs remained unchanged over the course of the study in herring stocks from Port Orchard (highly developed basin) and Squaxin (moderately developed basin). The year coefficient for 2 P 17 PCBs in Port Orchard herring was censored, because it was both indistinguishable from SRM variability and its partial correlation coefficient was \0.10; hence a 2 P 17 PCBs rate was reported as not significant for that stock (Table 2). Year was not significant for Squaxin in the multiple linear regression (p [ 0.05). Even though year was not a significant predictor of 2 P 17 PCBs for these two stocks, a predicted (dashed) line was calculated (shown in Fig. 4 for reference). Levels of 2 P 17 PCBs in both herring stocks from the northern, low-development basin (Semiahmoo and Cherry Point stocks) declined significantly, with annual rates of 4.0 and 2.3% ( Fig. 4; Table 2). MCL and lipids were either not significant, or weakly to moderately predictive in all stocks, with partial correlation coefficients ranging from 0.02 to 0.44. 2 P 17 PCBs in all predicted lines calculated for Fig. 4 were adjusted for significant covariates.
The P 11 PBDE levels declined in all herring stocks from all basins, with annual reductions ranging from 4.3 to 8.1%. Fish size (MCL) contributed significantly to explaining PBDE variation with partial r 2 ranging from 0.16 to 0.41), and lipids were weakly predictive for Squaxin (r 2 = 0.02, p \ 0.05) or insignificant (p [ 0.05) for the other three stocks (Table 2). Final PBDE regression models were moderately predictive, with r 2 ranging from 0.52 to 0.71; predicted lines were adjusted using grand mean values for significant coefficients (Fig. 4).
As with P 11 PBDEs, P 5 DDTs declined in all herring stocks, with rates ranging from 5.6% (Port Orchard) in the highly developed basin, 4.2% from Squaxin in the moderately developed basin, to 5.8 and 7.0% in the two stocks from the low-development basin (Cherry Point and Semiahmoo; Fig. 4). Final DDT models were moderately to strongly predictive, with r 2 ranging from 0.54 to 0.87 (Table 2). Along with year, MCL was a significant partial predictor for all four stocks, albeit most strongly for Squaxin herring, and grand mean values for significant covariates were used to adjust the predicted lines in Fig. 4.

English Sole Current Status
Concentrations of 2 P 17 PCBs were generally highest in English sole at highly developed sites; in the most current sampling year (2015) arithmetic mean concentrations Harbor) and sole from the low-development locations had the lowest concentrations; 20 ng/g wet wt from Anderson Island, and less than 10 ng/g wet wt from the remaining three. As with 2 P 17 PCBs, the greatest P 11 PBDE concentrations in 2015 occurred in three highly developed locations, although concentrations were relatively low overall: 4.3, 4.5, and 6.0 ng/g wet wt from Duwamish River, Seattle Waterfront, and Tacoma City Waterway (Table 1; Fig. 3). Concentrations of P 11 PBDEs in the seven locations with moderate to low development were less than 2 ng/g wet wt. Similar to P 11 PBDE, DDTs were low overall, with concentrations ranging from 0.32 to 4.0 ng/g wet wt across all locations. Greatest P 5 DDT concentrations occurred at the four highly developed sites, Duwamish River, Seattle Waterfront, Tacoma City Waterway, and Bremerton (4.0-1.3 ng/g wet wt). Although fish length, percent lipid, and sex ratio covariates (or their interactions) were sometimes significant in the ANCOVA models comparing locations in the current status analysis, they always contributed a minor amount (\3%) in explaining the variability in concentration among the ten sites for any of the POPs. Hence, to simplify the site locations analyses, all geometric means reported in Fig. 3 for English sole POPs were computed without covariate adjustments.

English Sole Time Trends
We did not observe a declining trend in 2 P 17 PCBs from any of the 10 index sites we monitored over the 18-year period from 1997 to 2015 (Fig. 5). Positive, statistically significant (p \ 0.05) year coefficients indicated increasing  Table 2). Year was weakly to moderately predictive for these sites (partial r 2 ranging from 0.16 to 0.41; Table 2), and the overall multiple linear regression models for these sites were at best moderately predictive, with final, adjusted r 2 slightly above 0.50 for only two sites, Tacoma City Waterway, and Everett. English sole from the remaining six index sites exhibited no significant change in PCB concentration over the 19 years. The year coefficient for four of these six sites was not significant (p [ 0.05), and it was censored as not significant because its partial r 2 was \0.10 at two sites (Anderson Island and Hood Canal). Sex ratio, lipids, or MCL contributed significantly to the final model for eight sites overall and predicted lines were adjusted using grand mean values of the coefficients for significant covariates (Fig. 5). Predicted lines for sites where there was no significant trend were calculated and are shown in Fig. 5 as dashed, for reference.
Although year was a weak predictor for P 11 PBDEs overall in English sole (r 2 ranged from 0.11 to 0.17), the concentrations of P 11 PBDEs with r 2 [ 0.10 declined significantly at two highly developed sites (Bremerton, 4.0% and Seattle Waterfront, 3.7%) and at one low-development site (Vendovi Island, 7.2%; Fig. 5; Table 2). The year coefficient was not significant for five sites (p [ 0.05) and it was censored as not significant, because its partial correlation coefficient was \0.10 at two sites (Everett and Hood Canal). Fish size, lipids, and sex ratio were weakly to moderately predictive (r 2 from 0.06 to 0.56) for nine sites, and predicted lines were adjusted using grand mean coefficient values for significant covariates (Fig. 5). Predicted lines for sites where the time trend was not significant were calculated and are shown in Fig. 5 as dashed, for reference. Fig. 2 Geometric mean concentration of 2 P 17 PCBs, P 11 PBDE, and P 5 DDTs (ng/g wet wt) in four stocks of Pacific herring from 2014, the most recent sampling period. The levels of basin development are illustrated by color; white for low, cyan for moderate and dark blue for high. For each POP class, sites with the same lower case letter were not significantly different from each other Fig. 3 Geometric mean concentration of 2 P 17 PCBs, P 11 PBDE, and P 5 DDTs (ng/g wet wt) in English sole from ten index sites in 2015, the most recent sampling effort. Levels of site development are illustrated by color; white for low, cyan for moderate and dark blue for high. For each POP class, sites with the same lower case letter were not significantly different from each other Table 2 Adjusted coefficients of determination (adj r 2 ) for forward stepwise multiple linear regression models and partial r 2 for each significant factor, for models testing the predictive value of year, sex ratio (% of male fish in a composite for English sole only), tissue lipid concentration (%), and fish size (mean composite length in mm) on log-transformed PCBs, PBDEs, and DDTs.
Mean  Concentrations of P 5 DDTs declined in English sole from two highly developed sites (Duwamish River, 4.2% and Bremerton, 3.6%) and from one low-development site (Strait of Georgia, 4.2%; Table 2; Fig. 5). One highly developed site (Tacoma City Waterway) exhibited an increase in P 5 DDTs, with an annual rate of 5.2%. Year was not significant (p [ 0.05) or censored as not significant (r 2 \ 0.10) for the remaining six sites ( Table 2). As with 2 P 17 PCBs and P 11 PBDE, full models for P 5 DDTs were weakly to moderately predictive, with r 2 ranging from 0.11 to 0.46. Predicted lines shown in Fig. 5 were calculated using grand mean values for significant covariates. Predicted lines for sites where the P 5 DDTs rate was not significant are shown as dashed, for reference.

Consistency in Monitoring Data
Monitoring programs are challenged by continually shifting and changing analytical methodologies. Two practices can mitigate loss of consistency when methods change: (1) conducting paired analyses by both old and new methods to generate correction factors, and (2) retaining archived tissues for reanalysis once new methods are developed, validated, and then adopted. Each of these practices requires additional cost for duplicate or reanalysis or raises concerns related to sample quality over time. In this study we employed both practices: PBDEs were always conducted using one methodology, with older archived samples analyzed along with current samples to build a time series.
Other long-term monitoring studies have employed this practice, including Ross et al. (2013) who used it to construct the Puget Sound harbor seal time trend line reported earlier herein. Correction factors generated from comparing outdated methods with newer methods appeared to have provided a good tool for correcting biases in our longterm monitoring program. Linear regressions of paired data (old vs. new methods) typically provided excellent fit, with strong coefficients of determination (usually [0.90), and sufficient homoscedasticity to accept significant correlation coefficients as correction factors.

POPs in Herring
The long-term trend analyses indicate PCB tissue residues declined in both stocks of Pacific herring from the lowdevelopment basin, however they remained unchanged in the highly and moderately developed basins. This pattern is consistent with PCB trends in harbor seals, a primary predator of herring, sampled from the south Puget Sound, the moderately developed basin. Although harbor seals exhibited a steep decline in blubber PCBs from the 1980s to 1990s (prior to this study), PCBs remained unchanged or declined only slightly in the most recent years (2003Ross et al. 2013). A number of studies have demonstrated exponential declines in PCBs in pelagic species in marine and large freshwater ecosystems in the 1980s, with concentrations leveling off over the past decade or two, including herring (Bignert et al. 1998;Miller et al. 2013) and sprat (Szlinder-Richert et al. 2009) from the Baltic Sea, and marine fishes, birds, and mammals from the Canadian Arctic (Braune et al. 2005). Coho and Chinook salmon in Fig. 4 Temporal trend in 2 P 17 PCBs, P 11 PBDEs, and P 5 DDTs in Pacific herring from four stocks in Puget Sound, WA, sampled from a highly developed basin (Port Orchard stock), moderately developed basin (Squaxin stock) and low-development basin (Semiahmoo and Cherry Point stocks). Dashed lines indicate regression models where year was not significant (p [ 0.05) or censored as trivial (partial r 2 \ 0.10) Arch Environ Contam Toxicol (2017) 73:207-229 221 the US Great Lakes followed a rapid first order kinetics PCB decline in the 1980s but subsequently leveled off to a non-zero steady state (Stow et al. 1994;Lamon et al. 1999). PCB concentrations in central and southern Puget Sound herring were considerably higher than levels reported in herring from the most highly developed areas of the southern Baltic Sea sampled during the same time period. Szlinder-Richert et al. (2009) reported mean ICES 7 PCBs in muscle tissue of Baltic herring ranging from 15 to 27 ng/ g wet wt and 360-820 ng/g lipid wt from 1997 to 2006. 17 PCBs, P 11 PBDEs, and P 5 DDTs in English sole from ten index sites in Puget Sound.
Dashed lines indicate regression models where year was not significant (p [ 0.05) or censored as trivial (partial r 2 \ 0.10) Because we sampled whole male fish of a relatively fixed size and reported an estimate of total PCBs (2 P 17 PCBs), a direct comparison of PCB wet wt results was inappropriate. Instead, we limited our comparison to concentrations of ICES 7 PCBs on lipid weight to mitigate the difference in tissue types and PCB summation methods. Whole Pacific herring from central and southern Puget Sound were roughly two times higher than those reported in southern Baltic herring muscle samples, with ICES 7 PCBs calculated in the Puget Sound herring ranging from 690 to 1800 ng/g lipid weight from a similar time period (1999-2006; Table 1). Southern Baltic herring PCB levels were more similar to concentrations we observed in the two stocks from the Strait of Georgia (least-developed basin), with concentrations ranging from 250 to 900 ng/g lipid. These comparisons were somewhat confounded by differences in fish size and sex; southern Baltic herring from Szlinder-Richert et al. (2009) were larger, with maximum sizes reported as 305 mm compared with mean composite lengths of Pacific herring ranging from 143 to 181 mm. Based on size, Pacific herring in the current study were likely younger than those used in the Baltic study, and so may have experienced a shorter duration of exposure to PCBs. Moreover, the selection of only male herring in the Puget Sound study likely would have increased PCB differences if female fish (represented in the Baltic samples) were depurating POPs during reproduction.
PCB concentrations in Puget Sound herring were within the range reported for two forage fish species, silversides (Menidia audens) and topsmelt (Atherinops affinis), in San Francisco Bay from 2010 (Greenfield and Allen 2013). These authors reported mean concentrations for sum of 209 PCBs (equivalent to our 2 P 17 PCBs) of 138 ng/g wet wt from ''ambient,'' randomly selected sites, and 441 ng/g wet wt from sites known to be highly PCB-contaminated. These species were selected by Greenfield and Allen (2013) for similar reasons herring were selected in the present study; herring, topsmelt, and silversides represent small, schooling, midwater prey upon which many higher predators rely and therefore can illustrate an important pathway for POPs like PCBs. Some significant differences between the species prevent an unqualified comparison: the San Francisco silversides and topsmelt were much smaller, roughly half the size of Puget Sound herring, and they exhibited somewhat lower lipids than two of the Puget Sound stocks. Additionally, they likely exhibit somewhat different feeding ecology; topsmelt and silversides were thought to represent local sediment conditions, by feeding on organisms near to the seafloor in shallow waters, whereas herring in Puget Sound likely forage more widely, in much deeper waters, on more pelagic prey.
These Puget Sound and other studies underscore the environmentally recalcitrant characteristics of PCBs, the importance of pelagic biota as a sink for these chemicals, and the role pelagic fishes play as vectors of PCBs to apex predators. Rice et al. (2012) reported Pacific herring as the second most frequently encountered forage fish in a study of Puget Sound surface waters, after juvenile Chinook salmon. These species were encountered in well over 50% of all surface trawls conducted in the same basins represented in the current study. Pacific herring are thought to exert bottom-up ecosystem control in Puget Sound, based on Ecopath/Ecosim modeling (Harvey et al. 2012), and indeed they are one of several midwater fishes that dominate the diet of harbor seals, Puget Sound's most abundant resident fish-eating marine mammal. Pacific herring, tomcod, hake and salmon were estimated to constitute over 60% of harbor seal diet (Cullon et al. 2005).
O'Neill and West (2009) posited residence time in Puget Sound (i.e., feeding on Puget Sound-resident prey) as a risk factor for PCB accumulation in migratory species that move through this system. Elevated levels of PCBs have been documented in species throughout Puget Sound's pelagic food web, from plankton (West et al. 2011a) to Puget Sound-resident pelagic fish predators (West et al. 2011b), migratory adult salmon (Cullon et al. 2009;O'Neill and West 2009), marine birds (Good et al. 2014), and marine mammals including resident harbor seals (Ross et al. 2004) and migratory fish-eating killer whales (Ross et al. 2000, Krahn et al. 2007. The monitoring program forming the basis of this paper compares POP tissue residues in fish to published critical body residues (CBR; Meador et al. 2011) to predict health impairment. Using this CBR concept, PCB and PBDE concentrations in some herring stocks are cause for concern. Current (2014) PCB concentrations in herring from Port Orchard (3400 ng/g lipid) and Squaxin (4800 ng/g lipid) stocks in this study (Table 1) exceeded a 2400 ng/g health effects threshold determined for juvenile salmon (Meador et al. 2002; the most appropriate CBR for herring from current literature). Because herring sampled in this study were male fish targeted at age 3 years, it is likely that older male fish could have substantially higher POP concentrations (and, conversely, lower levels in younger fish and females of spawning age). Juvenile salmon exposed to PCBs in the range of exposures we observed in Squaxin and Port Orchard herring exhibited a number of health effects including increased enzyme activity, altered thyroid hormones, and increased mortality. Trophic transfer of PCBs via herring also presents a health risk to their predators; 22% of adult Chinook salmon sampled from Puget Sound exceeded the Meador et al. (2002) CBR for PCBs (O'Neill and West 2009). Furthermore, adult Chinook salmon (Cullon et al. 2009;Mongillo et al. 2016) and other pelagic fishes (Cullon et al. 2005;West et al. 2011b) from the Salish Sea (including Puget Sound) were identified as key vectors of PCBs to killer whales and harbor seals, both important apex predators in the ecosystem. Using this exposure pathway, Hickie et al. (2007) predicted PCBs in Southern Resident killer whales feeding in Puget Sound would not fall below a 17-mg total PCBs/ kg blubber lipids health-effects threshold concentration until the year 2063. Additionally, the Washington Department of Health (WDOH) has recommended restricting human consumption of Chinook salmon residing in Puget Sound because of PCB contamination in that species (Washington Department of Health 2015).
Although PCBs are still of significant concern in Puget Sound's pelagic food web almost 40 years after their production and legal use was banned, decreasing trends of PBDEs and DDTs in all herring stocks point to the success of remediation and control efforts for POPs in Puget Sound. Washington State recommended reductions in PBDE use as early as 2006 (Washington Department of Ecology 2006), and its state legislature restricted further use of PBDEs in 2008. These actions were taken to reduce exposure of humans to PBDEs and to mitigate input of PBDEs to Puget Sound. PBDE reductions notwithstanding, current (2014) tissue levels may still be high enough to impair herring health. Methodological differences in measuring and reporting PBDE concentrations hinder direct interpretation of published PBDE CBRs with the concentrations we reported in Table 2. However, Tomy et al. (2004) and Arkoosh et al. (2017) reported variations in thyroid homeostasis related to dietary doses of PBDEs that produced tissue burdens similar to those we observed in wild herring. In particular, Arkoosh et al. (2017) observed depressed tri-iodothyronine (T3) in plasma of juvenile Chinook salmon with a PBDE body burden of 37 ng/g wet wt (for a mixture of two dominant congeners typically observed in fish samples, PBDE 47 and 99). We reported a similar range of P 11 PBDE concentrations; 30 and 44 ng/g wet wt from Squaxin and Port Orchard herring stocks in 2014, and PBDE 47 and 99 were also the most abundant congeners in our samples, accounting for approximately 67% of the P 11 PBDE concentration. Tomy et al. (2004) reported depressed thyroxine levels in plasma of lake trout fed ''low'' and ''high'' doses of a mixture of 13 PBDEs resulting in maximum lipid-adjusted tissue concentrations of 150 to 1600 ng PBDEs/g lipids (manual summation of concentrations in their Fig. 1) compared with 410-1000 ng/g lipids for the four herring stocks that we reported from Puget Sound. Arkoosh et al. (2010;Fig. 2a) reported increased disease susceptibility in juvenile Chinook salmon dosed with five commonly observed PBDE congeners (47, 99, 100, 153, and 154) resulting in total PBDE concentrations similar to what we observed in our field sampled herring.
Assuming that herring respond similarly to Chinook salmon, we would predict some impairment of herring health related to PBDE exposure.
Unlike PCBs and PBDEs, DDT concentrations in all herring were well below a CBR for whole fish of 600 ng/g wet wt (Beckvar et al. 2005); the greatest concentration of P 5 DDTs observed in the current study was 39 ng/g wet wt, suggesting DDTs alone posed less risk to herring than the other POPs we measured. However, a tissue residue approach designed to evaluate the toxicity of chemical mixtures (Dyer et al. 2011) for the 3 POP categories we measured may be appropriate for future analyses.

POPs in English Sole
Whereas Pacific herring (this report) and other pelagic species (West et al. 2011a, b) exhibited pervasive POP contamination of the pelagic food web in highly or moderately developed oceanographic basins (on the order of hundreds of square km) in Puget Sound, English sole illustrated POP conditions at a smaller spatial scale (probably a few square km), related to their proximity to contaminated sediments and limited home range (O'Neill et al. 2007). Results reported for English sole herein are congruent with sediment conditions reported by a companion Puget Sound sediment monitoring program, which measured contaminants in sediments collected near our English sole sampling sites. Long et al. (2005) characterized only approximately 1% of Puget Sound sediments as being significantly impaired by toxic chemicals, including PCBs and DDTs, all of which occurred in close proximity to the same four urban centers where we observed the most contaminated English sole (the bays and shorelines near Seattle, Tacoma, Bremerton, and Everett). English sole exhibited feeding-site fidelity and a small foraging area in an acoustic tagging study (O'Neill et al. 2007), further validating the use of this species as an indicator of more local, bay-scale conditions.
Relatively high PCB concentrations and their increase (or lack of decline) in English sole indicates efforts to control these chemicals in Puget Sound waters have not been effective enough to achieve declining trends. Moreover, mean ICES 7 PCB levels in English sole from four urban sites (55 ng/g from Seattle Waterfront, Tacoma City Waterway, Duwamish River and Bremerton sites combined) were approximately 12 times higher than the 4.6 ng/ g wet wt reported in confamilial European flounder (Platichthys flesus) from densely populated and highly urbanized areas in the Baltic Sea in the most recent year when both programs measured these species [ These observations are consistent with Brown et al. (1998), who documented increasing PCBs in English sole from 1984 to 1990 at two highly developed sites near the Seattle Waterfront and the Tacoma City Waterway sampled in the current study. They also reported no change in PCBs near our low-development Anderson Island site, where we also reported no declining PCB trend. Brown et al. (1998) concluded that redistribution of sediment particles from depositional areas in deeper parts of the contaminated urban bays or continued inputs from existing onshore sources likely accounted for the increase trends they observed in PCB concentration at these highly developed sites. Although it is difficult to track all in-water activities in the urban waterways where English sole were collected for this study, many of them experienced significant anthropogenic sediment-disturbing activities throughout our 19-year monitoring period, including sediment removal via dredging for navigation or for contaminant remediation. Other factors that may have prevented PCB declines in English sole include: (1) recirculation of persistent congeners in the benthic food web; (2) new PCB inputs that have not been adequately controlled; and (3) increasing wastewater and stormwater conveyance of PCBs related to increases in human population in the region that may offset source control efforts.
The P 11 PBDEs and P 5 DDTs declined or remained unchanged in all English sole with one exception, the 5.2% annual increase of P 5 DDTs in English sole from the Tacoma City Waterway. The levels of P 11 PBDEs and P 5 DDTs at many English sole sites were heading in a desirable direction (declining) and all concentrations were low (especially compared with 2 P 17 PCBs), with values that varied narrowly, near the limit of quantitation for many of the analytes used in the summations. In these cases, trends should be interpreted with caution because of the possibility of true trends being obscured by method performance near the limits of quantitation. In addition, recovery of DDT analytes in standard reference materials was particularly problematic because of interference with unknown chemicals, and inconsistent performance through time that was difficult to quantify. These cautions notwithstanding, P 5 DDTs in English sole from the Tacoma City Waterway is particularly notable, because it showed a relatively strong increasing trend. The Tacoma City Waterway index site is located near the mouth of the Puyallup River, which receives runoff from agricultural, residential, and urbanized lands. Johnson et al. (1988) noted recalcitrant DDTs from agricultural soils can be a significant source of DDTs in Washington rivers as a result of erosion, which may explain the source of DDTs at the Tacoma site. This site also has experienced a number of dredging and cleanup operations, as well as removal of contaminated shoreline soils, which may have exposed otherwise sequestered DDTs.
Similar to San Francisco Bay marine fish (Davis et al. 2007), a primary concern regarding PCBs in English sole is protection of human health from contaminated seafood. PCB concentrations in English sole were high enough to generate advice restricting consumption from five of the ten Puget Sound index sites in the current report; Tacoma City Waterway, Duwamish River, Eagle Harbor, Seattle Waterfront, and Bremerton (Washington Department of Health 2015). Furthermore, PCBs in English sole from these five sites, plus those from Anderson Island and Everett exceeded a human health screening value (10 ng/g wet wt) for marine fish in San Francisco Bay (Davis et al. 2007). PBDEs (0.28-6.0 ng/g wet wt) and DDTs (0.30-3.1 ng/g wet wt) in English sole were low relative to fish-health CBRs; that is, near to PBDE concentrations reported for experimental controls in Arkoosh et al. (2010) and more than 100 times lower than the Beckvar et al. (2005) DDT threshold of 600 ng/g wet wt.

Challenges to Reducing POPs
Evaluation of best management practices to mitigate POPs in the marine environment requires a full understanding of not only the sources of POPs to the system but also the fate and transport of POPs in the system (Arnot and Gobas 2004;Osterberg and Pelletier 2015). The POPs selected for tracking in this long-term monitoring program cover a wide range of hydrophobicity and half-lives, with most considered environmentally recalcitrant (Johnson et al. 2014). For example, Fisk et al. (1998) reported biomagnification factors (BMF) for the 16 PCB congeners we quantitated ranging from 1 to 16, and half-lives of 44-224 days, with slow or no metabolism of this group in fish (for a low dose scenario). Similarly, Tomy et al. (2004) reported BMFs for PBDEs in lake trout ranging from 2.1 to 18, and half-lives of 38-210 days, although accurate assessment of PBDE BMFs is hampered by biotransformation of highly brominated (e.g., deca-molecules) to lower brominated forms within fish. Although POPs can degrade in the environment by photolysis (e.g., PBDEs; Fang et al. 2008) the more likely degradation mechanisms for POPs in biota are biological, because light may not penetrate solid matrices, such as sediments and tissues (Johnson et al. 2014).
Metabolic breakdown of PCBs in marine sediments can occur by bacteria (Papale et al. 2016;Matturro et al. 2016); however, the extent of this degradation capacity is unknown.
POPs in our benthic indicator (English sole) clearly reflected sediment conditions where they foraged; however, POPs may enter pelagic food webs via more complex mechanisms. Others have shown POP uptake into pelagic food webs directly from surface waters via ad/absorption by bacteria, other microplankton, and phytoplankton (Larsson et al. 2000;Hudson et al. 2005). Sinking POPladen particles that may otherwise reach the seafloor in shallow systems can be intercepted by vertically migrating zooplankton, such as krill (Euphausia spp; Dilling et al. 1998) in deeper systems. These species are known to concentrate and feed on sinking particulates at the density gradient, where the sinking rate of particles is slowed (Dilling and Alldredge 2000) and which is a characteristic of deeper, vertically stratified systems, such as Puget Sound. The disparity in POP status and trends we observed between herring and English sole probably reflects these processes, and a deeper investigation of these processes may better inform the development of best management practices to mitigate POPs in systems like Puget Sound. Reduction of POPs in English sole may achieved most directly by reducing POPs in sediments from highly developed embayments. POP reductions in the pelagic food web may be more complicated and more responsive to reducing POPs in the water column, whether originating from stormwater or nonpoint sources, point-source inputs, sediment resuspension, atmospheric transport, or trophic recirculation within the pelagic food web.

Conclusions
Relatively high and static levels of PCBs in Pacific herring from the moderately and highly developed Puget Sound basins (this study), as well as from resident pelagic fishes in other large coastal estuaries, such as San Francisco Bay (Davis et al. 2007) and the Baltic Sea (Szlinder-Richert et al. 2009) underscore the persistence of PCBs in the pelagic food webs of nearshore, relatively enclosed marine ecosystems. Remediation efforts to date seem to have reduced PCBs from high levels in the 1970s and 1980s; however, levels of significant concern continue. PCBs in Puget Sound's benthic indicator species appear to be increasing in several urbanized areas, and not decreasing in any sampled location. Particular challenges identified by these studies may be to develop methods for reducing new PCB sources to the pelagic food web and to understand the causes for PCB increases in benthic species from urbanized locations. Puget Sound's pelagic prey base continues to be a hot-spot of PCBs in the northeastern Pacific Ocean region, resulting in longterm contamination of apex predators and other species important to healthy ecosystem function and to human use. Recovery of at least two species listed for protection and recovery by the U.S. Endangered Species Act and Canadian Species at Risk Act, southern resident killer whales (Krahn et al. 2002) and Chinook salmon (Myers et al. 1998) may be significantly hindered by their exposure to PCBs from their prey base in Puget Sound (Hickie et al. 2007;O'Neill and West 2009;Mongillo et al. 2016). PBDEs and DDTs appear to be on the decline in both benthic and pelagic species, except for a few notable exceptions, suggesting that efforts to reduce or remove these chemicals have been largely successful. Continued tracking of POPs in benthic species, such as English sole, will continue to provide a measure of efficacy of more local cleanup or remediation efforts.