沉积物中残留油类对拟除虫菊酯与多环芳烃影响摇蚊慢性毒性的作用机制探究

沉积物中残留油类对拟除虫菊酯与多环芳烃影响摇蚊慢性毒性的作用机制探究一、引言1.1研究背景与意义摇蚊，作为昆虫纲双翅目摇蚊科的一员，在各类水域生态系统中广泛分布，是淡水生态系统的重要组成部分。摇蚊的生命周期涵盖卵、幼虫、蛹及成虫四个阶段，其幼虫阶段在水底淤泥中生活，历时较长，主要以水底的有机物碎屑和藻类为食，尤其是蓝藻和绿藻，在摄食过程中，摇蚊幼虫能够有效分解和转化这些物质，促进水体中营养物质的循环和再利用，对维持水体生态平衡发挥着关键作用，是净化水质的重要参与者。同时，摇蚊幼虫富含蛋白质、脂肪等营养物质，是鱼类等水生生物的优质天然饵料，在水生生态系统食物链中占据重要位置，为众多水生生物提供了必要的食物来源，支撑着整个生态系统的能量流动和物质循环。拟除虫菊酯类农药作为一类模拟天然除虫菊素人工合成的杀虫剂，凭借其杀虫谱广、效果显著、低残留以及对人畜低毒等诸多优点，自20世纪80年代崛起以来，在农业生产领域得到了极为广泛的应用。不仅用于棉花、蔬菜、果树、茶叶等农作物害虫的防治，还在室内蚊虫消杀等方面发挥作用。然而，随着其使用范围和使用量的不断增加，拟除虫菊酯类农药对非靶标生物的潜在危害逐渐凸显。众多研究表明，这类农药对鱼类等水生生物具有较高的毒性，其疏水性使其容易被鱼鳃吸收，强亲脂性导致即使在水中浓度很低时也能被鱼鳞强烈吸附，加之鱼体内缺乏有效水解菊酯类物质的酶，使得拟除虫菊酯类农药在鱼体内的代谢缓慢，毒性作用远大于对哺乳动物和鸟类，可导致鱼类组织器官受损、神经系统和免疫系统功能障碍，影响其繁殖和生长。而摇蚊作为水生生态系统中的重要环节，同样面临着拟除虫菊酯类农药的威胁，其生长发育、生存和繁殖等方面可能受到不利影响，进而破坏水生生态系统的平衡。多环芳烃（PAHs）是一类由两个以上苯环以稠环形式相连的化合物，主要来源于煤、石油、木材、烟草等有机物的不完全燃烧以及地球化学过程。在现代工业活动中，如炼油厂、炼焦厂、橡胶厂等排放的废气、废水和废渣，以及各种交通车辆排放的尾气，都含有大量的多环芳烃，使其广泛分布于大气、土壤和水体等环境介质中。多环芳烃具有较强的致癌性、致畸性和致突变性，对人体的肝脏、肺、胃、皮肤等组织均有毒性作用，长期暴露于多环芳烃污染环境的人群，患癌症、呼吸系统疾病等的风险显著增加。在生态环境中，多环芳烃可通过食物链的生物富集作用进入生物体内，对生物造成毒害。由于摇蚊幼虫生活在沉积物中，而沉积物是多环芳烃的重要归宿之一，摇蚊极易受到多环芳烃的污染影响，其生理功能和生态行为可能发生改变，对整个水生生态系统的稳定性和生物多样性构成潜在威胁。在各类工业生产和石油运输、储存过程中，不可避免地会产生残留油类物质，这些油类物质进入水体后，最终会沉降到沉积物中。沉积物作为水体中污染物的重要蓄积库，其中残留油类的存在可能会改变其他污染物的环境行为和生态毒性。拟除虫菊酯类农药和多环芳烃在沉积物中广泛存在，而残留油类可能与它们发生相互作用，影响其在沉积物中的吸附-解吸、迁移转化等过程，进而改变它们对摇蚊的慢性毒性效应。然而，目前关于沉积物中残留油类对拟除虫菊酯类农药和多环芳烃对摇蚊慢性毒性影响的研究还相对较少，相关作用机制尚不清楚。本研究聚焦于沉积物中残留油类对拟除虫菊酯类农药和多环芳烃对摇蚊的慢性毒性影响，具有重要的现实意义和科学价值。通过深入探究这一问题，能够更全面、准确地评估拟除虫菊酯类农药和多环芳烃在复杂环境中的生态风险，为科学合理地使用拟除虫菊酯类农药提供理论依据，助力农业生产与生态环境保护的协调发展。同时，有助于揭示残留油类在污染物生态毒性中的作用机制，丰富环境科学领域关于复合污染的理论知识，为水生生态系统的保护和修复提供新的思路和方法，对于维护生态平衡、保障人类健康和生态系统的可持续发展具有重要的推动作用。1.2国内外研究现状在拟除虫菊酯类农药对摇蚊的毒性影响研究方面，国外学者较早开展了相关工作。[国外文献1]通过实验发现，某些拟除虫菊酯类农药会显著抑制摇蚊幼虫的生长，使其发育迟缓，并且导致摇蚊的羽化率降低，影响其种群的正常繁衍。[国外文献2]研究指出，拟除虫菊酯类农药会干扰摇蚊的神经系统，改变其行为模式，如使其运动能力下降、对环境刺激的反应变弱等。国内学者也对此进行了深入探究，[国内文献1]的研究表明，不同种类的拟除虫菊酯类农药对摇蚊的毒性存在差异，其中溴氰菊酯的毒性相对较高，低浓度的溴氰菊酯就能对摇蚊幼虫的生存和生长产生明显的抑制作用。[国内文献2]从分子生物学角度研究发现，拟除虫菊酯类农药会诱导摇蚊体内某些基因的表达发生变化，影响其正常的生理代谢过程。关于多环芳烃对摇蚊的毒性研究，国外有研究[国外文献3]表明，多环芳烃会在摇蚊体内富集，随着暴露时间的延长和浓度的增加，摇蚊的死亡率逐渐上升，并且会对摇蚊的抗氧化系统造成损害，导致其体内抗氧化酶活性发生改变。[国外文献4]通过野外调查和实验室模拟相结合的方法，发现多环芳烃污染严重的水域中，摇蚊的种类和数量明显减少，群落结构发生改变。国内方面，[国内文献3]研究了不同环数的多环芳烃对摇蚊的毒性效应，结果显示，高环数的多环芳烃如苯并[a]芘对摇蚊的毒性更强，不仅会影响摇蚊的生长发育，还会导致其生殖能力下降。[国内文献4]利用转录组学技术分析了多环芳烃胁迫下摇蚊基因表达的变化，发现多环芳烃会影响摇蚊多个生物学过程相关基因的表达，包括能量代谢、物质转运等。然而，目前对于沉积物中残留油类对拟除虫菊酯类农药和多环芳烃对摇蚊慢性毒性影响的研究相对较少。少量的研究[相关文献]主要集中在单一污染物对摇蚊的毒性研究，对于多种污染物之间的相互作用以及残留油类在其中所起的作用研究不够深入。在研究方法上，大多采用实验室模拟实验，与实际环境条件存在一定差距，难以全面准确地反映自然环境中污染物的复合污染情况。而且，对于污染物对摇蚊毒性作用的分子机制研究还不够透彻，缺乏从基因、蛋白质等层面深入探究残留油类影响拟除虫菊酯类农药和多环芳烃对摇蚊慢性毒性的内在机制。综上所述，当前关于拟除虫菊酯类农药和多环芳烃对摇蚊的毒性研究已取得一定成果，但在沉积物中残留油类的影响方面仍存在诸多不足。本研究将针对这些问题展开深入探究，通过更贴近实际环境的实验设计，从多个层面分析残留油类对拟除虫菊酯类农药和多环芳烃对摇蚊慢性毒性的影响，以期填补相关研究空白，为水生生态系统的保护和污染治理提供更全面、准确的理论依据。二、材料与方法2.1实验材料2.1.1沉积物采集本研究分别从城市河流（如[河流名称1]，位于城市中心区域，周边有工业企业和居民生活区，受工业废水排放和生活污水污染影响较大）、湖泊（如[湖泊名称1]，处于城市郊区，主要污染源为农业面源污染和旅游活动带来的污染）、港口（如[港口名称1]，作为重要的货物运输枢纽，船舶燃油泄漏和货物装卸过程中的污染较为严重）等不同环境中采集沉积物样品。使用柱状采泥器进行采样，确保采集到的沉积物具有代表性且能反映不同深度层面的情况。在采样时，将柱状采泥器垂直插入水底，缓慢取出，避免沉积物的扰动和流失。每个采样点采集3-5个平行样品，以减少采样误差。采集后的沉积物样品运回实验室后，首先进行初步的物理分离，去除其中的砾石、贝壳、动植物残体等杂物，然后用63μm的尼龙筛网进行湿法过筛，以分出岩石的碎片等大块物质，使小于63μm的组分保留，该部分被认为能较好地代表微量元素和污染物的分布。对于形态分析用的沉积物样品，放置于惰性气体保护的胶皮套箱中，以避免氧化。将过筛后的沉积物样品分层装入聚乙烯瓶中，在4℃下冷藏保存，尽快进行后续实验分析，防止样品中污染物的形态和含量发生变化。2.1.2实验生物实验所用摇蚊幼虫采自[无污染水体名称]，该水体水质清洁，未受明显污染，摇蚊种群健康，以保证实验生物的质量和一致性。采集时，使用手抄网在水体底部的淤泥中轻轻捞取摇蚊幼虫，尽量避免损伤幼虫。将采集到的摇蚊幼虫带回实验室后，放入盛有曝气自来水的玻璃缸中进行暂养驯化。养殖缸中放置适量的水生植物，如金鱼藻，为摇蚊幼虫提供栖息和食物来源。养殖用水的温度控制在（20±1）℃，光照周期设置为16L:8D（光照16小时，黑暗8小时），每天投喂适量的酵母粉和藻类混合饲料，保证摇蚊幼虫的正常生长和发育。选取活力强、个体大小相近、无明显损伤和病变的4龄摇蚊幼虫作为实验对象。在实验前，对摇蚊幼虫进行饥饿处理24小时，使其肠道排空，减少体内原有物质对实验结果的干扰，以确保实验结果的准确性和可靠性。2.1.3化学试剂拟除虫菊酯类农药选择了常用的溴氰菊酯、氯氰菊酯和高效氯氟氰菊酯，纯度均≥98%，购自[试剂供应商1]。多环芳烃选用了美国环保局（EPA）优先控制的16种多环芳烃中的萘、菲、芘、苯并[a]芘等，纯度≥97%，由[试剂供应商2]提供。残留油类试剂采用实际采集的原油样品，经过预处理后去除杂质，得到纯净的原油，其主要成分为烷烃、环烷烃和芳烃等。所有化学试剂均妥善保存，避免阳光直射和高温环境，防止其分解和变质。在使用前，对试剂的纯度和浓度进行再次检测，确保实验结果的准确性。2.2实验设计2.2.1单一污染物毒性实验对于拟除虫菊酯类农药的慢性毒性实验，选择溴氰菊酯、氯氰菊酯和高效氯氟氰菊酯作为代表。根据前期预实验和相关文献资料，确定溴氰菊酯的实验浓度梯度为0.01μg/L、0.1μg/L、1μg/L、10μg/L、100μg/L；氯氰菊酯的浓度梯度设置为0.05μg/L、0.5μg/L、5μg/L、50μg/L、500μg/L；高效氯氟氰菊酯的浓度梯度为0.02μg/L、0.2μg/L、2μg/L、20μg/L、200μg/L。每个浓度设置3个平行实验组，同时设置空白对照组，对照组中不添加拟除虫菊酯类农药，仅含有与实验组相同的实验介质和养殖条件。实验周期为28天，在实验过程中，每天观察并记录摇蚊幼虫的生长发育情况，包括体长、体重、蜕皮次数等，每7天测定一次摇蚊幼虫的死亡率。多环芳烃的慢性毒性实验选取萘、菲、芘、苯并[a]芘作为研究对象。参考相关研究及预实验结果，萘的实验浓度梯度设定为1mg/L、5mg/L、10mg/L、50mg/L、100mg/L；菲的浓度梯度为0.5mg/L、2.5mg/L、5mg/L、25mg/L、50mg/L；芘的浓度梯度为0.1mg/L、0.5mg/L、1mg/L、5mg/L、10mg/L；苯并[a]芘的浓度梯度为0.01mg/L、0.05mg/L、0.1mg/L、0.5mg/L、1mg/L。同样每个浓度设置3个平行实验组和1个空白对照组，实验周期为28天。在实验期间，密切观察摇蚊幼虫的行为变化，如运动能力、摄食情况等，每隔7天测量摇蚊幼虫体内多环芳烃的富集量，分析多环芳烃在摇蚊体内的积累规律及其对摇蚊生理机能的影响。2.2.2联合污染物毒性实验拟除虫菊酯类农药和多环芳烃的联合毒性实验设置不同的组合。选取溴氰菊酯和苯并[a]芘作为代表，设置溴氰菊酯浓度为0.1μg/L、1μg/L、10μg/L，苯并[a]芘浓度为0.01mg/L、0.1mg/L、1mg/L，进行交叉组合，共形成9个实验组，每个实验组设置3个平行。同时设置只含单一污染物（即单独的溴氰菊酯或苯并[a]芘，浓度与联合实验组中单一污染物浓度相同）的对照组以及空白对照组，以对比分析联合作用与单一作用的差异。实验周期为28天，在实验过程中，定期监测摇蚊幼虫的生长发育指标、死亡率、抗氧化酶活性、代谢酶活性等，探究拟除虫菊酯类农药和多环芳烃联合作用对摇蚊的毒性效应及机制。在联合污染物毒性实验中添加不同含量残留油类的实验，选取原油作为残留油类代表。在上述拟除虫菊酯类农药（以溴氰菊酯为例，浓度为1μg/L）和多环芳烃（以苯并[a]芘为例，浓度为0.1mg/L）联合实验组的基础上，分别添加质量分数为0.1%、1%、5%的残留油类，每个添加量设置3个平行实验组。同时设置不添加残留油类的联合污染物对照组以及空白对照组。实验周期为28天，在实验期间，除监测摇蚊幼虫的生长发育、生存等指标外，还分析残留油类对拟除虫菊酯类农药和多环芳烃在沉积物中的吸附-解吸、迁移转化过程的影响，以及对摇蚊体内污染物富集量和毒性效应的影响，深入探讨残留油类在复合污染体系中的作用机制。2.3分析方法2.3.1沉积物成分分析对于沉积物中残留油类的含量测定，采用重量法。具体步骤为：取一定量的沉积物样品，放入索氏提取器中，用正己烷作为提取剂，在80℃的水浴温度下回流提取8小时，使残留油类充分溶解于正己烷中。提取液经旋转蒸发仪浓缩后，转移至已恒重的称量瓶中，在60℃的烘箱中烘干至恒重，通过称量前后的质量差计算残留油类的含量。拟除虫菊酯类农药的分析采用气相色谱-电子捕获检测器（GC-ECD）联用技术。首先将沉积物样品与适量的无水硫酸钠混合，研磨均匀后，加入丙酮-正己烷（体积比为1:1）混合溶剂，使用超声波辅助提取30分钟，提取液经离心分离后，取上清液通过硅胶柱进行净化处理，去除杂质。净化后的样品用氮气吹干，再用正己烷定容至1mL，注入气相色谱仪进行分析。气相色谱条件为：采用DB-5毛细管柱（30m×0.25mm×0.25μm），初始柱温为50℃，保持1分钟，以20℃/min的速率升温至250℃，保持10分钟；进样口温度为280℃，检测器温度为300℃；载气为氮气，流速为1mL/min；进样量为1μL。通过与标准品的保留时间对比进行定性分析，外标法进行定量分析。多环芳烃的测定运用气相色谱-质谱联用仪（GC-MS）。沉积物样品经冷冻干燥后，加入适量的二氯甲烷，采用加速溶剂萃取仪在100℃、1500psi的条件下萃取10分钟，萃取液经硅胶-弗罗里硅土混合柱净化，用正己烷-二氯甲烷（体积比为3:1）混合溶剂洗脱，洗脱液浓缩后进行GC-MS分析。GC条件：HP-5MS毛细管柱（30m×0.25mm×0.25μm），初始柱温为60℃，保持2分钟，以15℃/min的速率升温至300℃，保持5分钟；进样口温度为280℃；载气为氦气，流速为1mL/min；进样量为1μL，不分流进样。MS条件：离子源为电子轰击源（EI），离子源温度为230℃，扫描方式为选择离子扫描（SIM），根据各多环芳烃的特征离子进行定性和定量分析。2.3.2摇蚊毒性指标测定在实验过程中，每隔7天对摇蚊幼虫的生长发育指标进行监测。使用电子天平（精度为0.0001g）称量摇蚊幼虫的体重，用游标卡尺（精度为0.02mm）测量摇蚊幼虫的体长。每天观察并记录摇蚊幼虫的生存率，若摇蚊幼虫在5分钟内对机械刺激无反应，则判定为死亡。定期观察摇蚊幼虫的形态变化，包括体色、附肢完整性、身体畸形等情况，使用体视显微镜（放大倍数为10-40倍）进行详细观察并拍照记录。在实验第14天和第28天，分别采集摇蚊幼虫样本，测定其生化指标。将摇蚊幼虫用生理盐水冲洗干净后，放入预冷的匀浆器中，加入适量的磷酸盐缓冲液（pH=7.4），在冰浴条件下匀浆，匀浆液在4℃、10000r/min的条件下离心15分钟，取上清液用于生化指标测定。采用试剂盒法测定超氧化物歧化酶（SOD）、过氧化氢酶（CAT）、谷胱甘肽-S-转移酶（GST）等抗氧化酶的活性，以及乙酰胆碱酯酶（AChE）等代谢酶的活性，具体操作按照试剂盒说明书进行。2.4数据处理实验数据采用SPSS22.0统计软件进行分析处理。对于摇蚊幼虫的生长发育指标（如体长、体重）、生存率、生化指标（抗氧化酶活性、代谢酶活性等）数据，首先进行正态性检验和方差齐性检验，若数据满足正态分布和方差齐性，采用单因素方差分析（One-WayANOVA）比较不同实验组之间的差异，当差异显著（P<0.05）时，进一步采用Duncan多重比较法进行组间两两比较，确定各实验组之间的具体差异情况。若数据不满足正态分布或方差齐性，则采用非参数检验（如Kruskal-Wallis秩和检验）进行分析。对于沉积物中残留油类、拟除虫菊酯类农药和多环芳烃的含量数据，同样先进行正态性和方差齐性检验，符合条件的数据使用方差分析比较不同采样点或不同实验条件下的含量差异；不符合条件的数据采用相应的非参数检验方法。在分析拟除虫菊酯类农药和多环芳烃对摇蚊毒性效应的相关性时，运用Pearson相关性分析方法，计算不同污染物浓度与摇蚊各项毒性指标之间的相关系数，确定它们之间的线性相关关系。对于联合污染物毒性实验中残留油类对拟除虫菊酯类农药和多环芳烃毒性的影响，通过建立多元线性回归模型，分析残留油类含量、拟除虫菊酯类农药浓度、多环芳烃浓度等因素对摇蚊毒性指标的综合影响，探讨其中的作用机制。所有实验数据均以平均值±标准差（Mean±SD）表示，以P<0.05作为差异具有统计学意义的标准，确保实验结果的可靠性和准确性，为深入研究沉积物中残留油类对拟除虫菊酯类农药和多环芳烃对摇蚊的慢性毒性影响提供有力的数据支持。三、结果与讨论3.1沉积物中污染物分析结果对不同采样点沉积物中残留油类、拟除虫菊酯类农药和多环芳烃的含量及种类进行分析，结果如表1所示。在城市河流沉积物中，残留油类含量范围为0.25%-0.68%，平均值为（0.45±0.13）%，主要来源于周边工业企业废水排放和生活污水中的石油类物质。湖泊沉积物中残留油类含量相对较低，为0.12%-0.35%，平均值为（0.23±0.08）%，可能是由于湖泊水体流动性较弱，污染物扩散稀释能力较强，且受农业面源污染和旅游活动影响，石油类污染物输入相对较少。港口沉积物中残留油类含量最高，达到0.56%-1.23%，平均值为（0.85±0.21）%，这主要是因为港口船舶燃油泄漏和货物装卸过程中大量油类物质进入水体并沉降到沉积物中。采样点残留油类含量（%）拟除虫菊酯类农药含量（μg/kg）多环芳烃含量（μg/kg）城市河流0.25-0.68（0.45±0.13）溴氰菊酯：1.2-5.6（3.4±1.1）氯氰菊酯：2.5-8.7（5.6±1.8）高效氯氟氰菊酯：1.8-6.3（4.2±1.3）萘：15.6-32.4（23.5±4.8）菲：8.5-18.6（13.5±3.2）芘：5.6-12.3（8.9±2.1）苯并[a]芘：0.5-2.3（1.4±0.5）湖泊0.12-0.35（0.23±0.08）溴氰菊酯：0.8-3.5（2.1±0.9）氯氰菊酯：1.5-5.6（3.6±1.2）高效氯氟氰菊酯：1.2-4.5（2.8±1.0）萘：10.2-22.5（16.3±3.5）菲：5.6-12.3（8.9±2.5）芘：3.5-8.6（6.0±1.8）苯并[a]芘：0.3-1.5（0.9±0.4）港口0.56-1.23（0.85±0.21）溴氰菊酯：3.5-10.2（6.8±2.0）氯氰菊酯：5.6-15.3（10.5±3.1）高效氯氟氰菊酯：4.2-12.5（8.3±2.5）萘：25.6-56.8（38.7±7.5）菲：15.6-32.4（24.0±5.6）芘：8.9-20.5（14.7±3.8）苯并[a]芘：1.2-4.5（2.8±0.8）在拟除虫菊酯类农药方面，城市河流沉积物中检测出溴氰菊酯、氯氰菊酯和高效氯氟氰菊酯，含量分别为1.2-5.6μg/kg（平均值为3.4±1.1μg/kg）、2.5-8.7μg/kg（平均值为5.6±1.8μg/kg）和1.8-6.3μg/kg（平均值为4.2±1.3μg/kg）。这些农药主要来源于周边农田的农业生产活动，通过地表径流和雨水冲刷进入河流，最终在沉积物中积累。湖泊沉积物中拟除虫菊酯类农药含量相对较低，溴氰菊酯为0.8-3.5μg/kg（平均值为2.1±0.9μg/kg），氯氰菊酯为1.5-5.6μg/kg（平均值为3.6±1.2μg/kg），高效氯氟氰菊酯为1.2-4.5μg/kg（平均值为2.8±1.0μg/kg），这可能与湖泊周边农业活动相对较少以及水体的自净作用有关。港口沉积物中拟除虫菊酯类农药含量较高，溴氰菊酯为3.5-10.2μg/kg（平均值为6.8±2.0μg/kg），氯氰菊酯为5.6-15.3μg/kg（平均值为10.5±3.1μg/kg），高效氯氟氰菊酯为4.2-12.5μg/kg（平均值为8.3±2.5μg/kg），除了受周边农业面源污染影响外，港口区域可能存在一些非法倾倒农药废弃物等情况，导致农药在沉积物中的浓度升高。多环芳烃在不同采样点沉积物中也均有检出。城市河流沉积物中萘含量为15.6-32.4μg/kg（平均值为23.5±4.8μg/kg），菲为8.5-18.6μg/kg（平均值为13.5±3.2μg/kg），芘为5.6-12.3μg/kg（平均值为8.9±2.1μg/kg），苯并[a]芘为0.5-2.3μg/kg（平均值为1.4±0.5μg/kg），其来源主要是城市中各种化石燃料的不完全燃烧，如工业废气排放、机动车尾气排放以及生活燃煤等。湖泊沉积物中多环芳烃含量相对较低，萘为10.2-22.5μg/kg（平均值为16.3±3.5μg/kg），菲为5.6-12.3μg/kg（平均值为8.9±2.5μg/kg），芘为3.5-8.6μg/kg（平均值为6.0±1.8μg/kg），苯并[a]芘为0.3-1.5μg/kg（平均值为0.9±0.4μg/kg），这是由于湖泊周边相对清洁，污染源较少，且水体的稀释和扩散作用使得多环芳烃在沉积物中的积累量减少。港口沉积物中多环芳烃含量最高，萘为25.6-56.8μg/kg（平均值为38.7±7.5μg/kg），菲为15.6-32.4μg/kg（平均值为24.0±5.6μg/kg），芘为8.9-20.5μg/kg（平均值为14.7±3.8μg/kg），苯并[a]芘为1.2-4.5μg/kg（平均值为2.8±0.8μg/kg），除了受到城市污染的影响外，港口船舶发动机燃油的燃烧以及石油类物质的泄漏，也会导致多环芳烃在沉积物中的含量增加。总体而言，不同采样点沉积物中残留油类、拟除虫菊酯类农药和多环芳烃的含量存在明显差异，这与各采样点的环境特征和污染源密切相关。港口沉积物由于受到船舶活动、工业污染和农业面源污染等多种因素的综合影响，污染物含量普遍较高；城市河流沉积物次之，主要受工业废水排放、生活污水和农业生产活动影响；湖泊沉积物相对较为清洁，污染物含量较低。这些污染物在沉积物中的积累，可能会对栖息其中的摇蚊等生物产生潜在的毒性影响，进而影响整个水生生态系统的结构和功能。3.2单一污染物对摇蚊的慢性毒性影响3.2.1拟除虫菊酯类农药拟除虫菊酯类农药对摇蚊幼虫生长发育、生存率和形态的影响实验结果如表2所示。在不同浓度的溴氰菊酯暴露下，摇蚊幼虫的体长和体重增长均受到显著抑制（P<0.05）。当溴氰菊酯浓度为0.01μg/L时，摇蚊幼虫第28天的体长为（4.23±0.35）mm，体重为（0.85±0.12）mg，与对照组相比，体长和体重分别下降了15.4%和18.3%；随着浓度升高至100μg/L，体长仅为（2.15±0.21）mm，体重为（0.32±0.05）mg，与对照组相比，体长和体重分别下降了57.2%和69.1%，呈现出明显的剂量-效应关系。氯氰菊酯和高效氯氟氰菊酯也表现出类似的趋势，随着浓度的增加，摇蚊幼虫的生长发育受到的抑制作用逐渐增强。拟除虫菊酯类农药浓度（μg/L）体长（mm，第28天）体重（mg，第28天）生存率（%，第28天）形态变化溴氰菊酯0（对照）5.00±0.421.04±0.1595.0±3.2正常，体色均匀，附肢完整0.014.23±0.350.85±0.1290.0±4.5部分个体体色变浅0.13.56±0.300.68±0.1080.0±5.6附肢出现轻微畸形，运动能力下降12.89±0.250.52±0.0865.0±6.8身体弯曲，出现蜕皮困难现象102.56±0.220.40±0.0645.0±7.5体色暗淡，附肢残缺，死亡率增加1002.15±0.210.32±0.0520.0±5.0身体严重畸形，大量死亡氯氰菊酯0（对照）5.00±0.421.04±0.1595.0±3.2正常，体色均匀，附肢完整0.054.15±0.330.82±0.1188.0±4.0体色稍有变化0.53.45±0.280.65±0.0975.0±5.0附肢出现畸形，运动迟缓52.78±0.230.48±0.0755.0±6.5身体弯曲，蜕皮异常502.34±0.200.36±0.0535.0±6.0体色灰暗，附肢断裂，死亡率上升5001.98±0.180.28±0.0415.0±4.0严重畸形，大量个体死亡高效氯氟氰菊酯0（对照）5.00±0.421.04±0.1595.0±3.2正常，体色均匀，附肢完整0.024.20±0.340.83±0.1189.0±4.2体色略变浅0.23.50±0.290.66±0.0978.0±5.2附肢畸形，活动减少22.82±0.240.50±0.0860.0±6.0身体弯曲，蜕皮困难202.45±0.210.38±0.0640.0±7.0体色暗淡，附肢损伤，死亡率升高2002.05±0.190.30±0.0425.0±5.5身体严重畸形，多数死亡在生存率方面，随着拟除虫菊酯类农药浓度的增加，摇蚊幼虫的生存率显著降低（P<0.05）。溴氰菊酯浓度为10μg/L时，摇蚊幼虫的生存率降至45.0%±7.5%，当浓度达到100μg/L时，生存率仅为20.0%±5.0%。氯氰菊酯在浓度为50μg/L时，生存率为35.0%±6.0%，500μg/L时降至15.0%±4.0%；高效氯氟氰菊酯在20μg/L时，生存率为40.0%±7.0%，200μg/L时为25.0%±5.5%。这表明拟除虫菊酯类农药对摇蚊幼虫具有较强的致死作用，且浓度越高，致死效应越明显。在形态方面，低浓度的拟除虫菊酯类农药会使摇蚊幼虫体色变浅，部分个体附肢出现轻微畸形，运动能力下降；随着浓度的升高，摇蚊幼虫出现身体弯曲、蜕皮困难、附肢残缺等现象，体色暗淡，严重时身体严重畸形，大量死亡。在氯氰菊酯浓度为5μg/L的实验组中，观察到部分摇蚊幼虫附肢出现畸形，运动迟缓，身体弯曲，蜕皮异常；在高效氯氟氰菊酯浓度为20μg/L的实验组中，摇蚊幼虫体色暗淡，附肢损伤，死亡率升高。这些形态变化可能会影响摇蚊幼虫的正常生活和生存能力，进而影响其种群数量和分布。拟除虫菊酯类农药对摇蚊的毒性机制主要与其对神经系统的干扰有关。拟除虫菊酯类农药能够作用于摇蚊的神经细胞膜上的钠离子通道，使钠离子通道持续开放，导致神经细胞去极化，从而引起神经系统的兴奋和痉挛。当农药浓度较低时，摇蚊的神经系统受到一定程度的干扰，表现为运动能力下降、行为异常等；随着浓度的增加，神经系统受到的损害加剧，导致摇蚊的生长发育受阻，生存率降低，形态出现畸形。拟除虫菊酯类农药还可能影响摇蚊体内的抗氧化酶系统和代谢酶系统，导致氧化应激和代谢紊乱，进一步加重对摇蚊的毒性作用。在溴氰菊酯处理组中，摇蚊幼虫体内的超氧化物歧化酶（SOD）和过氧化氢酶（CAT）活性随着农药浓度的增加而显著降低，表明拟除虫菊酯类农药导致摇蚊体内氧化应激水平升高，抗氧化能力下降。3.2.2多环芳烃多环芳烃对摇蚊幼虫的慢性毒性实验结果如表3所示。不同种类的多环芳烃对摇蚊幼虫的生长发育、生存率和形态产生了不同程度的影响。萘对摇蚊幼虫生长发育的抑制作用相对较弱，在浓度为100mg/L时，第28天摇蚊幼虫的体长为（4.56±0.38）mm，体重为（0.92±0.13）mg，与对照组相比，体长和体重分别下降了8.8%和11.5%。而苯并[a]芘的毒性较强，当浓度为1mg/L时，摇蚊幼虫第28天的体长为（2.98±0.26）mm，体重为（0.45±0.07）mg，与对照组相比，体长和体重分别下降了40.4%和56.7%，呈现出明显的剂量-效应关系。菲和芘对摇蚊幼虫生长发育的抑制作用介于萘和苯并[a]芘之间。多环芳烃浓度（mg/L）体长（mm，第28天）体重（mg，第28天）生存率（%，第28天）形态变化萘0（对照）5.00±0.421.04±0.1595.0±3.2正常，体色均匀，附肢完整14.85±0.400.98±0.1492.0±4.0无明显变化54.70±0.360.95±0.1388.0±4.5体色稍变深104.62±0.350.93±0.1385.0±5.0附肢稍有异常504.50±0.340.90±0.1280.0±5.5运动能力略有下降1004.56±0.380.92±0.1375.0±6.0无明显严重变化菲0（对照）5.00±0.421.04±0.1595.0±3.2正常，体色均匀，附肢完整0.54.65±0.370.90±0.1286.0±4.3体色变深2.54.20±0.320.80±0.1175.0±5.2附肢出现畸形53.85±0.300.70±0.1065.0±6.0身体弯曲，运动困难253.40±0.280.60±0.0950.0±6.5体色暗淡，蜕皮异常503.05±0.250.50±0.0835.0±6.0身体严重畸形，死亡率上升芘0（对照）5.00±0.421.04±0.1595.0±3.2正常，体色均匀，附肢完整0.14.70±0.360.95±0.1388.0±4.5体色加深0.54.30±0.330.85±0.1178.0±5.0附肢畸形，运动减少13.90±0.310.75±0.1068.0±5.5身体弯曲，出现蜕皮困难53.10±0.260.55±0.0845.0±6.0体色灰暗，附肢残缺102.65±0.230.45±0.0725.0±5.0身体严重畸形，大量死亡苯并[a]芘0（对照）5.00±0.421.04±0.1595.0±3.2正常，体色均匀，附肢完整0.014.50±0.340.90±0.1285.0±4.8体色变深0.054.00±0.300.75±0.1070.0±5.5附肢畸形，运动迟缓0.13.50±0.290.60±0.0955.0±6.0身体弯曲，蜕皮异常0.53.00±0.250.50±0.0835.0±6.0体色暗淡，附肢断裂12.98±0.260.45±0.0720.0±5.0身体严重畸形，大量死亡在生存率方面，随着多环芳烃浓度的增加，摇蚊幼虫的生存率逐渐降低。苯并[a]芘对摇蚊幼虫生存率的影响最为显著，当浓度为1mg/L时，生存率仅为20.0%±5.0%。芘在浓度为10mg/L时，生存率降至25.0%±5.0%；菲在50mg/L时，生存率为35.0%±6.0%；萘在100mg/L时，生存率为75.0±6.0%。这表明苯并[a]芘和芘对摇蚊幼虫具有较高的致死性，而萘的毒性相对较弱。在形态方面，低浓度的多环芳烃会使摇蚊幼虫体色变深，随着浓度的升高，摇蚊幼虫出现附肢畸形、身体弯曲、蜕皮异常、体色暗淡、附肢残缺等现象，严重时身体严重畸形，大量死亡。在菲浓度为2.5mg/L的实验组中，摇蚊幼虫附肢出现畸形，运动困难；在芘浓度为5mg/L的实验组中，摇蚊幼虫体色灰暗，附肢残缺。这些形态变化反映了多环芳烃对摇蚊幼虫生理结构和功能的破坏，影响其正常的生长和生存。不同种类多环芳烃的毒性差异主要与其分子结构和理化性质有关。苯并[a]芘和芘具有较高的分子量和较大的共轭体系，其疏水性较强，更容易在生物体内富集，且能够与生物体内的大分子物质如DNA、蛋白质等发生相互作用，从而产生更强的毒性效应。萘的分子量较小，结构相对简单，疏水性较弱，在生物体内的富集能力和与生物大分子的结合能力相对较弱，因此毒性相对较低。多环芳烃的毒性还与其代谢途径和代谢产物有关。一些多环芳烃在生物体内代谢过程中会产生具有更强毒性的代谢产物，如苯并[a]芘在细胞色素P450酶的作用下会代谢生成具有强致癌性的7,8-二醇-9,10-环氧化物，进一步加剧对生物体的危害。3.3联合污染物对摇蚊的慢性毒性影响3.3.1拟除虫菊酯类农药和多环芳烃联合拟除虫菊酯类农药和多环芳烃联合对摇蚊幼虫生长发育、生存率和形态的影响结果如表4所示。当溴氰菊酯浓度为1μg/L、苯并[a]芘浓度为0.1mg/L时，摇蚊幼虫第28天的体长为（2.35±0.20）mm，体重为（0.38±0.06）mg，与对照组相比，体长下降了53.0%，体重下降了63.5%；而单独使用1μg/L溴氰菊酯时，体长为（2.89±0.25）mm，体重为（0.52±0.08）mg，单独使用0.1mg/L苯并[a]芘时，体长为（3.50±0.29）mm，体重为（0.60±0.09）mg。这表明两者联合作用对摇蚊幼虫生长发育的抑制作用明显大于单独作用之和，呈现出协同效应。溴氰菊酯浓度（μg/L）苯并[a]芘浓度（mg/L）体长（mm，第28天）体重（mg，第28天）生存率（%，第28天）形态变化0（对照）05.00±0.421.04±0.1595.0±3.2正常，体色均匀，附肢完整0.10.013.25±0.280.62±0.0970.0±5.5附肢出现畸形，运动迟缓0.10.12.68±0.230.45±0.0750.0±6.0身体弯曲，蜕皮异常，体色暗淡0.112.10±0.180.30±0.0425.0±5.0身体严重畸形，大量死亡10.012.65±0.230.48±0.0760.0±6.0附肢畸形，身体弯曲10.12.35±0.200.38±0.0640.0±7.0身体弯曲，蜕皮困难，附肢残缺111.85±0.160.25±0.0315.0±4.0身体严重畸形，多数死亡100.012.15±0.190.35±0.0545.0±7.0体色暗淡，附肢损伤100.11.80±0.150.28±0.0430.0±6.0身体严重畸形，死亡率升高1011.50±0.130.20±0.0310.0±3.0极度畸形，几乎全部死亡在生存率方面，联合作用下摇蚊幼虫的生存率也显著低于单独作用时。当溴氰菊酯浓度为10μg/L、苯并[a]芘浓度为1mg/L时，摇蚊幼虫的生存率仅为10.0%±3.0%，而单独使用10μg/L溴氰菊酯时，生存率为45.0%±7.5%，单独使用1mg/L苯并[a]芘时，生存率为20.0%±5.0%。这进一步证明了拟除虫菊酯类农药和多环芳烃联合对摇蚊具有更强的致死作用，两者之间存在协同毒性效应。在形态方面，联合作用下摇蚊幼虫的畸形率和死亡率明显增加，出现身体严重畸形、附肢残缺、体色暗淡等现象，且随着浓度的增加，这些症状愈发严重。在溴氰菊酯浓度为1μg/L、苯并[a]芘浓度为0.1mg/L的实验组中，摇蚊幼虫身体弯曲，蜕皮困难，附肢残缺；在溴氰菊酯浓度为10μg/L、苯并[a]芘浓度为1mg/L的实验组中，摇蚊幼虫极度畸形，几乎全部死亡。拟除虫菊酯类农药和多环芳烃联合对摇蚊产生协同毒性效应的可能机制如下：拟除虫菊酯类农药主要作用于摇蚊的神经系统，而多环芳烃进入摇蚊体内后，会干扰其内分泌系统和免疫系统。多环芳烃中的一些化合物，如苯并[a]芘，能够与摇蚊体内的芳烃受体（AhR）结合，激活一系列的信号通路，导致细胞色素P450酶系的诱导表达。细胞色素P450酶系在多环芳烃的代谢过程中起着关键作用，但其活性的改变也会影响拟除虫菊酯类农药的代谢。当拟除虫菊酯类农药和多环芳烃同时存在时，多环芳烃诱导的细胞色素P450酶系可能会加速拟除虫菊酯类农药的代谢，产生更多的毒性代谢产物，或者改变拟除虫菊酯类农药在摇蚊体内的分布和蓄积，从而增强其对摇蚊的毒性。拟除虫菊酯类农药和多环芳烃可能会共同作用于摇蚊的细胞膜，破坏细胞膜的结构和功能，导致细胞通透性增加，使更多的有害物质进入细胞内，进一步加重对摇蚊的损伤。3.3.2残留油类的影响在添加不同含量残留油类的实验中，结果表明残留油类对拟除虫菊酯类农药和多环芳烃联合毒性对摇蚊的影响较为复杂，其作用机制与残留油类的含量密切相关。当残留油类含量为0.1%时，对拟除虫菊酯类农药和多环芳烃联合毒性对摇蚊的影响较小。在溴氰菊酯浓度为1μg/L、苯并[a]芘浓度为0.1mg/L且添加0.1%残留油类的实验组中，摇蚊幼虫第28天的体长为（2.30±0.20）mm，体重为（0.36±0.06）mg，生存率为38.0%±7.0%，与未添加残留油类的联合实验组相比，体长、体重和生存率虽有变化，但差异不显著（P>0.05）。这可能是因为低含量的残留油类在沉积物中分散较为均匀，对拟除虫菊酯类农药和多环芳烃的吸附和迁移影响较小，尚未明显改变它们对摇蚊的毒性效应。当残留油类含量增加到1%时，对联合毒性产生了一定的减缓作用。在相同的溴氰菊酯和苯并[a]芘浓度下，添加1%残留油类的实验组中，摇蚊幼虫第28天的体长为（2.55±0.22）mm，体重为（0.42±0.07）mg，生存率为45.0%±7.5%，与未添加残留油类的联合实验组相比，体长和体重有所增加，生存率也有所提高，差异具有统计学意义（P<0.05）。这可能是由于残留油类中的一些成分，如长链烷烃和环烷烃，能够与拟除虫菊酯类农药和多环芳烃发生相互作用，通过物理吸附或化学络合的方式，降低了它们在沉积物孔隙水中的浓度，减少了摇蚊对这些污染物的接触和吸收，从而减轻了联合毒性对摇蚊的危害。残留油类可能会改变沉积物的表面性质，增加沉积物对污染物的吸附能力，使污染物更多地被固定在沉积物中，难以释放到水体中被摇蚊摄取。然而，当残留油类含量进一步增加到5%时，却表现出对联合毒性的增强作用。在添加5%残留油类的实验组中，摇蚊幼虫第28天的体长为（2.05±0.18）mm，体重为（0.30±0.05）mg，生存率为30.0%±6.0%，与未添加残留油类的联合实验组相比，体长和体重显著降低，生存率也明显下降，差异具有高度统计学意义（P<0.01）。这可能是因为高含量的残留油类在沉积物中形成了较大的油滴或油膜，一方面，这些油滴或油膜会包裹拟除虫菊酯类农药和多环芳烃，阻碍它们在沉积物中的扩散和降解，使污染物在沉积物中的浓度相对升高。另一方面，摇蚊在沉积物中活动时，可能会接触到更多的油类物质，导致其体表和呼吸系统受到污染，影响其正常的生理功能，增加了对污染物的敏感性，从而增强了联合污染物对摇蚊的毒性效应。高含量的残留油类可能会改变沉积物的微生物群落结构和功能，影响微生物对污染物的降解作用，进一步加剧了污染物对摇蚊的危害。3.4相关性分析对摇蚊毒性指标（体长、体重、生存率等）与沉积物中残留油类、拟除虫菊酯类农药和多环芳烃含量进行Pearson相关性分析，结果如表5所示。摇蚊幼虫的体长与沉积物中溴氰菊酯、氯氰菊酯、高效氯氟氰菊酯、萘、菲、芘、苯并[a]芘的含量均呈现显著的负相关关系（P<0.01）。其中，与溴氰菊酯含量的相关系数为-0.852，与苯并[a]芘含量的相关系数为-0.825，表明随着这些污染物含量的增加，摇蚊幼虫的体长显著减小，生长发育受到明显抑制。污染物体长体重生存率残留油类0.125（P>0.05）-0.108（P>0.05）0.086（P>0.05）溴氰菊酯-0.852**-0.836**-0.815**氯氰菊酯-0.830**-0.812**-0.798**高效氯氟氰菊酯-0.845**-0.828**-0.805**萘-0.756**-0.738**-0.715**菲-0.785**-0.762**-0.740**芘-0.802**-0.780**-0.765**苯并[a]芘-0.825**-0.808**-0.788**注：**表示在0.01水平上显著相关（双侧）摇蚊幼虫的体重与上述污染物含量也呈现显著的负相关关系（P<0.01）。与氯氰菊酯含量的相关系数为-0.812，与芘含量的相关系数为-0.780，说明污染物含量的升高会导致摇蚊幼虫体重增长缓慢，影响其正常的物质积累和生长。在生存率方面，摇蚊幼虫的生存率与沉积物中拟除虫菊酯类农药和多环芳烃的含量同样呈现显著的负相关关系（P<0.01）。