基于自动探测方法的南海北部涡致海表温度锋面的特征研究

南海北部具有丰富的温度锋面和中尺度涡,它们调节着局地的热量和能量平衡。本文利用卫星海洋高度异常和海表温度数据,并基于自动探测方法,探究了2007年至2017年南海北部中尺度涡边缘的海表温度锋面（涡致锋面）特征。反气旋/气旋边缘出现锋面的概率可达20%。气旋涡在各个方向上出现锋面的概率比较均匀,反气旋涡的东北部和西南部出现锋面的概率大于西北部和东南部。中尺度涡致锋面的数量有明显的季节变化,而涡动能未表现出明显的季节变化。中尺度涡致锋区的总涡动能是中尺度涡内动能的3倍,并且反气旋涡致锋面的总涡动能明显强于气旋涡致锋面的总涡动能。中尺度涡致锋面的数量和涡动能的年际变化与厄尔尼诺南方涛动指数没有明显的相关性。本研究也讨论了中尺度涡致锋面的可能机制,但是中尺度涡对海表温度锋的贡献需要进一步定量研究。     

海洋中尺度涡旋源汇空间分布特征研究

借助卫星高度计数据,对中尺度涡进行识别和追踪;以16年内中尺度涡个数上的生消为判据,发现中尺度涡在除赤道外的全球大洋中生消频繁,但在海盆内区并没有明显的生成占优区或消亡占优区;而在中纬度近岸的狭窄东边界内中尺度涡生成居多,在另一侧近岸的狭窄西边界中尺度涡消亡居多。同时,我们以一阶斜压模态所对应的特征深度作为两层结构的内界面深度,并假定涡动能平均分配于正压模态和一阶斜压模态,计算得到了包含涡动能和涡有效重力位能的中尺度涡能量年平均净生成率和净耗散率,发现虽然海盆西边界是涡场能量耗散大于能量生成的区域,但强耗散过程实际上在海盆西侧内区的强流及其回流区均有发生。另外,中尺度涡生消个数差值的分布与中尺度涡能量净生成率和净耗散率的分布表明,虽然海盆东边界近岸区域内中尺度涡的生成居多,但该区域中尺度涡的能量偏弱,因此该区并非涡场能量的主要源区。     

西北太平洋中尺度涡半径与涡动能的统计关系

中尺度涡的半径与涡动能之间的关系及其内在机制的研究,对我们理解中尺度涡旋的三维结构特征及其变异规律有很大的帮助。本文利用AVISO卫星高度计融合数据,基于流场几何特征的识别方法对西北太平洋(10°～52°N, 120°～180°E)的中尺度涡进行了识别追踪,并由此分析了该区域内中尺度涡的半径与涡动能(Eddy Kinetic Energy, EKE)的统计关系。结果表明该区域中尺度涡的半径与EKE之间存在类似翻转高斯函数的递增关系,特别是在涡旋半径达到一定程度后,EKE将不随半径的变大而继续增强。而且该关系与纬度密切相关,相同半径条件下,涡旋的EKE近似与其所在纬度的第一斜压Rossby变形半径的平方呈反比。进一步分析表明该统计关系与中尺度涡所处的发展阶段没有明显的联系。该结果为建立一个实用性的参数化、归一化的中尺度涡模型提供了一定的理论基础。     

黑潮延伸体上游中尺度涡场的年代际振荡及其相关机制

黑潮延伸体上游区域的中尺度涡场的涡动能和涡特征尺度存在显著地年代际振荡,和黑潮延伸体路径的年代际变化有很好的相关性。当黑潮延伸体路径比较稳定时,其上游区域涡动能比较高,涡特征尺度比较大,反之相反。通过对黑潮延伸体上游区域的中尺度涡场进行集合分析发现:当黑潮延伸体处于稳定状态时,上游涡场几乎是各向均匀地,有轻微的径向伸长;而当黑潮延伸体处于不稳定状态时,上游的中尺度涡场有显著地纬向伸长。对与中尺度涡场的产生相关的线性斜压不稳定和正压不稳定进行了计算分析,结果显示,线性斜压不稳定不是控制中尺度涡场年代际变化的机制,而正压不稳定对中尺度涡场的年代际变化有积极的贡献。不稳定产生的中尺度涡之间存在非线性涡-涡相互作用。     

2000年8月南海中部与南部海洋温、盐与环流特征

根据2000年8—9月份南海中部与南部航次的温、盐资料,采用P—矢量诊断方法,结合ADCP测流资料和同期伪风应力资料以及TOPEX／Poseidon高度计资料,研究了2000年夏季风持续强迫之后南海大尺度环流与中尺度涡旋的空间结构。结果表明,南海夏季温度和盐度水平分布随深度有显著的变化：中层(250—400m左右)温、盐水平分布与其它各层的温度和盐度分布相比有很大的差异。用诊断方法计算得到的环流场与用TOPEX／Poseidon海面高度计资料计算得到的地转流场比较一致,即流场内部有多个中尺度的涡旋,主要有越南东南外海反气旋涡、中沙群岛东南反气旋涡以及南沙群岛东北角的气旋涡等,这说明南海中部与南部盛夏环流具有较强的地转分量和显著的多涡结构,并且这些中尺度涡在垂向上存在速度场的切变。     

南海中尺度涡边缘亚中尺度过程模式研究

基于高分辨率ROMS模式数据,对南海中尺度涡边缘的亚中尺度过程进行了研究。结果表明:在中尺度涡边缘处存在明显的温度细丝结构,涡度值要明显大于1,O~W参数趋于正值,即在涡旋边缘存在强速度剪切,加强了涡丝的形成。对亚中尺度过程动能空间分布分析发现,中尺度涡边缘区域动能所占的比例大约为88.4%,比中心区域动能所占的比例11.6%要大许多。从动能谱的分析来看,在中尺度涡边缘处动能谱的斜率趋近于k-2,与表层准地转理论(SQG)相吻合,而在中心处呈现k~(-3),与准地转理论(QG)相一致。     

Rhines效应对高度计观测的海洋中尺度涡空间分布特征的影响

Rhines效应是指Rossby波和大湍流（中尺度涡）相互作用,将涡动能量以波的形式传播出去,从而使中尺度涡发生形变,最终消亡的一种动力学机制。本文通过比较海洋里涡特征速度和Rossby长波波速的方法,研究了一种广义的Rhines效应对高度计观测的海洋中尺度涡空间分布特征的影响。结果显示,广义Rhines效应比只考虑行星涡度梯度的传统形式对中尺度涡的分布具有更显著的影响。大部分中尺度涡分布在涡特征速度（U_e）大于由广义Beta值计算的Rossby长波波速（U_cg）的区域。这些涡可以由动能反向串级过程获取能量,成长为振幅和空间尺度较大的涡。热带海域以外的“涡旋沙漠”区域,中尺度涡的数量稀少,强度很弱,大都分布于U_ecg的海域。广义Rhines效应可能是这些海域中尺度涡难以成长的动力学机制。     

基于四维变分同化方法的南海中尺度涡后报实验

海洋中尺度涡在本质上是属于满足准地转平衡的大尺度运动,因此理论上,其在短时间内的运动将主要受到准地转平衡关系的约束,而外部强迫场的影响在短期内不会明显改变其运动特征。基于上述思想,我们提出了一种基于四维变分同化初始场的中尺度涡旋预报方案。为了检验该方案的可行性,本文使用区域海洋模式(regional ocean modeling system, ROMS)和其内建的增量强约束四维变分同化(incremental strong constraint four dimensional variational, I4D-Var)模块,建立了一个南海海洋同化模拟系统。首先,通过I4D-Var方法将AVISO卫星高度计资料同化到海洋数值模拟中,获得了理想的中尺度涡同化模拟结果。同化、模式模拟和观测三者的中尺度涡统计结果表明,该同化系统模拟的南海中尺度涡的路径、半径、海表高度异常和振幅等特征信息与AVISO(Archiving ValidationandInterpolationofSatelliteOceanographicData)观测结果高度吻合,同时在深度上的分析表明,涡旋对应的温度、盐度和密度均得到有效的调整。然后,将该同化系统的模拟结果做为初始场,对某一特定时段的南海中尺度涡进行了后报模拟和结果的定量化分析。通过比较后报模拟与观测资料中对应涡旋的海表面高度异常(sea surface height anomalies, SSHA)相关系数、涡心差距和半径绝对误差,证明该方案的中尺度涡后报时效至少可达10 d以上。后报实验结果验证了该中尺度涡预报方案的可行性,从而为中尺度涡的预报提供一定的理论基础和可行性方案。     

基于HYCOM的南海中尺度涡数值模拟

结合AVISO(Archiving Validation and Interpolation of Satellite Oceanographic Data)高度计资料,利用改进的NERSCHYCOM(Nansen Environmental and Remote Sensing Center-Hybrid Coordinate Ocean Model)大洋环流模式,对南海中尺度涡进行数值模拟研究,主要包括中尺度涡的三维结构、南海EKE(Eddy Kinetic Energy,涡动动能)的垂向变化、黑潮中尺度涡的脱落以及涡旋近岸时的结构变化等。模式再现了2007年2月-3月菲律宾西侧海域的一次暖涡过程,探究了其生命期中各阶段的特征物理量的变化,对其成熟时期的涡旋结构研究表明,中尺度涡的结构呈现不对称性,涡旋两侧的流场空间范围和流场强度均不相同,涡旋的半径和中心位置随深度不断变化,并且由涡旋作用产生的升降流的中心与涡旋自身中心并不重合,二者之间有一定距离。初步探索EKE的垂向分布情况,认为南海年平均EKE在垂向变化上呈现三段式,主要部分分布在300m以浅深度,但同时垂向又能达到海洋深层。分析了一次黑潮中尺度涡脱落的模式模拟个例,推测黑潮中尺度涡脱落原因:黑潮流径西移、外海中尺度涡对黑潮的强迫、地形作用,并且结果表明从黑潮脱落的中尺度涡可以携带大量高温高盐水体进入南海,对南海的温盐性质产生很大的影响。初步探索涡旋近岸时的结构变化,涡旋靠近岸界时,受岸界挤压,流速在一段时间内会增大,继续靠近岸界,由于岸界的摩擦、海底的拖曳,导致能量耗散,流速减小,最终涡旋消亡。     

南海中尺度涡对深层流调控作用的个例研究

南海深层环流作为印太海域热盐环流的重要组成部分,研究其时间变异特征对于深入认识印太海域大洋环流具有重要意义。中尺度涡作为南海极为活跃的动力过程,有观测显示其影响深度可超过2000m。结合海表面高度计资料与潜标观测资料,识别南海区域中尺度涡,并分析了中尺度涡对深层流场以及温度场的影响。针对三个中尺度涡的潜标实测结果表明:在涡旋经过时,海洋深层流动有明显的变异,表现在速度、温度方面,速度增加量能达到5.5cm/s,温度变异可达到0.02℃。涡旋经过时,海洋上层、深层流向呈相反态势,第一斜压模动能显著增强。     

Numerical study on the eddy–mean flow interaction between a cyclonic eddy and Kuroshio

A three-dimensional primitive equations ocean model (POM) is employed to study the momentum and energy balance of a moving cyclonic eddy (CE) during eddy–mean flow interaction. The CE generated by an idealized typhoon forms to the east of the Philippine islands. A momentum balance analysis shows that the dynamics of the CE are generally dominated by the geostrophic current throughout the life cycle of the CE. An energy analysis suggests that the eddy kinetic energy (EKE) and the eddy potential energy (EPE) decay rapidly after generation. The maximum EPE initially appears at the surface of the eddy center and gradually appears in the subsurface layer. The largest baroclinic instability (BCI) initially occurs at the surface. For a CE moving along a trajectory, the conversion from mean potential energy (MPE) to the EPE is positive (negative) in the front (rear) part of the trajectory, and then the eddy transfers its EPE forward along its trajectory by means of the front (rear) part of the eddy obtaining (losing) EPE from (to) the mean flow. During the interaction stage, the northward flowing Kuroshio interacts with the southward flow on the western side of the eddy and the inverse velocity shear between the Kuroshio and the eddy causes the EKE to gradually develop east–west asymmetry. The largest barotropic instability (BTI) is found in the interaction zone. Advection term, pressure work, and friction term play the dominating role in eddy decay in the eddy zone, while BTI only dominates in the interaction zone.     

北太平洋北部的风能输入门户

风不仅驱动了上层海洋的环流,也是深层海洋运动的主要能量来源.本文主要研究了北太平洋北部的风能输入的季节性分布特征和年际变化趋势,包括风向表面波、表层地转流和表层非地转流的能量输入.基于SODA3数据的结果表明,风能输入门户随季节变化显著,其中黑潮延伸区是冬季门户,副极地流涡是春、秋季门户,大洋东边界则是夏季门户,能量输...     

Surface available gravitational potential energy in the world oceans

Satellite altimetry observations, including the upcoming Surface Water and Ocean Topography mission, provide snapshots of the global sea surface high anomaly field. The common practice in analyzing these surface elevation data is to convert them into surface velocity based on the geostrophic approximation. With increasing horizontal resolution in satellite observations, sea surface elevation data will contain many dynamical signals other than the geostrophic velocity. A new physical quantity, th...     

How faithfully will the geostrophic currents represent the existing ocean currents?

It is widely recognized that the geostrophic flows computed by the dynamic method of Bjerknes and collaborators represent the actual currents pretty faithfully. However, what would be the reason that a geostrophic current derived by only retaining the terms of Coriolis and the pressure gradient forces in the hydrodynamical equations agrees so closely with the actual ocean current of the same area? In this attempt was assumed an imaginative ocean of homogeneous water and uniform depth on a rotating earth but with neither continent nor islands. The average wind distribution observed along several meridians over the Pacific Ocean was assumed to prevail in this sea throughout with no variation in east-west direction. Taking the curvature of the earth surface, rotation of the earth, Coriolis forces, pressure gradients and the horizontal and vertical eddy viscosity into account, the equations of motion were solved and velocity components were derived for all latitudes. A comparison of the east-west components thus obtained with the corresponding components of the geostrophic flows, reveals that they agree well in higher latitudes but there appears a remarkable disagreement in lower latitudes. This means that a special care must be taken in replacing the existing currents with the geostrophic flows at lower latitudes.     

Subtleties in forcing eddy resolving ocean models with satellite wind data

Using new global satellite remote sensing data, we show that ignoring the ocean current dependence in the wind stress artificially increases global wind power input to the oceanic general circulation by about 32%, and more than doubles the input in the regions of strong ocean current systems. Scatterometer-derived wind stress naturally accounts for the moving ocean that is not included in traditional wind stress products. However, forcing an ocean model with a scatterometer-derived wind stress cannot actually account for the ocean current effect on the wind power input. The difference between the real and modeled surface eddy fields can reduce the damping associated with the ocean current dependence in wind stress, leading to a positive bias in global wind power input of about 23%. Most of this spurious energy flux goes directly to the fluctuation eddy field and is several times larger than the energy flux to real ocean eddies.     

Velocity Statistics Inferred from the TOPEX/Poseidon-Jason-1 Tandem Mission Data

Using a parallel-track approach to estimate geostrophic surface velocities, an estimate of the statistics of ocean geostrophic surface currents and momentum stresses is provided on a 10 km along-track resolution from the first 49 repeat cycles (16 months) of the Jason-TOPEX/Poseidon tandem altimetric sea surface height (SSH) data. Results are compared with estimates obtained in a traditional way from along-track SSH data at crossover points and with in situ, Acoustic Doppler Current Profiler (ADCP) measurements obtained on board the VOS Oleander along a nominal path connecting Bermuda with the U.S. mainland. Agreements with the Oleander data are reasonable when simultaneous (in space and time) sampling is available. However, amplitudes of parallel-track geostrophic velocity variances are about 25% lower as compared to Oleander measurements which represent geostrophic and ageostrophic flow components. Estimates of velocity variances show clear signs of an anisotropic eddy field in the vicinity of all major current systems. At the same time estimates of Reynolds stresses and eddy momentum fluxes show a convergence of eddy momentum in all those regions, suggesting a forcing of the mean flow by the eddy field there.     

Velocity Statistics Inferred from the TOPEX/Poseidon-Jason-1 Tandem Mission Data

Using a parallel-track approach to estimate geostrophic surface velocities, an estimate of the statistics of ocean geostrophic surface currents and momentum stresses is provided on a 10 km along-track resolution from the first 49 repeat cycles (16 months) of the Jason-TOPEX/Poseidon tandem altimetric sea surface height (SSH) data. Results are compared with estimates obtained in a traditional way from along-track SSH data at crossover points and with in situ, Acoustic Doppler Current Profiler (ADCP) measurements obtained on board the VOS Oleander along a nominal path connecting Bermuda with the U.S. mainland. Agreements with the Oleander data are reasonable when simultaneous (in space and time) sampling is available. However, amplitudes of parallel-track geostrophic velocity variances are about 25% lower as compared to Oleander measurements which represent geostrophic and ageostrophic flow components. Estimates of velocity variances show clear signs of an anisotropic eddy field in the vicinity of all major current systems. At the same time estimates of Reynolds stresses and eddy momentum fluxes show a convergence of eddy momentum in all those regions, suggesting a forcing of the mean flow by the eddy field there.     

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

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

Residual hydrodynamic fields after tsunami generation by an earthquake

The linear theory of long waves was applied to study horizontal motions of the water layer in a rotating ocean which appear after tsunami generation by an earthquake. The structures of residual potential and eddy fields are analyzed on the basis of the analytical solution of a model axisymmetric problem for an ocean of constant depth. The estimates of the horizontal displacements of water particles, velocity of the eddy current, and energy of the geostrophic eddy are calculated for typical conditions of the tsunami source. Particular features of the residual fields related to the existence of stable stratification are considered. Static and dynamic numerical models are described that allow us to calculate the residual potential field and its evolution related to the realistic events. The field of residual horizontal displacements of water particles for the catastrophic earthquake near the coasts of Japan on March 11, 2011, is calculated and analyzed.     

一种使用卫星高度计数据且基于密度聚类的新型海洋中尺度涡自动探测算法

中尺度涡旋是海洋中典型的中尺度现象,是海洋中能量传递的运输者,中尺度涡识别与提取是物理海洋学研究的重要内容之一,而中尺度涡自动发现算法是最基础的用于寻找与分析中尺度涡的工具。中尺度涡旋探测工作的数据来源主要为卫星高度计数据融合出的SLA数据,该数据可以客观的描述海洋表层高度状态。中尺度涡表示为SLA闭合等值线所包围的局部等值区域,涡旋识别需要从SLA数据中提取出稳定的闭合等值线结构。针对基于SLA数据中的中尺度涡探测的特点,本文提出了一种新的基于聚类方法的中尺度涡自动识别算法,通过对SLA数据集的分割与筛选将中尺度涡区域与背景区域分离,后建立区域内联系并将其映射到SLA地图上来提取中尺度涡结构。本文算法解决了传统探测算法中参数设定的敏感性问题,不需要进行稳定性测试,算法适应性增强。算法中加入了涡旋筛选机制,保证了结果的涡旋结构的稳定性,提高了识别准确率。在此基础上,本文选取了西北太平洋及中国南海地区进行了中尺度涡探测实验,实验结果展示出了本文算法在较传统算法提高算法效率的同时,也保持着较高的算法稳定性,可以在稳定识别各个单涡结构的同时识别稳定的多涡结构。