长江口是我国重要的三大河口之一。近年来由于人口增多,工农业迅速发展,大量的有机污染物,通过直接排放、雨水冲刷以及大气沉降进入水体,并由于其低溶解性和疏水性,富集在沉积物中。沉积物中这些有机污染物并不是单一存在的,而是以混合物的形式存在,污染物之间又可能存在相互作用,产生相加或协同等复合效应,对人类和底栖生物造成的危害可能远大于单一污染物。多环芳烃(polycyclic aromatic hydrocarbons, PAHs)是一类环境中普遍存在的持久性有机污染物(persistent organic pollutants, POPs),对人体具有“三致”效应(致癌、致畸和致突变)[1-2],美国环境保护局(United States Environmental Protection Agency, US EPA)于1976年将16种PAHs列为优先控制污染物。PAHs通过工业废水排放、海上运输以及石油泄漏等途径进入长江口及其相邻水域,并富集在沉积物[3-4]。然而,PAHs及其他未知污染物在复合污染中对生物的毒性作用贡献仍不明确。2019年长江流域全面实施“长江大保护”战略,长江经济带11省市强力推进长江生态系统保护修复。因此,明确表层沉积物中复合污染的生物毒性作用,评估PAHs的毒性贡献,并确定长江流域污染物削减控制重点,了解污染物浓度、风险和毒性的变化,对确定区域特征污染物变化,进一步有针对性地制定污染控制策略有重要意义。
由于有机化合物以混合物的形式存在,很难根据它们各自的浓度来评估沉积物中有机污染物的危害。此外,可能还存在一些由于方法或技术限制而未被检测到的污染物。因此,Escher等[5-6]建立了Iceberg模型,该模型基于生物分析当量浓度(BEQ),利用细胞毒性实验,综合评估环境样品中已检测到的或未检测到的化学物质对复合毒性效应的贡献程度。它可以评估沉积物中有机污染物的潜在危害,并广泛应用于地表水[7]、沉积物[8]和废水样品[9-11]中有机污染物的毒性贡献评估。通过比较生物分析的BEQ值(BEQbio)和化学分析的BEQ值(BEQchem),确定已检测到或未检测到的化学物质对混合物毒性效应的贡献程度[7],从而明确环境介质中的特征污染。
目前,对长江口沉积物中有机污染物的研究,主要集中在对其组成、分布及来源的分析[3-4,12-15],而对其复合毒性以及各个组分贡献程度的研究却尚未有报道。因此,本研究在长江口近海海域采集表层沉积物样品,运用气相色谱-质谱联用仪测定长江口沉积物中16种PAHs的浓度。在此基础上,分析其时空分布特征、发光细菌毒性当量以及PAHs毒性贡献,并就其环境风险进行评估,以期为控制和削减该区域PAHs的污染提供理论基础及科学依据。
仪器:冷冻干燥机(FreeZone 6,美国LABCONCO公司)、电子天平(ML104T/02,梅特勒-托利多仪器(上海)有限公司,中国)、微波消解仪(CEM MARS CLASSIC,美国CEM公司)、气相色谱-质谱联用仪(Agilent 8890A/5975C,安捷伦科技(中国)有限公司)、干式氮吹仪(QYN100-2,上海巧跃电子有限公司,中国)、涡旋振荡器(H-101,上海康禾光电仪器公司,中国)、立式高压灭菌器(LDZX-30KPS,上海申安医疗器械厂,中国)、马弗炉(KSL-1200X-H,合肥科晶材料技术有限公司,中国)、电热鼓风干燥箱(DHG-9070A,上海一恒科学仪器有限公司,中国)、低速大容量多管离心机(LXJ-IIB,上海安亭科学仪器厂,中国)、超声波清洗器(KQ3200DB,昆山市超声仪器有限公司,中国)、多功能微量移液器(Transferpette® S-12 Dig. 30~300 μL,普兰德(上海)贸易有限公司,中国)、纯水机(Milli-Q® Express 40,默克化工技术(上海)有限公司,中国)、医用低温保存箱(DW-86L578J,成都壹科医疗器械有限公司,中国)。
试剂:正己烷、丙酮(色谱纯,默克化工技术(上海)有限公司,中国);二氯甲烷(色谱纯,上海阿达玛斯试剂有限公司,中国);层析硅胶(60~100目,分析纯,上海阿拉丁生化科技股份有限公司,中国)置于马弗炉中650 ℃下灼烧4 h,冷却后转移至干燥器中待用;0.7 μm玻璃纤维滤膜(GF/F,英国Whatman公司);铜片(高纯试剂,5N,国药集团化学试剂有限公司,中国);高纯氮气(99.99%,上海利旦工业气体有限公司,中国);氯化钠(NaCl,分析纯,上海柯灵斯试剂有限公司,中国);六水氯化镁(MgCl2·6H2O,分析纯,上海阿拉丁生化科技股份有限公司,中国);氯化钾(KCl,分析纯,99.5%,上海阿拉丁生化科技股份有限公司,中国);3-N-吗啉丙烷磺酸(MOPS,上海阿拉丁生化科技股份有限公司,中国);氢氧化钠(NaOH,分析纯,上海阿拉丁生化科技股份有限公司,中国);盐酸(HCl,优级纯,上海柯灵斯试剂有限公司,中国);二甲基亚砜(DMSO,纯度>99.5%,西格玛奥德里奇(上海)贸易有限公司,中国);费氏弧菌冻干粉(北京海富达科技有限公司,中国);16种多环芳烃混合标准溶液(浓度为2 000 μg·mL-1,BePure)购自Sigma-Aldrich公司(NIST1647F),于4 ℃冰箱中保存待用。
于2019年和2020年12月对长江口海域沉积物进行了2次调查采集(图1)。其中2019年航次共调查采集17个站点(S1~S17),2020年航次共调查7个站点(L1~L7)。利用彼得森抓斗式采泥器(DXCN1/40,厦门登迅仪器设备有限公司)采集表层6 cm以上的沉积物,将其置于干净的铝箔袋中密封。置于加入冰袋的保温箱中,运回实验室后于-20 ℃冰箱保存。
图1 2019—2020年长江口采样点图
Fig. 1 Map of sampling points in the Yangtze
Estuary from 2019 to 2020
将沉积物样品冷冻干燥,用研钵研磨,过100目筛并存放于铝箔袋中,避光储藏待测。
1.3.1 沉积物中PAHs的萃取
称取3 g冻干并过筛的沉积物样品,用滤纸包好置于微波萃取管中,加入15 mL二氯甲烷,萃取完成后,将萃取溶液转移至10 mL玻璃离心管内,再用5 mL二氯甲烷溶液洗涤萃取管2次,转移至玻璃离心管内。再将2次萃取液浓缩并合并。每个离心管中加入1片铜片,浓缩至2 mL左右过滤。用硅胶柱净化上述溶液,二氯甲烷洗脱(10 mL),最后用氮气吹至近干,正己烷定容至0.5 mL,待测。
1.3.2 沉积物复合污染物的萃取
称取3 g冻干并过筛的沉积物样品,用滤纸包好置于微波萃取管中,加入10 mL二氯甲烷,萃取完成后,将萃取液转移至30 mL玻璃离心管内;再向微波萃取管中加入10 mL丙酮∶正己烷(V∶V=1∶1)混合液,萃取完成后,将萃取液转移至的玻璃离心管内,用5 mL丙酮∶正己烷(V∶V=1∶1)混合液洗涤萃取管2次,转移至玻璃离心管内。然后以3 000 r·min-1离心10 min;离心结束后,用0.7 μm玻璃纤维过滤器过滤并浓缩至2 mL左右合并,再氮吹至近干,用丙酮∶正己烷(V∶V=1∶1)混合液定容至1 mL,待用。
采用安捷伦公司气相色谱-质谱联用仪(Agilent 8890A/5975C)对样品中16种优先控制PAHs单体进行定量测定。色谱柱为HP-5-MS(30 m×0.25 mm×0.25 μm),载气为氦气,流速为1 mL·min-1,进样口温度260 ℃,进样量1 μL,进样方式不分流;柱温:70 ℃保持1 min,然后以10 ℃·min-1的速率升至160 ℃ ,再以5 ℃·min-1的速率升至280 ℃,保持5 min,最后以20 ℃·min-1的速率升至280 ℃,保持5 min。质谱条件:离子源类型为EI;离子源温度为280 ℃;溶剂延迟为5 min;扫描方式为选择离子;扫描类型为MS2全扫,MSM保持5 min,MS1 SIM保持36.5 min;扫描时间为107 ms。
16种US EPA规定的优先控制的PAHs分别是萘(NAP)、二氢苊(ANY)、芴(FLU)、苊(ANA)、菲(PHE)、蒽(ANT)、荧蒽(FLT)、芘(PYR)、(CHR)、苯并[a]蒽(BaA)、苯并[b]荧蒽(BbF)、苯并[k]荧蒽(BkF)、苯并[a]芘(BaP)、二苯并[a,h]蒽(DBA)、茚并[1,2,3-c,d]芘(IPY)和苯并[g,h,i]苝(BPE)。为保证实验结果的准确性和科学性,在实验开始前,对长江口沉积物样品进行了16种PAHs的加标回收率实验,其回收率在72%~116%之间,相对标准偏差均<15%,方法检出限为0.61~2.75 ng·g-1(以湿质量计)。所有样品的测定均设置空白实验以及2个平行,污染物浓度测定结果用平均值±标准偏差表示。
采用美国伯腾公司Synergy H1多功能微孔板检测仪对底泥样品总体毒性进行定量测定。向1 L纯水中加入0.3 g KCl、2.035 g MgCl2·6H2O、0.3 g MOPS、20 g NaCl和0.2 g NaOH配制缓冲液,将pH调节至7.0±0.2,高温灭菌后储存备用。取200 μL沉积物提取液转移到棕色小瓶中氮吹至干,将残留物重新溶解在DMSO中,并转移到缓冲液中。在透明96孔板中进行梯度稀释后,将100 μL缓冲液中的样品加入到白色96孔板中的50 μL的费氏弧菌(Vibrio fischeri)菌液中。在加入样品之前和培养30 min后,测量细菌的发光度。
按照ISO标准方法11348-3[16]和Escher等[17]计算生物发光抑制率的方法,用公式(1)计算效应中值浓度(EC50),其中s为浓度-效应曲线的斜率[18]。每种化学品的EC50值以mol·L-1为单位,而环境样品的效应浓度(EC)用相对富集因子(REF)表示[19-20]。样品REF(公式(2))是生物测定法(公式(3))的稀释因子与样品富集因子(公式(4))的乘积。沉积物提取物(kgsed, dw的富集因子为以提取物最终体积提取的沉积物质量,并根据公式(4)计算。每种生物测定的稀释因子均使用公式(3)计算[21-22]。
抑制率
(1)
REF=稀释因子生物测定 × 富集因子沉积物
(2)
稀释因
(3)
富集因
(4)
样品中化合物i的相对效应效价(REPi)用于计算不同化学物质之间的相互作用以及预测已知化学混合物的混合作用,由公式(5)获得。其中,EC50(参照化合物)指的是参照化合物效应浓度,EC50(i)指的是化合物i的效应浓度[6-7,11,23]。其中,本文中所用的参照化合物为苯酚溶液,化合物i为多环芳烃各单体,其EC50来自于相关参考文献[24]。
(5)
发光细菌毒性测试得到的生物分析当量BEQ值(BEQbio),可以由参照化合物的EC50除以样品的EC50计算得到(公式(6))[6-7,11,23]。
(6)
化学分析得到的化学分析当量BEQ值(BEQchem)可以计算为REPi与所有被测化合物i的浓度(Ci)的乘积之和(公式(7))[5-6,11,25]。
(7)
将BEQbio和BEQchem进行比较,确定检测到的化学物质对混合物效应的贡献程度(公式(8))[5, 10]。
贡献百分比
(8)
2019年和2020年长江口表层沉积物中16种PAHs的总浓度(ΣPAHs)的时空分布如图2所示。2019年长江口表层沉积物中16种PAHs的ΣPAHs范围为32.84~283.47 ng·g-1,平均值为123.16 ng·g-1。站点S3(283.47±29.94) ng·g-1浓度最高,其次是站点S7(279.07±16.75) ng·g-1和S16(234.15±0.80) ng·g-1,其原因可能是S3和S16这2个站点都靠近渔港,它们分别紧邻吕四渔港和沈家门渔港,船只过往频繁,汇集了船舶燃油及工业排放等多来源的污染物[26-27];而站点S7浓度较高是由于该站点ANY(30.38±1.04) ng·g-1、BkF(5.99±0.28) ng·g-1浓度较高。ANY与石油源污染有关,BkF与液体化石燃料燃烧有关[28],很可能是船舶燃油燃烧以及燃油泄漏导致该处的多环芳烃浓度较高。站点S6(43.86±3.54) ng·g-1和S9(32.84±1.30) ng·g-1附近海域PAHs浓度最低,主要由于该区域远离陆地,受陆源影响较小,并且没有突发污染事故发生。其他采样点沉积物中浓度范围为66.64~106.59 ng·g-1 (S1=(93.18±0.070) ng·g-1、S2=(83.33±5.47) ng·g-1、S4=(66.64±2.41) ng·g-1、S5=(106.59±0.35) ng·g-1、S8=(75.66±2.18) ng·g-1、S10=(102.62±6.89) ng·g-1、S14=(97.4±12.74) ng·g-1、S15=(95.96±9.03) ng·g-1和S17=(88.08±2.81) ng·g-1),空间分布均匀,点位之间无显著差异(OneWay ANOVA, P<0.05),说明研究区域的大部分点位受到长江口水流冲刷的影响,PAHs浓度得到一定程度的稀释。
2020年长江口表层沉积物中16种PAHs的ΣPAHs范围为66.93~132.64 ng·g-1,平均值为98.91 ng·g-1。2020年站点L6(132.64±9.95) ng·g-1和L7(130.14±8.55) ng·g-1的表层沉积物中PAHs浓度比较高,与2019年相比,2020年L6(2019年S15)附近区域PAHs浓度有所上升,可能存在某些污染源排放或突发原油泄漏。而在站点L3(68.01±5.8) ng·g-1和L5(66.93±6.22) ng·g-1处,沉积物中PAHs浓度有所降低。从时间变化来看,2020年PAHs的总体浓度有所降低。2019年的S1和S11站点分别与2020年的L1和L3站点重合。2020年的PAHs浓度在站点L1为(89.99±0.79) ng·g-1,在L3为(68.01±5.8) ng·g-1,比2019年数据S1=(93.18±0.070) ng·g-1和S11=(146.88±10.75) ng·g-1有所降低,可能与2020年新型冠状病毒疫情暴发,沿海港口受到影响,PAHs排放相对较少有关。将本研究与长江口其他年份调查研究对比,发现本研究PAHs浓度与2013年长江口及浙江省沿岸海域PAHs(31.8~384) ng·g-1 [3]浓度相当,并远低于2009—2010年的调查结果[29],呈现出一种逐年降低的趋势,可能与“长江大保护”的实施有关,长江经济带的生态环境得到逐年改善。
图2 2019—2020年长江口表层沉积物中PAHs含量时空分布
注:(a) 2019年,(b) 2020年;ΣPAHs为多环芳烃总浓度。
Fig. 2 Temporal and spatial distribution of PAHs in the sediments from the Yangtze River Estuary in 2019 and 2020
Note: (a) In 2019, (b) In 2020; ΣPAHs stands for the total concentration of polycyclic aromatic hydrocarbons.
根据沉积物质量基准(SQGs)提出的风险效应低值(effects range low, ERL)和风险效应高值(effects range median, ERM),本研究进一步评估了调查海域表层沉积物中PAHs的生态风险(表1)。结果表明,2019年部分点位的FLU浓度(S11~S17)和ANA浓度(S3、S10~S11、S14)以及2020年部分点位的FLU浓度(L1~L2、L6~L7)和ANA浓度(L3)处于ERL和ERM之间,说明以上站点表层沉积物中PAHs可以引发潜在生态风险。剩余站位的PAHs各单体浓度均低于ERL值,说明这些站点海域的沉积物中PAHs引发生态风险的可能性较低。2020年长江口沉积物的生态风险略低于2019年,这也与前文对PAHs总浓度的分析结果一致。对于BbF、BkF、IPY和BPF,由于没有最低安全值,因此这些物质只要存在对环境就会产生毒副作用[3, 30]。而在本次调查研究中,这些组分在各个站点(2019年S1~S17;2020年L4~L7)有部分检出,其中检出率最高的是BbF,其次是BkF,这可能是海上来往船只燃油燃烧导致的,需引起重视。
表1 长江口沉积物中PAHs的生态风险评估
Table 1 Ecological risk assessment of PAHs in sediments of the Yangtze Estuary
PAHs沉积物质量基准[31]Sediment quality benchmark[31]2019年In 20192020年In 2020ERL/(ng·g-1)ERM/(ng·g-1)浓度范围/(ng·g-1)Concentration range/(ng·g-1)风险程度Degree of risk浓度范围/(ng·g-1)Concentration range/(ng·g-1)风险程度Degree of riskNAP1602 100nd~60.07低风险 Low risknd~15.90低风险 Low riskANY44640nd~37.93低风险 Low risknd~5.80低风险 Low riskFLU19540nd~61.60部分潜在风险Some potential risksnd~66.27部分潜在风险Some potential risksANA16500nd~25.52部分潜在风险Some potential risksnd~29.41部分潜在风险Some potential risksPHE2401 500nd~2.29低风险 Low risknd~4.58低风险 Low riskANT85.31 100nd~3.82低风险 Low risknd~8.91低风险 Low riskFLT6005 100nd~47.14低风险 Low risknd~13.35低风险 Low riskPYR6652 600nd~213.38低风险 Low risk1.36~14.34低风险 Low riskCHR3842 800nd~233.51低风险 Low risk3.30~23.82低风险 Low riskBaA2611 600nd~23.58低风险 Low risknd~30.47低风险 Low riskBbFNANAnd~14.4部分潜在风险Some potential risksnd~31.96部分潜在风险Some potential risksBkFNANAnd~5.99部分潜在风险Some potential risksnd~8.85部分潜在风险Some potential risksBaP4301 600nd低风险 Low risknd~4.58低风险 Low riskDBA63.4260nd低风险 Low risknd~3.80低风险 Low riskIPYNANAnd低风险Low risknd~7.16部分潜在风险Some potential risksBPFNANAnd低风险Low risknd~2.90部分潜在风险Some potential risks
注:PAHs为多环芳烃;NAP为萘;ANY为苊烯;FLU为芴;ANA为苊;PHE为菲;ANT为蒽;FLT为荧蒽;PYR为芘;CHR为;BaA为苯并[a]蒽;BbF为苯并[b]荧蒽;BkF为苯并[k]荧蒽;BaP为苯并[a]芘;DBA为二苯并[a,h]蒽;IPY为茚并[1,2,3-cd]芘;BPF为苯并[g,h,i]苝;nd表示未检出;NA表示无效;ERL为风险效应低值(effects range low);ERM为风险效应高值(effects range median);若PAHs的含量小于ERL,则产生负面生态效应的可能性较小;若PAHs的含量在两者之间,则具有潜在的生态风险;若PAHs的含量大于ERM,则可能产生严重的生态风险[32];浓度单位(ng·g-1)以干质量计。
Notes: PAHs stands for polycyclic aromatic hydrocarbons; NAP stands for naphthalene; ANY stands for acenaphthylene; FLU stands for fluorene; ANA stands for acenaphthene; PHE stands for phenanthrene; ANT stands for anthracene; FLT stands for fluoranthene; PYR stands for pyrene; CHR stands for chrysene; BaA stands for benz[a]anthracene; BbF stands for benzo[b]fluoranthene; BkF stands for benzo[k]fluoranthene; BaP stands for benzo[a]pyrene; DBA stands for dibenz[a,h]anthracene; IPY stands for indeno[1,2,3-cd]pyrene; BPF stands for benzo[g,h,i]perylene; nd indicated not detected; NA indicated not available; ERL stands for effects range low; ERM stands for effects range median; if the content of PAHs is less than ERL, the possibility of negative ecological effects is less likely; if the content of PAHs is between the two, there are potential ecological risks; if the content of PAHs is greater than ERM, serious ecological risks may occur [32]; concentration units (ng·g-1) are based on dry mass.
本研究采用发光细菌毒性测试法评估长江口表层沉积物中有机污染物复合污染的细胞毒性效应,其量化指标为生物分析当量浓度(BEQbio)。长江口表层沉积物毒性效应的时空分布如图3所示,不同采样位点的BEQbio浓度存在一定的空间差异。在2019年采集的表层沉积物中,BEQbio值的范围在30.45~604.62 mg·kg-1之间,平均值为128.20 mg·kg-1。最大值出现在S6点,为604.62 mg·kg-1;其次是S5和S11,其BEQbio分别为427.13 mg·kg-1和267.15 mg·kg-1;最小值出现在S14处,其BEQbio为30.45 mg·kg-1;在站点S1~S4、S7~S10、S13、S15~S16处,沉积物中BEQbio含量水平相对较低,且较为接近,空间分布也较为均匀,表明该区域沉积物受外界有机污染影响较小。2019年S6处的BEQbio值最大,但其ΣPAHs浓度(43.86±3.54) ng·g-1却很低,表明PAHs并不是该站点海域主要的污染贡献者。该站点远离陆地,受陆源影响较小,污染物可能来源于洋流携带、大气沉降或某种污染物运输过程的突然泄漏,因此需要进一步筛选该区域的特征污染物。2020年采集的表层沉积物中,BEQbio值的范围在23.15~122.87 mg·kg-1之间,平均值为98.91 mg·kg-1。最大出现在L1点,为122.87 mg·kg-1;其次是L4,其BEQbio为81.06 mg·kg-1。最小值出现在L5处,其BEQbio为23.15 mg·kg-1。
2019年的S1和S11站点分别与2020年的L1和L3站点重合。二者相比,2020年的细胞毒性BEQbio在站点L1为(122.87±18.29) mg·kg-1,在L3为(37.16±3.77) mg·kg-1,与2019年数据S1=(45.30±2.62) mg·kg-1和S11=(267.15±4.60) mg·kg-1相比,L1处有所升高,L3处有所降低。这与PAHs浓度变化有所不同,说明PAHs总浓度与底泥的复合毒性不存在明显的相关性。
图3 长江口表层沉积物中有机提取物的毒性当量BEQbio的空间分布
注:(a) 2019年,(b) 2020年;BEQbio为生物分析当量浓度。
Fig. 3 Spatial distribution of toxic equivalent BEQbio of organic extracts in the surface sediments from the Yangtze River Estuary
Note: (a) In 2019, (b) In 2020; BEQbio stands for biological analysis equivalent concentration.
本研究使用浓度加和概念将检测到的化学物质的水平转换为BEQchem,并计算了沉积物提取物中PAHs的毒性当量及毒性贡献(表2)。结果发现,BEQchem在0.39~26.51 mg·kg-1之间。其中,2019年长江口表层沉积物中PAHs的毒性当量BEQchem在0.39~26.51 mg·kg-1之间,最大值出现在S2处,该处的PAHs毒性贡献占比最高,为37.59%,其中低环芳烃ANY(BEQchem=(23.75±31.43) mg·kg-1)的毒性贡献最多,约占89.60%,表明该站点受到石油源的污染较重;最小值出现在S9处,该处的PAHs毒性贡献占比相对较低,为1.00%,表明该站点受到PAHs和有机污染物影响较小。2020年长江口表层沉积物中PAHs的毒性当量BEQchem在0.56~7.47 mg·kg-1之间,最大值出现在L3处,该处的PAHs毒性贡献占比最高,为20.10%,其中芳烃BaA(BEQchem=(0.95±0.17) mg·kg-1)贡献最多,约占12.69%,该站点附近的海上油井内喷出的天然气发生自燃,导致该处的高环芳烃浓度较高;最小值出现在L5处,该处的PAHs毒性贡献占比相对较低,为2.42%,表明该站点受PAHs及有机污染物影响较小。
在2019年长江口沉积物的细胞毒性测试中,S6点BEQbio最高,但是通过对其PAHs的BEQchem计算,可知PAHs对其贡献极少,仅为0.081%,其他未知的有机污染物占毒性贡献的99.9%,说明PAHs不是该站点的特征污染物。2020年的L1点也存在这种情况。
在长江口沉积物毒性当量浓度中PAHs所占比例较小,将本研究中PAHs的毒性贡献数据与其他研究进行了比较。Hwang等[33]的研究表明,PAHs在沉积物中的毒性贡献在Yellow Sea(韩国管辖海域)为0.02%~9.4%(1.8%)、Yellow Sea(中国管辖海域)为0.09%~246%(18.9%)、中国渤海为0.30%~1 991%(93.3%);Jin等[34]的研究表明,PAHs在PM2.5中的毒性贡献为北京(26.5%)、广州(16.7%);而本研究2019年和2020年PAHs的毒性贡献分别为4.46%和4.25%,与以上研究结果相比,本研究PAHs的毒性贡献偏低。该结果也表明,本研究中检测到的PAHs仅能解释所观察到的复合毒性效应的一小部分,还需要进一步对其他未检测的化学物质进行测试分析。
表2 长江口表层沉积物中有机提取物的毒性当量BEQ的质量浓度
Table 2 The mass concentration of the toxic equivalent BEQ of organic extracts in
surface sediments at the Yangtze River Estuary (mg·kg-1)(以干质量计 Based on dry mass)
采样年份Year of sampling采样点位Point of sampling生物分析当量浓度(BEQbio)/(mg·kg-1)Biological analysis equivalent concentration (BEQbio)/(mg·kg-1)化学分析当量浓度(BEQchem)/(mg·kg-1)Chemical analysis equivalent concentration (BEQchem)/(mg·kg-1)贡献率/%Contribution rate/%2019S145.30±2.621.85±0.134.08S270.52±0.9226.51±34.4537.59S3100.19±8.372.30±0.272.30S462.55±0.681.72±0.0462.75S5427.13±42.292.02±0.180.47S6604.62±74.820.49±0.0360.081S759.01±4.722.05±0.163.47S852.09±4.261.72±0.0273.31S939.00±4.750.39±0.0241.00S1040.53±4.142.72±0.116.71S11267.15±4.601.91±0.0500.71S12133.97±6.441.42±0.0551.06S1347.68±1.442.07±0.414.34S1430.45±0.650.87±0.162.86S1591.26±7.431.64±0.0691.80S1668.21±9.131.35±0.0621.98S1739.82±2.150.50±0.0301.262020L1122.87±18.291.18±0.0980.96L263.53±2.391.78±0.0372.80L337.16±3.777.47±0.08120.10L481.06±10.000.96±0.0451.18L523.15±2.310.56±0.0252.42L680.37±9.020.91±0.00581.13L758.18±8.650.68±0.0221.17
综上所述,本研究表明:
(1) 长江口及其邻近海域表层沉积物中PAHs分析结果表明,2019年采集的表层沉积物样品中PAHs在渔港附近呈现较高浓度,2020年则在舟山岛附近呈现较高浓度;与2019年相比,2020年PAHs的总体浓度有所降低,且生态风险也略低于2019年。
(2) 2019年采集的表层沉积物中,BEQbio值在30.45~604.62 mg·kg-1之间(平均值为128.20 mg·kg-1)。其中,S6处毒性最强,S14处毒性最弱;而2020年的BEQbio值在23.15~122.87 mg·kg-1之间(平均值为66.62 mg·kg-1),其中,L1处毒性最强,L5处毒性最弱。
(3)在采集的表层沉积物中,2019年PAHs毒性贡献占比最高的是站点S2(37.59%),其主要贡献者是来自石油源的低环芳烃ANY,约占S2的89.60%;而2020年的是站点L3(20.10%),其主要贡献者是来自附近海上油井内喷出的天然气发生自燃产生的芳烃BaA,约占L3的12.69%。在长江口沉积物毒性当量浓度中,2019年和2020年的PAHs所占比例较小,其平均占比分别为4.46%和4.25%,表明检测到的PAHs仅能解释所观察到的复合毒性效应的一小部分,还需要进一步对其他未检测的化学物质进行测试分析。
通讯作者简介:王茜(1987—),女,博士,讲师,主要研究方向为有机污染物的环境行为、生物积累以及环境污染修复等。
[1] Hawliczek A, Nota B, Cenijn P, et al. Developmental toxicity and endocrine disrupting potency of 4-azapyrene, benzo[b]fluorene and retene in the zebrafish Danio rerio [J]. Reproductive Toxicology, 2012, 33(2): 213-223
[2] Botsou F, Hatzianestis I. Polycyclic aromatic hydrocarbons (PAHs) in marine sediments of the Hellenic coastal zone, eastern Mediterranean: Levels, sources and toxicological significance [J]. Journal of Soils and Sediments, 2012, 12(2): 265-277
[3] 母清林, 方杰, 邵君波, 等. 长江口及浙江近岸海域表层沉积物中多环芳烃分布、来源与风险评价[J]. 环境科学, 2015, 36(3): 839-846
Mu Q L, Fang J, Shao J B, et al. Distribution, sources and risk assessment of polycyclic aromatic hydrocarbons (PAHs) in surface sediments of Yangtze Estuary and Zhejiang coastal areas [J]. Environmental Science, 2015, 36(3): 839-846 (in Chinese)
[4] 欧冬妮, 刘敏, 许世远, 等. 长江口滨岸水和沉积物中多环芳烃分布特征与生态风险评价[J]. 环境科学, 2009, 30(10): 3043-3049
Ou D N, Liu M, Xu S Y, et al. Distribution and ecological risk assessment of polycyclic aromatic hydrocarbons in overlying waters and surface sediments from the Yangtze estuarine and coastal areas [J]. Environmental Science, 2009, 30(10): 3043-3049 (in Chinese)
[5] Escher B I, Stapleton H M, Schymanski E L. Tracking complex mixtures of chemicals in our changing environment [J]. Science, 2020, 367(6476): 388-392
[6] Escher B I, Neale P A, Villeneuve D L. The advantages of linear concentration-response curves for in vitro bioassays with environmental samples [J]. Environmental Toxicology and Chemistry, 2018, 37(9): 2273-2280
[7] Neale P A, Ait-Aissa S, Brack W, et al. Linking in vitro effects and detected organic micropollutants in surface water using mixture-toxicity modeling [J]. Environmental Science & Technology, 2015, 49(24): 14614-14624
[8] Hu X X, Shi W, Yu N Y, et al. Bioassay-directed identification of organic toxicants in water and sediment of Tai Lake, China [J]. Water Research, 2015, 73: 231-241
[9] Nivala J, Neale P A, Haasis T, et al. Application of cell-based bioassays to evaluate treatment efficacy of conventional and intensified treatment wetlands [J]. Environmental Science: Water Research & Technology, 2018, 4(2): 206-217
[10] Hebert A, Feliers C, Lecarpentier C, et al. Bioanalytical assessment of adaptive stress responses in drinking water: A predictive tool to differentiate between micropollutants and disinfection by-products [J]. Water Research, 2018, 132: 340-349
[11] König M, Escher B I, Neale P A, et al. Impact of untreated wastewater on a major European river evaluated with a combination of in vitro bioassays and chemical analysis [J]. Environmental Pollution, 2017, 220: 1220-1230
[12] 刘敏, 侯立军, 邹惠仙, 等. 长江口潮滩表层沉积物中多环芳烃分布特征[J]. 中国环境科学, 2001, 21(4): 343-346
Liu M, Hou L J, Zou H X, et al. Distribution characteristics of polycyclic aromatic hydrocarbons (PAHs) in surface sediments of tidalflats of the Yangtze Estuary [J]. China Environmental Science, 2001, 21(4): 343-346 (in Chinese)
[13] 尹方, 黄宏, 刘海玲, 等. 长江口及毗邻海域沉积物中多环芳烃分布、来源及风险评价[J]. 安全与环境学报, 2010, 10(3): 94-98
Yin F, Huang H, Liu H L, et al. PAHs in surface sediments from Yangtze River Estuary and its adjacent sea area: The distribution, sources and risk assessment [J]. Journal of Safety and Environment, 2010, 10(3): 94-98 (in Chinese)
[14] 沈小明, 吕爱娟, 沈加林, 等. 长江口启东-崇明岛航道沉积物中多环芳烃分布来源及生态风险评价[J]. 岩矿测试, 2014, 33(3): 374-380
Shen X M, Lv A J, Shen J L, et al. Distribution characteristics, sources and ecological risk assessment of polycyclic aromatic hydrocarbons in waterway sediments from Qidong and Chongming Island of Yangtze River Estuary [J]. Rock and Mineral Analysis, 2014, 33(3): 374-380 (in Chinese)
[15] 王波, 李正炎, 傅明珠, 等. 长江口及其邻近海域表层沉积物中多环芳烃的分布和生态风险评价[J]. 中国海洋大学学报: 自然科学版, 2007, 37(S1): 83-87
Wang B, Li Z Y, Fu M Z, et al. Distribution and ecological risk assessment of PAHs in surface sediments from the Yangtze Estuary and its adjacent areas [J]. Periodical of Ocean University of China, 2007, 37(S1): 83-87 (in Chinese)
[16] International Organization for Standardization (ISO). Water quality: Determination of the inhibitory effect of water samples on the light emission of Vibrio fischeri (luminescent bacteria test). Part 3: Method using freeze-dried bacteria [R]. Geneva: ISO, 2007
[17] Escher B I, Baumer A, Kai B, et al. General baseline toxicity QSAR for nonpolar, polar and ionisable chemicals and their mixtures in the bioluminescence inhibition assay with Aliivibrio fischeri [J]. Environmental Science: Processes and Impacts, 2017, 19: 414-428
[18] Tang J Y M, McCarty S, Glenn E, et al. Mixture effects of organic micropollutants present in water: Towards the development of effect-based water quality trigger values for baseline toxicity [J]. Water Research, 2013, 47(10): 3300-3314
[19] Escher B I, Allinson M, Altenburger R, et al. Benchmarking organic micropollutants in wastewater, recycled water and drinking water with in vitro bioassays [J]. Environmental Science & Technology, 2014, 48(3): 1940-1956
[20] Hashmi M A K, Escher B I, Krauss M, et al. Effect-directed analysis (EDA) of Danube River water sample receiving untreated municipal wastewater from Novi Sad, Serbia [J]. Science of the Total Environment, 2018, 624: 1072-1081
[21] Bräunig J, Tang J Y M, Warne M S J, et al. Bioanalytical effect-balance model to determine the bioavailability of organic contaminants in sediments affected by black and natural carbon [J]. Chemosphere, 2016, 156: 181-190
[22] Li J Y, Tang J Y M, Jin L, et al. Understanding bioavailability and toxicity of sediment-associated contaminants by combining passive sampling with in vitro bioassays in an urban river catchment [J]. Environmental Toxicology and Chemistry, 2013, 32(12): 2888-2896
[23] Neale P A, Altenburger R, Aït-Aïssa S, et al. Development of a bioanalytical test battery for water quality monitoring: Fingerprinting identified micropollutants and their contribution to effects in surface water [J]. Water Research, 2017, 123: 734-750
[24] Jacobs M W, Coates J A, Delfino J J, et al. Comparison of sediment extract microtox® toxicity with semi-volatile organic priority pollutant concentrations [J]. Archives of Environmental Contamination and Toxicology, 1993, 24(4): 461-468
[25] Neale P A, Munz N A, Aït-Aïssa S, et al. Integrating chemical analysis and bioanalysis to evaluate the contribution of wastewater effluent on the micropollutant burden in small streams [J]. The Science of the Total Environment, 2017, 576: 785-795
[26] 袁广旺, 郑江鹏, 花卫华, 等. 吕泗渔场表层沉积物中多环芳烃分布特征及生态风险评价[J]. 海洋开发与管理, 2015, 32(10): 78-83
[27] 万希鹏. 沈家门渔港海域环境质量评价的探讨[J]. 环境污染与防治, 1981, 3(2): 37-42
[28] 张玉凤, 吴金浩, 李楠, 等. 渤海北部表层沉积物中多环芳烃分布与来源分析[J]. 海洋环境科学, 2016, 35(1): 88-94, 122
Zhang Y F, Wu J H, Li N, et al. Distribution and source identification of polycyclic aromatic hydrocarbon of surface sediments from the North Bohai Sea [J]. Marine Environmental Science, 2016, 35(1): 88-94, 122 (in Chinese)
[29] 刘霞. 长江口及其邻近海域沉积物有机污染物分布特征及对有机质来源的指示意义[D]. 青岛: 中国海洋大学, 2012: 32-39
Liu X. Distribution of hydrophobic pollutants in sediment from Yangtze Estuary and adjacent area and itsimplication to sources of organic matter [D]. Qingdao: Ocean University of China, 2012: 32-39 (in Chinese)
[30] 张明, 唐访良, 吴志旭, 等. 千岛湖表层沉积物中多环芳烃污染特征及生态风险评价[J]. 中国环境科学, 2014, 34(1): 253-258
Zhang M, Tang F L, Wu Z X, et al. Pollution characteristics and ecological risk assessment of polycyclic aromatic hydrocarbons (PAHs) in surface sediments from Xin’anjiang Reservoir [J]. China Environmental Science, 2014, 34(1): 253-258 (in Chinese)
[31] Witt G. Polycyclic aromatic hydrocarbons in water and sediment of the Baltic Sea [J]. Marine Pollution Bulletin, 1995, 31(4-12): 237-248
[32] Long E R, Macdonald D D, Smith S L, et al. Incidence of adverse biological effects within ranges of chemical concentrations in marine and estuarine sediments [J]. Environmental Management, 1995, 19(1): 81-97
[33] Hwang K, Lee J, Kwon I, et al. Large-scale sediment toxicity assessment over the 15, 000 km of coastline in the Yellow and Bohai Seas, East Asia [J]. Science of the Total Environment, 2021, 792: 148371
[34] Jin L, Xie J W, Wong C K C, et al. Contributions of city-specific fine particulate matter (PM2.5) to differential in vitro oxidative stress and toxicity implications between Beijing and Guangzhou of China [J]. Environmental Science & Technology, 2019, 53(5): 2881-2891