内孤立波过陆架陆坡地形的数值模拟研究

内波是海洋中普遍存在的波动形式。内孤立波是典型的非线性内波,多发于陆架边缘海,如南海等海域,对陆架海域有重要影响。本文针对内孤立波在陆架地形上的传播问题,先基于弱非线性与全非线性数值模型,模拟了不同振幅、地形高度条件下内孤立波的演化的过程,探讨了动力系数对内孤立波演化过程的影响,对比了两模型的模拟结果在内孤立波演化过程、能量分配以及能量耗散的差异,后分析了南海的动力系数分布特征。结果表明,在内孤立波不发生破碎的情况下,弱非线性模型与全非线性模拟结果相近。当发生破碎过程时,弱非线性模型可准确模拟头波,但无法通过强非线性的破碎过程耗散能量,只能以裂变的方式辐射能量。在弱非线性模型中,随地形高度增加,频散系数减小到零,平方非线性系数由负转正,立方非线性系数绝对值增大一个量级,并主导陆架地形上内孤立波的演化过程。通过对比南海夏季与冬季非线性内波动力系数空间分布,发现内孤立波在传播过程由于夏季平方非线性效应、立方非线性效应与频散效应较强的影响,其在夏季更易发生陡化与裂变,波列发生频率高。     

基于X 波段雷达获取东沙群岛附近内波的传播速度

根据2009年6月24日15时40分至25日16时40分“科学一号”考察船在东沙岛东北部陆架上K106站进行的长达25 h的X波段雷达、温度链、声学多普勒流速剖面仪(ADCP)同步观测数据,在该海域利用 Radon 变换技术获取了本次观测到的内孤立波的传播速度.利用该孤立波到达前30 min 的ADCP 流速值,计算得到内波传播方向上的背景流流速为0.04 m/s.最终得出该孤立波的传播速度为3.04 m/s,传播方向约为297°.基于同样的处理方法,求解了本航次中获取的南海东北部陆架处其他站点部分内波的波速信息.南海东北部陆架处内波,主要向西或西北方向传播.本研究对增进东沙群岛附近内孤立波传播特征的认识,具有重要意义.     

拉格朗日描述下三层流体系统中界面内波的三阶渐近解

在流体力学中,描述流体运动有Lagrange方法和Euler方法.Euler方法是通过观测通过空间各固定位置点处流体质点的运动行为来描述流体运动规律,而Lagrange方法是跟踪各个流体质点,通过观测它们在时空运动中所走过的路径来描述流体的运动规律.在数学处理上,Euler方法较Lagrange方法简单,但Lagrange方法可以完全描述流体运动的整个流场的所有特性,而Euler方法却无法描述每个流体质点的运动轨迹.本文,我们研究具有刚性边界的三层流体系统中的界面内波,其中上层流体的密度比下层流体的密度大.通过在界面处引入朗格朗日匹配条件并使用微扰法得到了拉格朗日描述下的界面内波的一阶解、二阶解及三阶解,给出了质量输运速度、波频率、平均水平和质点运动轨迹的解.结果表明对于质量输运速度、波频率、平均水平和质点运动轨迹在界面处会有不连续性,但是我们发现在满足一定的三层流体水深比和密度比条件时这种不连续性将会消失.     

有限长正弦地形上的波浪传播问题研究

应用基于势流理论的频域边界元方法,研究了波浪在有限长正弦周期地形上的传播问题.开展了网格收敛性分析,并与物理模型试验结果做了对比,验证了数值方法的准确性.研究了正弦地形长度和幅度对水波反射率的影响,以及波浪在周期性地形上传播时的沿程波幅变化.研究发现,波浪反射系数随波数的变化及波浪幅值沿程的变化均呈包络状;当频率处于布拉格共振区域边缘且反射率为零时,周期地形上会出现极大的波幅;随海床波幅增大,第一布拉格共振带变宽且向低频移动.     

Is the small-scale turbulence an exclusive breaking product of oceanic internal waves

On the basis of the theoretical research results by the author and the literature published up to date, the analysis and the justification presented in this paper show that the breaking products of oceanic internal waves are not only turbulence, but also the fine-scale near-inertial internal waves (the oceanic reversible finestructure) for inertial waves and the internal solitary waves for internal tides respectively. It was found that the oceanic reversible finestructure may be induced by the effect of the horizontal component f (f = 2Ωcosφ) of the rotation vector on inertial waves. And a new instability of the theoretical shear and strain spectra due to the effect of f occurs at critical vertical wavenumber β c ≈ 0.1 cpm. It happens when the levels of shear and strain of the reversible finestructure are higher than those of inertial waves, which is induced by the effect of f along an "iso-potential-pycnal" of internal wave. If all breaking products of internal waves are taken into account, the average kinetic energy dissipation rate is an order of magnitude larger than the values of turbulence observed by microstructure measurements. The author’s theoretical research results are basically in agreement with those observed in IWEX, DRIFTER and PATCHEX experiments. An important impersonal fact is that on the mean temporal scale of thermohaline circulation these breaking products of internal waves exist simultaneously with turbulence. Because inertial waves are generated by winds at the surface, and internal tides are generated by strong tide-topography interactions, the analysis and justification in this paper support in principle the abyssal recipes Ⅱ:energetics of tidal and wind mixing by Munk & Wunsch in 1998, in despite of the results of microstructure measurements for the turbulent kinetic energy dissipation rate and the diapycnal turbulent eddy diffusivity.     

具有内波特性的渤海南部次重力波

次重力波(infragravity waves)是频率范围介于局地平均Brunt-V¾is¾l¾频率和0.05Hz(涌浪频率)之间的一种海洋波动,其能量在深水大洋中一般甚低,但在浅海中却很强.该波动的研究对于浅海海洋内部混合、近岸海洋工程、海洋生态学、海洋沉积学及军事海洋学均有重要的意义.Munk[1]很早就开始研究这类波动,然而主要由于观测技术的限制,对于该类波动的动力学机制的研究进展迟缓.     

考虑旋转向量水平分量的大洋内波方程的一个渐近解

对于考虑旋转向量水平分量的大洋内波方程,本文求出了满足该方程的波函数的一个渐近解.该渐近波函数的振荡特性得到研究,结果表明旋转向量水平分量对高模态(高波数)内波的影响是不可忽略的.     

海啸俘获波在弯曲海脊上的传播规律

越洋海啸会受大洋海脊的引导以俘获波的形式沿其传播上万千米,其携带的巨大能量会严重影响远场地区、威胁海岸地区的安全。本文基于MIKE21-BW模型,分别模拟0°(直海脊)至90°(直角弯曲海脊)不同弯曲角度海脊上俘获波的传播变形过程,并定量比较其能量分配。结果表明,海脊俘获波传至海脊转弯处,少部分能量会泄露出海脊重新以自由波的形式扩散至整个海域;部分能量会被反射回来形成与初始海啸波相反方向的俘获波沿海脊传播,反射的能量会随着海脊弯曲角度的增加而增加;还有一部分能量继续沿着弯曲的海脊向前传播,其随着海脊弯曲的角度增加而减小。     

Infragravity waves with internal wave characteristics in the south of the Bohai Sea of China

The measurements by using ADCP (500 KH) and CTD were made during August 2000 in the south (37°55''N, 120°25''E) of the Bohai Sea, where the water depth was about 16.5m. The data of horizontal velocity with sampling interval of 2 min in 7 layers were obtained. The power spec-trum analysis of these data indicates that there are very energetic infragrvity waves with a period of about 6 min. The coherence spetrum analysis and the analysis of temporal variation of shear show that these infragravity waves are mainly the free wave model (properties of edge waves), in the meantimethey possess some characteristics of internal waves, which are likely due to the distinctive marine environment in this area. It is speculated on that the instability processes (chiefly shear instability) of sheared stratified tidal flow owing to the effect of sea-floor slope in the coastal area might be the main mechanism generating these infragravity waves.     

南海西北陆坡区内潮与近惯性内波观测研究

通过对2006年南海西北部海域近3个月的全水深流速观测资料的分析,研究了该海区正压潮、内潮及近惯性内波的时空分布特征。结果表明,全日内潮明显强于半日内潮,且最大潮流均出现在海洋上层;内潮的主轴方向基本沿东南-西北方向,近似与局地等深线垂直;内潮能量显示出明显的时间长度约为半月的大小潮调制周期;全日内潮的coherent部分占全日内潮能量的70%,而半日内潮的coherent部分占半日内潮能量的53%;进一步研究发现半日内潮主要由第一模态主导,而全日内潮第三模态能量占总能量的比例仅次于第一模态且量值上与之相当;强风过程可激发出强的近惯性运动,暖涡使得近惯性内波能量更有效地向海洋深层传播,冷涡则不利于近惯性内波能量向下传播。     

On generation source sites of internal waves in the Luzon Strait

This effort aims to determine the generation source sites in the Luzon Strait for energetic, long-crest, transbasin internal waves （IW） observed in the northern South China Sea （NSCS）. The roles of islands distributed on eastern side of the strait, Kuroshio, submarine ridges, shoaling thennocline, and strait configuration played in the IW generation are examined using the cruise data analysis, satellite data interpretation, and dynamical analysis. The islands and channels on eastern side of the strait are excluded from a list of possible IW source sites owing to their unmatched horizontal dimensions to the scale of IW crest line length, and the relative low Reynolds number. The Kuroshio has a potential to be a radiator for the long-crest IW disturbances, meanwhile, the Kurosbio west （east） wing absorbs the eastward （westward） propagating IW disturbance. Namely, the Kuroshio blockades the outside west-east propagating IW disturbances. The 3-D configuration of the Luzon Strait is characterized by a sudden, more than one order widening of the cross-section areas at the outlets on both sides, providing a favorable condition for IW type initial disturbance formation. In the Luzon Strait, the thermocline is featured by a westward shoaling all the year around, providing the dynamical conditions for the amplitude growth （declination） to the westward （eastward） propagating IW type disturbance. Thus, the west slope of western submarine ridge at the western outlet of the Luzon Strait is a high possibility source sites for energetic, long-crest, transbasin IWs in the NSCS. The interpretation results of satellite SAR images during a 13 a period from 1995 to 2007 provide the convincing evidence for the conclusions.     

地形和层结变化对内潮生成影响的数值模拟研究

基于ROMS(Regional Ocean Modeling System)模式,设计了9组理想数值实验,研究了内潮的产生对地形和背景层结环境的敏感性。结果表明:高斯海脊两侧地形梯度较大的超临界坡面处是内潮的主要源区;地形和层结环境的变化对内潮生成有很大影响;如地形和层结保持不变,不同的网格水平分辨率和垂向层数条件下估算的内潮转换率存在差别;地形和跃层的相对空间关系对内潮生成及其引起的水体混合强度存在影响。高分辨率的地形和合适的背景层结是内潮模拟的关键。该结论对于内潮形成的数值模拟研究有一定参考价值。     

A study of SAR remote sensing of internal solitary waves in the north of the South China Sea: Ⅱ. Simulation of SAR signatures of internal solitary waves

A new model developed from the full-spectrum model of Lyzenga and Bennett （ 1988 ） is built up by us preliminarily through considering the impact of the ocean surface mixed layer turbulence on SAR remote sensing of internal solitary waves. In the partial differential equation of the action spectral density of the surface gravity-capillary waves the source function representing the contribution of the turbulence is added besides the usual source function representing the contribution of the wind. The source function is determined by applying the κ - ε model and adopting the Nasmyth spectrum of oceanic turbulence （ Nasmyth, 1970; Oakey, 1982; Fan, 2002） on the basis of the previous simulation results of internal tide transformation obtained in the companion paper （Fan et al. ,2008）. Either under relatively high wind speed, or under low wind speed, our model predicts significant large modulations of radar backscatter at all three bands （ L, C and X bands） for both VV and HH polarization. These results prove that considering the impact of ocean surface mixed layer turbulence on SAR remote sensing of internal solitary waves is reasonable and appropriate for settlement of the well-known problem of contemporary radar imaging models.     

背景流和混合层深度对内波表面信号的影响

主要利用Lamb 模式进行了四个数值敏感实验, 研究了背景流、混合层深度对南海东北部内波表面信号强度的影响, 并用MASCG 指标(表层流梯度绝对值的最大值)进行了量化。研究结果表明: 同向背景流能增强内波表面流梯度, 逆向背景流则反之。混合层深度对内波表面信号影响很大, 深的混合层能弱化内波表面信号强度。内波表面信号的量化依赖于海洋环境条件和内波动力参数, 本实验结论将加深对内波遥感观测的理解。     

两相流模型在南海内孤立波与沙波相互作用中的应用

为研究内孤立波与沙波的相互作用,本文对基于OpenFOAM的SedWaveFoam求解器进行改进,建立了内孤立波-泥沙运动欧拉两相流模型。在利用试验资料对模型进行验证的基础上,在南海北部典型代表性条件下,模拟分析了500 m水深位置沙波床面上内孤立波作用下的水动力变化和泥沙运动。结果表明,内孤立波逐渐离开沙波时,海底沙波背流面处出现与内孤立波背景流速反向的流速,在内孤立波导致的流场作用下,沙波床面上的泥沙悬起并运动到床面以上的水体中。振幅100 m的内孤立波可以导致床面以上14 m高的位置处出现约0.07 kg/m³的悬沙浓度。     

剪切流中界面波的二阶Stokes波解

考虑了均匀剪切流场中强非线性界面波,建立了基于任意水深处速度而不是通常所用的平均速度为速度变量的模型,分析了其色散关系,并求得了各层速度的二阶渐近解和界面内波波面位移的二阶Stokes解,揭示了波流之间的非线性相互作用和界面波解之间的非线性相互作用。     

南海北部孤立波对近岸波动的响应

On the basis of the time series observations from a temperature chain and an acoustic Doppler current profiler on the continental shelf of the northern South China Sea, a sequence of internal solitary waves （ISWs） and background waves （BWs, including internal tides and near-inertial waves） on the continental shelf were captured simultaneously after the transit of Typhoon Neast in October 2011. These measurementsprovided a unique opportunity to explore the influence of BWs on the ISWs. The BWs appeared a conversion on the current strength and vertical mode structure during the observational period. The BWs were dominated by weak and mode-one waves before October 2 and then turned to strong and high-mode waves after that time. Meanwhile, the ISWs displayed different wave structures before and after October 2, which was closely related to BWs＇ changes. According to the current profiles of BWs, the high-mode wave structure with strong current could significantly strengthen the vertical shear of ISWs in the near-surface layer and promote the breaking of ISWs, and thus it may play an important role in affecting the background current condition.     

Application of desingularized approach to water wave propagation over three-dimensional topography

A numerical approach based on desingularized boundary element method and mixed Eulerian–Lagrangian formulation [Zhang et al., 2006. Wave propagation in a fully nonlinear numerical wave tank: a desingularized method. Ocean Engineering 33, 2310–2331] is extended to solve the water wave propagation over arbitrary topography in a three-dimensional wave tank. A robust damping layer applicable for regular and irregular incident waves is employed to minimize the outgoing wave reflection back into the wave tank. Numerical results on the propagation of regular and irregular incident waves over the flat bottom and linear incident waves over an elliptical shoal show good concurrence with the corresponding analytical solutions and experimental data.     

An experimental study of stratified mixing caused by internal solitary waves in a two-layered fluid system over variable seabed topography

Stratified mixing is observed in a wave flume on an internal solitary wave (ISW) of depression or elevation type propagating over a submarine ridge. The submarine ridges, which comprise the seabed topography, are either semicircular or triangular. Tests are performed in a series of combinations of submarine ridges with different heights and ISW in different amplitudes within a two-layer fluid system. When the thickness of the top layer is less than that of the lower layer (i.e., H₁<H₂), a depression-type ISW may produce a strong hydraulic jump with downwards motion and continuous eddy diffusion. During diffusion, the leading profile of the ISW transforms a wrapped vortex on the front face of the ridge, and a vortex separation at the apex of the ridge. Meanwhile, an elevation-type ISW causes a vortex in the lee of a submarine ridge, which resembles a surface solitary wave in terms of wave transmission process. The degree of wave-obstacle interaction is determined by energy loss, which is induced by submarine ridge blockage. The experiment results suggest that degree of blocking can be applied to classify various degrees of ISW-obstacle encounter in the stratified two-layer fluid system.     

抛物型对称海脊引导波完整解析理论

越洋海啸能够被大洋海脊引导并沿海脊传播至远场地区,虽然传播速度较慢,但携带较大的能量,会对远场地区造成灾害影响,相关研究对于提高海啸传播特性本质的认识具有重要意义。本文基于线性长波方程,推导出了抛物型对称海脊引导波完整解析理论。研究表明对称型海脊同时存在对称与反对称引导波,其自由水面波动可以表示为虚宗量Bessel函数形式。利用海脊中心对称条件给出了描述其频率与波数的频散关系。基于所提理论进一步分析讨论了引导波的运动特性,包括频散关系、波速度、能量传播速度与波面空间分布等。本研究为揭示地形坡度由脊顶至两侧逐渐增加的海脊引导波运动特征,预测实际越洋海啸中最具威胁性海啸波的到达时间提供了理论依据。