模拟北大西洋公海核释放对定居及洄游鱼类的生态后果

《Fishes》：Modeling the Ecological Consequences of an Open-Sea Nuclear Release on Resident and Migratory Marine Fish in the North Atlantic

　　

摘要：核动力舰船或运载核武器的船舶发生意外放射性核素释放，可能使海洋生物暴露于放射性污染物中。此类事件的后果对于定居鱼类种群与经过受影响区域的洄游物种之间可能存在显著差异。尽管研究人员已针对沿海核事故评估了这些对比性响应，但公海情景下的类似评估仍然有限。研究人员模拟了北大西洋船舶放射性核素释放，以评估定居鱼类和洄游鱼类的暴露途径及潜在生态效应，因为若干运营核动力舰船的海军基地位于北美东海岸沿线。放射性核素的物理输运采用拉格朗日框架模拟，该框架纳入了洋流平流、三维湍流扩散、放射性衰变以及动态水-沉积物交换。该输运模型与一个四营养级食物网生物累积模型耦合，该食物网模型代表浮游植物、浮游动物、非食鱼性鱼类和食鱼性鱼类。该方法能够评估定居鱼类种群以及高流动性物种沿洄游路线的污染物摄取情况。蓝鳍金枪鱼（Thunnus thynnus）作为一种具有重要生态和商业价值的物种，被选为洄游案例研究对象。从电子标记数据重建的洄游路径被整合到模型中，以估算沿个体轨迹的放射性核素暴露。该组合建模框架为公海核释放如何对北大西洋定居鱼类群落和广域洄游物种产生差异化影响提供了新见解，对鱼类生态学、种群风险评估和渔业管理具有参考意义。

论文解读：北大西洋公海核释放对定居与洄游鱼类生态后果的模拟研究评述

一、研究背景与科学问题

海洋环境中放射性核素的人为释放可能来源于多种途径，其中核动力舰船或运载核武器的船舶发生事故是重要的潜在风险源。相较于陆基核设施事故（如沿海核电站）已有较多数值模拟研究，针对公海船舶核事故的生态后果评估十分有限。以往关于船舶核事故的模拟研究，如日本海K-431潜艇事故、库尔斯克号潜艇事故模拟、巴伦支海B-159潜艇事故分析等，大多未考虑放射性核素在水体与沉积物之间的相互作用，默认所有放射性核素始终溶解于水体中；对生物吸收的模拟也多采用简单的浓度因子平衡模型，缺乏动态食物网过程。部分针对陆基核设施事故的研究已将动态食物网模型纳入放射性核素输运模型，但此类模型尚未被系统应用于公海船舶核事故情景。

定居鱼类种群与经过受影响区域的洄游物种对核事故的响应可能显著不同。尽管研究人员此前已针对沿海核事故评估了这种对比性响应，但公海情景下缺乏可比照的评估。随着核动力海军舰队的扩张以及沉没或废弃核潜艇受到越来越多的地缘政治关注，此类公海情景的现实意义日益增强。无论是船体结构劣化还是打捞作业中的意外事件，均可能导致放射性核素释放，因此亟需建立能够同时评估定居鱼类和广域洄游物种生态后果的预测性建模框架。

本研究旨在模拟北大西洋公海船舶放射性核素释放，以评估其对定居鱼类和洄游鱼类的潜在生态效应。研究选用¹³⁷Cs作为目标放射性核素，因其化学性质与钾相似，易于在鱼肉中蓄积。大西洋蓝鳍金枪鱼（*Thunnus thynnus*）因其重要的生态和社会经济价值及高度洄游行为，被选为洄游案例研究对象。研究模拟了两种事故情景：公海瞬时释放和近岸释放，并比较了定居鱼类和洄游鱼类体内的放射性核素浓度。该论文发表在《Fishes》。

二、主要关键技术方法

研究人员采用耦合的拉格朗日放射性核素输运模型与四箱动态食物网生物累积模型。拉格朗日模型以粒子追踪方式表征放射性核素释放，纳入三维洋流平流、三维湍流扩散、放射性衰变及水-沉积物动态交换（以单一可逆反应动力学系数描述）。食物网模型包含浮游植物、浮游动物、非食鱼性鱼类和食鱼性鱼类四个营养级，各营养级通过水体直接吸收和摄食前一个营养级累积放射性核素，清除速率由各营养级特异性生物学半衰期控制。洋流数据来自HYCOM全球海洋环流模型（水平分辨率0.08°，32个垂向层次），模拟域覆盖西经25°–47°、北纬?82°至?62°。蓝鳍金枪鱼洄游轨迹来自斯坦福大学在美国和加拿大外海开展的电子标记释放实验（15条轨迹），仅选取在设定区域内停留至少90天的个体用于最坏情景分析。假设每次事故瞬时释放1 PBq的¹³⁷Cs。公海事故位置为长尾鲨号潜艇沉没点（41°46′ N, 65°3′ W），近岸事故位置距离佐治亚州金斯湾潜艇基地约100 km、水深50 m处。

三、主要研究结果

（一）公海事故情景：Th系列实验

以实验Th-15为例，模拟结果显示事故区域存在气旋式环流，放射性核素云团向海岸方向输运，该环流南侧受湾流约束。尽管释放点相对靠近模型边界，但局地环流驱动云团向模型内部移动。非食鱼性鱼类受¹³⁷Cs影响的区域大于水体受影响区域，这是因为鱼类摄入放射性核素后，由于清除速率缓慢，即使受污染水体已被洋流冲刷离开，鱼类体内仍保留放射性核素。

时间序列分析表明，蓝鳍金枪鱼在事故后第16天开始出现体内¹³⁷Cs浓度积累，随后在第40天和第60天因摄食更高度污染的饵料而出现更急剧的浓度升高，之后随着金枪鱼离开污染区域而缓慢清除。实验Th-15中金枪鱼肌肉最大¹³⁷Cs浓度约为1 Bq/kg。尽管金枪鱼经过高度污染的水体斑块，但其通过时间过短，不足以导致肌肉中大量蓄积¹³⁷Cs。沉积物作为放射性核素的延迟来源对金枪鱼肌肉浓度无明显影响。

传统方式计算的定居型食鱼鱼类（静态居民）最大¹³⁷Cs浓度约为30 Bq/kg，比洄游金枪鱼最大浓度高约30倍。深水沉积物未被放射性核素到达，仅部分大陆架受到污染，浓度可达10⁴ Bq/kg，甚至高于福岛事故后日本近海的实测值。这归因于密度跃层对垂直混合的屏障作用以及海洋垂直混合强度远弱于水平混合。

所有公海事故实验中，金枪鱼肌肉的最大浓度均处于相同数量级，表明公海位置事故对洄游金枪鱼几乎不产生显著影响。

（二）近岸事故情景：Co系列实验

以实验Co-14为例，模拟结果显示受湾流影响，污染水羽向东北方向延伸，随后放射性核素随弯曲和涡旋向南输运。尽管是瞬时释放，但污染羽从释放点扩展后，在事故位置东北方向数公里处出现明显浓度峰值，这归因于沉积物作为延迟放射性核素源的作用。由于释放发生在近海底，床底沉积物受到严重污染，放射性核素缓慢再溶解，使沉积物成为长期延迟源。

金枪鱼沿路径经过严重污染的沉积物区域，但整个模拟期间金枪鱼肌肉浓度始终低于1 Bq/kg。相比之下，定居食鱼鱼类的¹³⁷Cs浓度达到3×10⁴ Bq/kg。

四、讨论与结论

环境背景水平方面，北大西洋表层水典型¹³⁷Cs背景浓度约1.2 Bq/m³，为1950–1960年代大气核武器试验的遗留。非食鱼性鱼类背景浓度约为10^?2 Bq/kg湿重；食鱼性远洋鱼类通常低于10^?2 Bq/kg。蓝鳍金枪鱼全球平均¹³⁷Cs背景浓度约0.10 Bq/kg湿重。日本设定食品中¹³⁷Cs安全限值为100 Bq/kg，欧盟设定为1250 Bq/kg。

与以往观测和模拟研究比较，本研究预测的公海实验金枪鱼最大浓度约0.78 Bq/kg湿重，虽然高于环境背景值，但处于福岛事故后实测范围的较低区间，比国际食品安全限值低2–3个数量级。定居鱼类则表现出强烈得多的响应：非食鱼性鱼类达10³ Bq/kg量级，定居食鱼性鱼类约30 Bq/kg，远高于背景。结果与福岛事故后太平洋蓝鳍金枪鱼实测升高但低于安全限值的观测一致。

生态意义方面，核心结果不在于放射性核素浓度的绝对大小，而在于定居与高度洄游物种之间的显著差异。定居食鱼性鱼类浓度约为洄游金枪鱼最大值的30倍，比背景高两个数量级。将鱼类视为静止定居者的传统模型会显著高估洄游物种的放射性核素摄取量，说明在辐射生态风险评估中必须纳入真实的移动行为。放射性核素在鱼类中的积累受营养位置、摄食生态和代谢特征强烈影响，本研究中洄游金枪鱼极低的预测浓度表明其生态后果可忽略不计。污染沉积物局限于大陆架，不会进入深海，但动态水-沉积物交换模型捕捉到了其向水体缓慢再释放的过程，类似于爱尔兰海Sellafield核燃料后处理厂释放后的观测效应。

模型局限性包括：实验设计为最坏情景（事故日期与金枪鱼入域日期同步），预测浓度为保守上限值；浮游动物被湾流平流的过程未显式表示；可用洄游轨迹仅代表大西洋蓝鳍金枪鱼行为多样性的一部分；生物摄取参数和HYCOM洋流场存在不确定性；模型仅考虑直接海洋释放，未包括大气沉降输入。但所有这些局限性均不太可能改变主要结论，因为定居与洄游鱼类之间的显著差异主要由暴露持续时间控制，而非个别模型参数的不确定性。

研究结论：耦合拉格朗日放射性核素输运模型与四箱动态食物网模块的建模框架被用于模拟北大西洋船舶核事故后的¹³⁷Cs扩散与生物累积。公海事故情景下，各次模拟中洄游大西洋蓝鳍金枪鱼体内累积的¹³⁷Cs浓度均保持低水平，最高暴露个体最大值约0.78 Bq/kg湿重。放射性核素积累主要由金枪鱼在受污染水体中的短暂停留时间控制。定居食鱼性鱼类在事故点附近累积约30 Bq/kg，比任何洄游金枪鱼最大计算值高约30倍，比背景高两个数量级。污染沉积物局限于大陆架，但动态水-沉积物交换模型捕捉到其作为水体二次¹³⁷Cs源的作用。事故区域气旋式环流受湾流南侧约束，控制了大尺度扩散格局。近岸事故情景因位于大陆架上方，计算得到比公海案例更强的沉积物相互作用和更高的水体及定居生物浓度。总体而言，该研究将公海核舰船事故后果评估扩展到了静态和洄游鱼类种群，为该具有战略和生态意义的区域提供了新的定量估算，建模框架可迁移至其他海洋区域和污染物情景，为事故后环境评估和渔业管理规划提供了决策支持工具。