Contamination and Human Health Risk Assessment

0 downloads 0 Views 1MB Size Report
aromatic hydrocarbons in the aquatic environment. CRC Press,. Boca Raton, pp 151–184. Camilli R, Reddy CM, Yoerger DR, Van Mooy BAS, Jakuba MV,.
Arch Environ Contam Toxicol DOI 10.1007/s00244-017-0394-9

SPECIAL ISSUE: OCEAN SPILLS AND ACCIDENTS

Contamination and Human Health Risk Assessment of Polycyclic Aromatic Hydrocarbons (PAHs) in Oysters After the Wu Yi San Oil Spill in Korea Andrew Loh1,2 • Un Hyuk Yim1,2 Moonkoo Kim1,2



Sung Yong Ha1 • Joon Geon An1



Received: 30 November 2016 / Accepted: 14 March 2017 Ó Springer Science+Business Media New York 2017

Abstract After the collision of the Singapore-registered oil tanker M/V Wu Yi San into the oil terminal of Yeosu, Korea on January 31, 2014, approximately 900 m3 of oil and oil mixture were released from the ruptured pipelines. The oil affected more than 10 km of coastline along Gwangyang Bay. Emergency oil spill responses recovered bulk oil at sea and cleaned up the stranded oil on shore. As part of an emergency environmental impact assessment, region-wide monitoring of oil contamination in oyster had been conducted for 2 months. Highly elevated concentrations of polycyclic aromatic hydrocarbons (PAHs) were detected at most of the spill affected sites. Four days after the spill, the levels of PAHs in oysters increased dramatically to 627–81,000 ng/g, the average of which was 20 times higher than those found before the spill (321–4040 ng/g). The level of PAHs in these oysters increased until 10 days after the spill and then decreased. Due to the strong tidal current and easterly winter winds, the eastern part of the Bay—the Namhae region—was heavily contaminated compared with other regions. The accumulation and depuration of spilled oil in oyster corresponded with the duration and intensity of the cleanup activities, which is the first field observation in oil spill cases. Human health risk assessments showed that benzo[a]pyrene equivalent concentrations exceeded levels of

& Un Hyuk Yim [email protected] 1

Oil and POPs Research Group, South Sea Research Institute, KIOST, Geoje 53201, Republic of Korea

2

Marine Environmental Science Major, Korea University of Science and Technology, Daejeon 34113, Republic of Korea

concern in the highly contaminated sites, even 60 days after the spill.

On January 31, 2014, the Singapore-registered oil tanker Wu Yi San collided into the oil terminal of a petrochemical complex in Yeosu, Korea. The collision with the vessel impacted a bridge connecting a pier to the petrochemical facility, rupturing three floating oil pipelines. From the ruptured pipelines, approximately 900 m3 of Basrah Light crude oil, naphtha, and an oil–water mixture were spilled into the sea. The spilled oil mixture spread rapidly over 10 km of coastline along Gwangyang Bay. Directly after the spill, emergency oil spill responses removed the bulk of the oil at sea over 4 days, and stranded oil on the shoreline was cleaned extensively for a month after the spill. Cleanup using high-pressure sprays, flushing, absorption pads, snares, skimmers, and booms were applied to the affected shoreline areas. Stranded oil can release toxic oil components in the intertidal zone and can affect the surrounding environment, including oyster and mussel beds (Carls et al. 2001). Toxicity of petroleum related products often is caused by the presence of organic compounds, such as saturate, aromatics, and carboxylic acids (Neff 2004). Among them, polycyclic aromatic hydrocarbons (PAHs) have highly varying toxic potencies, and groups of parent and alkylated PAHs in multi-media samples are routinely monitored after oil spills. Main sources of PAHs can be divided as petrogenic or pyrogenic. Petrogenic PAHs originate from petroleum and their refined products, and pyrogenic PAHs are produced from oxygen-depleted, high-temperature combustion of fossil fuels and biomass. Petrogenic or pyrogenic sources of PAHs can be distinguished through interpretation of PAHs profiles. Generally, petrogenic PAH

123

Arch Environ Contam Toxicol

exhibits a bell-shaped pattern of alkylated PAHs and depletion of HMW PAHs, and pyrogenic PAH exhibits a skewed pattern with predominance of unsubstituted compounds over their alkylated homologues (Douglas et al. 2007; Lake et al. 1979; Youngblood and Blumer 1975). Bivalves are known for their ability to accumulate significant amounts of contaminants through direct or indirect exposure from the surrounding environment (Farrington 1991). Their low metabolic capacity and high body burden allow them to accumulate wide ranges and high levels of contaminants (Cossa 1989). Since the first proposal of Goldberg’s ‘‘Mussel Watch’’ program in 1975, regional or international monitoring programs using mussel or oysters have been continuously conducted (Farrington et al. 2016). In the same context, bivalves have been extensively used as biomonitors for petroleum hydrocarbons and to measure the rate of recovery after oil spill accidents (Elmgren et al. 1983; MLTM 2013; Short and Heintz 1997; Tunnell et al. 1981). Right after the Wu Yi San oil spill, multi-media environmental monitoring, including oil, seawater, sediment, pore water, and bivalves, were conducted to evaluate acute environmental damages by the spilled oil. This study mainly focused on the oyster monitoring to identify (1) the spatial distribution and temporal changes of petroleum derived PAHs, and (2) the potential human health risks from consuming the contaminated oysters.

Experimental Methods Study Area and Sample Collection Gwangyang Bay is a semi-enclosed embayment on the southern coast of the Korean Peninsula (Fig. 1). The western part of the bay is relatively hydro-dynamically static with poor water exchange and sluggish circulation; the water depth is less than 5 m and the tidal exchange volume averages 4%. The Eastern part of the bay is the main channel for tidal exchange; water depth is more than 20 m and tidal exchange volume is on the average of 14%. Tidal exchange volume ranges from 20% at spring tide to 8% at neap tide, with an average of 14% (POSCO 1983). Gwangyang Bay is a highly industrialized area with steel and petrochemical complexes as the dominant industry; the area is thus prone to intentional or unintentional contamination of petroleum hydrocarbons into the surrounding environment. This area was designated as one of the Special Management Area (SMA), and as part of the SMA program, some of the priority pollutants have been routinely monitored in bivalves and sediments (Hong et al. 2011; KORDI 2003).

123

One of the main types of spilled oil (Basrah Light crude oil) was collected from the tanks leading to the ruptured pipelines. Stranded oil was collected from the shorelines a day after the spill. Oyster samples (Crassostrea gigas) were collected from a total of 21 stations along the spill regions in Gwangyang Bay, which includes 8 stations in Yeosu, 12 stations in Namhae, and 1 station in Hadong (Fig. 1). Sampling was performed at intervals of 4, 10, 30, and 60 days after the spill for environmental impact assessments. Pre-spill monitoring had been conducted at seven stations (St.1, St.2, St.4, St.7, St.8, St.18, and St.21) in April 2012. More than 20 oysters were randomly collected from each station and were shucked from their shells then stored at -20 °C before analysis. Analytical Methods Our analysis of the petroleum hydrocarbons in the oil was based on the method described by Wang et al. (1994), with slight modifications (Yim et al. 2011). For the oil samples, approximately 20 mg of oil was dissolved in dichloromethane (GC2 grade, Burdick and Jackson) and fractionated using activated silica gel columns (100–200 mesh, Sigma-Aldrich). Saturate fractions (F1) were collected using hexane (GC2 grade, Burdick and Jackson) and aromatic fractions (F2) were eluted with a mixture of dichloromethane and hexane (1:1 volume). Equal volumes of fractions F1 and F2 were mixed to form a fraction denoted as F3, which were then used for the analysis of the total amount of petroleum hydrocarbons. Surrogate and GC-internal standards were added into each fraction for quantitation and quality controls. For the oyster samples, extraction and clean up procedures were based on the method described by Jin et al. (2008). Approximately 20 g of homogenized oyster tissue sample material was mixed with anhydrous sodium sulfate and then extracted for 16 h with a Soxhlet apparatus through 200-mL aliquots of GC-grade dichloromethane. Surrogate standards were added before the extraction. The extracts were then sequentially cleaned with silica gelalumina columns followed by high-pressure liquid chromatography equipped with a gel permeation column to remove interferences such as lipids. Cleanup columns were prepared by slurry packing of 10 g of alumina (deactivated with 1% water), followed by 20 g silica gel (deactivated with 5% water), and topped with 1 cm of anhydrous sodium sulfate. The eluates were concentrated to 0.5 mL and then spiked with internal standards and finally adjusted to an accurate preinjection volume of 1 mL for gas chromatography–mass spectrometry analysis. The saturated compounds and the total petroleum hydrocarbon (TPH) within the oil samples were analyzed using an Agilent 7890 GC equipped with a flame-

Arch Environ Contam Toxicol Fig. 1 Location of sampling site in Gwangyang Bay. The oil spill site is represented by red star. Post-spill monitoring sites are represented by black circle. Pre-spill ? post-spill monitoring sites are represented by blue square

Hadong

Gwangyang

St.21 St.20

St.1

St.19 St.18 St.17 St.16

Spill Site St.3

St.2

St.15 St.14 St.4

Namhae

St.13

Yeosu

St.12 St.11

St.5 St.6

St.10 St.9 St.7

South Korea

St.8

ionization detector (FID) and a DB-5 capillary column (30 m 9 0.32 mm 9 0.25 lm; Agilent J&W, Santa Clara, CA). The oven temperature was programmed to start at an initial temperature of 50 °C (2 min) and then ramped to 300 °C at 10 °C/min and held for 16 min. The injector and detector temperatures were set at 300 °C. GC-detectable TPH, resolved peaks, and the unresolved complex mixture (UCM) were analyzed. The GC-resolved peaks ranged from n-C8 to n-C40 including the selected isoprenoids, pristine, and phytane. The PAHs were analyzed using an HP 5890 GC, equipped with an HP 5972 mass selective detector (MS); a DB-5MS (30 m 9 0.25 mm 9 0.25 lm; Agilent J&W) was used for analysis of the PAHs in both the oil and bivalve samples. The mass spectrometer was operated in a selected ion-monitoring (SIM) mode based on the molecular ions of the target PAHs. The oven temperature was

programmed to start at an initial temperature of 60 °C (1.5 min), then ramped up to 300 °C at 4 °C/min, and held for 10 min. The injector temperature was set at 300 °C and the transfer line was set at 280 °C. We analyzed 16 United States Environmental Protection Agency (US EPA) priority PAHs (16 PAHs) and alkylated PAHs. The alkylated PAHs included C1- to C4- naphthalene, C1- to C3-fluorene, C1to C4-phenanthrene, C1- to C3-dibenzothiophene, and C1to C3-chrysene. Total PAHs refer to the sum of the 16 PAHs, the alkylated PAHs, and an additional four other PAHs (dibenzothiophene, retene, benzo[e]pyrene, and perylene). Abbreviations of the target PAHs are as follows: naphthalene (Nap), acenaphthylene (Acl), acenaphthene (Act), fluorene (Flu), phenanthrene (Phe), anthracene (Ant), dibenzothiophene (DBT), fluoranthene (Flr), pyrene (Pyr), benzo[a]anthracene (BaA), chrysene (Chr), benzo[b]fluoranthene (BbF), benzo[k]fluoranthene (BkF),

123

123

1.84 ± 0.04 1.46 ± 0.03 2.42 ± 0.03 1.88 ± 0.045 7.0 ± 2.0 17.0 ± 4.7 Values were based on the average of six analyzed stranded oil samples

107.5 ± 22.1 140.1 ± 31.5 19.4 ± 4.2 33.5 ± 3.3 22.2 ± 1.4 Strandeda

24.9 ± 5.8

21.8 Basrah Light

a

1.46 2.51 1.90 13.3 43 141 356 10.2

Asphatene (%) Resin (%) Aromatic (%) Saturate (%) Oil type

Table 1 Chemical properties of source and stranded oil

Of the three types of oil spilled, the highest amount released was Basrah Light crude oil (BLC), followed by naphtha and then the oil mixture, with a ratio of 4:2:1, respectively. BLC is categorized as a Group III oil according to the persistence classification guidelines, which have specific gravities between 0.85 and 0.95 and have tendency to lose up to 40% of volume through evaporation and, thus, have high potential for long-term environmental impact (ITOPF 2014). The percentage of aromatics (58.5%) in BLC was much higher than that of the saturate (21.8%), while percentages of resins and asphaltenes were relatively low (9.5 and 10.2%, respectively; Table 1). Chromatograms showed that n-alkanes ranged from C8 to C40 with C8 being most abundant. The PAHs and their alkyl homologues ranged from naphthalene to benzo(ghi)perylene, where alkyl homologues of naphthalene (C1- to C4-) and dibenzothiophene (C1- to C3-) were the most dominant (Fig. 2a). Six stranded oil samples collected one day after the spill showed compositions indicative of a slightly weathered state. Compared with spill source oil, the percentage of aromatics (24.9%) decreased, and those of resins (33.5%) and asphaltenes (19.4%) increased correspondingly, reflecting the initial weathering effects (Table 1). GC-FID chromatograms showed n-alkanes ranging from C11 to C40

TPH lg/mg

Chemical Properties of Source and Stranded Oil

9.5

UCM lg/mg

Results and Discussion

58.5

R Alkane lg/ mg

R PAHs lg/ mg

2-/3-mD/1mD

4-mD/1mD

C31S/R hopane

C29/C30 hopane

benzo[e]pyrene (BeP), benzo[a]pyrene (BaP), perylene (Per), indeno(1,2,3-cd)pyrene (IP), dibenz(a,h)anthracene (DA), benzo(ghi)perylene (BP), and retene (Ret). Quality control samples were processed in a manner identical to the actual samples. A minimum of one procedural blank was used for each batch of 14 samples. Blank levels were no more than three times higher than the detection limit. The method detection limit for target PAHs ranged from 0.46 to 0.92 ng/g, except for naphthalene, which had a detection limit of 3.36 ng/g. Surrogate standards were used to determine the recovery of each sample and quantify the analytes. The surrogate standard recoveries satisfied acceptable ranges (40–120%). The certified reference materials (National Institute of Standards and Technology 1974b) were analyzed to monitor the performance of the analytical methods. The recoveries of the PAHs in the certified materials were within an acceptable range (certified values ±20%). The results reported in this study were based on dry weight. The wet:dry ratio was determined by calculating the weight loss after drying of approximately 2 g of homogenized sample at 105 ± 5 °C for more than 30 min.

1.84

Arch Environ Contam Toxicol

Arch Environ Contam Toxicol

(a) 35

Composition (%)

30 25 20 15 10 5

C0Nap C1Nap C2Nap C3Nap C4Nap Acl Act C0Flu C1Flu C2Flu C3Flu C0DBT C1DBT C2DBT C3DBT C0Phe Ant C1Phe C2Phe C3Phe C4Phe Ret Flr Pyr BaA C0Chr C1Chr C2Chr C3Chr BbF BkF BeP BaP Per IP DA BP

0

Compounds

(b) 35

Composition (%)

30 25 20 15 10 5

C0Nap C1Nap C2Nap C3Nap C4Nap Acl Act C0Flu C1Flu C2Flu C3Flu C0DBT C1DBT C2DBT C3DBT C0Phe Ant C1Phe C2Phe C3Phe C4Phe Ret Flr Pyr BaA C0Chr C1Chr C2Chr C3Chr BbF BkF BeP BaP Per IP DA BP

0

Compounds Fig. 2 PAHs composition profiles of spilled oil (a) and stranded oil (b). Refer to ‘‘Analytical Methods’’ section for abbreviation of PAHs

with C15 being most abundant. In contrast to the loss of low molecular weight alkyl homologues (C0- and C1-) of naphthalene, an increase of alkyl homologues (C2- and C3-) of dibenzothiophene was observed in the stranded oil (Fig. 2b). Representative source fingerprinting indices of the PAHs and biomarkers were well matched between the source oil and the stranded oil (Table 1). Temporal Changes and Spatial Distribution of the PAHs Concentration of the total PAHs in the oysters collected before the spill ranged from 210 to 4040 ng/g with an average of 1200 ng/g (Fig. 3a). Stations in close proximity to the petrochemical and steel complexes showed relatively higher concentrations compared with other stations. PAH levels drastically increased (by nearly 20 times at highest contaminated site) 4 days post-spill (627–81,000 ng/g with an average of 33,800 ng/g), and increased even higher over

the first 10 days post-spill (2000–110,000 ng/g with an average of 35,000 ng/g; Fig. 3b, c). PAHs levels then decreased over the following 30 days (923–63,000 ng/g with an average of 15,200 ng/g) and then to even lower levels over the following 60 days after the spill (515–41,000 ng/g with an average of 5900 ng/g; Fig. 3d, e). These temporal trends were similar with those previously found after the Hebei Spirit oil spill, where total PAHs in oyster right after the spill were approximately 40–500 times higher than pre-spill and decreased exponentially over time (Yim et al. 2012). The movement and stranding of spilled oil contaminated the entire Gwangyang Bay coastline to varying degrees. Notably, massive amounts of oil moved into the Sindeok beach (St. 4) near the spill site, and caused severe contamination in intertidal areas. Station 4 showed the highest PAH concentration in oysters (110,000 ng/g), and there were prolonged effects from the residual oil for 2 months after the spill. Due to strong tidal movements and easterly

123

Arch Environ Contam Toxicol

Total PAHs (ng/g dw.)

Total PAHs (ng/g dw.)

Hadong

80000

Namhae

100000

Yeosu

Total PAHs (ng/g dw.)

(a) 120000

60000 40000 20000

NS

0

120000

1

2

1

2

3

4

NS

NS

5

6

7

8

NS

NS

NS

NS

NS

NS

NS NS NS

9

10

11

12

13

14

15

NS

NS

NS

12

13

14

18

NS

NS

19

20

21

16

17

NS

NS

NS

NS

NS

16

17

18

19

20

21

(b)

100000 80000 60000 40000 20000

NS

0

120000

NS

NS

NS

NS

3

4

5

6

7

8

NS 9

10

11

15

(c)

100000 80000 60000 40000 20000

NS NS

0

NS

NS

NS

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

Total PAHs (ng/g dw.)

(d) 120000 100000 80000 60000 40000 20000 0

Total PAHs (ng/g dw.)

(e) 120000 100000 80000 60000 40000 20000 0

Fig. 3 Spatiotemporal changes of PAHs in bivalves after the Wu Yi San oil spill. Stations marked with ‘‘NS’’ represents stations that were not sampled. Each station number corresponds to the station shown in Fig. 1. a Pre-spill, b 4 days, c 10 days, d 30 days, e 60 days

123

Arch Environ Contam Toxicol

Temporal Changes in PAH Composition The characteristics of a bell-shaped pattern within the alkyl homologue series of PAHs (3- to 4-ring PAHs) and some pyrogenic 4- to 6-ring PAHs suggests a mixed source of pyrogenic and petrogenic in the pre-spill oysters (Fig. 5a;

(a)

40000

20000

10000

0

0

Depuraon phase

30000

Accumulaon phase

Alkylated PAHs (ng/g dw.)

50000

10

20

30

40

50

60

70

60

70

Duration after spill (days)

(b) 35000 30000

20000 15000 10000 5000 0

0

Depuraon phase

25000

Accumulaon phase

Target PAHs (ng/g dw.)

winter winds, the spilled oil spread rapidly, mainly into the eastern part of the bay. Compared with the western part of the bay, the eastern coastline was more highly contaminated. Twelve sites in Namhae showed highly elevated PAH concentrations; seven of these sites exceeded 40,000 ng/g. Other regions in Gwangyang Bay also exhibited significant oiling. Seven sites in Yeosu and one site in Hadong showed elevated PAHs in bivalves and similar temporal changes at the heavily contaminated sites. Station 8 which is located at the southern most of Gwangyang Bay showed no significant changes in PAH concentrations and is regarded as a reference site. The temporal trends of PAHs in bivalves after the spill corresponded with the timeline of the emergency oil spill responses. Initial cleanup activities to recover bulk oil at sea ended 4 days after the spill. Because most of the oil slick was stranded at shorelines, extensive shoreline cleanup had been performed for 10 days after the spill. The oiling status of the spilled oil on the affected areas had been monitored daily by a shoreline cleanup assessment team (SCAT), and the amount of visible oil on the shoreline was reported to be significantly decreased 10 days post spill. Secondary cleanup to remove residual oil continued for 30 days after the spill (KCG 2014). During extensive oil spill cleanups such as high- or lowpressure water flushing, surf washing, and tilling, stranded oil is resuspended, thus facilitating dispersion and dissolution (Lee et al. 1997; Lunel et al. 1996). Filter-feeding bivalves such as oysters and mussels are susceptible to water-column disturbances caused by cleanup events (Baumard et al. 1999). In this case, it could be inferred that the accumulation and depuration of petrogenic PAHs in oysters were largely governed by the duration and extent of the shoreline cleanup activities. Thus, the initial increasing trend of PAHs in oysters for 10 days post spill could be interpreted to be the accumulation phase, and decreasing trends for 2 months after the spill might represent the depuration phase. This was well demonstrated in temporal changes of PAHs at St. 4 and St. 18 (Fig. 4). In the accumulation phase, PAH concentrations drastically increased to as high as 45,000 and 80,000 ng/g, respectively, followed by a further increase to 65,000 and 110,000 ng/g, respectively in the 10 days post-spill. After that, PAHs were depurated exponentially to 7000 and 40,000 ng/g, respectively over 2 months, and they are expected to further decrease to pre-spill levels over time.

10

20

30

40

50

Duration after spill (days)

Fig. 4 Temporal changes of total PAHs accumulated in oysters compared with pre-spill data. a St. 4, b St. 18

Douglas et al. 2007; Stout 2007). Compared with pre-spill data, elevated proportions of naphthalene and dibenzothiophene were found in oysters collected four days after the spill (Fig. 5b). The bell-shaped dibenzothiophene and phenanthrene alkyl homologue series and the elevated proportions of naphthalene suggest a petrogenic source of contamination (Stout et al. 2002; Wang and Fingas 2003). The proportion of naphthalene subsequently decreased by an order of magnitude and that of dibenzothiophene increased significantly 10 days after the spill (Fig. 5c). Concentrations of the less persistent C1- and C2- alkyl homologues of dibenzothiophene and phenanthrene were substantially lower than the C3- and C4- homologues after 30 days (Fig. 5d). These temporal changes were more apparent in 60 days after the spill (Fig. 5e). Diagnostic ratios often are used as indicators for source identifications after oil spills. For example, alkylated dibenzothiophene over alkylated phenanthrene are degradation-resistant, which make them useful source-specific

123

C0Nap C1Nap C2Nap C3Nap C4Nap Acl Act C0Flu C1Flu C2Flu C3Flu C0DBT C1DBT C2DBT C3DBT C0Phe Ant C1Phe C2Phe C3Phe C4Phe Ret Flr Pyr BaA C0Chr C1Chr C2Chr C3Chr BbF BkF BeP BaP Per IP DA BP

Composition (%)

C0Nap C1Nap C2Nap C3Nap C4Nap Acl Act C0Flu C1Flu C2Flu C3Flu C0DBT C1DBT C2DBT C3DBT C0Phe Ant C1Phe C2Phe C3Phe C4Phe Ret Flr Pyr BaA C0Chr C1Chr C2Chr C3Chr BbF BkF BeP BaP Per IP DA BP

Composition (%)

C0Nap C1Nap C2Nap C3Nap C4Nap Acl Act C0Flu C1Flu C2Flu C3Flu C0DBT C1DBT C2DBT C3DBT C0Phe Ant C1Phe C2Phe C3Phe C4Phe Ret Flr Pyr BaA C0Chr C1Chr C2Chr C3Chr BbF BkF BeP BaP Per IP DA BP

Composition (%)

C0Nap C1Nap C2Nap C3Nap C4Nap Acl Act C0Flu C1Flu C2Flu C3Flu C0DBT C1DBT C2DBT C3DBT C0Phe Ant C1Phe C2Phe C3Phe C4Phe Ret Flr Pyr BaA C0Chr C1Chr C2Chr C3Chr BbF BkF BeP BaP Per IP DA BP

Composition (%)

C0Nap C1Nap C2Nap C3Nap C4Nap Acl Act C0Flu C1Flu C2Flu C3Flu C0DBT C1DBT C2DBT C3DBT C0Phe Ant C1Phe C2Phe C3Phe C4Phe Ret Flr Pyr BaA C0Chr C1Chr C2Chr C3Chr BbF BkF BeP BaP Per IP DA BP

Fig. 5 Temporal changes of PAHs compositions in oysters within the spill region. Refer to ‘‘Analytical Methods’’ section for abbreviation of PAHs. a Prespill, b 4 days, c 10 days, d 30 days, e 60 days

123 Composition (%)

Arch Environ Contam Toxicol

(a) 35

30

25

20

15

10 5

0

(b) Compounds

35

30

25

20

15

10

5

0

(c) Compounds

35

30

25

20

15

10

5

0

(d) Compounds

35

30

25

20

15

10

5

0

(e) 35

Compounds

30

25

20

15

10

5

0

Compounds

Arch Environ Contam Toxicol

1978). Just like other marine organisms, oysters undergo both phase I and II enzymatic biotransformation of xenobiotics (Livingstone 1985; Moore 1980; Buhler and Williams 1989; Foureman 1989). During phase I enzymatic biotransoformation, xenobiotics that are insufficiently polar cannot be excreted and therefore need to be conjugated before excretion. The presence of sulfur in dibenzothiophene allows it to be excreted more easily from the organism without the need for conjugation processes in phase II of enzymatic biotransformation (Campo et al. 2010).

8

C3DBT/C3Phe ratio

7 6 5 4 3

Stranded oil Pre-spill 4 days 10 days 30 days 60 days

2 1 0 0

1

2

3

4

5

6

C2DBT/C2Phe ratio

Fig. 6 Ratios showing source recognition indices

markers in abiotic samples (Douglas et al. 1996; Wang et al. 1999; Yim et al. 2011). These double ratios were plotted to see their applicability in oyster samples (Fig. 6). Double ratios (C2D/C2P vs C3D/C3P) gradually decrease with time. Interestingly, the temporal changes of these ratios lie on a line connecting two points, which start from stranded oil to pre-spill. This suggests the vulnerability of dibenzothiophene to depuration and/or metabolization compared with phenanthrenes in the oyster, which will be further discussed. Selective Depuration of Dibenzothiophene Although the fate of dibenzothiophene and phenanthrene is generally known to be similar in the abiotic environment, once they are accumulated in biota, their metabolic pathways are unclear. Accumulation rate of these two compounds were reported similar in previous studies (Eastmond et al. 1984; Herbes and Risi 1978), but the rates of depuration of these two compounds were different (Eastmond et al. 1984). Similar to reports by Eastmond et al. (1984), this study observed that dibenzothiophene was depurated more rapidly than phenanthrene from the organism. The more rapid depuration of dibenzothiophene than phenanthrene in this study might be due to (1) slight difference in polarity between dibenzothiophene and phenanthrene, and (2) difference in metabolic pathways in the organism. Parent and their respective alkyl homologues of dibenzothiophene and phenanthrene are known to have slight polarity differences (Meador et al. 1995), but this slight difference can significantly influence the enzymatic biotransformation processes in an organism (Moore 1980). Although oysters are known to have little metabolic capacities, these can significantly influence the rates of elimination of xenobiotics from the organism (Anderson

Depuration Rate and Half-Life of PAHs in Oysters The accumulation of organic contaminants is governed by type and duration of exposure, as well as by the metabolic activity of the organism (Djomo et al. 1996; McKim and Goeden 1982). The levels of contamination in an organism are driven by its ability to bioaccumulate contaminants through all routes of exposure from the ambient environment (Arnot and Gobas, 2006). When the organism is not consistently exposed to a source, the accumulated contaminants can be released, thus resulting in decreased levels of contamination. This process, also known as depuration, occurs when organisms resume normal activities, such as filtration, respiration, and reproduction (Cusson et al. 2005; Lee et al. 2008). The apparent first-order depuration of 16 PAHs, alkylated PAHs, and total PAHs were calculated by fitting the data into a linear least square analysis (Branson et al. 1975; Djomo et al. 1996; Neely 1979). In this study, depuration rate and half-life was determined using only substantially detected priority PAHs, such as naphthalene, fluorene, phenanthrene, fluoranthene, pyrene, benz[a]anthracene and chrysene for PAHs, and alkylated homologues of naphthalene, fluorene, phenanthrene, and chrysene for alkylated PAHs. ½Concentration at any given time ¼ ½Concentration at beginning of depuration  Kd t ½lnCa ¼ ½lnCt  Kd t Ca = concentration of contaminants at any given time, Ct = concentration of contaminants at the beginning of depuration, Kd = depuration rate constant, and t = time in terms of days. The depuration half-life (t1/2) also was computed. t1=2 ¼ ln 2=Kd Average depuration-rate constants for the target compound groups at all of the stations varied from -0.015 to -0.038; 16 PAHs showed the lowest rates, and the alkylated PAHs showed the highest rates (Table 2). As alkyl groups of

123

Arch Environ Contam Toxicol Table 2 Depuration-rates and half-lives of PAHs in oysters

Apparent depuration-rate constant (K) PAHs

-0.015

47.54

Alkylated PAHs

-0.038

18.23

Target ? alkylated PAHs

-0.037

18.72

naphthalenes include almost 50% of the initial volume of the alkylated PAHs, loss of naphthalene alkyl group series accounted for the high depuration rates of alkylated PAHs compared with the 16 PAHs. The half-life estimates for all of the stations within the spill region ranged from 47.5 days for the 16 PAHs, 18.7 days for the TPAHs, and 18.2 days for the alkylated PAHs. Because alkylated PAHs comprised more than 90% of the TPAHs, the TPAHs depuration rate and half-life was highly dependent on those of the alkylated PAHs. Human Health Risk Assessment As a basis for human health risk assessment, PAHs are separated into two subclasses based on their carcinogenicity: carcinogenic and noncarcinogenic (Reilly and York 2001). In application of the cancer slope factor, the toxic equivalency factor (TEF) approach often is used to estimate the carcinogenic activity of PAHs relative to BaP. Of the wide range of PAHs, only BaP has been well characterized in terms of toxicology and has been commonly used to determine ‘‘Levels of Concern’’ through the use of BaP equivalents (BaPE) (Yender et al. 2002). The following equation was used to determine levels of concern for public health (LoC: in ng/g lipid wt.) for carcinogenic PAH compounds (BaPE) potentially found in oysters (Reilly and York 2001; Yender et al. 2002). LoC ðBaPEÞ ¼ ðRL  BW  ATÞ=ðCSF  CR  EDÞ RL is the risk level. BW is the average consumer body weight. AT is the averaging time (i.e., life expectancy). CSF is the cancer slope factor of BaP. CR is the consumption rate (daily amount of seafood consumed). ED is the assumed exposure duration. Using regional statistics on the Korean lifestyle, LoC was calculated to be 3.35 ng/g lipid wt. Levels of TEF for the PAHs range from 0.001 to 5; BaP represents 1 as recommended by Nisbet and Lagoy (1992). Tissue concentrations of PAHs other than BaP were multiplied by their respective TEF and added to the tissue concentrations of BaP for the determination of BaPE. Cancer risk interpretations were performed by comparing LoC and BaPE values from each station (Fig. 7). BaPE concentrations exceeded the LoC for oysters in five stations (St. 4, 16, 17, 18, and 19) over the 10 days directly after the spill, which highlights the potential for human health risks. This exceedance was prolonged even after 60 days postspill. Korea has the most productive shellfish aquaculture

123

Half-life (t1/2; days)

farming areas in the world and supplies large amounts of the shellfish products to local and nonlocal consumers (Kwon et al. 2004; Yoon et al. 2003), and hence persistence of contaminants in shellfish suggests a need for further long-term investigation and health risk assessments in this region. Comparison with Previous Oil Spills To briefly compare the levels of contamination in oysters, we selectively compared several components between the Wu Yi San oil spill and previous large spills (Table 3). It is common to find a rapid increase followed by a decrease in petroleum derived PAHs in bivalves after an oil spill. However, this is exceptional for the Deepwater Horizon oil spill in the sense that it was an offshore deep-sea blowout. Unlike other spills, the oil released from the Deepwater Horizon blowout allows it to undergo more weathering processes before being stranded on shorelines. The increase in petroleum derived PAHs in bivalves after an oil spill is mainly due to dissolved and dispersed oil to which they were exposed, but in the Deepwater Horizon oil spill, no significant differences of petroleum derived PAHs were observed in bivalves that were exposed to the highly weathered emulsions. After the largest oil spill in Korea—the Hebei Spirit oil spill in 2007—advances in the fields of environmental health, ecological risk assessment, and development of new technologies provided better preparedness for responses to similar accidents that occur. Compared with other large spills, emergency oil spill responses were deployed within 24 h after the Wu Yi San oil spill. Right after deployment, multimedia samples were collected and analyzed for short- and long-term effects. The rapid emergency environmental impact assessment and consecutive post-spill monitoring along with prespill data provided better information of the fate and effects of the spill. This also allows the identification of accumulation and depuration phases in oyster which is the first field observation in oil spill cases, and contributes to protect human health from consuming contaminated seafood.

Conclusions After the Hebei Spirit oil spill on December 7, 2007, advances in national preparedness have been elevated, including prevention, response, and restoration of marine

Arch Environ Contam Toxicol

(a) 5

BaP equivalants

BaP equivalancy = 3.35 4

3

2

1

0 1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

(b) 5

BaP equivalants

BaP equivalancy = 3.35 4

3

2

1

0 1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

(c) 5

BaP equivalancy = 3.35 BaP equivalants

4

3

2

1

0 1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

(d) 5

BaP equivalancy = 3.35 BaP equivalants

4

3

2

1

0 1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

123

Arch Environ Contam Toxicol b Fig. 7 Temporal changes of benzo[a]pyrene equivalent concentra-

tions. Level of concern (LoC) was calculated based on statistics of Korean lifestyle. Each station number corresponds to the station shown in Fig. 1. a Pre-spill, b 4 days, c 10 days, d 30 days, e 60 days

Table 3 Comparison of the Wu Yi San oil spill with previous large spills Spill case

Wu Yi San

Date of spill

31st Jan 2014

24th March 1989

7th Dec 2007

20th April 2010

Spill type

Coastal pipeline leakage

Grounding of oil tankerm,r

Oil tanker collisionw,v

Deep-sea blowoutb,c,e

Type of oil spilled

Medium crude oil (Basrah Light) 900 m3

Medium crude oil (Alaskan North Slope)m,r 42,000 m3;m,r

Mixture of middle and heavy crude oil (Kuwait Export, Iranian Heavy, and Upper Zakum)v,w 12,547 m3;w

Light crude oil (Louisiana Sweet)b 500,000 m3;b

Application of cleanup

Mechanical only

Mechanical and chemicall,m

Mechanical and chemicalw,x

Mainly chemicalb,u

State of oil stranded on shorelines right after spill

Fresh to slightly weathered Right after spill

Slightly to highly weatheredo,p

Fresh to slightly weatheredv,w,

Highly weathereds,t

Approximately 1–3 dayso,r

Right after spillw,x

Approximately 25–40 daysj

Duration of cleanup for stranded oil

10–30 days

2 monthsh

1–2 yearsf,w

1–2 yearsi,k,n

Form of oil potentially exposed to bivalves (along/stranded on shorelines) Bivalve species

Dissolved and dispersed

Dissolved and dispersedo,p

Dissolved and dispersedf,x

Mainly emulsiona,d,g,t

Crassostrea gigas

Mytilus trossuluso,p

Crassostrea gigas

Crassostrea virginicab,d,q

2000–110,000

5000–234,000o

8000–100,000w

No significant changesd,q

Amount of oil spilled

Time for oil to reach shorelines (stranded)

Level of contamination in oysters post-spill (PAHs; ng/g dw)

Exxon Valdez

Hebei Spirit

Deepwater Horizon

a

Aeppli et al. (2012); bBeyer et al. (2016); cCamilli et al. (2010); dCarmichael et al. (2012); eDiercks et al. (2010); fHong et al. (2014); gLiu et al. (2012); hMearns (1996); iMichel et al. (2013); jOSAT (2011); kOwens et al. (2011); lPeterson (2001); mPetterson et al. (2003); nSantner et al. (2011); oShort and Babcock (1996); pShort and Harris (1996); qSoniat et al. (2011); rSpies et al. (1996); sSteffy et al. (2013); tStout et al. (2016); u Thibodeaux et al. (2011); vYim et al. (2011); wYim et al. (2012); xYim and Shim (2017)

oil spills. Despite these national efforts, negligence with regard to safety caused the Wu Yi San oil spill. This spill has reawakened awareness of the importance of emergency responses for marine disasters, and this study was conducted as a part of such emergency environmental responses. Multimedia, environmental monitoring has been conducted for 2 months; oysters clearly demonstrated that the temporal changes of petroleum hydrocarbons contamination in the intertidal environment. Levels of contamination in oysters corresponded with the time and duration of removal of residual oils during cleanup processes. These temporal changes of contaminants were identified as accumulation and depuration phases, which is the first field observation in oil spill cases.

123

Human health risk assessment demonstrated actual risks of consuming contaminated oysters produced in these regions. Risk levels decreased over time, but one site exceeded the LoC 2 months after the spill, suggesting the necessity of further, long-term investigations of shellfish products in these regions. Temporal changes of bioaccumulated PAHs in oyster and their depuration rates all indicated that there remained persistent effects from the residual oil after the Wu Yi San oil spill, thus emphasizing the importance of a long-term monitoring study. Acknowledgements This work was supported by the Project entitled ‘‘Oil Spill Environmental Impact Assessment and Environmental Restoration (PM59291)’’ funded by the Ministry of Oceans and Fisheries of Korea.

Arch Environ Contam Toxicol

References Aeppli C, Carmichael CA, Nelson RK, Lemkau KL, Graham WM, Redmond MC, Valentine DL, Reddy CM (2012) Oil weathering after the Deepwater Horizon disaster led to the formation of oxygenated residues. Environ Sci Technol 46:8799–8807 Anderson RS (1978) Benzo[a]pyrene metabolism in the American oyster Crassostrea virginica. EPA Ecol Res Ser Monogr (EPA600/3-78-009) Arnot JA, Gobas FAPC (2006) A review of bioconcentration factor (BCF) and bioaccumulation factor (BAF) assessments for organic chemicals in aquatic organisms. NRC Can Environ Rev 14:257–297 Baumard P, Budzinski H, Garrigues P, Narbonne JF, Burgeot T, Michel X, Bellocq J (1999) Polycyclic aromatic hydrocarbon (PAH) burden of mussels (Mytilus sp) in different marine environments in relation with sediment PAH contamination, and bioavailability. Mar Environ Res 47:415–439 Beyer J, Trannum HC, Bakke T, Hodson PV, Collier TK (2016) Environmental effects of the Deepwater Horizon oil spill: a review. Mar Pollut Bull 110:28–51 Branson DR, Blau GE, Alexander HC, Neely WB (1975) Bioconcentration of 2, 20 , 4, 40 -tetrachlorodiphenyl in rainbow trout as measured by an accelerated test. Trans Am Fish Soc 4:785–792 Buhler DR, Williams DE (1989) Enzymes involved in metabolism of PAH by fishes and other aquatic animals: oxidative enzymes (or phase I enzymes). In: Varanasi U (ed) Metabolism of polycyclic aromatic hydrocarbons in the aquatic environment. CRC Press, Boca Raton, pp 151–184 Camilli R, Reddy CM, Yoerger DR, Van Mooy BAS, Jakuba MV, Kinsey JC, McIntyre CP, Sylva SP, Maloney JV (2010) Tracking hydrocarbon plume transport and biodegradation at deepwater horizon. Science 330:201–204 Campo L, Rossella F, Pavanello S, Mielzynska D, Siwinska E, Kapka L, Bertazzi PA, Fustinon S (2010) Urinary profiles to assess polycyclic aromatic hydrocarbons exposure in coke-oven workers. Toxicol Lett 492(1):72–78 Carls MG, Babcock MM, Harris PM, Irvine GV, Cusick JA, Rice SD (2001) Persistence of oiling in mussel beds after the Exxon Valdez oil spill. Mar Environ Res 51(2):167–190 Carmichael RH, Jones AL, Patterson HK, Walton WC, Perez-Huerta A, Overton EB, Dailey M, Willett KL (2012) Assimilation of oilderived elements by oysters due to the Deepwater Horizon oil spill. Environ Sci Tehnol 46:12787–12795 Cossa D (1989) A review of the Mytilus sp. As a quantitative indicator of cadmium and mercury contamination in coastal waters. Oceanol Acta 12:417–432 Cusson M, Tremblay R, Daigle G, Roussy M (2005) Modeling the depuration potential of blue mussels (Mytilus sp.) in response to thermal shock. Aquaculture 250(1):183–193 Diercks A, Highsmith RC, Asper VL, Joung D, Zhou Z, Guo L, Shiller AM, Joye SB, Teske AP, Guinasso N, Wade TL, Lohrenz SE (2010) Characterization of sub-surface polycyclic aromatic hydrocarbons at the Deepwater Horizon site. Geophys Res Lett 37:6 Djomo JE, Garrigues P, Narbonne JF (1996) Uptake and depuration of polycyclic aromatic hydrocarbons from sediments by the zebrafish (Brachydanio Rerio). Environ Toxicol Chem 15(7):1177–1181 Douglas GS, Bence AE, Prince RC, McMillen SJ, Butler EL (1996) Environmental stability of selected petroleum hydrocarbon source and weathering ratios. Environ Sci Technol 30:2332–2339 Douglas GS, Emsbo-Mattingly SD, Stout SA, Uhler AD, McCarthy KJ (2007) Chemical fingerprinting methods. In: Murphy BL, Morisson RD (eds) Introduction to environmental forensics, 2nd edn. Academic, New York, pp 311–454

Eastmond DA, Booth GM, Lee ML (1984) Toxicity, accumulation, and elimination of polycyclic aromatic sulfur heterocycles in Daphnia magna. Arch Environ Contam Toxicol 13:105–111 Elmgren R, Hansson S, Larsson U, Sundelin B, Boehm PD (1983) The Tsesis oil spill: acute and long-term impact on the benthos. Mar Biol 73:51–65 Farrington JW (1991) Biogeochemical processes governing exposure and uptake of organic pollutant compounds in aquatic organisms. Environ Health Perspect 90:75–84 Farrington JW, Tripp BW, Tanabe S, Subramanian A, Sericano JL, Wade TL, Knap AH (2016) Edward D. Goldberg’s proposal of ‘‘the Mussel Watch’’: reflections after 40 years. Mar Pollut Bull 110(1):501–510 Foureman GL (1989) Enzymes involved in metabolism of PAH by fishes and other aquatic animals: hydrolysis and conjugation enzymes (or phase II) enzymes. In: Varanasi U (ed) Metabolism of polycyclic aromatic hydrocarbons in the aquatic environment. CRC Press, Boca Raton, pp 185–202 Herbes SE, Risi GF (1978) Metabolic alteration and excretion of anthracene by Daphnia pulex. Bull Environ Contam Toxicol 19:147 Hong SH, Kannan N, Yim UH, Choi J, Shim WJ (2011) Polychlorinated biphenyls (PCBs) in a benthic ecosystem in Gwangyang Bay, South Korea. Mar Pollut Bull 62:2863–2868 Hong S, Khim JS, Ryu J, Kang S, Shim WJ, Yim UH (2014) Environmental and ecological effects and recoveries after five years of the Hebei Spirit oil spill, Taean, Korea. Ocean Coast Manage 102:522–532 ITOPF (2014) Fate of marine oil spills. Technical information paper series 02. The International Tanker Owners Pollution Federation Limited, London, UK Jin Y, Hong SH, Li D, Shim WJ, Lee SS (2008) Distribution of persistent organic pollutants in bivalves from the northeast coast of China. Mar Pollut Bull 57(6):775–781 KCG (2014) Daily report of emergency oil spill response for the Wu Yi San oil spill, Korea Coast Guard, Republic of Korea KORDI (2003) A study on the management model for environmental pollution of special management area in Namhae: (I) Gwangyang Bay study. 2001–2003. Korea Ocean Research and Development Institute Report # BSPE 836-001577-7 Kwon H, Lee C, Jun B, Yun J, Weon S, Koopman B (2004) Recycling waste oyster shells for eutrophication control. Resour Conserv Recy 41:75–82 Lake JL, Norwood C, Dimock C, Bowen R (1979) Origins of polycyclic aromatic hydrocarbons in estuarine sediments. Geochim Cosmochim Ac 43:1847–1854 Lee K, Lunel T, Wood P, Swannell R, Stoffyn-Egli P (1997) Shoreline cleanup acceleration of clay-oil flocculation process. In: The 1997 international oil spill conference, pp 235–240 Lee RL, Lovatelli A, Ababouch L (2008) Bivalve depuration: fundamental and practical aspects. Food and Agriculture Organization. Fisheries Technical Paper Number 511. FAO of UN, Rome 139 Liu Z, Liu J, Zhu Q, Wu W (2012) The weathering of oil after the Deepwater Horizon oil spill: insights from the chemical compositions of the oil from the sea surface, salt marshes and sediments. Environ Res Lett 7(3):035302 Livingstone DR (1985) Responses of the detoxification/toxication enzyme systems of molluscs to organic pollutants and xenobiotics. Mar Pollut Bull 16(4):158–164 Lunel T, Lee K, Swannell R, Wood P, Rusin J, Bailey N, Halliwell C, Davies L, Sommerville M, Dobie A, Mitchell D, McDonagh M (1996) Shoreline cleanup during the Sea Empress incident: the role of surf washing (clay-oil flocculation), dispersants and bioremediation. In: Proceedings of the 19th arctic and marine oilspill program (AMOP) technical seminar, environment Canada, Ottawa ON, pp 1521–1540

123

Arch Environ Contam Toxicol McKim JM, Goeden HM (1982) A direct measure of the uptake efficiency of a xenobiotic chemical across the gills of brook trout (Salvelinus fontinalis) under normoxic and hypoxic conditions. Comp Biochem Phys C72:65–74 Meador J, Stein J, Reichert W, Varanasi U (1995) Bioaccumulation of polycyclic aromatic hydrocarbons by marine organisms. Rev Environ Contam Toxicol 143:79–165 Mearns AJ (1996) Exxon Valdez shoreline treatment and operations: Implications for responses, assessment, monitoring, and research. In: Proceedings of the Exxon Valdez oil spill symposium, American Fisheries Society Symposium 18, Maryland, USA Michel J, Owens EH, Zengel S, Graham A, Nixon Z, Allard T, Holton W, Reimer PD, Lamarche A, White M, Rutherford N, Childs C, Mauseth G, Challenger G, Taylor E (2013) Extent and degree of shoreline oiling: deepwater Horizon oil spill, Gulf of Mexico, USA. Plos ONE 8:9 MLTM (2013) Environmental impact assessment and environmental restoration of the Hebei Spirit oil spill. Technical Report 2013. Ministry of Land, Transport and Maritime Affairs, Republic of Korea Moore MN (1980) Cytochemical determination of cellular responses to environmental stressors in marine organisms. In: McIntyre AD, Pearce JB (eds) Biological effects of marine pollution and the problems of monitoring, vol 179. Rapp P V Reuns Cons Int Explor Mer, Copenhague, pp 7–15 Neely WB (1979) Estimating the rate constants for the uptake and clearance of chemicals by fish. Environ Sci Technol 13:1506–1510 Neff JM (2004) Bioaccumulation in marine organisms: effect of contaminants from oil well produced water. Elsevier Science, Amsterdam Nisbet ICT, Lagoy PK (1992) Toxic equivalency factors (TEFs) for polycyclic aromatic hydrocarbons (PAHs). Regul Toxicol Pharm 16:290–300 OSAT (2011) Summary report for fate and effects of remnant oil in the beach environment. Report to federal on-scene coordinator, Operational Science Advisory Team, p 64 Owens E, Taylor E, Graham A, Castle R (2011) Sand beach treatment studies and field trials conducted during the Deepwater HorizonMacondo response operation. In: Proceedings of the 2011 international oil spill conference, American Petroleum Institute, Washington DC Peterson CH (2001) The ‘‘Exxon Valdez’’ oil spill in Alaska: acute, indirect and chronic effects on the ecosystem. Adv Mar Biol 39:1–103 Petterson CH, Rice SD, Short JW, Esler D, Bodkin JL, Brenda EB, Irons BI (2003) Long-term ecosystem response to the Exxon Valdez oil spill. Science 302(5653):2082–2086 POSCO (1983) Evaluation of the environmental effect of Gwangyang steel mill. POSCO, Pohang, pp 83–1046 Reilly TI, York RK (2001) Guidance on sensory testing and monitoring of seafood for presence of petroleum taint following an oil spill. NOAA technical memorandum NOS OR & R 9. Seattle: Office of Response and Restoration, National Oceanic and Atmospheric Administration, p 109 Santner RM, Cocklan-Vendl J, Michel E, Owens, Stong B, Taylor E (2011) Shoreline treatment during the Deepwater HorizonMacondo response. In: Proceedings of the 2011 international oil spill conference, American Petroleum Institute, Washington DC, p9 Short JW, Babcock MM (1996) Prespill and postspill concentrations of hydrocarbons in mussels and sediments in Prince William Sound. In: Proceedings of the Exxon Valdez oil spill symposium, American Fisheries Society Symposium 18, Maryland, USA

123

Short JW, Harris PM (1996) Petroleum hydrocarbons in caged mussels deployed in Prince William Sound after the Exxon Valdez oil spill. In: Proceedings of the Exxon Valdez oil spill symposium, American Fisheries Society symposium 18, Maryland, USA Short JW, Heintz RA (1997) Identification of Exxon Valdez oil in sediments and tissues from Prince William Sound and the northwestern Gulf of Alaska based on a PAH weathering model. Environ Sci Technol 31:2375–2384 Soniat TM, King SM, Tarr MA, Thorne MA (2011) Chemical and physiological measures on oysters (Crassostrea virginica) from oil-exposed sites in Louisiana. J Shellfish Res 30:713–717 Spies RB, Rice SD, Wolfe DA, Wright BA (1996) The effects of the Exxon Valdez oil spill on the Alaskan coastal environment. In: Proceedings of the Exxon Valdez oil spill symposium, American Fisheries Society Symposium 18, Maryland, USA Steffy D, Nichols A, Morgan LJ, Gibbs R (2013) Evidence that the Deepwater Horizon oil spill caused a change in the nickel, chromium, and lead average seasonal concentrations occurring in sea bottom sediment collected from the eastern Gulf of Mexico continental shelf between the years 2009 and 2011. Water Air Soil Pollut 224:1–11 Stout SA (2007) Characterization and source of unknown ‘‘tar-like material’’ and ‘‘slag’’ in a former oil field in Compton, California. Environ Forens 8:265–282 Stout SA, Magar VX, Uhler RM, McCarthy KJ, Emsbo-Mattingly S (2002) Chemical fingerprinting of hydrocarbons. In: Murphy BL, Morrison RD (eds) Introduction to environmental forensics, 1st edn. Academic, New York, pp 137–260 Stout SA, Payne JR, Emsbo-Mattingly SD, Baker G (2016) Weathering of field-collected floating and stranded Macondo oils during and shortly after the Deepwater Horizon oil spill. Mar Pollut Bull 105:7–22 Thibodeaux LJ, Valsaraj KT, John VT, Papadopoulos KD, Pratt LR, Pesika NS (2011) Marine oil fate: knowledge gaps, basic research, and development needs; a perspective based on the Deepwater Horizon spill. Front Microbiol 5:12 Tunnell JW, Jr Dokken QR, Kindinger ME, Thebeau LC (1981) Effects of the Ixtoc I oil spill on the intertidal and subtidal infaunal populations along lower Texas coast barrier island beaches. In: Proceedings of the 1981 oil spill conference. American Petroleum Institute, Washington, DC, pp 467–475 Wang ZD, Fingas M (2003) Development of oil hydrocarbon fingerprinting and identification techniques. Mar Pollut Bull 47:423–452 Wang ZD, Fingas M, Li K (1994) Fractionation of a light crude-oil and identification and quantitation of aliphatic, aromatic, and biomarker compounds by GC-FID and GC-MS 1. J Chromatogr Sci 32:361–366 Wang ZD, Fingas M, Landriault M, Sigouin L, Grenon S, Zhang D (1999) Source identification and tracking of an unknown spilled oil from Quebec (1998) by unique biomarkers and diagnostic ratios of ‘‘source-specific-marker’’ compounds. In: Proceedings of the 22nd Arctic and Marine oilspill program technical seminar. Calgary, Alberta, pp 107–128 Yender R, Michel J, Lord C (2002) Managing seafood safety after an oil spill. Seattle: hazardous materials response division, office of response and restoration, National Oceanic and Atmospheric Administration, p 72 Yim UH, Ha SY, An JG, Won JH, Han GM, Hong SH, Kim M, Jung JH, Shim WJ (2011) Fingerprint and weathering characteristics of stranded oils after the Hebei Spirit oil spill. J Hazard Mater 197:60–69 Yim UH, Kim M, Ha SY, Kim S, Shim WJ (2012) Oil spill environmental forensics: the Hebei Spirit oil spill case. Environ Sci Technol 46(12):6431–6437

Arch Environ Contam Toxicol Yim UH, Shim WJ (2017) The Hebei Spirit oil spill, 2007. In: Fingas M (ed) Oil spill science and technology, 2nd edn. Edmonton, AB, Canada, pp 925–937 Yoon G, Kim B, Kim B, Han S (2003) Chemical-mechanical characteristics of crushed oyster-shell. Waste Manage 23: 825–834

Youngblood WW, Blumer M (1975) Polycyclic aromatic hydrocarbons in the environment: homologous series in soils and recent marine sediments. Geochim Cosmochim Ac 39:1303– 1314

123