海洋学报

西太平洋黑洞涡旋物质输运能力分析

PDF下载

田丰林 ¹^,² , 杨笑坤 ¹ , 刘潇 ¹ , 何遒 ¹ , 陈戈 ¹^,²^,^*

海洋学报 | 论文 2021,43(12): 1-14

收起

海洋学报 | 论文 2021, 43(12): 1-14

西太平洋黑洞涡旋物质输运能力分析

全屏

田丰林¹^,², 杨笑坤¹, 刘潇¹, 何遒¹, 陈戈¹^,²^,^*

作者信息

¹中国海洋大学信息科学与工程学部海洋技术学院，山东青岛 266100

²青岛海洋科学与技术试点国家实验室区域海洋动力学与数值模拟功能实验室，山东青岛 266237

田丰林（1978—），男，山东省青岛市人，博士，副教授，主要从事海洋信息可视化与海洋遥感工作。E-mail：tianfenglin@ouc.edu.cn

通讯作者:

陈戈（1965—），男，博士，教授，博士生导师，主要从事海洋遥感与大数据工作。E-mail：gechen@ouc.edu.cn

Analysis of black-hole eddy on material transport in the western Pacific

Fenglin Tian¹^,², Xiaokun Yang¹, Xiao Liu¹, Qiu He¹, Ge Chen¹^,²^,^*

Affiliations

¹Department of Marine Technology, College of Information Science and Engineering, Ocean University of China, Qingdao 266100, China

²Regional Ocean Dynamics and Numerical Simulation Functional Laboratory, Pilot National Laboratory of Marine Science and Technology (Qingdao), Qingdao 266237, China

出版时间: 2021-12-30 doi: 10.12284/hyxb2021163

文章导航

摘要

收起

基于椭圆型拉格朗日拟序结构（elliptic Lagrangian Coherent Structures, eLCSs）可以提取出黑洞涡旋，它具有较强的输运能力并且随流场长时间运动后仍保持连贯，类似于海洋中的“黑洞”，故被称为黑洞涡旋。本文基于地转流速度场数据，针对西太平洋中一个典型的黑洞涡旋（Eddy A）使用eLCSs的方法提取涡旋边界，并进行分析研究。利用海表温度、海表盐度和叶绿素浓度数据分析验证Eddy A水平方向上物质输运的连贯性；使用Argo浮标得到的不同深度的温度、盐度和溶解氧数据对Eddy A垂直方向的物质相干性进行分析验证。通过分析证明，在较长的时间尺度上，黑洞涡旋边界较欧拉涡旋边界的连贯性更强，能更加客观地描述物质输运。

关键词

黑洞涡旋 / 连贯性 / 物质输运 / 三维结构

Abstract

收起

Black-hole eddy can be extracted based on elliptic Lagrangian Coherent Structures (eLCSs), which can transport material and maintain coherent under geostrophic for a long time. It was similar to the black hole in the ocean, so it was called black-hole eddy. In this paper, the boundary of black-hole eddy was extracted based on the data of the geostrophic flow velocity field. Using the method of eLCSs and choosing a targeted eddy (Eddy A) in the western Pacific to analyze. Sea surface temperature, sea surface salinity and chlorophyll concentration data are used to verify that the Eddy A is coherent in horizontal material transport. The temperature, salinity and dissolved oxygen data obtained by Argo in different depths are used to prove the coherence in vertical of Eddy A. The result demonstrates that the boundary of black-hole eddy is more coherent than Eulerian boundary in a long time scale and can describe the transport of material more objectively.

Key words

black-hole eddy / coherence / material transport / three-dimensional structures

引用本文

田丰林, 杨笑坤, 刘潇, 何遒, 陈戈. 西太平洋黑洞涡旋物质输运能力分析. 海洋学报, 2021 , 43 (12) : 1 -14 . DOI: 10.12284/hyxb2021163

Fenglin Tian, Xiaokun Yang, Xiao Liu, Qiu He, Ge Chen. Analysis of black-hole eddy on material transport in the western Pacific[J]. Haiyang Xuebao, 2021 , 43 (12) : 1 -14 . DOI: 10.12284/hyxb2021163

正文

收起

1 引言

收起

中尺度涡旋是一种海洋中广泛存在，空间尺度从数十千米至数百千米不等，寿命从数十天至数百天不等的旋转水体。根据旋转方向，中尺度涡旋可以分为气旋型（北半球为逆时针）和反气旋型（北半球为顺时针）^[1-2]，它的运动影响着海洋环流、大规模水体和生物分布，对盐度、热量和其他海水化学元素的混合、输运也非常重要^[3-6]。

中尺度涡旋的识别方法可分为欧拉方法和拉格朗日方法。欧拉方法根据瞬时海表面高度或海流瞬时特征来识别涡旋，可以利用卫星高度计数据实现自动化识别全球涡旋^[7]。近期的欧拉涡旋识别和追踪方法采用并行计算的方法^{[3, 8-9]}，大大提高了对全球欧拉涡旋的识别效率。欧拉方法的识别结果会随着参考系的旋转或线性加速而变化，使研究结果的不确定性增加^[10]。欧拉方法高估了中尺度涡旋的输运能力^[7]，在非稳定流场中，不依赖于参考系的拉格朗日拟序结构（Lagrangian Coherent Structures，LCSs）方法能更客观地描述涡旋对水体的输运能力^[11]。欧拉涡旋描述短期内物质的输运很有效，但在非稳定流场中，欧拉涡旋边界随流场运动会很快失去连贯性，沿涡旋移动轨迹出现拉伸和成丝的现象；而拉格朗日涡旋可以描述一定时间尺度内涡旋的边界，在时间尺度内，涡旋的边界作为物质屏障可以携带内部的水体运动，并且不会产生大量的泄露^[7]。黑洞涡旋是一种连贯性很强的拉格朗日涡旋，Beron-Vera等^[7]认为黑洞涡旋的定义是涡旋边界所处的物质带符合条件：涡旋在随着流场变形后，边界所处的物质带的变形应该与边界的变形一致，不应该出现显著的拉伸或者剪切变形。黑洞涡旋边界的理论是寻找流场平流输运下具有最小变形的边界，在计算中，求得嵌套的椭圆拉格朗日拟序结构（eLCSs）的最外层闭合曲线，即为黑洞涡旋的边界，由于该边界的数学解和宇宙中黑洞边界的数学解等价，且物质进入涡旋内部就很难出来，故被命名为黑洞涡旋^[12]。与其他拉格朗日涡旋相比，黑洞涡旋的识别过程比较复杂，对涡旋的连贯性要求很高，要求在给定时间尺度内没有任何水体的丢失，相比于其他的拉格朗日涡旋能更好地评估涡旋的输运能力^[13]，它具有很好的包裹性和连贯性^[14]。

在LCSs的早期研究中，有限时间李雅普诺夫指数方法（Finite-time Lyapunov Exponent, FTLE）^[15]和有限空间李雅普诺夫指数方法（Finite-size Lyapunov Exponent, FSLE）^[16-17]被用来评估海洋流场的扰动增长^[18]。虽然FTLE在无剪切LCSs提取中发挥了重要作用并得到了广泛应用，但其连贯性比较差，会在特定情况下产生错误^[19]。因此，一种新的拉格朗日变分测地线理论被提出^[20]，它能更准确地描述物质输运并且可以通过流体变形有效地提取LCSs，这一方法也经过了理论流场的验证^[3]。Haller和Beron-Vera^[12]提出黑洞涡旋的识别方法，并比较了拉格朗日涡旋和欧拉涡旋的物质输运能力，Beron-Vera等^[14]通过二维遥感数据和Drifter数据证明了拉格朗日涡旋确实具有强连贯性。由于黑洞涡旋的识别方法对涡旋的连贯性要求较高并且计算过程比较复杂，Haller等^[21]以及Abernathey和Haller^[22]又提出通过计算拉格朗日平均涡度偏差（Lagrange Averaged Vorticity Deviation, LAVD）的方法识别亚中尺度拉格朗日涡旋。相比于LAVD涡旋，黑洞涡旋的连贯性更强，也能更好地反映拉格朗日涡旋对物质的输运能力。Abernathey和Haller^[22]，Rypina等^[23]都使用地转流数据识别涡旋边界。

中尺度涡旋三维结构的研究主要基于卫星高度计数据和观测数据的综合分析，可以使用遥感数据和浮标数据来评估涡旋对物质的输运能力^{[5, 24]}。Dong等^[25]利用卫星高度计数据和Argo数据分析涡旋在水平和垂直方向上热量和盐度的输运能力，并比较了不同区域涡旋对热量和淡水的输运能力。目前关于欧拉涡旋的三维结构研究较多^[26-27]。Hu等^[28]基于卫星高度计和Argo数据研究一个气旋涡的三维结构和物理性质。Xu等^[26]利用Argo剖面数据观察不同深度上中尺度涡旋对模态水俯冲和输运的影响，并计算了北太平洋涡旋扩散系数。Dong等^[29]和Lin等^[30]在南海和南加州湾采用了三维欧拉涡旋识别方法获得三维涡旋数据集，证明了大多数涡旋的垂直结构不超过400 m。Blazevski和Haller^[31]研究了非稳定流场中eLCSs的三维结构，在二维流场的基础上证明了三维流场中拉格朗日运输边界的稳定性和有效性。总体来说，与欧拉涡旋相比，基于三维观测数据的拉格朗日涡旋研究较少。

为了补充拉格朗日涡旋的研究，验证黑洞涡旋的物质输运能力，本文使用卫星高度计数据和Argo浮标数据来分析涡旋在水平和垂直方向上输运物质的能力。首先，识别出不同时间尺度（T=30 d，60 d，90 d和120 d）的黑洞涡旋边界，同时使用虚拟粒子来验证黑洞涡旋的强连贯性，并与欧拉边界进行对比。然后，将海表盐度（Sea Surface Salinity，SSS）、海表温度（Sea Surface Temperature，SST）和叶绿素浓度（Chlotophyll Concentration，CHL）数据与边界进行叠加分析，研究黑洞涡旋在水平方向上对物质的输运能力。最后，利用Argo浮标提供的三维采样数据绘制Eddy A的剖面图，验证黑洞涡旋垂直方向上作为输运边界的有效性，以补充拉格朗日涡旋的三维验证。我们利用Argo剖面数据（盐度、温度和溶解氧（Dissolved Oxygen，DO））分析不同深度黑洞涡旋的物质输运能力，还通过剖面图发现了黑洞涡旋盐度异常的三明治结构。

2 数据和方法

收起

2.1 数据

2.1.1 海平面高度异常数据和速度场数据

本文使用法国AVISO（Archiving Validation and Interpretation of Satellite Oceanographic）数据中心的海平面高度异常（MSLA-H）数据和通过地转流公式计算得到的速度场数据（MSLA-UV, https://www.aviso.altimetry.fr/index.php?id=1271），空间分辨率为25 km×25 km，时间分辨率为1 d。

2.1.2 海表面温度数据

本文使用美国国家航空航天局（NASA）（https://podaac.jpl.nasa.gov）的高分辨率红外辐射计（AVHRR）提供的海表温度数据，空间分辨率为25 km×25 km，时间分辨率为1 d。

2.1.3 海表面盐度数据

本文使用哥白尼海洋环境监测中心（CMEMS）（http://marine.copernicus.eu/）提供的周平均SSS数据，空间分辨率为25 km×25 km。

2.1.4 叶绿素浓度数据

叶绿素浓度数据由CMEMS （http://marine.copernicus.eu/）提供。本文使用基于SeaWiFS、MODIS、MERIS和OLCI-S3A&S3B等传感器的融合四级地图产品，空间分辨率为4 km×4 km，时间分辨率为1 d。

2.1.5 Argo浮标数据

本文使用2014年3月底投放的17个Argo浮标（世界气象组织的IDs 2901550−2901566）具有增强时空分辨率的采样数据（温度、盐度和溶解氧数据，垂直分辨率为2 m）^[26]，时间分辨率为1 d，初始位置如图1所示。其中一个ID为2901552的Argo大部分盐度数据标记较差，我们在研究过程中删除了该浮标的数据。

2.1.6 澳大利亚联邦科学与工业研究组织区域海洋地图数据

本文使用澳大利亚联邦科学与工业研究组织（CSIRO）（http://www.marine.csiro.au/～dunn/）提供的区域海洋地图数据（CARSs）计算温度异常和盐度异常。在现代海洋测量的时间尺度上，CARSs数据提供了空间分辨率为50 km×50 km的海洋气候态平均网格数据。

2.1.7 研究范围

基于Argo浮标与Eddy A整个生命周期的轨迹，我们选取西太平洋（20°～50°N，130°～160°E）作为研究区域。

2.2 黑洞涡旋识别方法

在二维速度场$ v(x,t) $作用下，一个区域内的粒子经过时间$ t $后的运动轨迹为对速度场求偏导

(1)${{x}}={{v}}\left( {{{x}},{{t}}} \right).$

定义流场的轨迹场如下：

(2)$F_{{{{t}}_0}}^t({x_0}):=x(t;{t_0},{x_0}).$

这表示粒子在初始时刻t₀的初始位置$ {x}_{0} $在流场中经过时间$ t $后到达新位置的运动轨迹。

柯西−格林应变张量$ {C}_{{t}_{0}}^{t}\left({x}_{0}\right) $定义如下：

(3)$C_{{t_0}}^t\left( {{x_0}} \right)={\left[ {\nabla F_{{t_0}}^t\left( {{x_0}} \right)} \right]^{\rm{T}}}\nabla F_{{t_0}}^t\left( {{x_0}} \right),$

式中，$ \nabla$表示梯度计算。

物质曲线$ \eta $是在时间区间$ \left[{t}_{0},{t}_{1}\right] $内同一拉伸系数$ \lambda $的封闭环。在$ {t}_{0} $时刻，物质曲线是由柯西格林应变张量的两个线性无关特征值和对应的特征向量组成

(4)$\eta _ \pm ^\lambda \left( r \right)=\sqrt {\frac{{{\lambda _2}\left( {{x_0}} \right) - {\lambda ^2}}}{{{\lambda _2}\left( {{x_0}} \right) - {\lambda _1}\left( {{x_0}} \right)}}} {\xi _1}\left( {{x_0}} \right) \pm \sqrt {\frac{{{\lambda ^2} - {\lambda _1}\left( {{x_0}} \right)}}{{{\lambda _2}\left( {{x_0}} \right) - {\lambda _1}\left( {{x_0}} \right)}}} {\xi _2}\left( {{x_0}} \right),$

式中，λ₁和λ₂为特征值；$ {\xi }_{1} $和$ {\xi }_{2} $是对应的特征向量，$ {\lambda }_{1}\left({x}_{0}\right) < {\lambda }^{2} < {\lambda }_{2}\left({x}_{0}\right) $。

闭合曲线$ r\left(s\right) $满足下列微分公式

(5)$r'\left( s \right)=\eta _ \pm ^\lambda \left( r \right).$

选取适当庞加莱截面上的点由式（5）进行积分计算得到的封闭曲线即为λ-lines^[32]。为了适应不同的涡旋，从庞加莱截面定义的λ-lines的数量（${\lambda }_{{\rm{number}}}$）为

(6)${\lambda _{{\rm{number}}}}={L_{{\rm{PS}}}}/{l_{\rm{Res}}},$

式中，${L}_{{\rm{PS}}}$是庞加莱截面的长度；$\rm{{l}_{Res}}$是速度场辅助网格的长度。

图2详细介绍了黑洞涡旋边界的提取。选取涡旋内涡度极值点作为涡心^[12]。基于张量场，选取一个明显在涡旋外部的点，两点相连即为庞加莱截面，进行积分计算得到不同的λ-lines，如图2所示。如果存在外层的闭合λ-line，该涡旋就是黑洞涡旋，其边界就是最外层的闭合λ-line，其λ值不一定等于1。如果存在λ=1的闭合λ-line，那么涡旋表现出极强的连贯性。当涡旋只包含λ≠1的最外层边界时，该涡旋为次连贯拉格朗日涡旋，λ<1连贯性增强，λ>1连贯性减弱^[12]。基于本文的识别结果，Eddy A存在封闭的λ=1的λ-line，也有一个λ=0.9（连贯性增强）的最外层边界，因此，Eddy A是一个连贯性很强的黑洞涡旋^[12]。

3 结果与讨论

收起

3.1 遥感数据验证

为了分析不同时间尺度的黑洞涡旋的物质输运能力，我们识别初始时刻为2014年4月10日，时间尺度为T=30 d，60 d，90 d和120 d的黑洞涡旋边界。研究中使用的欧拉边界是一种基于海平面高度（Sea Surface Height, SSH）得到的有效边界^[33]。本文通过虚拟粒子随流场的运动轨迹来比较黑洞涡旋边界和欧拉边界的连贯性^[14]（图3）。通过观察初始位置（图3a），我们发现欧拉边界的面积大于所有时间尺度的黑洞涡旋边界，时间尺度越长的黑洞涡旋边界的面积越小。边界随流场运动到第30 d时，欧拉边界出现了严重的成丝现象，时间尺度为30 d的黑洞涡旋边界保持着很好的连贯性（图3b）。时间尺度为60 d、90 d和120 d的黑洞涡旋边界在随地转流运动的过程中，位置发生了改变，但是没有出现明显的丝状结构，仍然包裹着边界内部的粒子移动（图3c，图3d）。由此可见，黑洞涡旋边界在时间尺度内保持着连贯性，边界阻碍着粒子与周围水体混合。对于一些黑洞涡旋边界，即使随流场运动的时间超过时间尺度，边界出现丝状结构，但粒子还是被包围在边界内，没有分散。这体现出黑洞涡旋的边界随流场运动时具有稳定性，并且对其内部的物质具有较强的包裹能力。

本文将海表面遥感数据（SST异常、SSS异常和CHL异常数据）与黑洞涡旋边界进行叠加，分析黑洞涡旋在水平方向上对物质的输运能力（图4至图6）。图4a展示了初始时刻（2014年4月10日）Eddy A的欧拉边界，以及时间尺度分别为30 d，60 d，90 d和120 d的黑洞涡旋边界。通过观察图4，我们发现黑洞涡旋边界内存在高温度异常。随着地转流运动，欧拉边界分散严重，而黑洞涡旋的位置发生变化，但在时间尺度内仍保持连贯性，并且图4c和图4d的黑洞涡旋边界内存在明显的高温度异常。通过SSS异常与Eddy A边界的叠加分析可以发现，黑洞涡旋内部对应低盐度异常特征，并且这一特征在黑洞涡旋边界和地转流运动的整个周期中都存在（图5）。图6a中，黑洞涡旋内叶绿素浓度较低，图6b、图6c和图6d中Eddy A的黑洞涡旋边界内也对应低叶绿素浓度。综上所述，黑洞涡旋边界被认为是物质输运的屏障，与海表面物质分布密切相关，证明了黑洞涡旋边界可以描述海表面属性的边界。Eddy A作为一个反气旋涡使海表面水体向下运动，表现为高海表面温度异常、低盐度异常和低叶绿素浓度异常^[34]，这与我们的研究结果一致。图4a中海表面温度异常与涡旋边界没有严格对应，这是由于我们使用的是海表面温度异常数据，但涡旋是三维水体，所引起的是三维上的温度异常，Eddy A的暖核可能在海表面以下，在海表面体现的温度异常不够典型。图4b右侧的子图中，由于欧拉边界随地转流变形扭曲，在黑洞涡旋旁边出现旋转结构可能是没有在积分时间内保持连贯性的涡旋或涡流结构。图6b中叶绿素浓度异常与边界的对应不如图6a和图6c严谨，这是因为叶绿素浓度受光照、降水和生物因素等一系列海洋参量的影响，影响原因比较复杂，但是根据图6a和图6c，我们认为该涡旋在大部分时间内可以很好的对应。另外，时间尺度90 d和120 d的黑洞涡旋边界会在切向上出现细小的成丝现象，这是切向上的剪切力导致的，不影响边界的连贯性。总体来说，尽管欧拉边界可以描述瞬时的涡旋位置和范围，但是高估了欧拉涡旋对物质的输运能力，而黑洞涡旋能更客观地描述在一定时间内的物质的输运边界。

3.2 Argo数据验证

本文利用Argo浮标得到的温度、盐度和溶解氧数据，对Eddy A进行不同深度的物质输运能力分析。将Argo数据和CARSs的气候态数据进行处理得到不同深度的温度异常、盐度异常和溶解氧数据，用于追踪涡旋内部的物质变化。通过图7可以看到在90 d内有多个Argo浮标位于Eddy A内部，然而，Argo浮标本身的密度和物理特性与海水不同，所以浮标轨迹不能精确地反映海流方向，但Argo的观测数据可以反映涡旋内部盐度、温度和溶解氧的属性。

欧拉方法描述的是涡旋的瞬时边界，随流场运动会很快失去连贯性，因此本文只研究黑洞涡旋边界在一定时间内的剖面结构。图8为Eddy A边界内外温度和盐度异常随深度变化图。温度异常的峰值出现在深度为400 m附近（图8a、图8d、图8g和图8j）。盐度异常在深度为200 m时出现峰值后迅速下降，然后从500 m又开始上升，在700 m附近再次出现峰值（图8b、图8e、图8h和图8k），形成很明显的三明治结构。水体的垂直分布是决定盐度异常剖面的重要因素，Eddy A的上升运动使盐度较低的北太平洋中间层水（North Pacific Intermediate Water）上升，降低北太平洋热带水（North Pacific Tropical Water）的盐度，所以深度为200～500 m时，深度越小盐度异常越低；而盐度较高的北太平洋热带水被抬升后使浅层海水（0～200 m）呈现正盐度异常。在更深层的海水，由于北太平洋深层水（the North Pacific Deep Water）的盐度相对较高，Eddy A的运动将其抬升，使得深度为500～800 m时低盐度异常减弱^[35]。溶解氧的浓度随深度不断增加，在深度为0～400 m时变化最快。由于受到各种海面因素的影响，海表面的信号可能与深层信号不一致。不同时间，不同深度涡旋内外物质分布不同，但涡旋内物质随深度变化基本保持一致，特别是温度和盐度，这充分说明黑洞涡旋在垂直方向上也能包裹物质运动。

通过建立Argo数据点坐标系，可以观察到不同深度的温度和盐度异常分布与黑洞涡旋边界的关系（图9）。通过计算边界经过流场运动后与涡旋中心的距离，得到不同时间积分的平均拉格朗日边界，根据Eddy A的空间范围，我们忽略与涡心相对距离超过500 km的Argo数据点。根据温度和盐度异常数据得到拟合曲线并计算拟合曲线二阶导数为0的点，该点坐标表示拟合曲线变化最快的位置。虽然温度异常和盐度异常的拟合曲线在浅层海水比较平缓（图9a，图9b），但还是和较深层一样呈现涡旋中心高异常值，涡旋边界低异常值的趋势。随着深度的增加，相较于图9a、图9b、图9i和图9j，Argo数据点的温盐异常值距涡心越远下降趋势越明显（图9c至图9h），这一趋势在深度为900 m时减弱（图9i，图9j）。结合不同深度来看，二阶导数为0的点几乎总是在积分时间大于等于90 d的平均黑洞涡旋边界附近，即长时间积分的黑洞涡旋边界可以反映垂直方向上物质和热量变化最快的位置。

从图10中可以观察到，随着海表面水从辐合中下降，不同深度的DO值变低。并且在较深层的水体中，涡旋中心附近的DO值高于边缘处。随着深度的增加，拟合曲线的下降趋势更加明显，说明在深水中，靠近涡旋中心的水体相对于边缘水体富氧程度更高。和温度异常和盐度异常相同，DO梯度变化最大的点总是在长时间积分的拉格朗日边界附近。这表明Eddy A具有将上层富氧水汇聚并将其输送到底层的能力。这些发现证明了随着流场变化，黑洞涡旋在垂直方向也具有对内部物质的包裹能力，并且能将上层热量和物质输运到深层。

4 结语

收起

本文使用基于eLCSs的黑洞涡旋识别方法，确定西太平洋一个典型的黑洞涡旋边界。通过虚拟粒子随着流场的变化来验证黑洞涡旋比欧拉涡旋更具有连贯性。利用遥感数据SST、SSS和CHL数据分析黑洞涡旋在水平方向上的物质输运能力，随着时间的变化，黑洞涡旋边界的移动与海表面温度、海表面盐度和叶绿素浓度的移动相吻合，证明了黑洞涡旋在海表面方向上作为输运边界的客观性。首次利用Argo浮标获取的采样数据（温度、盐度和DO）形成黑洞涡旋剖面图，证明了不同时间尺度的黑洞涡旋在垂直方向上可以裹挟涡旋内部的水体运动，而不同深度的温度、盐度和DO变化最快的点与长时间尺度的黑洞涡旋边界的一致性体现了其作为物质输运边界的可靠性。

LCS方法强调了物质在流场中的连贯性，因此拉格朗日涡旋具有非常卓越的客观性。我们将着重于利用拉格朗日方法的优点来推进涡旋的识别和验证工作。随着更高时空分辨率遥感数据的出现，我们能更精细地分析海洋结构，更准确地评估物质输运情况，本文中黑洞涡旋不能准确地反映海表面物质分布的问题也可以得到解决。此外，随着模式数据和再分析数据的不断升级，对海洋涡旋二维和三维结构的细致刻画值得进一步研究，也需要具有更精确的空间和时间分辨率的观测数据来精细地评价中尺度涡旋的三维物质输运能力。在未来工作中，我们会对黑洞涡旋边界内存在的高温度异常、低盐度异常和低叶绿素浓度异常进行定量刻画，并进一步定量分析光照、降水和生物因素等对涡旋物质包裹性的影响。

基金

收起

国家自然科学基金重点项目（42030406）；国家重点研发计划（2016YFC1401008，2016YFC1402608）；山东省海洋科技基金（2018SDKJ0102-08）；国家重点研发计划（2017YFA0603203）；自然资源部重点研究项目（KF202005085）；中欧龙计划5科技合作项目（58393）

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Chen G, Hou Y, Chu X. Mesoscale eddies in the South China Sea: Mean properties, spatiotemporal variability, and impact on thermohaline structure[J]. Journal of Geophysical Research: Oceans, 2011, 116(C6): C06018.

Faghmous J H, Styles L, Mithal V, et al. Eddyscan: A physically consistent ocean eddy monitoring application[C]//2012 Conference on Intelligent Data Understanding. Washington, DC: IEEE, 2012: 96−103.

Liu Y, Chen G, Sun M, et al. A parallel SLA-based algorithm for global mesoscale eddy identification[J]. Journal of Atmospheric and Oceanic Technology, 2016, 33(12): 2743−2754.

Faghmous J H, Le M, Uluyol M, et al. A parameter-free spatio-temporal pattern mining model to catalog global ocean dynamics[C]//2013 IEEE 13th International Conference on Data Mining. Dallas: IEEE, 2013: 151−160.

Chelton D B, Schlax M G, Samelson R M. Global observations of nonlinear mesoscale eddies[J]. Progress in Oceanography, 2011, 91(2): 167−216.

Chaigneau A, Le Texier M, Eldin G, et al. Vertical structure of mesoscale eddies in the eastern South Pacific Ocean: A composite analysis from altimetry and Argo profiling floats[J]. Journal of Geophysical Research: Oceans, 2011, 116(C11): C11025.

Beron-Vera F J, Wang Y, Olascoaga M J, et al. Objective detection of oceanic eddies and the Agulhas leakage[J]. Journal of Physical Oceanography, 2013, 43(7): 1426−1438.

Sun M, Tian F, Liu Y, et al. An improved automatic algorithm for global eddy tracking using satellite altimeter data[J]. Remote Sensing, 2017, 9(3): 206.

Tian Fenglin, Wu Di, Yuan Liming, et al. Impacts of the efficiencies of identification and tracking algorithms on the statistical properties of global mesoscale eddies using merged altimeter data.[J] International Journal of Remote Sensing, 2020,41(8): 2835−2860.

夏琼. 典型海洋中尺度涡旋特征分析研究[D]. 青岛: 中国科学院大学(中国科学院海洋研究所), 2018.

Xia Qiong. Study on the characteristics of typical oceanic mesoscale eddies[D]. Qingdao: University of Chinese Academy of Sciences (Institute of Oceanology, Chinese Academy of Sciences), 2018.

Wang Y, Olascoaga M J, Beron-Vera F J. Coherent water transport across the South Atlantic[J]. Geophysical Research Letters, 2015, 42(10): 4072−4079.

Haller G, Beron-Vera F J. Coherent Lagrangian vortices: The black holes of turbulence[J]. Journal of Fluid Mechanics, 2013(731): R4-1−R4-10.

Haller G. Lagrangian coherent structures[J]. Annual Review of Fluid Mechanics, 2015, 47: 137−162.

Beron-Vera F J, Hadjighasem A, Xia Q, et al. Coherent Lagrangian swirls among submesoscale motions[J]. Proceedings of the National Academy of Sciences, 2019, 116(37): 18251−18256.

Haller G, Yuan G. Lagrangian coherent structures and mixing in two-dimensional turbulence[J]. Physica D: Nonlinear Phenomena, 2000, 147(3/4): 352−370.

V Artale, Boffetta G, Celani A, et al. Dispersion of passive tracers in closed basins: Beyond the diffusion coefficient[J]. Physics of Fluids, 1997, 9(11): 3162−3171.

Aurell E, Boffetta G, Crisanti A, et al. Predictability in the large: an extension of the concept of Lyapunov exponent[J]. Journal of Physics A: Mathematical and General, 1997, 30(1): 1−26.

d'Ovidio F, Fernández V, Hernández-García E, et al. Mixing structures in the Mediterranean Sea from finite-size Lyapunov exponents[J]. Geophysical Research Letters, 2004, 31(17): L17203.

Haller G. A variational theory of hyperbolic Lagrangian coherent structures[J]. Physica D: Nonlinear Phenomena, 2011, 240(7): 574−598.

Haller G, Beron-Vera F J. Geodesic theory of transport barriers in two-dimensional flows[J]. Physica D: Nonlinear Phenomena, 2012, 241(20): 1680−1702.

Haller G, Hadjighasem A, Farazmand M, et al. Defining coherent vortices objectively from the vorticity[J]. Journal of Fluid Mechanics, 2016, 795: 136−173.

Abernathey R, Haller G. Transport by Lagrangian Vortices in the Eastern Pacific[J]. Journal of Physical Oceanography, 2018, 48(3): 667−685.

Rypina I I, Kamenkovich I, Berloff P, et al. Eddy-induced particle dispersion in the near-surface north Atlantic[J]. Journal of Physical Oceanography, 2012, 42(12): 2206−2228.

Beron-Vera F J, Olascoaga M J, Wang Y, et al. Enduring Lagrangian coherence of a Loop Current ring assessed using independent observations[J]. Scientific Reports, 2018, 8(1): 1−12.

Dong C, McWilliams J C, Liu Y, et al. Global heat and salt transports by eddy movement[J]. Nature Communications, 2014, 5: 3294.

Xu Lixiao, Li Peiliang, Xie ShangPing, et al. Observing mesoscale eddy effects on mode-water subduction and transport in the North Pacific[J]. Nature Communications, 2016, 7: 10505.

Zhang Zhengguang, Wang Wei, Qiu Bo. Oceanic mass transport by mesoscale eddies[J]. Science, 2014, 345(6194): 322−324.

Hu Jianyu, Gan Jianping, Sun Zhenyu, et al. Observed three-dimensional structure of a cold eddy in the southwestern South China Sea[J]. Journal of Geophysical Research:Oceans, 2011, 116(C5): C05016.

Dong Changming, Lin Xiayan, Liu Yu, et al. Three-dimensional oceanic eddy analysis in the Southern California Bight from a numerical product[J]. Journal of Geophysical Research: Oceans, 2012, 117(C7): C00H14.

Lin Xiayan, Dong Changming, Chen Dake, et al. Three-dimensional properties of mesoscale eddies in the South China Sea based on eddy-resolving model output[J]. Deep-Sea Research Part I: Oceanographic Research Papers, 2015, 99: 46−64.

Blazevski D, Haller G. Hyperbolic and elliptic transport barriers in three-dimensional unsteady flows[J]. Physica D: Nonlinear Phenomena, 2014, 273−274: 46−62.

Onu K, Huhn F, Haller G. LCS Tool: A computational platform for Lagrangian coherent structures[J]. Journal of Computational Science, 2015, 7: 26−36.

Mason E, Pascual A, McWilliams J C. A new sea surface height-based code for oceanic mesoscale eddy tracking[J]. Journal of Atmospheric and Oceanic Technology, 2014, 31(5): 1181−1188.

Leterme S C, Pingree R D. The gulf stream, rings and north atlantic eddy structures from remote sensing (altimeter and seawifs)[J]. Journal of Marine Systems, 2008, 69(3/4): 177−190.

Yang Guang, Wang Fan, Li Yuanlong, et al. Mesoscale eddies in the northwestern subtropical Pacific Ocean: Statistical characteristics and three-dimensional structures[J]. Journal of Geophysical Research Oceans, 2013, 118(4): 1906−1925.

2021年第43卷第12期

PDF下载

128

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2021163

接收时间：2020-08-06
首发时间：2026-02-26
出版时间：2021-12-30

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2020-08-06
修回日期：2020-12-28

基金

国家自然科学基金重点项目（42030406）；国家重点研发计划（2016YFC1401008，2016YFC1402608）；山东省海洋科技基金（2018SDKJ0102-08）；国家重点研发计划（2017YFA0603203）；自然资源部重点研究项目（KF202005085）；中欧龙计划5科技合作项目（58393）

作者信息

¹中国海洋大学信息科学与工程学部海洋技术学院，山东青岛 266100

²青岛海洋科学与技术试点国家实验室区域海洋动力学与数值模拟功能实验室，山东青岛 266237

通讯作者:

陈戈（1965—），男，博士，教授，博士生导师，主要从事海洋遥感与大数据工作。E-mail：gechen@ouc.edu.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2021163

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT