基于漂流浮标观测分析日本福岛以东表层海洋物质的散播轨迹*

卢 锡 胡石建, 4 王 凡, 4 胡敦欣, 4 刘海龙 刁新源, 4 于仁成, 4

基于漂流浮标观测分析日本福岛以东表层海洋物质的散播轨迹*

卢 锡1, 2, 3胡石建1, 2, 3, 4①王 凡1, 2, 3, 4胡敦欣1, 2, 3, 4刘海龙5, 3刁新源1, 2, 3, 4于仁成1, 2, 3, 4

(1.中国科学院海洋研究所 海洋环流与波动重点实验室 山东青岛 266071; 2. 中国科学院海洋大科学研究中心 山东青岛 266071; 3. 中国科学院大学 北京 100049; 4. 青岛海洋科学与技术试点国家实验室 山东青岛 266237; 5. 中国科学院大气物理研究所 大气科学和地球流体力学数值模拟国家重点实验室 北京 100029)

认识海洋中的物质如何散播对于理解海洋环境变化和人类活动污染在海洋中的扩散过程具有非常重要的意义。利用历史海表漂流浮标观测数据, 对日本福岛以东海域的表层物质散播轨迹进行了拉格朗日示踪分析和观测模拟试验研究。结果发现, 福岛以东海域海表浮标的散播路径主要分为东、南两支, 其中速度较快的东支为主要通道, 沿黑潮延伸体汇入北太平洋流, 最短用时大约22个月即可到达北美西海岸; 南支则沿黑潮延伸体以南的大范围南向流向西南方向运移, 速度较慢且明显受涡旋活动影响, 最快大约5个月即可到达吕宋海峡和中国台湾以东海域, 进而进入南海和东海等中国近海海域。通过开展观测模拟试验, 发现海表浮标散播的概率密度分布呈现以福岛附近海域为核心、向西南和正东方向递减扩展的形态, 其中, 到达中国近海的浮标主要通过吕宋海峡进入。文章详细讨论了研究结果的局限性、不足之处, 以及因基于大量现场观测而具备的重要参考价值。

福岛; 核污染; 西太平洋; 扩散轨迹; 漂流浮标

海洋中的水团和物质在海洋环流等海洋动力过程的作用下时刻发生着运移、扩散和循环, 塑造着海洋环境, 人类活动输入到海洋中的各类物质也在复杂的海洋动力和生物地球化学过程调整下发生着扩散和演变。因此, 研究物质在海洋动力过程调整下如何进行散播和扩散对于理解海洋环境变化和人类活动产生的污染物在海洋中的扩散过程具有非常重要的意义。

2011年3月11日, 日本东北部太平洋海域发生的强烈地震引发巨大海啸, 导致日本福岛第一核电站(Fukushima DAIICHI nuclear power plant, FDNPP)发生核泄漏。该事故是人类历史上最严重的核泄漏事故之一, 不仅对日本造成严重的经济损失和人员伤亡(Buesseler, 2011)、威胁周边国家的国民健康和生态环境(Fukuda, 2017), 也引发了全球性的大气和海洋污染(Aoyama, 2016)。研究表明该事故对海底沉积物的污染将持续数十年(Buesseler, 2014), 并将通过底栖鱼类的生物富集持续威胁周边国家的食物安全(Maderich, 2014)。此外, 受福岛附近早春盛行的强西风影响, 释放到大气的核物质往东扩散和沉积(Morino, 2011), 其中70%～80%扩散并沉积于西北太平洋(Yoshida, 2012), 约18%沉积在日本岛, 1.9%沉积在其他陆地(Stohl, 2012)。

事故发生后, 受损的核电站产生了严峻的进一步核污染问题。由于事故产生的大量核废水所含放射性物质成分极其复杂, 难以对其进行循环利用, 因而目前主要存于储水罐中(李小龙等, 2016), 到2021年4月15日已存核废水120.1577万t, 且每天新增核废水近110 t。2021年4月13日, 日本政府正式决定将福岛第一核电站的核废水排入海洋, 这导致西北太平洋可能将长期面临核废水输入和严重污染的威胁。放射性核素从福岛以东海域排放后在海洋中如何扩散, 成为国际社会广泛关注和迫切需要解决的重大问题之一(Chen, 2021; Zhao, 2021)。

物质在海洋中的扩散过程受到复杂的海洋动力过程调控, 其中海流引起的平流和扩散是最主要的海洋动力过程之一。福岛以东海域是世界上最强的洋流之一——黑潮及其延伸体所在海域, 其表层流速可达1.2 m/s, 流核位于100 m以浅(图1)。黑潮经其延伸体汇入北太平洋流后, 经加利福尼亚流等进入北赤道流, 北赤道流西行回到西太平洋后分叉形成黑潮和棉兰老流等西边界流, 对西太平洋和东亚、东南亚国家沿岸产生影响(图1, Hu, 2015, 2016)。2011年核事故产生的核污水在黑潮及其下游海流的平流作用下传播到中太平洋, 进而散布到整个太平洋海盆(Kawamura, 2011; Aoyama, 2013a)。

图1 北太平洋海盆流场分布(a)以及福岛核电站附近洋流以及生命周期达到12个月的浮标轨迹起点位置分布(b)

注: 颜色和箭头表示LICOM (LASG/IAP climate ocean model)数值模式输出的50 m以浅气候态平均速度, 对速度不低于0.05 m/s的流场区域, 使用黑色矢量图显示, 不足0.05 m/s的流场区域, 则使用蓝色矢量图显示; b中五角星表示福岛第一核电站位置, 黑色方框为所选浮标轨迹起点采样区域, 三角形表示生命周期达到12个月的浮标起点位置, 虚线分割区域为基于漂流浮标观测数据展开的连续排放模拟试验取样区域

针对2011年福岛核事故产生的放射性核素对海洋的影响, 前人通过数值模拟(Behrens, 2012; Smith, 2017; Rong, 2018)和船舶走航水体采样数据分析(Leon, 2011; Aoyama, 2013b; Kaeriyama, 2013; Deng, 2020)等方法研究了放射性核素的传播路径、覆盖范围和传播机制。部分研究还通过布放漂流浮标对福岛放射性核素的早期传播进行了分析(Buesseler, 2012; Rypina, 2013)。前人研究表明, 2011年福岛核事故产生的放射性核素在海洋中的传播过程主要有3个阶段:

第1阶段, 事故发生后, 放射性核素沿福岛海岸的微弱南向流汇入黑潮, 随后沿黑潮及其延伸体往东输送, 黑潮成为该阶段核素传播的南部边界(Kawamura, 2011), 该过程大约需要7～9个月。在这个过程中, 由于涡旋的存在, 一方面形成局地高浓度放射性核素区(Buesseler, 2012; Rypina,2013); 另一方面也增强放射性核素的垂向传播(Budyansky, 2015; Garraffo, 2016), 促进放射性核素随副热带模态水的潜沉(Chen, 2021)。

第2阶段, 放射性核素传播到中太平洋以后, 一方面沿北太平洋流继续往东传播, 但是速度明显下降(Aoyama, 2013b, 2015; 刘广山, 2015); 另一方面通过副热带模态水的形成和潜沉向南迁移(Kumamoto, 2014; Budyansky, 2015; Garraffo, 2016; 许丽晓等, 2021), 同时还存在沿北太平洋流的向北扩散(Kaeriyama, 2013), 该过程大约在事故发生后1～3 a完成。

第3阶段, 放射性核素随北太平洋流传播到北美西岸(Leon, 2011; Behrens, 2012; Zhao,2015), 进而通过太平洋海洋环流系统散播到整个北太平洋海盆(Behrens, 2012; Povinec, 2013), 该过程大约在事故发生后3～6 a完成。最终, 这些放射性核素将随海洋环流系统传播到印度洋和大西洋(Nakano, 2012; Rossi, 2013)。

2011年福岛核事故对中国近海也具有显著影响。事故发生后, 2012年冬季在中国南海、台湾海峡等海域观测到来自2011年福岛核事故的放射性核素(Deng, 2020; 林武辉等, 2020); Aoyama等(2013b)通过船舶走航水体采样, 发现2011年福岛核事故产生的放射性核素到2013年已经传播到中国海域; 到2015年已经遍布中国海域(赵云霞等, 2017)。这些研究表明, 部分放射性核素可以穿越黑潮阻碍, 到达中国近海, 但是福岛核污水影响中国近海的具体路径和扩散过程却尚不清楚。

Chen等(2021)指出, 现存的120多万t核废水排入海洋之后可能存在3条路径: 受涡旋、黑潮南部再循环和湍流扩散过程影响而从福岛周边海域向西进入中国东海的路径1; 随黑潮向东并在副热带中东太平洋潜沉后进入吕宋海峡的路径2; 沿副热带环流圈到达北美西岸后, 随北赤道流到达西边界, 重新循环的路径3。Chen等(2021)提出的3条路径与2011年福岛核事故产生的核素传播过程吻合, 也与林武辉等(2021)指出的受多尺度海洋水文动力过程影响的“被动”迁移路径和Zhao等(2021)通过数值模拟得到的放射性核素传播通道一致。除了海洋动力过程导致的污染物质扩散以外, 生物载体驱动的“主动”迁移过程也具有重要作用, 且响应更快、可能对人类健康产生更大影响(林武辉等, 2021)。

前人开展的大量研究对于理解2011年福岛核事故泄漏的放射性核素扩散过程、探究未来福岛以东污染物质如何散播具有非常重要的意义, 但也还存在一些缺陷和不足之处。船舶走航水体采样数据站位较少, 难以准确刻画核素传播路径。数值模拟方面, 现有数值模式的模拟能力还有限, 模式模拟的海洋流场与实际观测之间存在显著差异, 对放射性核素扩散过程的模拟结果和实际观测之间存在明显不同(Periáeez, 2015)。以传播到北美西岸的时间为例, 不同数值模拟试验的结果差别很大, 广泛分布于3～6 a范围内(Behrens, 2012; Nakano, 2012; Rossi, 2013; Zhao, 2015, 2021), 但是仅在事故发生后2 a多, 就已经在加拿大附近海域观测到了来自福岛核电站的放射性核素(Smith, 2015)。此外, 数值模拟试验初始场设置对大气沉积的设定也可能影响分析结果(Tsumune, 2012)。因此, 从更多的观测数据中寻找放射性核素散播路径的证据具有很重要的意义。

为此, 文章将使用过去数10 a以来获取的海表漂流浮标观测数据进行轨迹示踪分析, 并开展基于观测的模拟试验, 以研究福岛以东海域表层海洋物质的传播路径和扩散过程。本文的第一部分将介绍所使用的数据和研究方法, 第二部分将重点介绍主要的研究结果, 第三部分讨论研究结果的局限性和主要影响因素, 最后在第四部分总结主要结论和意义。

1 数据和方法

1.1 观测数据及其处理

本文采用的主要数据是全球漂流浮标项目漂流浮标集(Global Drifter Program drifting buoy collection)。该数据集包括自1979年2月14日至2020年9月2日布放的所有漂流浮标的全球漂流浮标项目编号(identification, ID)、世界气象组织浮标识别编号(World Meteorological Organization, WMO)、浮标布放和回收的日期、位置, 以及浮标轨迹时间序列。该数据集对海表浮标观测进行了质量控制和6 h插值等处理(Centurioni, 2019), 浮标经纬坐标值平均误差小于0.01°, 能够比较准确地展示真实的表层海流。考虑到核废水的污染物质可能主要在表层快速扩散(Macdonald, 2020), 使用海表漂流浮标数据开展轨迹示踪分析具有合理性。

参考前人研究, 结合福岛核电站的位置(38.3°N, 142.4°E)及其附近海陆分布和洋流情况, 选取所有流经福岛以东海域(140°～150°E, 30°～45°N, 图1黑色方框)的海表漂流浮标作为分析对象。其中, 初始布放坐标位于所选区域内的浮标, 以该坐标作为其轨迹起点; 初始布放坐标位于该区域以外的浮标, 则以其进入该区域后的第一个坐标作为其轨迹起点。由于实际布放过程中存在对同一浮标的多次布放, 本文将存在多次布放的浮标的每一段轨迹视作独立的时间序列, 并分别确定其轨迹起点。考虑到周边海域表层海流流速情况, 将6 h内浮标轨迹变化超过1°的时刻视为突变点(即6 h内平均速度超过约5 m/s), 仅保留突变点之前的轨迹时间序列, 并将突变后的浮标轨迹视作在突变点的一次重新布放。

统计结果显示, 轨迹起点位于福岛以东海域的可用浮标共计692个(表1), 包括直接布放于该海域的219个浮标和从其他海域进入的473个浮标。这692个浮标布放于1989～2020年间, 最长生命周期超过10 a, 其中生命周期大于12个月的浮标有309个(由于数据集截止到2020年9月, 因而2020年布放的浮标生命周期均尚未超过12个月)。

表1 轨迹起点位于福岛以东海域(140°～150°E, 30°～45°N)的不同生命周期漂流浮标数, 以及进入吕宋海峡和中国台湾以东海域(120°～125°E, 12°～30°N)、北美西海岸海域(120°～130°W, 25°～50°N)的时间和数量

Tab.1 Numbers of drifting buoys with different durations launched from the region of eastern Fukushima (140°～150°E, 30°～45°N) and reached the region of Luzon Strait and eastern Taiwan (120°～125°E, 12°～30°N) and the western coast of North America (120°～130°W, 25°～50°N)

注: 首次有浮标进入两海域的月份用加粗带下划线标出, 生命周期月份统计为向上取整

此外, 为讨论流场的重要作用, 使用基于Argo浮标观测的格点化温盐数据产品RG Argo Climatology(Roemmich-Gilson Argo Climatology) (Roemmich,2009)计算了参考深度为2 000 m的地转流速度(Hu, 2021), 以及LICOM3 (LASG/IAP climate ocean model, version 3)输出了上层50 m气候态平均流场。RG Argo数据集使用加权最小二乘法拟合了2004～2020年的Argo浮标观测数据廓线, 其水平分辨率为1°×1°, 垂直方向2 000 m分58层。LICOM3数值模式较好地再现了海洋环流和温盐的变化, 其水平分辨率为0.1°×0.1°, 上层50 m分10层。

1.2 观测模拟试验

为模拟福岛以东海域表层海洋物质在不断排入的情况下的扩散过程, 基于漂流浮标观测数据设计开展了一系列连续排放的模拟试验。

生命周期达6、12、18和24个月的浮标数分别为478、309、214和157个。为模拟福岛附近海表物质的连续输入, 对生命周期达到特定月数(6、12、18和24个月)的浮标分别进行随机取样, 取样次数为相应月数。每次随机抽取的浮标, 视为1次在轨迹起点的重新布放, 将自抽取月份从轨迹起点开始运动, 完成所有取样后1个月所有浮标停止运动。由于大量浮标轨迹起点(359个)集中在采样海域西侧边界的黑潮流入区域(图1), 每月重新布放的浮标, 按照1︰24的比例(浮标轨迹起点散布于福岛以东海域的25个网格内)选取, 即选取其他24个网格内浮标总数的90%, 并将该数量的1/24作为在西侧边界浮标起点密集区域的浮标选取数。

对某一个1°×1°的网格, 浮标流经其位置的概率密度(单位: %)定义为

其中,为试验中流经该网格点的总浮标数,为试验中选取的总浮标数。

2 结果

2.1 福岛以东表层海洋物质在太平洋的散播轨迹

图2展示了浮标从日本福岛以东海域开始运动后的散播轨迹。结果显示, 浮标的散播轨迹大致可分为向东和向南散播的两个分支: 大部分浮标(超过70%)沿黑潮及其延伸体往东传播, 以下称为东支; 少部分浮标往西南运动, 比东向传播的浮标速度慢, 以下称为南支。东支浮标主要沿Chen等(2021)提出的路径2、路径3运动, 南支浮标则主要沿路径1运动。浮标轨迹散播方向和范围与Zhao等(2021)对3H浓度数值模拟的结果较为一致, 但是浮标轨迹数据显示的散播速度明显比模拟的传播速度更快。

为保证充足的漂流浮标观测数据, 所选692个浮标轨迹起点所在福岛以东海域范围(140°～150°E, 30°～45°N)大于前人研究中选取的范围(140°～145°E, 35°～39°S) (Buesseler, 2012; Rypina, 2013)。但是数值模拟(Behrens, 2012)和船舶走航水体采样数据分析(Aoyama, 2013b)结果均显示, 2011年福岛核事故产生的放射性核素仅需不足1个月即可遍布所选采样区域, 图2b中浮标开始运动1个月后的轨迹也已遍布采样海域。

在前人研究中, 黑潮及其延伸体被看作核素传播第1阶段(约7～9个月)的南边界(Kameník, 2013; Garraffo, 2016), 但是浮标在福岛以东海域开始运动5个月后, 即可到达吕宋海峡和中国台湾以东海域(图2c), 这比数值模拟的结果更快(Behrens, 2012); 而东支浮标轨迹则明显脱离黑潮延伸体, 越过日界线、进入广阔而流速较为缓慢的北太平洋流。Aoyama等2013年在日界线附近海域采样观测到源于2011年福岛核事故的核素(Aoyama, 2013b), 但图2c表明这可能并不是首批到达日界线的核素。

从日本福岛以东海域开始漂移大约11个月后, 大量浮标抵达吕宋海峡和中国台湾以东海域(图2d), 值得注意的是, 此时已有1个浮标进入南海并最终抵达我国福建沿海。分析表明, 该浮标生命周期不足11个月, 从福岛以东海域出发9个月后进入吕宋海峡, 从吕宋海峡到达福建沿海仅用了69 d。浮标从福岛以东海域出发后11个月到3 a期间, 吕宋海峡和中国台湾以东海域均存在比较密集的浮标轨迹, 且持续有浮标到达中国近海(图2d～2f)。

东支路径上的浮标轨迹抵达日界线后进一步往东延伸, 但是其北部范围无明显扩张。浮标从福岛以东海域出发后最快21个月即可到达加拿大西海岸, 这早于Smith等(2015)的研究结果(大约27个月); 22个月后, 已有多个浮标进入美国西岸海域(图2e); 33个月后的浮标运动轨迹已经遍布整个北美西岸(图2f)。在较早到达美国西岸海域的5个浮标中(截止时间为24个月), 有2个位于采样海域西侧边界黑潮流入区域, 3个位于黑潮延伸体以北。

总体上看, 福岛以东海域的浮标在开始运动大约3 a后, 其轨迹已广泛散布于北太平洋海盆, 这要比Zhao等(2021)数值模拟的3H传播结果更快。其中, 东支浮标在到达北美西海岸后, 一部分沿加利福尼亚流往南进入北赤道流, 进而向太平洋西边界海域运移; 另一部分则往北沿阿拉斯加流进入北太平洋副极地环流。进入北赤道流的浮标鲜有直接跨越北赤道流南边界, 大多数随北赤道流回到太平洋西边界海域, 这与Chen等(2021)提出的路径3一致。3 a后的浮标运动轨迹(图2f)分布范围已经无明显变化。但是, 由于生命周期达到3、4和5 a的浮标数从73锐减到15和3 (表1), 无法对更长时间的散播轨迹进行进一步的分析。

2.2 概率密度分布

为研究福岛以东海域表层海洋物质在不断输入情况下的密度分布, 基于浮标观测数据开展了一系列连续排放的模拟试验。图3表明, 连续排放观测模拟试验得到的概率密度分布与浮标散播轨迹的总体形态一致。不同截止时间的概率密度分布相互吻合, 整体模态为以福岛以东海域为核心、往正东和西南方向延伸和衰减, 这与Zhao等(2021)模拟的海表3H浓度传播模态一致, 也与Chen等(2021)提出的传播路径一致。

图2 692个浮标的轨迹起点(a)和从福岛以东海域开始运动的生命周期超过1个月(b)、5个月(c)、11个月(d)、22个月(e)、33个月(f)的浮标的运动轨迹

注: 黑色方框为福岛以东浮标起点采样区域; 黄色三角形和玫红色圆形分别代表浮标轨迹的起点和终点, 到达中国近海的浮标轨迹已经用黄色加粗; 红色曲线表示浮标运动轨迹; b～f中的浮标数分别为692、511、325、151和89

南支浮标概率密度表现为吕宋海峡和中国台湾以东海域持续增加, 该海域截止时间为6个月的试验的概率密度为0.1% (图3a), 而截止时间为24个月的试验的概率密度为2%。东支浮标概率密度表现为持续地向东扩散, 截止时间为24个月的观测模拟试验(图3d)结果显示, 北美沿岸已经有概率密度信号, 早于Zhao等(2021)模拟的结果。

值得注意的是, 截止时间为12个月的试验结果显示, 南海有概率密度信号(图3b), 但截止时间为24个月的试验结果显示, 南海没有概率密度信号(图3d), 这是由于进入南海的浮标生命周期均未达到24个月, 因而未能体现在图3d (截止生命周期为24个月)中, 但这并不代表布放后24个月内没有浮标进入南海。

1.池塘准备 3月初，池塘施用75kg/亩生石灰进行干法清塘，3天后排干池水晒塘。放苗前10天左右，池塘水位加至10cm，以4～5m的行距、2～3m的株距种植轮叶黑藻。种植一周后提高水位至50cm左右，根据池塘肥力，放苗前7天左右使用一次生物肥，培育生物饵料。池塘进水用100目长筒状筛绢网过滤，防止敌害进入。前期保持较低水位，有利于轮叶黑藻生长和水体升温。

2.3 东支路径分析

轨迹起点位于福岛以东海域且最终到达北美西岸海域的浮标共有27个(表1, 布放于1993～2016年), 其中最短用时22个月, 最长用时63个月, 这意味着物质在进入福岛以东表层海洋后, 将在2～5 a内持续往北美输送。这27个浮标的轨迹起点主要位于黑潮及其延伸体范围内(图4a), 黑潮延伸体以南的197个浮标中有19个生命周期超过22个月, 但这19个浮标均没有流经北美西岸。较早到达美国西岸海域的5个浮标中(截止时间为24个月), 有2个源于采样海域西侧边界黑潮流入区域, 3个源于黑潮延伸体以北。在轨迹起点位于黑潮流入区域且生命周期超过22个月的67个浮标中, 有12个流经美国西岸海域(最快23个月, 最慢63个月); 在轨迹起点位于黑潮延伸体以北且生命周期超过22个月的65个浮标中, 有15个流经北美西岸海域(最快22个月, 最慢41个月)。北太平洋流在东太平洋海域分叉并形成向北的阿拉斯加流和向南的加利福尼亚流(Cheng, 2014), 部分浮标随阿拉斯加流往北运移, 但是这部分浮标随后并没有完全汇入阿拉斯加流, 一部分在近岸处向南转向, 最终沿加利福尼亚流进入北赤道流, 最终仅有少量浮标进入阿拉斯加流。

为探究北太平洋流分叉后浮标的再分配过程, 专门针对北太平洋流分叉区域(145°～155°W, 30°～48°N, 图4b)和北美西岸海域(123°55′～130°W, 30°～50°N, 图4c)进行了浮标轨迹分析。轨迹起点位于北太平洋流分叉区域的46个浮标, 只有1个在到达北美西岸后往北汇入阿拉斯加流(图4b)。轨迹起点位于北美西岸的335个浮标(图4c), 绝大多数往南或西南运动, 仅有14个浮标往北运动, 其中1个最终到达阿留申群岛。进入阿拉斯加流的浮标, 绝大部分是沿北美西海岸向北而来, 仅有少量直接来自北太平洋流在东太平洋的北分支(图4c)。东支路径上的浮标, 大部分到达北美西岸后往南汇入加利福尼亚流, 最终通过北赤道流到达西边界, 沿副热带环流圈运动(图4a), 亦即沿Chen等(2021)提出的路径3运动。

图3 生命周期达到6个月(a)、12个月(b)、18个月(c)和24个月(d)的漂流浮标分别进行连续排放观测模拟试验的概率密度分布图

注: 黑色方框为所选浮标布放区域; a～d的浮标数分别为478、309、214和139

图4 进入北美西岸海域(120°～130°W, 25°～50°N)的27个浮标的运动轨迹(a), 起点位于北太平洋流分叉海域(145°～155°W, 30°～48°N, 黑色方框)的46个浮标的运动轨迹(b), 起点位于北美西岸海域(123°55′～130°W, 30°～50°N, 黑色方框)的335个浮标的运动轨迹(c)

注: 黑色方框为所选浮标轨迹起点所在区域, 黄色三角形和玫红色圆形分别代表浮标轨迹的起点和终点; 红色曲线表示浮标运动轨迹

2.4 南支路径分析

南支路径上的浮标比东支更早进入中国近海, 因此非常重要。轨迹起点位于福岛以东海域并最终进入吕宋海峡和中国台湾以东海域的19个浮标(均于2002年以后开始运动), 最短用时5个月, 最长用时32个月(表1), 轨迹起点和运动轨迹如图5a所示。到达吕宋海峡和中国台湾以东海域后, 这19个浮标的漂移轨迹存在3种路径(其中4个浮标的生命周期止于吕宋海峡): 往南进入中国南海(4个)、往北经台湾暖流穿越台湾海峡进入中国东海(2个)以及通过黑潮流套进出吕宋海峡(9个), 这意味着到达吕宋海峡和中国台湾以东海域的浮标, 约有40%进入了我国南海和东海。黑潮以东存在强烈的中尺度涡旋活动和显著的涡致经向输运(Xu, 2016, 2017), Chen等(2021)也指出黑潮以南的再循环和涡旋是调控核素往我国近海传播的重要因素, 南支的浮标在穿越黑潮延伸体前后, 运动路径很弯曲、复杂(图5), 因此, 中尺度涡旋活动可能发挥了非常重要的作用。

为进一步探究浮标进入吕宋海峡后的运动, 重新选取了轨迹起点位于吕宋海峡附近海域(120°～123°E, 18°～23°N)的810个浮标并进行轨迹分析(图5b)。结果显示, 源于吕宋海峡附近海域的浮标轨迹遍布我国东南沿海, 往北逼近长江入海口, 往南经卡里马塔海峡进入印尼海。其中167个浮标进入中国南海和台湾海峡; 进入东海的20多个浮标, 一部分从南海穿越台湾海峡直接进入, 另一部分从中国台湾东北海域脱离黑潮主轴而进入。

图5 进入吕宋海峡和中国台湾以东海域(120°～125°E, 18°～30°N)的19个浮标的运动轨迹(a); 起点位于吕宋海峡(120°～123°E, 18°～23°N)的810个浮标的运动轨迹(b); 起点位于中国台湾东部沿海区域(122.5°～124°E, 24°～26.5°N)的584个浮标运动轨迹(c)

注: 黑色方框为所选浮标轨迹起点所在区域, 黄色三角形和玫红色圆形分别代表浮标轨迹的起点和终点; 红色曲线表示浮标运动轨迹; a中进入中国近海的浮标轨迹用深黄色标出

为研究浮标进入中国台湾以东海域后的具体运动轨迹, 对轨迹起点位于中国台湾东部海域(122°～ 123.5°E, 24°～26.5°N)的584个浮标的轨迹进行分析(图5c)。源于该海域的浮标运动轨迹主要受黑潮及其延伸体影响, 但值得注意的是, 有一部分可从东海南部绕过中国台湾岛、经台湾海峡进入南海。

从吕宋海峡开始运动的浮标, 其概率密度分布模态为以吕宋海峡为核心, 分别沿黑潮和南海贯穿流流向东北-西南方向递减(图6)。东北方向沿黑潮及其延伸体运动, 伴随黑潮延伸体以东的缓慢散播。西南方向主要在南海内部散播, 扩散到中国南海沿岸和越南近岸等海域。仅1个月后, 南海北部有大片概率密度信号(超过10%), 且3、6、12、18、24个月后的浮标概率分布模态(图6b～6f)渐趋一致。

图7展示了浮标从中国台湾东北海域开始运动1、3、6、12、18和24个月后的概率密度分布。结果显示, 在开始运动后不同月份, 其概率密度均以黑潮及其延伸体为核心, 并沿黑潮流向衰减(图7)。黑潮区域内浮标概率密度可达90%以上, 黑潮以外区域的浮标概率密度较小——不到10%。浮标开始运动3个月后, 我国南海(1%)、东海(2%)和黄海南部区域(2%)均有浮标概率密度信号(图7b)。从图7中浮标概率密度分布的演变可以看出, 浮标到达中国台湾东北海域后, 主要有3条路径: (1)沿黑潮流轴扩散; (2) 随黑潮以南的再循环扩散并到达吕宋海峡、进入南海; (3) 跨越东海陆架并沿台湾海峡扩散到东海和南海。

值得注意的是, 在截止时间超过12个月的试验中, 从吕宋海峡(图6d～6f)和中国台湾以东海域(图7d～7f)开始运动的浮标概率密度分布主要模态渐趋稳定, 但是在我国近海范围概率密度逐渐减小甚至消失, 这是因为从两区域出发的浮标, 一旦进入我国近海, 其生命周期会迅速结束, 这就导致大部分浮标生命周期短于截止时间, 因而未能在概率密度分布中体现, 但是这并不代表截止时间内没有浮标进入我国近海。

图6 轨迹起点位于吕宋海峡附近海域(120°～123°E, 18°～23°N)且生命周期到达1个月(a)、3个月(b)、6个月(c)、12个月(d)、18个月(e)、24个月(f)的漂流浮标分别进行连续布放模拟试验的浮标概率密度分布图

注: 黑色方框为所选浮标轨迹起点所在区域; a～f中的浮标个数分别为810、493、232、102、63和44

图7 轨迹起点位于中国台湾以东海域(122°～123.5°E, 24°～26.5°N)且生命周期到达1个月(a)、3个月(b)、6个月(c)、12个月(d)、18个月(e)、24个月(f)的漂流浮标分别进行连续布放模拟试验的浮标概率密度分布图

注: 黑色方框为所选浮标轨迹起点所在区域; a～f中的浮标个数分别为584、372、219、117、78和51

2.5 流场的作用: 重要个例分析

图8a展示了进入中国近海的浮标运移轨迹以及2014年7月50 m以浅的平均地转流场(参考深度2 000 m), 其中深黄色粗线是2014年7月开始运动(35°32′N, 145°38′E)且最终到达吕宋海峡的浮标(记为A)的运动轨迹, 图8b展示了LICOM数值模式输出的福岛附近海域50 m以浅气候态平均流场。福岛以东海域存在明显的表层南向流(图8), 且南向速度可达0.4 m/s (图8b), 成为福岛以东海域表层海洋物质穿越黑潮并向西南运移的重要通道。图8b显示, 143°～ 145°E处黑潮延伸体向东南弯曲, 在其以南区域, 尽管海表流速较小, 但是仍有明显南向流速分量。在这些南向流场的作用下, 福岛以东海洋表层的物质可以向西南方向运动, 直抵吕宋海峡和中国台湾以东海域(图8a)。

从前述结论可以看出, 最终到达吕宋海峡和中国台湾以东海域的浮标, 其轨迹起点多位于黑潮延伸体以南, 因此, 难以确定黑潮延伸体内是否存在浮标向南穿越的通道。但是, 浮标A (轨迹起点位于35°32’N, 145°38’E, 布放于2014年7月)的轨迹和海表平均地转流场(2014年7月, 参考深度2 000 m)分布显示, 浮标A轨迹起点位于黑潮延伸体范围内, 并在穿越黑潮延伸体后呈涡旋性运动(图8a), 最终到达吕宋海峡(2017年2月), 在此之前经历了31个月的曲折运动。浮标A穿越黑潮后的涡旋性运动, 可能受反气旋涡对北太平洋副热带模态水的西南向运输影响(Liu, 2007; Xu, 2017)。由于浮标A的轨迹与最终进入中国近海的6个浮标轨迹存在重叠(图8a), 福岛附近表层海洋物质可以通过类似浮标A的路径穿越黑潮, 经黑潮及其延伸体以南较弱南向流到达吕宋海峡, 在黑潮(Chen, 2011a)或中尺度涡(Chen, 2011b; 许丽晓等, 2021)等影响下最终进入中国近海, 这表明福岛以东海表物质可以直接穿越黑潮延伸体并经南向散播路径进入中国近海。

图8 进入中国近海的浮标运移轨迹和基于RG_Argo数据计算的2014年7月50 m以浅的平均地转流速度分布(a)和LICOM数值模式输出的福岛附近海域50 m以浅气候态平均流场(b)

注: a中地转流计算的参考深度为2 000 m, 黑色方框为所选浮标轨迹起点所在区域, 黄色三角形和玫红色圆形分别代表浮标轨迹的起点和终点, 黄色曲线表示浮标A的运动轨迹, 红色曲线表示进入中国近海的浮标运移轨迹; b中颜色表示气候态平均经向流, 向北为正; 速度超过0.05 m/s的流场区域, 使用黑色矢量图显示, 不足0.05 m/s的流场区域, 使用蓝色矢量图显示

3 讨论

研究结果还存在缺点和不足之处, 以下从流场变化和浮标轨迹起点区域选择的影响以及结果的不确定性等方面进行讨论。

3.1 流场变化的影响

浮标运动主要由流场驱动, 因此对单个浮标来说, 流场的变化会对该浮标的轨迹产生影响。为讨论流场的季节变化对浮标轨迹的影响, 图9a～9d展示了福岛以东海域生命周期超过12个月的309个浮标在不同季节开始运动的轨迹。从结果可见, 不同季节开始运动的浮标轨迹存在显著的差异(注意, 夏秋季开始运动的浮标比冬春季多), 主要体现在轨迹路径和终点的具体位置以及少数远距离散播的浮标位置。但是, 从浮标轨迹的总体分布来看, 不同季节的轨迹分布范围大体一致, 东支和南支路径在4个季节均存在。这表明在不同时间释放粒子, 尽管具体的流程可能有较大差异、对粒子个体的轨迹或有影响, 但对粒子轨迹的统计特征影响不大。

太平洋的年际和年代际变化也很显著, 比如厄尔尼诺-南方涛动(El Niño-Southern Oscillation)和太平洋年代际振荡(Pacific Decadal Oscillation)都有可能导致太平洋流场的年际和年代际变化(Hu, 2014, 2017, 2020a, 2020b), 不同年份释放的粒子可能会有不同的散播路径。为探究流场的低频变化对浮标轨迹的影响, 我们对比分析了2000～2009年(125个)和2010～2019年(140个)从福岛以东海域开始运动且生命周期超过12个月的浮标轨迹(图9e～9f)。结果类似于不同季节的粒子释放情景, 在两个十年期间开始运动的浮标轨迹范围和传播路径统计特征较为一致。这表明在大量粒子释放的情况下, 其散播轨迹路径是具有较高可信度的。

图9 轨迹起点位于福岛以东海域且生命周期超过12个月的浮标分别于春季(3～5月, a)、夏季(6～8月, b)、秋季(9～11月, c)和冬季(12月～次年2月, d), 以及2000～2009年(e)和2010～2019年(f)开始运动的轨迹

注: 黑色方框为所选浮标轨迹起点所在区域; 黄色三角形和玫红色圆形分别代表浮标轨迹的起点和终点; 红色曲线表示浮标运动轨迹; a～f中的浮标数分别为56、102、93、58、125和140

事实上, 按照日本政府公布的计划, 核废水排放过程将会持续20～30 a, 所选浮标从福岛以东海域开始运动的时间也涵盖了近20多年, 可以认为粒子在散播过程中受到的流场影响可能是随机分布的, 其统计特征是稳定的。尽管北太平洋海洋环流存在多种时间尺度的变化, 利用所选能够较好地刻画核素实际的传播路径和范围, 因此具有一定的参考意义。

3.2 浮标轨迹起点区域范围的影响

前人研究中选取的用以探究2011年放射性核素散播路径的漂流浮标布放范围为140°～145°E, 35°～39°S区域(Buesseler, 2012; Rypina, 2013), 由于流经该海域的浮标数较少(共422个, 其中生命周期超过12个月的浮标184个), 本文参考福岛以东海陆分布和流场特征, 并在分析浮标运动轨迹的基础上, 选择流经140°～150°E, 30°～45°N区域的692个浮标(生命周期超过12个月的浮标309个)进行分析。

图10a展示了轨迹起点位于前人采样区域(140°～ 145°E, 35°～39°S)且生命周期超过12个月的184个浮标的运动轨迹, 整体也大致分为向东和向南两个分支, 在北美西岸和中国近海均有浮标到达, 就统计特征而言, 其整体轨迹分布与源于本文采样区域(140°～150°E, 30°～45°N)的浮标比较一致。图10a中前人采样区域内的浮标开始运动后1个月, 其轨迹明显覆盖并超出本文浮标采样海域, 这意味着采样海域范围的扩大对浮标轨迹散播速度和范围影响较小。

轨迹起点位于本文所选海域北部、中部和南部的浮标轨迹如图10b～10d所示, 从不同海域开始运动的浮标轨迹表现出较大不同, 浮标轨迹的具体路径以及个别浮标的终点位置有较大差异, 而且从南部海域(190个)出发的浮标远多于从北部(66个)、中部海域(53个)出发的浮标, 这也导致图10d中的浮标轨迹密集得多。整体来看, 东支路径在从不同海域出发的浮标轨迹中均比较明显, 遍布整个北太平洋海盆。南支路径在从中部和南部海域出发的浮标轨迹中比较清晰(图10c～10d), 但是图10b中展示的从北部海域出发的浮标, 仅有1个明显经黑潮延伸体南部微弱南向流往南传播, 而且最终没有到达太平洋西边界。考虑到北部海域出发浮标数较少, 而且到达西边界的浮标生命周期多数不足12个月, 北部海域的浮标较难通过南向路径散播到太平洋西边界, 较难通过Chen等(2021)提出的路径1进入我国近海, 但并不意味着路径1对北部海域完全封闭。

图10 轨迹起点分别位于福岛附近前人采样海域(140°～145°E, 35°～39°N, a)和本文采样海域北部(140°～150°E, 40°～45°N, b)、中部(140°～150°E, 35°～40°N, c)和南部(140°～150°E, 30°～45°N, d)且生命周期超过12个月的浮标运动轨迹

注: 黑色方框为所选浮标轨迹起点所在区域; 黄色三角形和玫红色圆形分别代表浮标轨迹的起点和终点; 红色曲线表示浮标运动轨迹; a～d中的浮标数分别为184、66、53和190

在浮标筛选过程中受黑潮影响, 浮标轨迹起点在所选海域西南边界黑潮流入区域富集, 尽管在观测模拟试验中控制了该区域抽取浮标数, 黑潮的影响仍有可能被放大。前人研究所选区域(140°～145°E, 35°～39°S)主要位于黑潮延伸体以北, 虽然其南侧仍然受到黑潮延伸体的影响, 但是该区域较小且贴近福岛, 黑潮对该区域浮标概率密度分布影响也较小。轨迹起点位于前人采样海域的浮标截止时间为6个月和12个月的概率密度分布范围(图11)只比源于本文所选海域的浮标概率密度分布范围中略小(图2a～2b), 但是二者量级和模态较为一致, 黑潮对浮标散播轨迹的影响可能并没有被放大。

3.3 不确定性与缺陷

对于分析福岛以东表层海洋物质扩散、尤其是核素扩散这种复杂过程而言, 仅用海表浮标来开展分析显然存在不确定性和缺陷。

一方面, 浮标是刚体, 主要受表层海洋平流动力过程的影响, 而表层海洋物质的扩散还受到混合、生物地球化学等过程的影响。比如, 林武辉等(2021)指出福岛核废水中的物质迁移存在水文动力驱动的“随波逐流”和以生物为载体的“搭乘便车”两种方式, 而本文仅考虑了其中一部分物理过程。此外, 可用浮标仅692个, 且生命周期超过12个月的浮标仅309个, 浮标出发点的时间跨度过长、空间分布也不均匀, 也可导致潜在的不确定性。

另一方面, 浮标观测数据自身也存在缺陷。现有的历史浮标观测数据库中, 可用浮标寿命长短不一, 具有较大随机性, 只能通过对浮标轨迹分析后进行“接力”取样, 以获得更长的浮标轨迹时间序列。通过在吕宋海峡、中国台湾以东海域、中东太平洋和北美西岸海域进行“接力”取样, 虽然能够推测后续的浮标运动, 但是“接力”前后浮标轨迹并不连续, 存在时空误差, 且“接力”采样范围较大, 具有不确定性。生命周期超过3 a的浮标过少(73个), 限制了对跨海盆输运等长周期过程的分析。不同年份开始运动的浮标数差别很大, 2011年最多(73个), 12个年份不足10个。由此开展的长截止时间观测模拟试验中, 近岸浮标生命周期迅速结束, 导致岸界附近的概率密度被低估, 但是按照Zhao等(2021)数值模拟的结果显示, 长期来看(5～10 a后),3H浓度在日本以东海岸附近最高。

图11 轨迹起点位于前人采样海域(140°～145°E, 35°～39°S)且生命周期到达6个月(a)和12个月(b)的漂流浮标分别进行连续布放试验的概率密度分布图

注: 黑色方框为浮标起点采样区域; a和b中的浮标数分别为289和184

源于福岛以东海域的浮标, 截止时间为12个月的不同季节连续布放观测模拟试验中概率密度分布如图12a～12d所示, 不同季节的概率分布模态和范围比较类似, 但是高值区(>20%)差别较大, 主要集中在黑潮延伸体区域(图12a, 12c, 12d), 仅在夏季(图12b)高值区蔓延到整个福岛以东海域, 2000～2009年和2010～2019年的高值区分布也存在类似差异(图12e～12f)。Zhao等(2021)年基于数值模拟预测的3H浓度变化模态与图12e中模态更为接近, 但是浓度变化范围扩张速度要比浮标观测模拟试验的概率密度慢得多。浮标长短不一的生命周期限制了对平均态的分析, 概率密度分布能够揭示表层海洋物质传播的一般规律, 但是无法体现放射性核素浓度变化, 对周边海洋环境和人类健康的影响需要更多研究。

图12 起点位于福岛以东海域且生命周期到达12个月的漂流浮标分别于春季(3～5月, a)、夏季(6～8月, b)、秋季(9～11月, c)和冬季(12月～次年2月, d), 以及2000～2009年(e)和2010～2019年(f)进行连续布放模拟试验的浮标概率密度分布图

注: 黑色方框为所选浮标轨迹起点所在区域; a～f中的浮标数分别为56、102、93、58、125和140

4 小结

4.1 结论

利用全球漂流浮标项目的海表漂流浮标观测数据, 对1989～2020年间流经福岛附近海域(140°～150°E, 30°～45°N)的692个漂流浮标进行了轨迹诊断和示踪分析, 并基于这些历史浮标轨迹数据开展了浮标连续布放试验、计算了概率密度分布, 研究了福岛附近海域表层海洋物质可能的散播轨迹和主要路径。

漂流浮标轨迹分析显示, 福岛附近海域表层海洋物质的传播路径可大致分为向东和向南两个分支, 这与Chen等(2021)和林武辉等(2021)的路径分析结果大体一致。东支速度较快且为主要传播路径, 浮标从福岛以东海域开始运动后沿黑潮延伸体迅速东传, 随后沿北太平洋副热带环流运动, 大约5个月后可越过日界线、22个月后可到达北美西岸, 其中大部分浮标沿加利福尼亚流往南进入北赤道流后向西边界输运, 少部分沿阿拉斯加流往北抵达阿留申岛链。

南支运动速度较慢、路径复杂曲折, 且存在明显涡旋特征, 但距离中国近海近, 只需5个月后即可抵达吕宋海峡, 进入中国近海。福岛以东的大范围海表南向流, 成为该海域表层海洋物质穿越黑潮直接向西南方向运移并抵达吕宋海峡和中国台湾以东海域的重要通道。在穿越黑潮延伸体前后, 南支浮标路径很弯曲、复杂, 说明黑潮以东存在的强烈的中尺度涡旋活动可能具有重要作用。值得注意的是, 尽管南支浮标数相对东支较少, 但是由于南支路程较短且直抵西太平洋沿岸各国, 因此更加值得关注。

由于漂流浮标到达陆架附近后生命周期会迅速结束, 较难体现在长截止时间模拟试验中。但是, 基于从福岛以东海域出发的浮标轨迹分析和连续排放观测模拟实验结果, 浮标自福岛以东海域开始运动5个月后即可到达吕宋海峡和中国台湾以东海域(图2c), 其概率密度分布模态呈现以福岛以东海域为核心、往正东和西南方向延伸和衰减。从吕宋海峡开始运动的浮标, 有大约21%会进入我国南海和台湾海峡等近海海域(图5b), 其概率密度分布模态呈现以吕宋海峡为核心, 分别沿黑潮和南海贯穿流往东北-西南方向递减。

4.2 意义与展望

尽管研究结果仍存在一定缺陷和不确定性, 但是利用历史浮标观测数据开展分析具有其独特的优势和重要的意义。

一方面, 相较于前人研究中大量使用的数值模拟数据, 浮标观测数据具有扎实的观测基础, 反映了真实海洋中发生的表层刚性粒子扩散过程。海表浮标在散播过程中自然且完整地受到了真实海洋中各种宏观时空尺度的物理过程的影响, 相比之下, 数值模拟的海洋流场与实际的海洋流场存在各种差异和不确定性, 许多物理过程甚至无法模拟。因此, 浮标的观测基础具有不可替代的优势。

另一方面, 该研究使用了目前已知的最大规模的浮标观测数据集, 代表了目前在该海域的最丰富的海洋表层浮标轨迹观测, 因此, 基于此得到的浮标轨迹的统计特征具有重要价值。虽然前人已有少量研究使用了漂流浮标数据来分析2011年福岛核事故的核素传播路径(Buesseler, 2012; Rypina, 2013), 但是这些研究使用的漂流浮标数据存在布放时间单一、浮标观测时间序列过短(不足4个月)等问题, 不足以揭示海洋物质的长期散播过程。文章选择的浮标数远多于前人研究, 且所用浮标的生命周期相对较长, 能够基本涵盖完整的表层物质散播过程, 因此具有重要的意义。

海洋物质循环是海洋环境安全的重要内容。当前背景下, 海洋将长期面临人类活动污染带来的威胁。未来, 有必要开展如下三个方面工作: (1) 有计划地开展长期、大范围、多参数的海洋观测。目前现场观测资料非常缺乏, 导致对相关海域的水文动力过程(环流、潜沉、中尺度涡输运等过程)的具体影响、化学属性的特征和演变规律的认识不清楚; (2) 进一步提升高分辨率数值模式对物质循环和扩散的模拟水平, 将现场观测与高分辨率海洋数值模拟充分结合起来, 增强数值模式对西太平洋物质扩散的模拟能力; (3) 研究海洋物质扩散过程对海洋环境和生态系统的影响。

刘广山, 2015. 日本福岛核电站事故后的海洋放射化学[J]. 核化学与放射化学, 37(5): 341-354.

许丽晓, 刘秦玉, 2021. 海洋涡旋在模态水形成与输运中的作用[J]. 地球科学进展, 36(9): 883-898.

李小龙, 徐春艳, 刘新华, 等, 2016. 福岛事故废液滞留和循环处理进程分析[J]. 核科学与工程, 36(5): 701-708.

林武辉, 何建华, 余克服, 等, 2020. 海洋中90Sr: 日本周边海域与南海的对比[J]. 海洋学报, 42(10): 47-58.

林武辉, 余克服, 杜金秋, 等, 2021. 日本福岛核废水排海情景下海洋生态环境影响与应对[J]. 科学通报, 65(55): 4500-4509.

赵云霞, 韩磊, 曲大鹏, 等, 2017. 福岛核事故释放的137Cs在西北太平洋海水中输运的模拟与预测[J]. 海洋科学进展, 35(2): 221-233.

AOYAMA M, HAMAJIMA Y, HULT M,, 2016.134Cs and137Cs in the North Pacific Ocean derived from the March 2011 TEPCO Fukushima Dai-ichi Nuclear Power Plant accident, Japan. Part one: surface pathway and vertical distributions [J]. Journal of Oceanography, 72(1): 53-65, doi: 10.1007/s10872-015-0335-z.

AOYAMA M, KAJINO M, TANAKA T Y,, 2015.134Cs and137Cs in the North Pacific Ocean derived from the March 2011 TEPCO Fukushima Dai-ichi Nuclear Power Plant accident, Japan. Part two: estimation of134Cs and137Cs inventories in the North Pacific Ocean [J]. Journal of Oceanography, 72(1): 67-76, doi: 10.1007/s10872-015- 0332-2.

AOYAMA M, TSUMUNE D, HAMAJIMA Y, 2013a. Distribution of137Cs and134Cs in the North Pacific Ocean: impacts of the TEPCO Fukushima-Daiichi NPP accident [J]. Journal of Radioanalytical and Nuclear Chemistry, 296(1): 535-539, doi: 10.1007/s10967-012-2033-2.

AOYAMA M, UEMATSU M, TSUMUNE D,, 2013b. Surface pathway of radioactive plume of TEPCO Fukushima NPP1 released134Cs and137Cs [J]. Biogeosciences, 10(5): 3067-3078, doi: 10.5194/bg-10-3067-2013.

BEHRENS E, SCHWARZKOPF F U, LÜBBECKE J F,, 2012. Model simulations on the long-term dispersal of137Cs released into the Pacific Ocean off Fukushima [J]. Environmental Research Letters, 7(3): 034004.

BUDYANSKY M V, GORYACHEV V A, KAPLUNENKO D D,, 2015. Role of mesoscale eddies in transport of Fukushima-derived cesium isotopes in the ocean [J]. Deep Sea Research Part I: Oceanographic Research Papers, 96: 15-27, doi: 10.1016/j.dsr.2014.09.007.

BUESSELER K, 2014. Fukushima and ocean radioactivity [J]. Oceanography, 27(1): 92-105, doi: 10.5670/oceanog.2014.02.

BUESSELER K, AOYAMA M, FUKASAWA M, 2011. Impacts of the Fukushima nuclear power plants on marine radioactivity [J]. Environmental Science & Technology, 45(23): 9931-9935, doi: 10.1021/es202816c.

BUESSELER K O, JAYNE S R, FISHER N S,, 2012. Fukushima-derived radionuclides in the ocean and biota off Japan [J]. Proceedings of the National Academy of Sciences of the United States of America, 109(16): 5984-5988, doi: 10.1073/pnas.1120794109.

CENTURIONI L R, TURTON J, LUMPKIN R,, 2019. Globalobservations of essential climate and ocean variables at the air–sea interface [J]. Frontiers in Marine Science, 6: 419, doi: 10.3389/fmars.2019.00419.

CHEN G X, HOU Y J, CHU X Q, 2011b. Mesoscale eddies in the South China Sea: Mean properties, spatiotemporal variability, and impact on thermohaline structure [J]. Journal of Geophysical Research, 116(C6): C06018, doi: 10.1029/ 2010jc006716.

CHEN G X, HU P, HOU Y J,, 2011a. Intrusion of the Kuroshio into the South China Sea, in September 2008 [J]. Journal of Oceanography, 67(4): 439-448, doi: 10.1007/ s10872-011-0047-y.

CHEN G X, WANG Q, CHU X Q, 2021. Accelerated spread of Fukushima’s waste water by ocean circulation [J]. The Innovation, 2(2): 100119, doi: 10.1016/j.xinn.2021.100119.

CHENG Y H, HO C R, ZHENG Q A,, 2014. Statistical Characteristics of Mesoscale Eddies in the North Pacific Derived from Satellite Altimetry [J]. Remote Sensing, 6(6): 5164-5183, doi: 10.3390/rs6065164.

DENG F F, LIN F, YU W,, 2020. The distributions of134Cs,137Cs and90Sr in the northwest Pacific seawater in the winter of 2012 [J]. Marine Pollution Bulletin, 152: 110900, doi: 10.1016/j.marpolbul.2020.110900.

FUKUDA M, AONO T, YAMAZAKI S,, 2017. Dissolved radiocaesium in seawater off the coast of Fukushima during 2013-2015 [J]. Journal of Radioanalytical and Nuclear Chemistry, 311(2): 1479-1484, doi: 10.1007/s10967-016- 5009-9.

GARRAFFO Z D, KIM H C, MEHRA A,, 2016. Modeling of137Cs as a tracer in a regional model for the Western Pacific, after the Fukushima–Daiichi Nuclear Power Plant accident of March 2011 [J]. Weather and Forecasting, 31(2): 553-579, doi: 10.1175/waf-d-13-00101.1.

HU S J, HU D X, 2014. Variability of the Pacific North equatorial current from repeated shipboard acoustic Doppler current profiler measurements [J]. Journal of Oceanography, 70(6): 559-571, doi: 10.1007/s10872-014-0253-5.

HU S J, HU D X, GUAN C,, 2016. Interannual variability of the mindanao current/undercurrent in direct observations and numerical simulations [J]. Journal of Physical Oceanography, 46(2): 483-499, doi: 10.1175/jpo-d-15- 0092.1.

HU S J, LU X, LI S H,, 2021. Multi-decadal trends in the tropical Pacific western boundary currents retrieved from historical hydrological observations [J]. Science China Earth Sciences, 64(4): 600-610, doi: 10.1007/s11430-020-9703-4.

HU S J, SPRINTALL J, 2017. Observed Strengthening of Interbasin Exchange via the Indonesian Seas Due to Rainfall Intensification [J]. Geophysical Research Letters, 44(3): 1448-1456, doi: 10.1002/2016GL072494.

HU S J, SPRINTALL J, GUAN C,, 2020a. Deep-reaching acceleration of global mean ocean circulation over the past two decades [J]. Science Advances, 6(6): eaax7727, doi: 10.1126/sciadv.aax7727.

HU S J, SPRINTALL J, GUAN C,, 2020b. Observed triple mode of salinity variability in the thermocline of tropical Pacific Ocean [J]. Journal of Geophysical Research: Oceans, 125(9): e2020JC016210, doi: 10.1029/2020JC016210.

HU D X, WU L X, CAI W J,, 2015. Pacific western boundary currents and their roles in climate [J]. Nature, 522(7556): 299-308, doi: 10.1038/nature14504.

KAERIYAMA H, AMBE D, SHIMIZU Y,, 2013. Direct observation of134Cs and137Cs in surface seawater in the western and central North Pacific after the Fukushima Dai-ichi nuclear power plant accident [J]. Biogeosciences, 10(6): 4287-4295, doi: 10.5194/bg-10-4287-2013.

KAMENíK J, DULAIOVA H, BUESSELER K O,, 2013. Cesium-134 and 137 activities in the central North Pacific Ocean after the Fukushima Dai-ichi Nuclear Power Plant accident [J]. Biogeosciences, 10(9): 6045-6052, doi: 10.5194/bg-10-6045-2013.

KAWAMURA H, KOBAYASHI T, FURUNO A,, 2011. Preliminary numerical experiments on oceanic dispersion of131I and137Cs discharged into the ocean because of the Fukushima Daiichi nuclear power plant disaster [J]. Journal of Nuclear Science and Technology, 48(11): 1349-1356, doi: 10.1080/18811248.2011.9711826.

KUMAMOTO Y, AOYAMA M, HAMAJIMA Y,, 2014. Southward spreading of the Fukushima-derived radiocesium across the Kuroshio Extension in the North Pacific [J]. Scientific Reports, 4: 4276, doi: 10.1038/srep04276.

LEON J D, JAFFE D A, KASPAR J,, 2011. Arrival time and magnitude of airborne fission products from the Fukushima, Japan, reactor incident as measured in Seattle, WA, USA [J]. Journal of Environmental Radioactivity, 102(11): 1032-1038, doi: 10.1016/j.jenvrad.2011.06.005.

LIU Q Y, HU H B, 2007. A subsurface pathway for low potential vorticity transport from the central North Pacific toward Taiwan Island [J]. Geophysical Research Letters, 34(12): L12710, doi: 10.1029/2007gl029510.

MACDONALD A M, YOSHIDA S, PIKE S M,, 2020. A Fukushima tracer perspective on four years of North Pacific mode water evolution [J]. Deep Sea Research Part I: Oceanographic Research Papers, 166: 103379, doi: 10.1016/j.dsr.2020.103379.

MADERICH V, JUNG K T, BEZHENAR R,, 2014. Dispersion and fate of90Sr in the Northwestern Pacific and adjacent seas: global fallout and the Fukushima Dai-ichi accident [J]. Science of the Total Environment, 494-495: 261-271, doi: 10.1016/j.scitotenv.2014.06.136.

MORINO Y, OHARA T, NISHIZAWA M, 2011. Atmospheric behavior, deposition, and budget of radioactive materials from the Fukushima Daiichi nuclear power plant in March 2011 [J]. Geophysical Research Letters, 38(7): L00G11, doi: 10.1029/2011gl048689.

NAKANO M, POVINEC P P, 2012. Long-term simulations of the137Cs dispersion from the Fukushima accident in the world ocean [J]. Journal of Environmental Radioactivity, 111: 109-115, doi: 10.1016/j.jenvrad.2011.12.001.

PERIÁÑEZ R, BROVCHENKO I, DUFFA C,, 2015. A new comparison of marine dispersion model performances for Fukushima Dai-ichi releases in the frame of IAEA MODARIA program [J]. Journal of Environmental Radioactivity, 150: 247-269, doi: 10.1016/j.jenvrad.2015. 09.003.

POVINEC P P, AOYAMA M, BIDDULPH D,, 2013. Cesium, iodine and tritium in NW Pacific waters – a comparison of the Fukushima impact with global fallout [J]. Biogeosciences, 10(8): 5481-5496, doi: 10.5194/bg-10- 5481-2013.

ROEMMICH D, GILSON J, 2009. The 2004-2008 mean and annual cycle of temperature, salinity, and steric height in the global ocean from the Argo Program [J]. Progress in Oceanography, 82(2): 81-100, doi: 10.1016/j.pocean.2009. 03.004.

RONG Y N, LIANG X S, 2018. A Study of the impact of the Fukushima nuclear leak on east China coastal regions [J]. Atmosphere-Ocean, 56(4): 254-267, doi: 10.1080/07055900. 2017.1421139.

ROSSI V, VAN SEBILLE E, SEN GUPTA A,, 2013. Multi-decadal projections of surface and interior pathways of the Fukushima Cesium-137 radioactive plume [J]. Deep Sea Research Part I: Oceanographic Research Papers, 80: 37-46, doi: 10.1016/j.dsr.2013.05.015.

RYPINA I I, JAYNE S R, YOSHIDA S,, 2013. Short-term dispersal of Fukushima-derived radionuclides off Japan: modeling efforts and model-data intercomparison [J]. Biogeosciences, 10(7): 4973-4990, doi: 10.5194/bg-10-4973-2013.

SMITH J N, BROWN R M, WILLIAMS W J,, 2015. Arrival of the Fukushima radioactivity plume in North American continental waters [J]. Proceedings of the National Academy of Sciences of the United States of America, 112(5): 1310-1315, doi: 10.1073/pnas.1412814112.

SMITH J N, ROSSI V, BUESSELER K O,, 2017. Recent transport history of Fukushima radioactivity in the northeast pacific ocean [J]. Environmental Science & Technology, 51(18): 10494-10502, doi: 10.1021/acs.est.7b02712.

STOHL A, SEIBERT P, WOTAWA G,, 2012. Xenon-133 and caesium-137 releases into the atmosphere from the Fukushima Dai-ichi nuclear power plant: determination of the source term, atmospheric dispersion, and deposition [J]. Atmospheric Chemistry and Physics, 12(5): 2313-2343, doi: 10.5194/acp-12-2313-2012.

TSUMUNE D, TSUBONO T, AOYAMA M,, 2012. Distribution of oceanic137Cs from the Fukushima Dai-ichi nuclear power plant simulated numerically by a regional ocean model [J]. Journal of Environmental Radioactivity, 111: 100-108, doi: 10.1016/j.jenvrad.2011.10.007.

XU L X, LI P L, XIE S P,, 2016. Observing mesoscale eddy effects on mode-water subduction and transport in the North Pacific [J]. Nature Communications, 7: 10505, doi: 10.1038/ncomms10505.

XU L X, XIE S P, LIU Q Y,, 2017. Evolution of the North Pacific subtropical mode water in anticyclonic eddies [J]. Journal of Geophysical Research: Oceans, 122(12): 10118-10130, doi: 10.1002/2017JC013450.

YOSHIDA N, KANDA J, 2012. Tracking the Fukushima radionuclides [J]. Science, 336(6085): 1115-1116.

ZHAO C, WANG G, QIAO F L,, 2015. A numerical investigation into the long-term behaviors of Fukushima- derived137Cs in the ocean [J]. Acta Oceanologica Sinica, 34(12): 37-43, doi: 10.1007/s13131-015-0775-8.

ZHAO C, WANG G, ZHANG M,, 2021. Transport and dispersion of tritium from the radioactive water of the Fukushima Daiichi nuclear plant [J]. Marine Pollution Bulletin, 169: 112515, doi: 10.1016/j.marpolbul.2021. 112515.

ANALYSIS OF SPREADING TRAJECTORY OF SURFACE OCEAN MATERIALS FROM EAST OF FUKUSHIMA OF JAPAN BASED ON DRIFTING BUOY OBSERVATIONS

LU Xi1, 2, 3, HU Shi-Jian1, 2, 3, 4, WANG Fan1, 2, 3, 4, HU Dun-Xin1, 2, 3, 4, LIU Hai-Long5, 3, DIAO Xin-Yuan1, 2, 3, 4, YU Ren-Cheng1, 2, 3, 4

(1. CAS Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Qingdao 266071, China; 2. Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China; 4. Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China; 5. State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Chinese Academy of Sciences Institute of Atmospheric Physics, Beijing 100029, China)

Knowing how materials disperse at the ocean surface is very important for understanding the process of marine environmental change and anthropogenic pollution diffusion in the ocean. Lagrangian trajectory analysis and observation-based simulation experiments were performed to examine the spreading trajectory of surface materials from the east of Fukushima of Japan using historical sea surface drifting buoys observation data. The surface spreading path of drifting buoys originated from east of Fukushima could be divided into mainly two branches: eastern and southern ones. The eastern branch is the main channel in higher speed and merges into the North Pacific Current along the Kuroshio extension, with an earliest arrival to North America in about 22 months. The southern branch in lower speed migrates to the southwest along a widespread southward currents south of the Kuroshio extension, and are obviously affected by eddy activities. The southern branch reaches the Luzon Strait and east of Taiwan in about 5 months at the earliest, and then enters the South China Sea and the East China Sea. The probability density distribution of surface drifting buoys in simulation experiments based on observation spreads with the Fukushima vicinity as the core, and gradually expands to the southwest and east. The drifting buoys arriving in coastal China are mainly originated from the Luzon Strait. Confined to the specialties of drifting buoy observations, the conclusions have certain limitations from omitting processes such as mixing, diffusion, and biogeochemistry. However, the conclusions based on observation are of great importance for reference.

Fukushima; nuclear pollution; western Pacifica Ocean; diffusion tracks; drifting buoy

P733

10.11693/hyhz20211100292

*中国科学院重点部署项目, ZDBS-SSW-DQC002号; 中国科学院战略性先导科技专项项目, XDB42010403号; 国家自然科学基金项目, 42022040号; 山东省自然科学基金资助项目, ZR2020JQ18号; 中国科学院创新交叉团队项目, JCTD-2020-12号; 中国科学院青年创新促进会会员项目, 2018240号。卢 锡, 博士研究生, E-mail: luxi@qdio.ac.cn

胡石建, 研究员, 博士生导师, E-mail: sjhu@qdio.ac.cn

2021-11-23,

2022-01-11