林萬(wàn)濤 林一驊 石蘭芳 莫嘉琪

1)(中國(guó)科學(xué)院大氣物理研究所,大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100029)

2)(南京信息工程大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,南京 210044)

3)(安徽師范大學(xué)數(shù)學(xué)計(jì)算機(jī)科學(xué)學(xué)院,蕪湖 241003)

(2013年2月14日收到;2013年3月22日收到修改稿)

1 引言

厄爾尼諾和南方濤動(dòng)(ENSO)分別是發(fā)生在熱帶大氣和海洋中的異常事件,它嚴(yán)重地影響全球各地和區(qū)預(yù)氣防為候當(dāng)和前生學(xué)態(tài)術(shù)等界方所面關(guān)的注變化[1-,9對(duì)].它海的洋規(guī)和律大的氣研之究間的交互作用影響著氣候的波動(dòng),這對(duì)種交互作用的研究已經(jīng)集中在厄爾尼諾-南方濤動(dòng)現(xiàn)象上,它每隔3到4年周期地出現(xiàn),并主要發(fā)生在熱帶太平洋區(qū)域.厄爾尼諾和南方濤動(dòng)的振蕩性態(tài)是海-氣流動(dòng)的正負(fù)兩種反饋的結(jié)果,這兩種反饋決定了海表溫度(SST)的變化,并與南方濤動(dòng)流動(dòng)強(qiáng)度有關(guān),導(dǎo)致了弱信風(fēng)沿著赤道行進(jìn).弱信風(fēng)驅(qū)動(dòng)著海洋流動(dòng)變化,加強(qiáng)了SST的異常.海洋-大氣的正反饋與耦需合研的究不東、穩(wěn)定西性太平導(dǎo)致洋了振子赤模道太型平的異洋溫常度關(guān)變系.化許,因多此學(xué)者已使用不同的方法對(duì)ENSO局部和整體的性態(tài)做了多方位的討論[1-6].對(duì)于復(fù)雜的全球海-氣耦合模型,通過(guò)簡(jiǎn)化的海-氣非線性相互作用物理過(guò)程得到的振子概念模型能更容易描述海-氣耦合過(guò)程的本質(zhì)和物理機(jī)理,從而能描述ENSO的某些重要物理現(xiàn)象.因此許多學(xué)者提出了研究ENSO的各種振子形式的動(dòng)力系統(tǒng)模型[7-12],如“時(shí)滯振子”,“平流-反射振子”,“西太平洋振子”,“充電-放電振子”等理論的研究.本文研究一類(lèi)ENSO海-氣耦合振子震動(dòng)動(dòng)力學(xué)模型[11].近來(lái),許多學(xué)者研究了非線性奇攝動(dòng)問(wèn)題,一些近似方法被優(yōu)化,包括邊界層法、平均化法、匹配近似展開(kāi)方法和多重尺度法等[13-17].作者等也利用攝動(dòng)理論等方法研究了一類(lèi)非線性問(wèn)題[18-30].本文是利用攝動(dòng)理論來(lái)討論一類(lèi)ENSO海-氣耦合振子.

2 ENSO耦合振子模型

考慮如下ENSO耦合振子動(dòng)力學(xué)模型[7-9]：

其中T表示赤道東太平洋的SST距平,h表示赤道西太平洋的溫躍層厚度距平,C,D,E,Rh表示正參數(shù),它們的詳細(xì)定義和物理意義參見(jiàn)文獻(xiàn)[7,8],ε為小的正參數(shù).本文是求ENSO耦合振子動(dòng)力學(xué)模型(1),(2)的震蕩近似解.

首先對(duì)方程(1)兩邊進(jìn)行d/d t+Rh運(yùn)算,考慮到(2)式,當(dāng)C＞Rh,DE＞CRh時(shí),可得

其中

由(3)式,對(duì)應(yīng)的方程：

這就是van der Pol方程.由van der Pol方程的性態(tài)知,方程(4)當(dāng)T＜Tc為負(fù)阻尼方程,當(dāng)T＞Tc為正阻尼方程.此外,方程(4)可改寫(xiě)為如下等價(jià)的平面系統(tǒng)：

不難看出,系統(tǒng)(5)有一個(gè)平衡態(tài)(T?,U?)=(0,0).其對(duì)應(yīng)的特征方程和特征根為

由此可知[11],當(dāng)μ2＜4ω02時(shí),即在弱阻尼的情況下,平衡態(tài)(T?,U?)=(0,0)是一個(gè)不平衡的焦點(diǎn);當(dāng)μ2＞4ω02時(shí),即在強(qiáng)阻尼的情況下,平衡態(tài)(T?,U?)=(0,0)是一個(gè)不平衡的結(jié)點(diǎn).因此,這時(shí)海-氣耦合振子當(dāng)氣候稍有偏離平衡態(tài)時(shí),將會(huì)處于不穩(wěn)定狀態(tài),并且還不難看出,當(dāng)μ=0(無(wú)阻尼的情形),van der Pol方程(3)是一個(gè)Hamilton系統(tǒng).當(dāng)ε很小時(shí),在無(wú)阻尼情形時(shí),海表溫度距將做周期振蕩.當(dāng)μ/=0的情形,在弱阻尼的條件下,存在兩個(gè)平衡態(tài)：一個(gè)平衡態(tài)就是(T?,U?)=(0,0),便是不穩(wěn)定的焦點(diǎn);另一個(gè)平衡態(tài)就是在相平面(T,U)上,是一個(gè)閉合的極限環(huán).它是van der Pol方程在弱阻尼的情況下求得的孤立波周期解(周期吸引子),即當(dāng)運(yùn)動(dòng)脫離平衡態(tài)后,有可能進(jìn)入一個(gè)周期的運(yùn)動(dòng),其周期為2π/ω0的氣候震蕩狀態(tài).

3 ENSO模型的振蕩近似解

對(duì)于小的參數(shù)ε,在弱阻尼的條件下,我們用Krylov-Boglinbov方法來(lái)求弱阻尼方程(3)振蕩近似解.方程(3)可改寫(xiě)為

當(dāng)參數(shù)ε=0時(shí),方程(6)的解為

其中a,φ為常數(shù).對(duì)于ε/=0時(shí),我們采用廣義的“常數(shù)變易法”.將方程的解看作仍由方程(6)的形式,但a,φ為時(shí)間t的函數(shù),且滿(mǎn)足條件

現(xiàn)對(duì)(7)式求導(dǎo)

比較(8),(9)式得

再對(duì)(8)式求導(dǎo),有

將上式代入方程(6),并考慮到(7)式,得

從(10),(11)式解出 d a/d t,dφ/d t,得到

為了解(12),(13)式,注意到這些方程右邊關(guān)于變量ψ是周期函數(shù),所以有

因?yàn)棣攀切×?a和φ是時(shí)間t的慢變函數(shù),所以在時(shí)間周期2π/ω0區(qū)間內(nèi)它們的變化是很小的.把這些方程右邊的a和φ當(dāng)作常數(shù),在一個(gè)周期[t,t+(2π/ω0)]上求 (12),(13)式的平均值,得到

其中

即

于是由(14)—(16)式知

積分(17),(18)式,得到

其中a0和φ0為常數(shù).將(19),(20)式代入(7)式,我們便得到弱阻尼van der Pol方程(4)的震蕩近似解：

再將(21)式代入方程(2),得

由方程(22),有解

故(21),(23)式便是ENSO耦合振子動(dòng)力學(xué)模型(1),(2)的震蕩近似解.用同樣的方法,還可繼續(xù)得到ENSO耦合振子動(dòng)力學(xué)模型更高次的震蕩近似解.

4 舉例

現(xiàn)舉例說(shuō)明ENSO耦合振子動(dòng)力學(xué)模型解的震蕩性態(tài).為了方便起見(jiàn),我們假設(shè)如下無(wú)量綱方程：

比較模型(1),(2)知C=2,D=Rh=1,E=3,μ=1,ω0=1.現(xiàn)求模型(24),(25)式的近似解.由(24),(25)式,對(duì)應(yīng)于方程(5)為

由(8)式,設(shè)方程(26)的解為

由(14),(15)式

其中

即有

于是由(14)—(16)式知

積分(31),(32)式,得到

其中a0和φ0為常數(shù).將(33),(34)式代入(27)式,我們便得到方程(26)的震蕩近似解：

將(35)式代入方程(25),得

由方程(36),有解

其中a0,b0為常數(shù).故(35),(37)式便是ENSO耦合振子動(dòng)力學(xué)模型(24),(25)的震蕩近似解.現(xiàn)設(shè)T(0)=h(0)=1,這時(shí)由解(35),(37)式得

由方程(36),有解

當(dāng)在初始條件T(0)=h(0)=1,ε=0.5下,震蕩近似解(38),(39)曲線圖見(jiàn)圖1和圖2所示.由圖1和圖2可以看出,ENSO耦合振子動(dòng)力學(xué)模型(24),(25)具有震蕩型的近似解.

5 結(jié)論

由于大氣物理中的復(fù)雜性,需要建立它的基本模型系統(tǒng)的方程,并去求解它.參數(shù)變值法和平均法是一個(gè)簡(jiǎn)單而有效的方法.它類(lèi)同于變系數(shù)非齊次常微分方程的常數(shù)變易法,但比常數(shù)變易法更深入,加上了取變量函數(shù)平均值的過(guò)程,以達(dá)到求得近似解的目的.這種方法對(duì)消去解的長(zhǎng)期項(xiàng)的有界解用相應(yīng)的近似解來(lái)替代很有效,同樣,這種近似解保留了解析表達(dá)式的特點(diǎn),以便可以繼續(xù)對(duì)它使用解析運(yùn)算.因此可以利用它對(duì)模型的相應(yīng)的定性性質(zhì)做進(jìn)一步的探討,這對(duì)進(jìn)一步了解和預(yù)測(cè)模型的有關(guān)性質(zhì)更為方便和實(shí)用,特別是對(duì)大氣物理方面的氣象預(yù)報(bào)等的描述.

圖1 ENSO耦合振子海表溫度T(t)近似震蕩曲線圖

圖2 ENSO耦合振子溫躍層厚度h(t)近似震蕩曲線圖

[1]Feng GL,Dong WJ,Jia X J,Cao H X 2002 Acta Phys.Sin.51 1181(in Chinese)[封國(guó)林,董文杰,賈曉靜,曹鴻興2002物理學(xué)報(bào)51 1181]

[2]Guan X P,He Y H,Fan Z P 2003 Acta Phys.Sin.52 276(in Chinese)[關(guān)新平,何宴輝,范正平2003物理學(xué)報(bào)52 276]

[3]Li CG 2003 Acta Phys.Sin.52 2114(in Chinese)[李春貴2003物理學(xué)報(bào)52 2114]

[4]Li Z,Han CZ 2002 Chin.Phys.11 9

[5]Lin J,Xu Y S2003 Chin.Phys.12 1049

[6]Hu Y Q 2002 Introductory to Atmospheric Thermodynamics and Dynamics—Linear and Nonlinear Thermodynamicsof the Atmospheric Non-equilibrium State(Beijng：Geology Pub.)(in Chinese)[胡隱樵2002大氣熱力動(dòng)力學(xué)導(dǎo)論——大氣非平衡態(tài)線性和非線性熱力學(xué)(北京：地質(zhì)出版社)]

[7]Jin FF 1997 J.Armos.Sci.54 811

[8]Jin FF 1997 J.Armos.Sci.54 830

[9]Wang C 2001 1999 J.Adv.Atmosph.Sci.18 674

[10]Feng G L,Dong W J 2005 Acta Meteo.Sin.63 864(in Chinese)[封國(guó)林,董文杰2005氣象學(xué)報(bào)63 864]

[11]Zhao Q,Liu SK,Liu SD 2012 Acta Phys.Sin.61 220201(in Chinese)[趙強(qiáng),劉式適,劉式達(dá)2012物理學(xué)報(bào)61 220201]

[12]Liu S K,Liu S D 2012 Nonlinear Equations in Physics(Beijing：Peking University Press)(in Chinese)[劉式適,劉式達(dá) 2012物理學(xué)中的非線性方程(北京：北京大學(xué)出版社)]

[13]de Jager E M,Jiang F R 1996 The Theory of Singular Perturbation(Amsterdam：North-Holland Publishing Co.)

[14]Barbu L,Morosanu G 2007 Singularly Perturbed Boundary-Value Problems(Basel：Birkhauserm Verlag AG)

[15]Ramos M 2009 J.Math.Anal.Appl.352 246

[16]D’Aprile T,Pistoia A 2010 J.Diff.Eqs.248 556

[17]Faye L,Frenod E,Seck D 2011 Discrete Contin.Dyn.Sys.29 1001

[18]Lin WT,Ji ZZ,Wang B 2002 Prog.Nat.Sci.12 102(in Chinese)[林萬(wàn)濤,季仲貞,王斌2002自然科學(xué)進(jìn)展12 102]

[19]Lin W T,Mo JQ 2004 Chin.Sci.Bull.48(suppl II)5

[20]Lin W T,Lin Y H,Mo JQ 2012 Chin.Phys.B 21 010204

[21]Shi L F,Ouyang C,Chen L H,Mo J Q 2012 Acta Phys.Sin.61 050203(in Chinese)[石蘭芳,歐陽(yáng)成,陳麗華,莫嘉琪2012物理學(xué)報(bào)61 050203]

[22]Shi L F,Lin W T,Lin Y H,Mo JQ 2013 Acta Phys.Sin.62 010201(in Chinese)石蘭芳,林萬(wàn)濤,林一驊,莫嘉琪2013物理學(xué)報(bào)62 010201]

[23]Mo JQ,Lin Y H,Lin W T,Chen L H 2012 Chin.Geog.Sci.22 42

[24]Mo JQ 2009 Sciencein China G 39 568

[25]Mo JQ,Lin Y H,Lin WT 2010 Acta Phys.Sin.59 6701(in Chinese)[莫嘉琪,林一驊,林萬(wàn)濤2010物理學(xué)報(bào)59 6701]

[26]Mo JQ,Lin W T,Lin Y H 2011 Acta Phys.Sin.60 080202(in Chinese)[莫嘉琪,林萬(wàn)濤,林一驊2011物理學(xué)報(bào)60 080202]

[27]Mo JQ 2011 Acta Phys.Sin.60 090203(in Chinese)[莫嘉琪2011物理學(xué)報(bào)60 090203]

[28]Mo JQ 2010 Commun.Theor.Phys.53 440

[29]Mo JQ,Lin Y H,Lin WT 2010 Chin.Phys.B 19 030202

[30]Mo JQ,Lin WT,Lin Y H 2011 Chin.Phys.B 20 070205