摘要:在考虑到海洋内波的二维效应时,含有横向微扰的二维模型可以较好地模拟海洋内波的传播,因此为更好地描述海洋内波这一复杂非线性现象的传播特性,利用二维柱KadomtsevPetviashvili模型研究了海洋内波非线性相互作用的传播特征。首先,利用解析法获得该二维模型的双抛物线型孤子解,同时模拟研究了南海东沙群岛的海洋内波的弹性碰撞及随时间变化的传播特性;其次,基于不同的海域背景参数分析研究了上升型和下降型抛物线海洋内波弹性碰撞引起的海洋表层流速的变化。经比较分析,本文模拟结果与已有的遥感卫星图像较一致,因此可为进一步深入研究海洋内波的传播特性、动力学机制解及内波参数反演等提供一定的借鉴。
关键词:双抛物线型海洋内孤立波;柱KadomtsevPetviashvili模型;解析解
中图分类号:P731.24文献标志码:A
利用水平一维动力学模型如Kortewegde Vries(KdV)方程,科学家们已经对世界多个海域的海洋内波的传播特征进行模拟研究,包括安德曼海、苏禄海、巴士海峡、Sulu海、New York湾和中国南海海域、中国台湾东北海域[19],通过实测资料分析、实验室实验、数值模拟、遥感图像分析和解析研究等方法,分析了各海域的海洋内波空间分布特征及季节变化、内波的起源与激发、极性转换等现象。另一方面,在考虑海洋内波的二维效应时,为更好地描述海洋内波这一复杂的非线性现象的传播特性,研究人员逐渐利用相应的二维模型对海洋内波的传播进行模拟和海洋动力学参数反演研究。通过众多研究表明利用二维KadomtsevPetviashvili(KP)方程可较好地描述海洋内波的传播特征[1012]。
1柱KP模型及孤子解
一般地,描述分层水体结构中海洋内波的传播特征二维柱KP模型可表示为
(ut+αuux+βuxxx+c0ux)x+12tux+c0t2uyy=0,(1)
其中:α=3c0(h1-h2)/(2h2/h2)为非线性系数;β=cuh2/h2/6为频散系数;c0=[gΔρ/ρ0h1h2/(h1+h2)]1/2 为线性波速;g是重力加速度;ρ0是海水的平均密度;Δρ是上下分层的密度差。
由式(1)可知,频散系数β>0,故非线性系数α的正负将直接影响海洋内波的波形。在深水域h1
通过构建变换u=6α/βW[X,Y,T],X=x-c0t,Y=y/c0α和T=αt,则可将式(1)变换为经典的柱KP模型:(WT+6WWX+WXXX)X+12TWX+1T2WYY=0,如文献[13]所示,经典柱KP模型是完全可积的,具有多孤子解、Lax对和无穷守恒律等众多可积性质,因此利用上述变换则可推出式(1)的双孤波解:
(1)的双孤波进行模拟研究,同时分析其在SAR图像上的信号特征变化。
2海洋内波仿真模型
海洋内波的SAR成像物理机制主要包括[1417]:海洋内波的传播引起海洋内部流场及表层流速的变化,形成海面辐聚和辐散;变化的表层流场对海洋表面微尺度波的调制作用改变了海洋表面粗糙度;表面微尺度波对SAR雷达波的调制作用通过上述3个过程的综合调制可获得包含丰富海洋内波信息的微波后向散射图像。在密度稳定层化的海洋中,式(1)的上升型海洋内波所引起的海洋表面波的水平流速为
其中i=1,2表示内波数,而j=1,2则分别对应层化流体上层和下层的水平流速。当入射角为20°~70°时,海面的雷达后向散射以Bragg散射为主,因此Δσ0和σ00的关系可表示为
其中:σ00为归一化的雷达波海面后向散射截面;Δσ0 是内波流场对雷达波后向散射截面的贡献;μ和γ分别为张弛时间、Bragg波的群速度与相位速度比值,而是雷达飞行方向与内波传播方向的夹角。若表示为dB的形式,则式(4)可改写为
即为进行海洋内波SAR成像模拟的公式。
3抛物线型海洋内波仿真模拟
据大量的现场观测与遥感数据证实,在中国南海海域存在丰富海洋内波,且该海域的内波存在较明显的时空特征,如文献[14]所示,在南海北部主要有3个区域存在海洋内波现象,包括吕宋海峡、东沙群岛和海南东北部,传播方向以西向为主,每年的4月-9月为海洋内波出现的高峰期,1月-3月和11月-12月为低峰期。图1为1998年4月28日在中国南海海域观测到的海洋内波现象[18],由图1可看出,在海洋内波经过东沙群岛时,其被分成2列,分别从岛的两侧通过,而在通过该岛后,当两者相遇时却发生了非线性相互作用[17]。根据该区域的条件密度氏和浮力频率的分布图水深图及海底地形变化(见图1 b)-图1 d)),平均海水深度可设为h=472 m,上层混合层深h1=60 m,底层深h2=412 m,海洋上下分层密度差Δρ/ρ=0.003 6,在海岸附近,最终的底层深为h2=40 m。
图1中国南海海域观测到的海洋内波现象及海洋背景参数
Fig.1Internal waves in the Dongsha Islands of South China Sea and ocean background parameters
在较深的海域,图1 d)中A1处的海域深度h1=60 m和h2=412 m具有下降型海洋内波,其在SAR图像上为明暗相间条纹特征。图2为下降型双海洋内波在某一时刻的弹性碰撞过程,与图1 a)相比,可以发现两者具有较好的一致性,2个双抛物线型海洋内波在弹性碰撞时除有微小相位移动外在交换角色后仍能保持原有的特性不变。图3为下降型抛物线海洋内波非线性相互作用引起的海洋表层流速的变化特征(暗亮相间的条纹),与在SAR图像上的特征(亮暗相间条纹)相反,图3 c)和图3 d)为海洋内波弹性碰撞引起的海洋表层流速变化分别在y=0和y=10时的信号特征。
随着时间的推移,抛物线型海洋内波的弯曲程度也在不断发生变化,如图4所示。
在较浅的海域,如图1 d)中A2处的海域水深约300 m,混合层深h1=160 m,底层深度h2=h-h1=140 m,故存在上升型海洋内波,其在SAR图像上为暗亮相间条纹特征。图5为双抛物线上升型海洋内波弹性碰撞在某一时刻的传播特征及在y=-5时截面特征。图6为上升型抛物线海洋内波弹性碰撞引起的海洋表层流速的变化特征(明暗相间的条纹),与在SAR图像上的特征相反,图6 b)和图6 c)为上升型双抛物线海洋内波弹性碰撞引起的海洋表层流速变化分别在y=-5和y=21时的信号特征。
4结论
本文主要是利用二维柱KP模型模拟了南海东沙群岛的双抛物线型海洋内波的传播特征,同时分析了上升型和下降型抛物线海洋内波非线性相互作用所引起的海洋表层流速的变化。由于缺乏有效历史观测数据和海洋环境参数,故在本文中只是初步的定性分析讨论了双抛物线型海洋内波弹性碰撞的传播特征,而有关的基于二维柱KP模型的海洋内波振幅、传播速度等动力学参数的反演还需要在下一步工作中进行深入研究,只希望本文的研究可为海洋内波进一步的研究提供一定的参考。
参考文献/References:
[1]LIU A K, APEL J R, HOLBROOK J R. Nonlinear internal wave evolution in the Sulu Sea [J]. Journal of Physical Oceanography, 1985, 15: 16131624.
[2]LIU A K. Analysis of nonlinear internal waves in the New York Bight[J].Journal of Geophysical Research, 1988, 93: 1231712329.
[3]LIU A K, CHANG Y S, HSU M K, et al. Evolution of nonlinear internal waves in China Seas [J]. Journal of Geophysical Research, 1998, 103: 79958008.
[4]ZHAO Z X, KLEMAS V, ZHENG Q A, et al. Remote sensing evidence for baroclinic tide origin of internal solitary waves in the northeastern South China Sea [J]. Geophysical Research Letters, 2004, 31: 14.
[5]CAI S Q, XIE J S, XU J X, et al. Monthly variation of some parameters about internal solitary waves in the South China Sea [J]. DeepSea Research I, 2014, 84:7385.
[6]吕海滨, 申辉, 何宜军. 背景流和混合层深度对内波表面信号的影响 [J]. 海洋科学, 2013, 37(9):2023.
LYU Haibin, SHEN Hui, HE Yijun. The effects of background current and pycnocline depth on surface signals of internal waves [J]. Marine Sciences, 2013, 37(9):2023.
[7]刘永军, 宋金宝, 黄小峰. 非传统近似下海洋内波的一类WKB近似解 [J]. 海洋科学, 2011, 35(1):8287.
LIU Yongjun,SONG Jinbao,HUANG Xiaofeng. A type of WKB solutions for internal waves under the nontraditional approximation [J]. Marine Sciences, 2011, 35(1):8287.
[8]GUO C, CHEN X. A review of internal solitary wave dynamics in the northern South China Sea [J]. Progress in Oceanography, 2014, 121: 723.
[9]KUMAR S, SINGH K, GUPTA R K. Dynamics of internal waves in a stratified ocean modeled by the extended Gardner equation with timedependent coefficients [J]. Ocean Engineering, 2013, 70:8187.
[10]李群, 孙丽, 徐肇廷. 内孤立波波波相互作用的数值模拟 [J].中国海洋大学学报, 2007, 37(6):868872.
LI Qun, SUN Li, XU Zhaoting. Modeling the wavewave interaction between internal solitary waves[J]. Periodical of Ocean University of China, 2007, 37(6):868872.
[11]MENG J M, ZHANG Z L, ZHAO J S. The simulation of the SAR image of internal solitary waves in Alboran Sea [J]. Journal of Hydrodynamics Series B, 2001(3):8892.
[12]PIERINI S. A model for the alboran sea internal solitary waves [J]. Journal of Physical Oceanography,1989,19:722755.
[13]LI J, ZHANG H Q, XU T, et al. Symbolic computation on the multisolitonlike solutions of the cylindrical KadomtsevPetviashvili equation from dusty plasmas [J]. Journal of Physics Mathematical and Theoretical,2007, 40:76437657.
[14]蔡树群, 甘子钧, 龙小敏. 南海北部孤立子内波的一些特征和演变 [J]. 科学通报,2001,46: 12451250.
CAI Shuqun, GAN Zijun, LONG Xiaomin.Some characteristics and evolution of the internal soliton in the northern South China Sea [J].Chinese Science Bulletin, 2001, 46: 12451250.
[15]李海艳.利用合成孔径雷达研究海洋内波 [D].青岛:中国海洋大学, 2004.
LI Haiyan.Studying Ocean Internal Waves with SAR [D].Qingdao:Ocean University of China, 2004.
[16]杨劲松.合成孔径雷达海面风场、海浪和内波遥感技术 [M]. 北京:海洋出版社,2005.
YANG Jinsong. Ocean Surface Wind Wave and Internal Wave Remote Sensing Technology Based on Synthetic Aperture Radar [M]. Beijing: Ocean Press, 2005.
[17]申辉.海洋内波的遥感与数值模拟研究[D].青岛:中国科学院海洋研究所, 2005.
SHEN Hui. Remote Sensing and Numerical Modeling of Oceanic Internal Waves [D]. Qingdao:Institute of Oceanology of Chinese Academy of Sciences, 2005.
[18]HSU M K, LIU A K. Nonlinear internal waves in the South China Sea [J]. Canadian Journal of Remote Sensing, 2000, 26: 72