黑潮源区海流季节内变化观测分析
段静, 陈朝晖, 吴立新
中国海洋大学物理海洋教育部重点实验室,山东 青岛 26003

作者简介:段静(1990-),女,甘肃南皮人,硕士研究生,主要从事海洋环流及气候方面研究.E-mail:duanjingalice@hotmail.com

摘要

通过分析2010年11月至2011年7月黑潮源区的单点潜标观测数据,得出在上层北向黑潮之下400~700 m的南向深层海流存在季节内变化(冬季>夏季>春季)。结合卫星高度计数据、SODA 2.2.4数据集、WOD09数据模式资料对观测结果进行分析,得到深层南向流季节内变化的一个可能影响机制为:观测点气候态背景流与中尺度涡旋运动扰动的叠加。进一步,结合HYCOM模式数据计算源地黑潮流主轴,得到在观测期间潜标放置点同流轴的纬向距离存在季节内变化(冬季>夏季>春季),得出深层南向流季节内变化的另一个可能影响机制为:潜标放置点与源地黑潮主轴的纬向距离会直接影响观测结果,潜标越靠近主轴测得流速越强,相反则越弱。

关键词: 潜标; 吕宋潜流; 季节内变化; 中尺度涡旋; 黑潮流轴
中图分类号:P731.27 文献标志码:A 文章编号:1001-8166(2014)04-0523-08
Study on Intraseasonal Variation of Current at the Source Region of Kuroshio by Analyzing the Buoy Observation Data
Duan Jing, Chen Zhaohui, Wu Lixin
Key Laboratory of Physical Oceanography, Ministry of Education, Ocean University of China, Qingdao 266003, China
Abstract

Using subsurface buoy data during November 2010 to July 2011, the south flow, which is under the northward Kuroshio Current between 400 m to 700 m, was studied. According to the time series available from subsurface buoy, the seasonality of the south flow strength is winter > summer > spring. In addition to the Acoustic Doppler Current Profiler (ADCP) data, satellite altimeter data, the SODA 2.2.4 and the World Ocean Database 2009 (WOD09) were further analyzed and a possible influence mechanism is: the intraseasonal variation (winter>summer>spring) of the south flow is mostly influenced by mesoscale eddies east of Luzon, indicating that the observed flow is a combination of background flow and eddy induced flow. With the consideration of oscillating of the source of the Kuroshio which affected the on observation, we calculated the axis of the source of the Kuroshio based on HYCOM data by using weighted average method. A number of features of the Kuroshio axis were identified. During November 2010 to July 2011, significant intraseasonal variability of the distance between Kuroshio axis and the longitude of buoy (122.7°E) was exposed. The whole layer had influence on the axis from surface to 700 m. Therefore, another possible influence mechanism is: The observed flow is also influenced by the horizontal oscillation of the Kuroshio axis, the observed meridional velocity is higher when the distance is smaller, and vice versa.

Keyword: Buoy; Luzon Undercurrent; Intraseasonal variation; Mesoscale eddy; Kuroshio axis
1 引言

黑潮(Kuroshio)作为北赤道流(North Equatorial Current,NEC)分叉后沿北太平洋西边界北向流动的强大暖流分支[ 1],对西北太平洋及我国近海的海洋生态环境、海洋工程以及国防安全等都有巨大影响,与我国东部天气系统和气候也有密切关联。其中,位于吕宋岛东部的黑潮源头在很大程度上控制着黑潮与南海的能量及水交换,对南海的动力环境有显著影响,并控制着东海黑潮的强弱[ 2]。研究表明,黑潮源头的流速分布以及流量强弱是黑潮在吕宋海峡形变的关键因素[ 3],因此对黑潮源区海流的直接观测能够为认识黑潮源头的时空变化特征提供直接的依据。

由于黑潮与之下的南向潜流关系紧密,因此,近年来,该潜流也得到了越来越多的关注。Hu等[ 4]利用1986-1988年间70个盐度/电导率—温度—深度剖面测量仪(Conductivity-Temperature-Depth, CTD)数据,发现黑潮源地可能存在流速为4 cm/s、流量约为2 Sv(1 Sv=1x106 m3/s)的潜流。Qu等[ 5]进一步采用1986到1991年的水文资料,确认潜流的存在并将其命名为吕宋潜流(Luzon Undercurrent)。他们发现吕宋潜流的平均地转输运量为3.6 Sv且28%的水体输运来源于低盐的北太平洋中层水;吕宋潜流最大速度位于700 m,大小约为7 cm/s。目前,对于吕宋潜流的研究多关注于其地转特征及源地[ 6, 7, 8, 9, 10],而对于它的直接观测仍非常缺乏。

2010年,国际组织CLIVAR/WCRP (Climate Variability and Predictability/World Climate Research Programme) 批准的西北太平洋海洋环流与气候实验(Northwestern Pacific Ocean Circulation and Climate Experiment, NPOCE)计划[ 11],首次利用次表层潜标与声学多普勒测流计(Acoustic Doppler Current Profilers, ADCP)对黑潮源区海流进行了长时间观测。Hu等[ 12]利用该次观测资料,通过分析2010年11月至2011年7月期间的海流数据,发现吕宋潜流存在有季节内变化:冬季 > 夏季 > 春季,并且是由西向传播的中尺度涡旋运动导致,同时指出该变化存在70~80天的周期。

本文采用与 Hu等[ 12]相同的潜标观测海流数据,进一步结合卫星高度计数据、SODA (Simple Ocean Data Assimilation)2.2.4、WOD09数据以及HYCOM数据资料,除了考虑中尺度涡旋的作用,还引入了对源地黑潮流轴摆动影响的分析,阐明了源地黑潮流轴的摆动亦会影响潜标观测结果的观点,从而进一步补充和完善了Hu等的结论。

2 数据及方法介绍
2.1 潜标观测数据

本文使用的观测数据来自NPOCE项目投放于18°1.491′N/122°38.044′E位置处的次表层潜标。潜标所装备的相关仪器的参数介绍以及初步数据处理,可以参考文献[12]。

2.2 SODA 2.2.4海洋环流同化数据和WOD09数据

本文选取了1960—2008年SODA 2.2.4版本[ 13]的水平流速( , )数据。SODA 2.2.4提供超过一百年时间尺度的同化资料,为长期的气候研究提供了丰富可靠的海洋数据(详细介绍可参考http://iridl.ldeo.columbia.edu/SOURCES/.CARTON-GIESE/SODA/.v2p0p2-4/)。针对本文对黑潮源区气候态海流流动特征的总结需求,SODA 2.2.4数据集为潜标数据研究提供了气候态背景场。

此外,本文使用了美国国家海洋和大气管理局(National Oceanic and Atmospheric Administration, NOAA)海洋气候实验室(Ocean Climate Laboratory, OCL)在大量真实实测和历史数据基础上出产的第三个数据集和分析产品WOD09(World Ocean Database 2009)(详细介绍可参考http://www.nodc.noaa.gov/OC5/WOD/datageo.html)。本文采用1960—2012年0°~30°N、120°~160°E的温度、盐度资料,并使用Cressman插值[ 14, 15, 16]和热成风关系及动力计算的处理方法[ 17],将WOD09的温盐数据转换为三维的地转流速数据,其中地转流计算的参考面选择为1500m。

2.3 卫星高度计资料和HYCOM数据

本文使用的卫星高度计资料来源于AVISO(Archivings Validation and Interpretation of Satellite Oceanographic Data)提供的网格化高度计数据。该数据融合了Jason-1、Topex/Poseidon、Envisat、ERS-1/2、GEOSAT等卫星资料。本文选用的海表面高度异常(Sea Level Anomaly,SLA)及相应的地转流经向、纬向流速异常( )数据的空间范围为120 °~130°E、15°~25°N,时间范围为1992年11月至2012年12月。由于我们主要关注季节内时间尺度以上的变化,因此对SLA、相应地转流流速异常数据进行了3周低通滤波处理以剔除其高频变化。

此外,美国海军研究实验室(Naval Research Laboratory,NRL)的混合坐标模式(HYbird Coordinate Ocean Model,HYCOM)数据,它是具有涡分辨能力的高精度资料(具体介绍可以借鉴本文参考文献[9])。本文利用观测点处2010年11月至2011年7月期间上层700 m的经向流速分量( )数据,分析了源地黑潮流轴的摆动及其垂向的变化。

3 潜标观测结果分析

通过对潜标观测数据的处理,本文得到与Hu等[ 12]一致的结果。结果显示,仅在2011年3月20日至2011年5月20日期间, 200~700 m的潜标观测经向流速 ( ) 表现为一致北向,而其下方的南向流动消失。基于对观测和模式数据的分析,Hu等[ 12]表明,观测中下层南向流动为吕宋潜流且流速有明显的季节内变化(冬季 > 夏季 > 春季)。然而,由于缺乏长期连续的观测,黑潮源区流动的平均态无法从较短的观测序列中直接获得,对吕宋潜流的可能季节内变化的刻画无法从观测获得充足的支持,因此本文将结合多种资料分析总结黑潮源区流动的气候态特征以辅助对潜标观测结果的分析。

为了解黑潮源区海流的气候态特征,需对20年以上的长时间海洋数据进行研究,尽管HYCOM非同化数据集提供有1979—2003年的高分辨率北太平洋海洋流速数据,但非同化数据相对于同化数据集其模拟的准确性更低。因此,本文利用SODA 2.2.4流速数据,描绘了黑潮源区附近的水平流速分布及观测点处 的垂向分布(图1)。其中,18°N断面气候态平均的经向流速分量 ( ) 曲线表明:在上层较强的北向流动之下存在很弱的南向潜流,且潜流流核位于700 m深处(图1d)。由气候态不同深度的水平流场分布可知,NEC分叉点随着深度的增加不断向北移动会导致在黑潮源区的深层存在南向流动(图1a,b, c)。同时,WOD09与SODA 2.2.4资料的结果也都显示随着深度的增加经向流速零线(NEC分叉点)位置不断北偏(图1e, f)。因此,基于上述分析我们认同潜标在600 m深度处观测到的南向流动为吕宋潜流的结论。然而,由于缺乏长期观测和多点观测的辅助,单点潜标的观测环境将影响观测结果,因此本文认为观测中深层南向流动呈现出的存在-消失-存在的变化特征还不能确定为吕宋潜流的季节内变化,针对该变化的可能影响机制将在下章探讨。

图1 黑潮源区气候态海流结构
黑潮源地附近 (a) 5 m, (b) 200 m, (c) 600 m 的水平流场( us, vs),星号代表NEC分叉点的位置(d) 气候态18°N断面 vs分量,实线代表北向流速,间隔为5 cm/s;虚线为南向流,间隔为0.8 cm/s;菲律宾以外(东)5°平均 vs分量(e) WOD09, (f) SODA 2.2.4
Fig.1 Climatology current structure in the source region of KuroshioSODA
2.2.4 current velocity at (a) 5 m, (b) 200 m, (c) 600 m, the asterisk is the bifurcation location of NEC; (d) Climatology vs at 18°N; solid linenorthward current, counter interval 5 cm/s; Dash line-southward current, counter interval 0.8 cm/s; vs averaged over 5° east of Philippines (e) from WOD09, (f) from SODA 2.2.4

4 影响机制讨论
4.1 源地黑潮区域涡旋运动的影响

西北太平洋是全球大洋涡旋运动最丰富的地区之一,涡旋的影响深度可至深层[ 18]。为了研究涡旋运动对观测的影响,本文对 逐层距平,结果表明:观测流速整层变化一致,且存在约60天的周期(图2a);经过20~90天带通滤波后该周期更为显著(图2b)。结合1400 m处海流计观测的流矢量变化(图2c),分析可知潜标观测经向流速分量整层受到一致扰动的影响,且影响深度可达1400 m。进一步,比较200~700 m垂向平均 的矩平值与 可以发现: 异常同卫星观测 的变化一致(图3a)。进而,同潜标观测位置处 异常和 所反映的气候态月平均季节变化进行对比(图3b)可知: 异常和 同潜标观测反映的流速异常变化一致,三者均在3~5月表现为较强的北向流动异常,而7月为较强的南向流动异常。同时,分析不同观测月份 水平结构(图3d)可知:2010年12月至2011年7月观测位置处有强反气旋环流异常,2011年3~5月则几乎没有反气旋环流异常甚至表现为气旋式环流异常。鉴于SLA可以反映涡旋运动,因此推断观测所反映的周期为60天的变化可能是由吕宋岛以东的涡旋运动引起的。

图2 潜标观测海流流速
(a) vb逐层异常值;(b)20~90天带通滤波 vb异常值;(c) 1 400 m海流计观测流失量
Fig.2 Current velocity of subsurface buoy observation
(a) vb anomaly at different depth; (b) vb anomaly with a band pass filter(20~90 days); (c) Buoy measured velocity vectors at 1 400 m

图3 观测点流速季节变化及水平分布
(a) 观测点处200~700 m垂向平均的 vb异常(红线)与 vSLA(黑线);(b)1960—2008年观测点处气候态 vs异常的季节变化;(c) 1992—2012年观测点处气候态平均 vSLA′的季节变化;(d) 2010年12月至2011年7月期间 uSLA′和 vSLA′在120 °~130°E,15°~25°N的水平分布;黑色圆点代表潜标放置位置
Fig.3 Seasonal variation and horizontal distribution of current velocity in the source region of Kuroshio
(a) Anomaly of vertically averaged (200~700 m) vb (red) and vSLA (black) at the research site; (b) Intraseasonal variation of climatology v s anomaly near the research site during 1960-2008; (c) Intraseasonal variation of vSLAnear the research site during 1992-2012; (d) uSLA′ and vSLA field of 120 °~130°E,15°~25°N during December 2010 to July 2011; Black dot is the location of the subsurface buoy

进一步,结合观测点气候态平均背景流的垂向分布特征,影响观测中南向逆流变化的可能形成机制为:观测点处气候态平均背景流在0~600 m表现为较强的北向海流,深层为较弱南向流动(图4a),当涡旋运动引起北向流动异常扰动时,上层北向流动进一步增加,深层南向流减弱(潜标观测2011年3~5月);而当引起南向流动异常扰动时,上层北向流动减弱,深层南向流则加强(潜标观测2010年12月)。从而,观测点处气候态平均背景流与涡旋扰动线性叠加可以由示意图 (图4b)表示。

图4 (a) 观测点处气候态 随深度的变化。(b) 观测点处经向流速受到涡旋影响机制示意图。Fig.4 (a) Vertical structure of at the research site; (b) A schematic diagram showing the mechanism of eddy influenced meridional current velocity structure.

因此,观测点处背景流与涡旋运动扰动的叠加会使深层南向流呈现出在3~5月流速加强,而在7月流速减小甚至反向的特点,很好的对应了潜标观测的结果。然而,由于黑潮是一支北向流动的强流,其流轴的摆动也会影响单点潜标的观测,因而,我们将在下一节通过描述黑潮流轴摆动来进一步完善对潜标观测的分析。

4.2 源地黑潮流主轴的变化及影响

鉴于黑潮源区海流与观测仪器的空间位置关系会影响潜标观测,为进一步分析观测结果,本节将对黑潮流轴摆动进行分析。由于源地黑潮紧邻菲律宾岛且纬向范围较小,因此本文选用了1/12° HYCOM资料提供的2010年11月-2011年6月经向流速分量数据来分析源地黑潮流轴。本文定义源地黑潮流轴为:观测点以东2度内经向流速分量与经度的加权平均:

(1)

其中, 代表逐点的经向流速、 代表逐点经度。

图5a表明,上200 m平均和200~700 m平均的源地黑潮流轴在2011年3~5月(MAM)期间最靠近潜标所在经度(122.7°E)。由于前文中潜标观测主要集中在200 m之下,因此本节以200~700 m平均为主要分析对象,得到源地黑潮流轴与潜标观测点之间的纬向距离存在季节内变化:冬季 > 夏季 > 春季。由于黑潮流轴摆动与潜标观测所呈现的深层南向逆流都存在季节内时间尺度上的变化,因此在相同时间尺度下流轴摆动可以直接影响到潜标观测。进一步结合潜标观测结果进行定量分析可知,在2011年3~5月(MAM)流主轴位置距潜标观测位置的距离平均值为59.41 km,相应的经向流速观测值平均为2.136 m/s;而2010年12月平均距离为156.9 km,相应的经向流速平均为-0.186 m/s。因此,观测点处200 m-700 m垂向平均流轴与垂向平均 存在对应关系,但是流轴摆动的影响深度和逐层影响力度也与单点观测结果密切相关,因而通过计算刻画有源地黑潮流轴随深度的变化(图5b),结果表明:源地黑潮流轴摆动呈现树状,在350 m以浅摆动幅度大而深层幅度较小,并且同一深度最大摆动幅度可达约200 km。从而,根据流轴摆动与观测流速变化之间存在的显著对应关系,其影响机制可总结为:潜标放置点与源地黑潮主轴的相对位置关系会直接影响观测结果,越靠近主轴流速越强;而影响主轴变化的因素除了上文提到的涡旋运动以外,西边界流的强非线性变化也是造成主轴摆动的因素之一,对于这一点我们将在以后的研究中进一步阐述。

图5 HYCOM模式中源地黑潮流轴
(a) 2010年11月至2011年6月期间HYCOM资料中上层700 m源地黑潮流轴的变化;(b) 2010年11月至2011年6月期间HYCOM资料中源地黑潮流轴的垂向分布;虚线代表潜标放置处的经度
Fig.5 Axis of Kuroshio in the source region in HYCOM data
(a) Variations of Kuroshio axis during November 2010 to June 2011 based on HYCOM output in upper 700 m; (b) Vertical structure of Kuroshio axis from November 2010 to June 2011 based on HYCOM output;Dash line is the longitude of the subsurface buoy

综合上述分析,本文认为潜标能够对吕宋潜流进行高质量观测,且观测中南向逆流季节内变化可以由涡旋运动引起,也可由源地黑潮流轴摆动造成。因此,潜标观测呈现的南向逆流季节内变化的两种影响机制为:①涡旋运动与背景流叠加;②源地黑潮流轴摆动在季节内时间尺度上的影响。总结对潜标观测结果的分析,本文认为由于短期的潜标单点观测不能呈现流动的完整三维信息,进而无法明确观测到的是吕宋潜流本身的季节内变化特征。

5 结论与讨论

黑潮源区的海流结构不仅对黑潮与南海的能量及水交换以及南海的动力环境有显著的影响,而且对东海黑潮的强弱有调控作用。本文利用黑潮源地的潜标观测数据、卫星数据、WOD09数据集、SODA 2.2.4再分析数据和HYCOM模式数据对黑潮源区的海流结构进行分析,得到以下结论:

(1) 分析SODA 2.2.4和WOD09资料于18°N断面的气候态平均经向流速分量,结果显示上层600 m以浅为北向强流(黑潮),下层为南向逆流(吕宋潜流)。与潜标观测结果对比可知,观测中的深层南向流动是吕宋潜流,这与Hu等[ 12]的观点达成共识。

(2) 2011年3~5月观测中深层南向逆流消失,表明南向流动存在季节内变化:冬季 > 夏季 > 春季。对潜标观测结果逐层计算距平并结合1400 m处的海流计数据,证实观测整层受到统一的扰动影响。进一步,比较垂向平均 距平值和 ,表明二者变化一致且存在约60天的显著周期。同时,分析 水平结构,可知南向流动季节内变化可能是由黑潮源地的涡旋运动引起。不同于Hu等[ 12]的是,根据观测地气候态背景流与扰动季节变化的叠加示意图,本文认为观测中深层南向流动的季节内变化不一定是吕宋潜流本身的季节内变化特征,而是涡旋运动在背景流下的一种表现。

(3) 利用HYCOM资料计算源地黑潮流主轴,发现流轴有正压特点,其摆动存在季节内变化:冬季 > 夏季 > 春季,因此认为流轴摆动可以影响观测结果。结合观测流速定量分析可知,流轴摆动同潜标观测的经向流速分量变化一致,流轴与观测点纬向距离的大小与观测值成反比,距离越小观测值越大。另外,分析源地黑潮流轴的垂向分布,得到树状结构,证明了黑潮源区上层700 m 的流轴摆动都会对潜标产生影响。因此,需要有更长时间的观测序列才能剔除由于涡旋以及强非线性导致的流轴摆动对单点潜标观测带来的干扰,从而描绘源地黑潮流动的多时间尺度特征。

The authors have declared that no competing interests exist.

参考文献
[1] Wu Lixin, Chen Zhaohui. Progresses and challenges in observational studies of physical oceanography[J]. Advances in Earth Science, 2013, 28(5): 542-551.
[吴立新, 陈朝晖. 物理海洋观测研究的进展与挑战[J]. 地球科学进展, 2013, 28(5): 542-551. ] [本文引用:1]
[2] Wu Lixin, Liu Qinyu, Hu Dunxin, et al. Variability of the subtropical gyre in North Pacific and its impacts on dynamic environment of China marginal seas[J]. Advances in Earth Science, 2007, 22(12): 1224-1230.
[吴立新, 刘秦玉, 胡敦欣, . 北太平洋副热带环流变异及其对我国近海动力环境的影响[J]. 地球科学进展, 2007, 22(12): 1224-1230. ] [本文引用:1] [CJCR: 1.388]
[3] Yuan D L, Han W Q, Hu D X. Surface Kuroshio path in the Luzon Strait area derived from satellite remote sensing data[J]. Journal of Geophysical Research, 2006, 111: C11007, doi: 10.1029/2005JC003412. [本文引用:1]
[4] Hu D, Cui M. The western boundary current in the far-western Pacific Ocean[C]∥Picaut J, Lukas R, Delcroix T, eds. Proceedings of Western International Meeting and Workship on TOGA COARE. New Caledonia, 1989: 123-134. [本文引用:1]
[5] Qu T, Kagimoto T, Yamagata T. A subsurface countercurrent along the east coast of Luzon[J]. Deep Sea Research, 1997, 44: 413-423. [本文引用:1]
[6] Qu T, Mitsudera H, Yamagata T. On the western boundary currents in the Philippine Sea[J]. Journal of Geophysical Research, 1998, 103(C4): 7537-7548. [本文引用:1]
[7] Qu T, Mitsudera H, Yamagata T. A climatology of the circulation and water mass distribution near the Philippine Coast[J]. Journal of Physical Oceanography, 1999, 29(7): 1488-1505. [本文引用:1] [JCR: 3.179]
[8] Xie L, Tian J, Hu D, et al. A quasi-synoptic interpretation of water mass distribution and circulation in the western North Pacific II: Circulation[J]. Chinese Journal of Oceanology Limnology, 2009, 27(4): 955-965. [本文引用:1] [JCR: 0.577] [CJCR: 0.3746]
[9] Gao S, Qu T, Hu D. Origin and pathway of the Luzon undercurrent identified by a simulated adjoint tracer[J]. Journal of Geophysical Research, 2012, 117(C5): C05011, doi: 10.1029/2011JC007748. [本文引用:1]
[10] Wang Q, Hu D. Origin of the Luzon undercurrent[J]. Bulletin of Marine Science, 2012, 88(1): 51-60. [本文引用:1] [JCR: 1.297]
[11] Hu D, Wang F, Wu L, et al. Northwestern Pacific Ocean Circulation and Climate Experiment (NPOCE) Science/Implementation Plan[M]. Beijing: China Ocean Press, 2011. [本文引用:1]
[12] Hu D X, Hu S, Wu L, et al. Direct measurements of the Luzon undercurrent[J]. Journal of Physical Oceanography, 2013, 43: 1417-1425. [本文引用:6] [JCR: 3.179]
[13] Carton J A, Giese B S. A reanalysis of ocean climate using Simple Ocean Data Assimilation (SODA)[J]. American Meteorological Society, 2008, 136: 2999-3017. [本文引用:1] [JCR: 6.591]
[14] Barnes S L. A technique for maximizing details in numerical weather map analysis[J]. Jounal of Applied Meteorology, 1964, 3: 396-409. [本文引用:1]
[15] Cressman G P. An operational objective analysis system[J]. Monthly Weather Review, 1959, 87: 367-374. [本文引用:1] [JCR: 2.758]
[16] Gand in L S. The Objective Analysis of meteorological fields. translated from the Russian. Jerusalem (Israel Program for Scientific Translations), 1965. Pp. vi, 242: 53 Figures; 28 Tables[J]. Quarterly Journal of the Royal Meteorological Society, 1996, 92(393): 447. [本文引用:1] [JCR: 3.327]
[17] Wang Guihua, Li Rongfeng. Advance in the study of deducing ocean circulation from hydrographic data[J]. Advances in Earth Science, 2004, 19(1): 100-106.
[王桂华, 李荣凤. 利用水文资料推测海洋流场的研究进展[J]. 地球科学进展, 2004, 19(1): 100-106. ] [本文引用:1] [CJCR: 1.388]
[18] Adams D K, McGillicuddy Jr D J, Luis Zamudio, et al. Surface-generated mesoscale Eddies transport deep-sea products from hydrothermal vents[J]. Science, 2011, 332(6029): 580-583. [本文引用:1]