地球物理学进展  2014, Vol. 29 Issue (1): 57-60   PDF    
Coriolis力作用下的β效应与层结效应的Rossby孤立波
刘全生1, 宋健1,2, 杨联贵1     
1. 内蒙古大学数学科学学院, 呼和浩特 010021;
2. 内蒙古工业大学理学院, 呼和浩特 010051
摘要:从包含完整Coriolis力的大气运动方程组出发, 利用半地转近似导出了β效应、层结效应和地球旋转水平分量fh共同作用下的非线性Rossby波满足KdV方程以及KdV-mKdV方程.结果表明β效应、层结效应和地球旋转水平分量对Rossby的作用.
关键词Coriolis力     β效应     层结效应     KdV和KdV-mKdV方程    
Solitary Rossby waves with beta effect and stratification effect in the Coriolis force
LIU Quan-sheng1, Song Jian1,2, YANG Lian-gui1    
1. School of Mathematical Science, Inner Mongolia University, Hohhot 010021, China;
2. College of Sciences, Inner Mongolia University of Technology, Hohhot 010051, China
Abstract: The nonlinear Rossby waves with bate effect and stratification effect in a horizontal component of the earth rotation are studied by using the semi-geostrophic approximate, which have a constant of the Coriolis parameter. The solitary Rossby waves described by KdV and KdV-mKdV equation are derived. The results show the horizontal component, beta effect and stratification effect can be important to the Rossby waves.
Key words: Coriolis force     β effect     stratification effect     KdV and KdV-mKdV equation    

0 引 言

地球流体中,地球旋转对许多现象都有重要的影响,它是通过Coriolis力作用而产生,传统近似是采用准地转近似,即在控制方程中省略了地球旋转作用的水平分量,但从动力学角度而言,一直有争议(Philips 19661968Veronis 1968Wangsness 1970Burger 1991White et al.,1995)对纬向动量平衡的尺度分析说明,对于行星大尺度的大气运动保留地球旋转水平分量的作用是合理的. 现在Coriolis力参数水平分量在地球流体中的重要性得到了广泛的重视(Sun, 1995; Leibovich, 1985; Draghici, 1987; Raymond, 2001; Kasahara, 2003; Durran, 2004; Grimshaw, 1975; Hoskins, 1975; 刘式适等, 1987; 赵强等, 2004). 赵强等人(2004,2008)采用地转近似(Hoskins, 1975, 刘式适等, 1987)给出了完整Coriolis力作用下非线性Rossby波的精确解,杨洁等 (2010)在Boussinesq近似下考虑了无摩擦的Burger模式给出了完整Corilois力和热源影响下超长波的解析解. 他们虽然考虑了Coriolis参数的水平分量对地球物理流体的影响,但忽视了β平面近似模式下β参数随纬度变量y的变化,同时也没有考虑层结效应对Rossby波的影响. 赵强等 (2006)讨论了雅可比椭圆函数在大气和海洋动力学的应用,巢纪平等 (2008)解释了热带扰动在大尺度经圈中的行为,张立凤等 (2011a,2011b)讨论了真实基流中Rossby演变特点,这些说明了大尺度非线性运动的大气特征.随着人类文明的发展加剧了全球变暖和全球海平面大幅度的升高,这都说明地形会随时间发生变化.张永垂等 (2010)利用15年(1993-2007年)月平均的海表面高度(SSH)异常资料,分析了Rossby波的传播作用在中纬度海域的副热带环流中西部和夏威夷岛以东起着重要作用, 刘秦玉等 (2007)根据对卫星观测的海平面高度资料的分析,这两支向东的逆流区Rossby波的特性不同:副热带逆流区70~210天周期振荡对应的Rossby波西传过程中增幅,在台湾以东振幅达到最大;而在夏威夷背风逆流区,70~210天周期Rossby波在西传过程中不出现增幅现象.本文在完整Coriolis作用下,采用半地转方法给出了β效应与层结效应的Rossby孤立波的演变.

1 控制方程

在Boussinesq近似下,绝热无摩擦的完整Coriolis力大气运动方程组为

其中0为位温的典型值,N(z)是Brunt-Vaisala频率,它是表征层结的作用, 是位温,是势函数,f是Coriolis力参数,g是重力加速度,fhCoriolis参数水平分量,fh=2Ωcos是纬度.

Grimshaw(1975)指出方程组(1)中fh是常数,f随纬度变化,并且假定Coriolis参数随纬度变化较小,即f=f0+β(y)y, |β(y)y| <0,其中f0是常数,β(y)是纬度变量y的函数.

由文献(赵强等,2008)方程组(1)化为

方程组(2)中包括了重力惯性波和Rossby波两类波动,为了克服准地转模式下波特解不能突出非线性作用的缺点,采用半地转近似. 由于在半地转近似中,非线性平流项的平流风为实际风,被平流风为地转风,因此半地转近似部分考虑了非地转平流效应,它滤去了高频快波,同时又能反映Rossby波的非线性特征,方程组(2)化为

其中,,它们分别是地转风的x,y方向的分量.

从方程组(3)可以得到了大气中地转水平分量、β效应与层结效应下的非线性Rossby波的方程组为

方程组(4)中第一式为地转水平分量与β效应下的涡度方程,第二式是地转水平分量、β效应下的散度方程,对于大尺度运动已将涡度方程的 项中的v和散度方程的 项中的u分别用地转风vg和ug代替(赵强等,2008).

将方程组(4)中的涡度方程化为

方程(5)中算子.

与文献(赵强等,2008)中的方程组(5)相比,在方程组(4)中的涡度方程与散度方程考虑了β效应,即β是纬度变量y的函数与层结效应N(z). 方程组(4)式是完整Coriolis力作用下β效应与层结效应的非线性Rossby波的基本数学模式.

2 KdV、KdV-mKdV方程的导出

设方程组(4)的波动解为(赵强等,2008)

这里,θ是位相函数,k,l,n分别为x,y,z方向的波数,ω是圆频率. 将(6)式代入方程组(4)得

将(7)式中的前三式积分一次(积分常数取为零)得

将(8)式中第三式代入(8)中的第一式和(7)中的第四式,得到

消去(9)式中的Φ,并且在nW+ω≠0的条件下得

其中,F(W)是W的非线性函数,W=0是方程(10)的平衡点,将F(W)在W=0处作Taylor级数展开得

2.1 线性Rossby波

如果(11)式右端只取第一项,得到

(12)式是描述线性Rossby波的常微分方程,它若表示波动,必须要求

(13)式为大气中线性Rossby波的频散关系,它与β效应、层结效应和水平分量fh有关.如果P′(θ),I(θ)为常数,(13)式是大气中三维线性Rossby波的频散关系,它只与地球旋转水平分量有关(赵强等,2008;杨洁等,2010).

2.2 带有强迫项的KdV方程

(11)式的右端取前两项,并对θ求导一次,得

这里对θ的一阶导数. 方程(14)是带有强迫项的KdV方程,它是描述非线性Rossby波在β效应、层结效应和完整Coriolis力作用下满足带有强迫项的KdV方程. 如果P′(θ),I(θ)为常数,则强迫项Q1(θ)=0,此时方程(14)为KdV方程,其孤立波解为(赵强等,2008;杨洁等, 2010)

2.3 带有强迫项的KdV-mKdV方程

对(11)式右端取前三项,并对θ求导一次,得

对θ的一阶导数. 方程(16)是带有强迫项的KdV-mKdV方程,它刻画了非线性Rossby波在β效应、层结效应和完整Coriolis力作用下满足非线性KdV-mKdV方程. 如果P′(θ),I(θ)为常数,则强迫项Q2(θ)=0,此时方程(16)为KdV-mKdV方程,其类孤立子解为(赵强等,2008;杨洁等,2010)

其中

3 结 论

从包含完整Coriolis力和β效应、层结效应的绝热无摩擦的大气运动方程组出发,在半地转近似下利用Taylor级数展开法,得到了非线性Rossby波动满足带有强迫项的KdV方程和KdV-mKdV方程,并给出了没有强迫项时的孤立波解及其类孤立子解. 分析表明,Coriolis参数分量fh、β效应和层结效应对非线性Rossby有影响.

参考文献
[1] Burger A P. 1991. The potential vorticity equation: from planetary to small scale[J]. Tellus A, 43(3): 191-197, doi:10.1034/j.1600-0870. 1991. t01-2-00002. x.
[2] Chao J P, Xu Z. 2008. The behavior of tropical disturbance in large-scale meridional flow[J]. Chinese J. Geophys.   (in Chinese), 51(6): 1657-1662.
[3] Draghici I. 1987. Non-hydrostatic Coriolis effects in an isentropic coordinate frame[J]. Meteor. Hydrol., 17: 45-54.
[4] Durran D, Bretherton C. 2004. Comments on “the roles of the horizontal component of the Earth’s angular velocity in nonhydrostatic linear models”[J]. Journal of the Atmospheric Sciences, 61(15):1982-1986,   doi:10.1175/1520-0469(2004)061<1982:COTROT>20.CO;2.
[5] Grimshaw R H J. 1975. A note on the β-plane approximation[J]. Tellus, 27(4): 351-357,   doi:10.1111/j.2153-3490.1975.tb01685.x.
[6] Hoskins B J. 1975. The geostrophic momentum approximation and the semi-geostrophic equations[J]. Journal of Atmospheric Sciences, 32(2): 233-242,    doi:10.1175/1520-0469(1975)032<0233: TGMAAT> 2. 0.CO;2.
[7] Kasahara A. 2003. The roles of the horizontal component of the earth’s angular velocity in nonhydrostatic linear models[J]. J. Atmos. Sci., 60(8): 1085-1095,    doi:10.1175/1520-0469(2003)60<1085: TROTHC>2.0.CO;2.
[8] Leibovich S, Lele S K. 1985. The influence of the horizontal component of Earth’s angular velocity on the instability of the Ekman layer[J]. J. Fluid Mech., 150: 41-87, doi:10.  1017/S0022112085000039.
[9] Liu Q Y, Li L J. 2007. Baroclinic stability of oceanic long Rossby wave in the North Pacific Subtropical Countercurrent[J]. Chinese J. Geophys.   (in Chinese), 50(1): 83-91.
[10] Liu S K, Liu S D. 1987. Nonlinear waves with semi-geostrophic flow[J]. Acta Meteor. Sin.   (in Chinese), 45(3): 257-266.
[11] Philips N A. 1966. The equations of motion for a shallow rotating atmosphere and the traditional approximation[J]. J. Atmos. Sci., 23(5): 626-628,    doi:10.1175/1520-0469(1966)023<0626: TEOMFA>2. 0.CO;2.
[12] Philips N A, Reply to G. 1968. Veronis’s comments on Phillips(1966)[J]. J. Atmos. Sci., 25(6): 1155-1157,    doi:10.1175/1520-0469(1968)025<1155: R>2. 0.CO;2.
[13] Raymond W H. 2001. Kelvin waves: rotationally induced circulations[J]. Dyn. Atmos. Oceans., 34(1): 23-43, doi:10.  1016/S0377-0265(01)00058-6.
[14] Sun W Y. 1995. Unsymmetrical symmetric instability[J]. Quarterly Journal of the Royal Meteorological Society, 121(522): 419-431, doi:10.1002/qj.   49712152209.
[15] Veronis G. 1968. Comments on Phillip’s (1966) proposed simplification of the equations of motion for a shallow rotating atmosphere[J]. J. Atmos. Sci., 25(6): 1154-1155,    doi:10.1175/1520-0469(1968)025<1154: COPPSO>2. 0.CO;2.
[16] Wangsness R K. 1970. Comments on “The equations of motion for a shallow rotating atmosphere and the traditional approximation”[J]. J. Atmos. Sci., 27(3): 504-506,    doi:10.1175/1520-0469(1970)027<0504: COEOMF>2. 0.CO;2.
[17] White A A, Bromley R A. 1995. Dynamically consistent, quasi-hydrostatic equations for global models with a complete representation of the Coriolis force[J]. Quart. J. Roy. Meteor. Soc., 121(522): 339-418, doi:10.1002/qj.   49712152208.
[18] Yang J, Zhao Q. 2010. Analytical solutions to ultra-long waves with the complete Coriolis force and heating[J]. Acta Phys. Sin.   (in Chinese), 59(2): 750-753.
[19] Zhao Q, Liu S K. 2004. Effects of the horizontal component of the Earth’s rotation on Rossby waves[J]. Chinese J. Atmos. Sci. (in Chinese), 28(2): 146-150,    doi:10.3878/j.issn. 1006-9895. 2004. 01.14.
[20] Zhang Y C, Zhang L F, Wang Y G. 2010.Interannual sea level variability in the North Pacific Ocean and its mechanisms[J]. Chinese J. Geophys. (in Chinese), 53(2): 247-255,    doi:10.3969/j.issn. 0001-5733. 2010.02.003.
[21] Zhao Q, Yu X. 2008. Exact solutions to the nonlinear Rossby waves with a complete representation of the Coriolis force[J]. Chinese J.   Geophys (in Chinese), 51(5): 1304-1308.
[22] Zhang L F, Zhang L, Li C Y. 2011. The evolution characteristic of nonlinear Rossby wave in the real basic flowⅠ: Numerical model design[J]. Chinese J. Geophys. (in Chinese), 54(3): 629-633,    doi:10.3969/j.issn. 0001-5733. 2011. 03.001.
[23] Zhang L F, Zhang L, Li C Y. 2011. The evolution characteristic of nonlinear Rossby wave in the real basic flow Ⅱ: The energy and structure evolution of Rossby wave and the effect of the initial field on it[J]. Chinese J. Geophys. (in Chinese), 54(4): 886-895,    doi:10.3969/j.issn. 0001-5733. 2011. 04.002.
[24] Zhao Q, Liu S K. 2006. Application of Jacobi elliptic functions in the atmospheric and oceanic dynamics: studies on two-dimensional nonlinear Rossby waves[J]. Chinese J. Geophys.   (in Chinese), 49(4): 965-970.
[25] 巢纪平, 徐昭. 2008. 热带扰动在大尺度经圈中的行为[J].   地球物理学报, 51(6): 1657-1662.
[26] 刘秦玉, 李丽娟. 2007. 北太平洋副热带向东逆流区Rossby波斜压稳定性[J].   地球物理学报, 50(1): 83-91.
[27] 刘式适, 刘式达. 1987. 半地转近似下的非线性波[J].   气象学报, 45(3): 257-266.
[28] 杨洁, 赵强. 2010. 有完整Coriolis 力和热源影响下超长波的解析解[J].   物理学报, 59(2): 750-753.
[29] 赵强, 刘式适. 2004. 地球旋转水平分量对Rossby波的影响[J]. 大气科学, 28(2): 146-150,    doi:10.3878/j.issn. 1006-9895. 2004. 01.14.
[30] 张永垂, 张立凤, 王业桂. 2010. 北太平洋海表面高度的年际变化及其机制[J]. 地球物理学报, 53(2): 247-255,    doi:10.3969/j.issn. 0001-5733. 2010.02.003.
[31] 赵强, 于鑫. 2008. 完整Coriolis力作用下非线性Rossby波的精确解[J].   地球物理学报, 51(5): 1304-1308.
[32] 张立凤, 张亮, 李崇银. 2011a. 真实基流中非线性Rossby波演变特征(一): 数值模式设计[J]. 地球物理学, 54(3): 629-633,    doi:10.3969/j.issn. 0001-5733. 2011. 03.001.
[33] 张立凤, 张亮, 李崇银. 2011b. 真实基流中非线性Rossby波演变特征(二): 能量、结构演变及初始场影响[J]. 地球物理学报, 54(4):886-895,  doi:10.3969/j.issn.0001-5733.2011.04.002.
[34] 赵强, 刘式适. 2006. 雅可比椭圆函数在大气和海洋动力学中的应用: 二维非线性Rossby波研究[J].  地球物理学报, 49(4): 965-970.