張翔,鄧冰,張銘

背景流與地形對海洋內(nèi)波影響初探

張翔1,鄧冰2,張銘3

（1.海軍海洋水文氣象中心,北京100161;2.北京應(yīng)用氣象研究所，北京100029;3.解放軍理工大學(xué)氣象學(xué)院大氣環(huán)流與短期氣候預(yù)測實(shí)驗(yàn)室,江蘇南京211101）

采用Boussinesq近似，研制了準(zhǔn)三維、非線性、非靜力的斜壓數(shù)值模式。利用該模式做了背景流與海底地形對高頻海洋內(nèi)波影響的數(shù)值試驗(yàn)。本文得到的主要結(jié)論有：各數(shù)值試驗(yàn)初始1小時(shí)內(nèi)，有適應(yīng)過程，并最終形成垂直方向?yàn)閱稳Νh(huán)流的海洋內(nèi)波波包；以后進(jìn)入海洋內(nèi)波的演變過程，并在非線性作用下形成橢圓余弦波波包，并向孤立內(nèi)波演變。當(dāng)無背景流時(shí)，地形會激發(fā)出兩支傳播方向相反的海洋內(nèi)波波包，水深大處波包傳播較慢。當(dāng)存在定常背景流時(shí)，在其與地形共同作用下，在地形處會激發(fā)出較強(qiáng)的內(nèi)波波包，若地形足夠陡峭，則內(nèi)波波包截陷于該陡峭處，造成該處強(qiáng)烈的垂直運(yùn)動。背景流隨時(shí)間的變化會對在臺階地形躍階處的內(nèi)波波包產(chǎn)生不同的影響。

海洋內(nèi)波；地形；背景流

1 引言

海洋內(nèi)波是海洋中普遍存在的一種現(xiàn)象，它對海洋熱鹽的細(xì)微結(jié)構(gòu)和海水能量的垂直交換起著非常重要的作用。自被觀測以來人們對它的研究不斷深入。關(guān)于海洋內(nèi)波生成機(jī)制（生成源）、演變規(guī)律的研究一直倍受關(guān)注，各國學(xué)者對此做了大量的研究工作[1-6]?？偟膩碚f，目前對海洋內(nèi)波的生成源已有了一定的認(rèn)識，但內(nèi)波在每一種生成源中的具體生成、演變機(jī)制仍不十分清楚，需做進(jìn)一步的理論分析和數(shù)值研究。對海洋內(nèi)波的生成機(jī)制和演變規(guī)律進(jìn)行數(shù)值試驗(yàn)是海洋內(nèi)波研究的一個(gè)重要手段。海洋內(nèi)波屬海洋中尺度現(xiàn)象，其空間尺度分布范圍較大。短波長的海洋內(nèi)波其水平尺度約幾公里，甚至更短；波長較長的海洋內(nèi)波可達(dá)50 km以上[7]。關(guān)于內(nèi)波生成和演變的數(shù)值研究,許多作者采用靜力近似或弱非線性約化模式[8-12],其無疑對海洋內(nèi)波研究起到了積極作用，但也存在某些局限。我們利用海洋無粘絕熱的Boussinesq方程，建立了一個(gè)準(zhǔn)三維、非線性、非靜力的斜壓數(shù)值海洋模式。該模式可在微機(jī)平臺上進(jìn)行數(shù)值積分。本文采用該模式初步探討了地形及背景流對高頻海洋內(nèi)波的演變、傳播和結(jié)構(gòu)的影響，以便為高頻海洋內(nèi)波的深入研究打下基礎(chǔ)。

2 數(shù)值模式和計(jì)算方案

2.1 數(shù)值模式

考慮無粘，非絕熱無摩擦的Boussinesq海洋方程組，引入熱力學(xué)變量（密度ρ、溫度T、鹽度S）的基本場（背景場）ρˉ、Tˉ、Sˉ，并將熱力學(xué)變量分為基本場和對基本場的偏差，即引入密度偏差ρ′、溫度偏差T′、鹽度偏差S′；再引入隨z, t變化的潮流基本流uˉ(z,t)、vˉ(z,t)，這里z方向由海底垂直向上指向海面；若再設(shè)擾動沿x方向傳播且在y方向分布是均勻的，即取水平準(zhǔn)三維近似，則對于擾動量有?( )/?y=0；在引入流函數(shù)Ψ后，可得以下該數(shù)值模式的控制方程組：

在此α、γ是海水狀態(tài)方程ρ=ρ0[1-α(T-T0)+γ(S-S0)]中的常數(shù)（ρ0、T0、S0為海水密度、溫度、鹽度的典型值）；因通常有?ρˉ/?z＜0，故N2＞0，N2為層結(jié)參數(shù)。將ρˉ在水平方向展開取線性項(xiàng)后，則N2可近似看成僅為z

對式（1—3）中的水平渦度ζ有：

求出ζ后則解泊松方程（8）可得Ψ，再用公式u′=-?Ψ/?z，w′=?Ψ/?x則能求得u′, w′。

因本文僅研究內(nèi)波，在該問題中可取以下的邊條件。在海面，采用剛蓋近似，即有：

這里H為常數(shù)，其即為最大水深，注意到這里設(shè)z=H處為海面。在海底則考慮了海底地形，并采用階梯地形來逼近實(shí)際地形，這樣就有：

這里zs(x,y)為階梯地形函數(shù)。

2.2 數(shù)值計(jì)算方案

模式水平方向取400個(gè)格點(diǎn)，格距為100 m，其水平范圍為40 km；垂直方向取24層，層距為50 m，垂直范圍為1.2 km。為方便，變量在水平和垂直方向均定義在同一網(wǎng)格點(diǎn)上?？臻g差分取中央差。上邊界采用剛蓋邊界，下邊界考慮階梯地形，水平側(cè)邊界取海綿邊界條件，地形側(cè)壁則用剛壁邊界條件。時(shí)間積分取前差疊代方案，疊代3次，積分時(shí)間步長取為6秒。在求解預(yù)報(bào)方程組（1）的過程中，必須求解泊松方程（8），此時(shí)可用超松弛迭代法來求解，并取松弛因子的值為0.6.。

參考內(nèi)波易發(fā)地區(qū)的情況，本文取兩種類型的海底地形，圖1a、b分別給出了其廓線圖：一種在模式水平范圍的右側(cè)有高0.5 km的海底臺地，另一種則在該范圍中部存在高0.45 km的海底山脈，其底部寬度為2 km。以下分別就有關(guān)試驗(yàn)做出說明并給出試驗(yàn)結(jié)果。

圖1 數(shù)值模擬地形圖

3 無背景流時(shí)數(shù)值試驗(yàn)結(jié)果

3.1 地形I的結(jié)果

取背景流uˉ=0，此時(shí)M2=0，S2=0。層結(jié)參數(shù)采用典型值，取地轉(zhuǎn)參數(shù)f=7×10-5s-1，取初始擾動流函數(shù)場如圖2中0時(shí)所示，擾動垂直波長為1.2 km，水平波長為5 km（流函數(shù)環(huán)流圈水平尺度為2.5 km）；而取初始v′, σ′場均為0。地形廓線取圖1a。以下各幅圖的橫坐標(biāo)和縱坐標(biāo)與圖1相同,不再贅述。

圖2給出了流函數(shù)演變（單位：小時(shí)，初始為0小時(shí)，下同），隨著積分時(shí)間增加，流函數(shù)場由初始時(shí)刻在垂直方向有兩個(gè)環(huán)流中心簡諧波的分布形式在一小時(shí)內(nèi)就演變?yōu)樵谠诖怪狈较蛑挥幸粋€(gè)環(huán)流中心的內(nèi)波波包形式，其中心出現(xiàn)在中上層，在水平方向該波包中心（波包振幅最大處，下同）出現(xiàn)在地形躍變處的前后，在離該處較遠(yuǎn)的地方，波動已十分微弱，此時(shí)波包中波動垂直環(huán)流圈的水平尺度約3—4 km，較初始擾動流函數(shù)環(huán)流圈的水平尺度要長。以后該波包分為兩支，向左、右兩方向傳播；在地形躍變處的左方因水較深，其環(huán)流中心也有下沉。第5小時(shí)在地形躍變處的波動已近乎消失，向左、右傳播的波包已離開地形躍變處；且因前者所處水深要較后者深，故其波包傳播速度也較后者慢。以后則兩支波包的結(jié)構(gòu)和傳播均無大的變化，只是強(qiáng)度減弱，波包中心傳播至邊界附近，而在地形躍變處波動已基本消失。

3.2 地形II的結(jié)果

圖2 無背景流時(shí)流函數(shù)演變圖（地形I）（流函數(shù)單位：m2s-1）

圖3 無背景流時(shí)流函數(shù)演變圖（地形II）（流函數(shù)單位：m2s-1）

我們還采用地形II，計(jì)算了此時(shí)內(nèi)波波包的演變。計(jì)算中除地形取地形II外，其它試驗(yàn)條件均不變。圖3給出了流函數(shù)演變，擾動在第1小時(shí)，與上類似，在地形附近產(chǎn)生了流函數(shù)垂向?yàn)閱稳Νh(huán)流的內(nèi)波波包。該環(huán)流中心也在中上層。到第2小時(shí)，就產(chǎn)生了向左、右傳播的兩支波包，波包中波動流函數(shù)垂直環(huán)流圈的水平尺度約3 km，隨著時(shí)間的增加，兩支波包繼續(xù)向左右兩方向傳播。因地形兩側(cè)水深相同，故這兩支波包傳播的速度相同，而在山脈地形處則波動消失。

為了說明以上數(shù)值試驗(yàn)中波動的性質(zhì)，計(jì)算了不同時(shí)次散度和渦度的比值。結(jié)果表明，對以上兩種地形的情況，在波動區(qū)，流場散度和渦度的比值（散渦比）均遠(yuǎn)大于1；這說明該內(nèi)波的性質(zhì)是非地轉(zhuǎn)的，屬重力慣性波。

4 有背景流時(shí)數(shù)值試驗(yàn)結(jié)果

4.1 試驗(yàn)1

在該試驗(yàn)中我們考慮了背景流的影響，取背景流僅隨z變化，背景流流動方向?yàn)閤方向（見圖4的右方），其垂直于地形方向（y方向，且地形在該方向不變），即取uˉ=uˉ(z), vˉ=0。為方便，在500 m以下格點(diǎn)取背景流為0，在500 m處取0.75 m/s，在該處以上格點(diǎn)則取其為1.5 m/s，地形取廓線I的分布形式，所取其它條件則與3.1小節(jié)中的相同，此時(shí)有S2=0，而M2=-f?uˉ/?z。圖4給出了流函數(shù)的演變（初始時(shí)流函數(shù)參見圖2中的0時(shí)次）。

由圖可見，此時(shí)與無背景流時(shí)地形I的結(jié)果（參見3.1小節(jié)）差別較大。在第1小時(shí)，內(nèi)波波包主要集中在背景流上游的地形躍變處附近，也即該處的左側(cè)附近；其在垂直方向也只有一個(gè)環(huán)流中心，但水平尺度減小。第2小時(shí)，該波包疑受臺階地形的反射，逆背景流向左傳播；而在地形躍變處右側(cè)附近的臺地上方，出現(xiàn)了一個(gè)流函數(shù)的正環(huán)流圈。以后隨著積分時(shí)間增加，該正環(huán)流圈在繼續(xù)增強(qiáng)，到第5、6小時(shí)，在該正環(huán)流圈的背景流下游方向，出現(xiàn)了多個(gè)正環(huán)流圈中心。第6小時(shí)，該正環(huán)流圈中心即內(nèi)波波包中心強(qiáng)度達(dá)最大；以后其強(qiáng)度略有減弱，但范圍繼續(xù)向右側(cè)擴(kuò)大?？偟目磥碓谂R近地形躍變處的右方，有強(qiáng)的

正環(huán)流中心存在，且其中心在水深約400 m處，并為該地形躍變所俘獲。

圖4 有背景流時(shí)流函數(shù)演變圖（地形I，試驗(yàn)1）（流函數(shù)單位：m2s-1）

圖5 有背景流時(shí)流函數(shù)演變圖（地形II，試驗(yàn)2）（流函數(shù)單位：m2s-1)

4.2 試驗(yàn)2

該試驗(yàn)采用地形II，其為大洋中常見的海底地形。圖5給出了數(shù)值計(jì)算的流函數(shù)演變（初始時(shí)分布參見圖3中的0時(shí)），除地形外各參數(shù)和背景流的取法均同4.1節(jié)，即試驗(yàn)1。由圖可見，在1小時(shí)即有了明顯變化。在臨近該地形的背景流下游處出現(xiàn)了很強(qiáng)的流函數(shù)正環(huán)流圈，在垂直方向也由初始的2圈環(huán)流變?yōu)?圈環(huán)流，強(qiáng)環(huán)流中心處于中層；在臨近該地形的背景流上游處，出現(xiàn)了較強(qiáng)的垂直方向?yàn)?圈的流函數(shù)負(fù)環(huán)流圈，其中心也在中層；而這些環(huán)流圈的水平尺度則變短。隨著積分時(shí)間增加，該地形附近的環(huán)流圈進(jìn)一步增強(qiáng)，即內(nèi)波波包增強(qiáng)；到第3小時(shí)，在相鄰該地形的下游處出現(xiàn)了很強(qiáng)的正環(huán)流圈，到第4小時(shí)該正環(huán)流圈分裂成2個(gè)環(huán)流中心，并順背景流向下游傳播；在臨近該地形的背景流上游處，則有多個(gè)負(fù)中心環(huán)流圈，并疑受該地形的反射，其逆背景流傳播。隨著積分時(shí)間的進(jìn)一步延長，在臨近該地形的環(huán)流中心繼續(xù)向背景流的上、下游傳播；到第10小時(shí)在該地形的下坡又出現(xiàn)了一個(gè)正環(huán)流圈，并順背景流傳播。而在該地形的上游，仍然是多個(gè)負(fù)中心的環(huán)流圈。

4.3 試驗(yàn)3

為了探討背景流大小的影響，仍取地形II，但減小了背景流強(qiáng)度。在此取背景流在山脈地形以上為uˉ=1 m/s,在500 m以下（最高地形以下）為0 m/s，500 m處（地形山尖）取0.5 m/s；該地形的背景流強(qiáng)度僅是4.2節(jié)試驗(yàn)2中背景流的2/3。除背景流外，其它參數(shù)和初始擾動均取得與4.2節(jié)試驗(yàn)2中相同。

由圖6所示，當(dāng)背景流減小時(shí)，在積分的第1小時(shí)，在臨近山脈附近產(chǎn)生呈垂直分布的單圈環(huán)流，即內(nèi)波波包，環(huán)流中心主要在500 m左右。隨著積分時(shí)間增加，環(huán)流中心上移，并在海底山脈上坡處存在負(fù)環(huán)流圈，而在下坡處則存在正環(huán)流圈，環(huán)流中心強(qiáng)度即內(nèi)波波包中心也隨積分時(shí)間的增加而增強(qiáng)?？偟恼f來，至積分6小時(shí)，波包的演變與4.2節(jié)試驗(yàn)2中大致相似，只不過在此因背景流較弱，故波包中心的強(qiáng)度較弱，且逆背景流傳播波包的強(qiáng)度也較弱。

在積分到第22小時(shí)，在離該山脈較遠(yuǎn)的下游區(qū)域，出現(xiàn)了多個(gè)正環(huán)流圈中心。到32小時(shí)在臨近山脈的左、右側(cè)均為較強(qiáng)的負(fù)環(huán)流圈，臨近山脈右側(cè)負(fù)環(huán)流圈的下游，則有強(qiáng)正環(huán)流圈，即波包中心；而在右側(cè)負(fù)環(huán)流圈的上游則幾乎無波動。在32小時(shí)后，擾動的分布大體仍保持這種態(tài)勢，不過在山脈的下游，僅出現(xiàn)了單個(gè)的強(qiáng)正環(huán)流圈，即內(nèi)波波包中心，而此時(shí)正、負(fù)環(huán)流圈的分布出現(xiàn)了明顯的不對稱性。

4.4 試驗(yàn)4

在以上試驗(yàn)中，背景流均取為定常；而在某些海域海水流動常以潮流為主，其呈周期性的變化。為了探討潮流對海洋內(nèi)波的作用，設(shè)定背景流隨時(shí)間作周期變化,取類似正規(guī)全日潮變化的時(shí)變背景流。在本小節(jié)試驗(yàn)中，除時(shí)變背景流取周期為24小時(shí)的簡諧波變化外，其它環(huán)境條件和初始擾動均與模擬試驗(yàn)2相同（地形也取地形II），這里取M2=-f?uˉ(z,t)/?z，S2=0時(shí)變背景流在12小時(shí)前是正向（向x正方向）流動的，在第6小時(shí)達(dá)最強(qiáng)；12小時(shí)后，背景流反向，呈反向（向x反方向）流動，到18小時(shí)流速達(dá)最強(qiáng)。背景流的振幅分布則取4.1小節(jié)試驗(yàn)1中的情況（定常背景流）。

圖7給出了該時(shí)變背景流下數(shù)值計(jì)算的流函數(shù)演變（初始情況參見圖2的0小時(shí)）。在第6小時(shí)背景流呈最大正（右）向流，在地形躍變處存在垂直單圈正環(huán)流，即海洋內(nèi)波波包中心，該正環(huán)流中心高度約900 m。以后隨著積分時(shí)間增加，正向背景流開始減弱，但該內(nèi)波波包中心還在繼續(xù)增強(qiáng)，到第10小時(shí)該中心在地形躍變處達(dá)到最強(qiáng)。到12小時(shí)，正向背景流已減弱為0，該中心也減弱，因受地形反射，反向傳播（左移）。繼續(xù)對時(shí)間積分，則背景流改變方向?yàn)榉聪蛄?，該中心繼續(xù)反向傳播（左移），到16小時(shí)該反向流進(jìn)一步增強(qiáng)，該中心也繼續(xù)反向傳播（左移）并增強(qiáng)，并在該中心正環(huán)流圈右側(cè)出現(xiàn)了較強(qiáng)的負(fù)環(huán)流圈，負(fù)環(huán)流圈的強(qiáng)度與水平尺度均要較正環(huán)流圈小，也即出現(xiàn)了不對稱性。至19小時(shí)，反向背景流又轉(zhuǎn)為正向背景流，而此時(shí)該正、負(fù)環(huán)流圈已移至邊界。第12小時(shí)及以后，在臺階地形躍變處，則基本沒有波動。

圖6 有背景流時(shí)流函數(shù)演變圖（地形II，試驗(yàn)3）（流函數(shù)單位：m2s-1）

圖7 時(shí)變背景流下流函數(shù)演變圖（地形I，試驗(yàn)4）（流函數(shù)單位：m2s-1）

5 討論

5.1 內(nèi)波的垂直運(yùn)動

由w′=?Ψ/?x知，明顯的垂直運(yùn)動發(fā)生在流函數(shù)正、負(fù)環(huán)流圈之間，如右側(cè)環(huán)流圈為正，左側(cè)為負(fù)，則兩者之間有上升運(yùn)動，否則反之，有下沉運(yùn)動。由以上各試驗(yàn)結(jié)果看，在1小時(shí)以后，最強(qiáng)的垂直運(yùn)動均發(fā)生在海洋的中層，在海面和海底附近垂直運(yùn)動均很小，而水平運(yùn)動在此處均較強(qiáng)（由u′=-?Ψ/?z可知）。這種運(yùn)動的分布體現(xiàn)了海洋內(nèi)波的特點(diǎn)。

5.2 運(yùn)動的適應(yīng)和演變

從以上各試驗(yàn)中可見，因模式的初始場由人為給定，而非實(shí)際觀測，故具有較強(qiáng)的不協(xié)調(diào)性。這樣在各試驗(yàn)的開始，必然要發(fā)生各物理量間的調(diào)整，以便達(dá)到協(xié)調(diào)。該過程可稱為適應(yīng)過程，其時(shí)間尺度很短，在1小時(shí)之內(nèi)。各試驗(yàn)共同的特點(diǎn)是，由初始時(shí)垂直方向的2個(gè)環(huán)流圈調(diào)整為1個(gè)環(huán)流圈。調(diào)整完畢后，垂直運(yùn)動和水平運(yùn)動具有5.1小節(jié)中的海洋內(nèi)波特點(diǎn)。此后運(yùn)動進(jìn)入演變過程，這時(shí)因背景流與地形的不同，各試驗(yàn)的演變情況各有不同，但共同的特點(diǎn)是演變過程時(shí)間尺度較長，要遠(yuǎn)大于適應(yīng)過程的1小時(shí)。隨著時(shí)間的推移，若運(yùn)動是穩(wěn)定的且又有耗散的話，則運(yùn)動的能量將被耗散掉，運(yùn)動將趨于靜止。因海洋內(nèi)波是海洋中尺度現(xiàn)象，其性質(zhì)與準(zhǔn)地轉(zhuǎn)的海洋Rossby波明顯不同，其是重力波或重力慣性波，是非地轉(zhuǎn)的，故該演變過程屬海洋重力慣性波的演變，而非準(zhǔn)地轉(zhuǎn)的Rossby波的演變。葉篤正,巢紀(jì)平曾討論過在大氣中關(guān)于中尺度運(yùn)動的適應(yīng)、演變等問題[13]，海洋與大氣同屬地球流體，本文的試驗(yàn)也驗(yàn)證了他們的結(jié)論。

5.3 單純地形的影響

本文采用了兩種地形，即海底臺階地形（地形I）與海底山脈地形（地形II）；前者與巴士海峽中的情形類似。從以上試驗(yàn)結(jié)果看，地形觸發(fā)內(nèi)波是顯而易見的。當(dāng)不考慮背景流時(shí)，對于以上兩種地形，均在適應(yīng)過程結(jié)束后，在地形（指海底臺階地形躍變處或海底山脈處）附近激發(fā)出海洋內(nèi)波，以后內(nèi)波波包分別從此處向左、右方向傳播，在經(jīng)過一段時(shí)間后，地形處則無波動；而內(nèi)波波包的傳播速度則與海水深度有關(guān)，深度大處傳播速度慢。

5.4 背景流與地形的共同影響

在背景流與地形共同作用下，內(nèi)波波包的演變與單純地形的結(jié)果有明顯不同。此時(shí)背景流不但影響波包的強(qiáng)度及其傳播速度，而且影響環(huán)流結(jié)構(gòu)。

在海底臺階地形與定常背景流共同作用下，波包發(fā)展的強(qiáng)中心位于該地形的躍階處，其不隨積分時(shí)間增加而向下游傳播；該波包可認(rèn)為被該地形躍階所俘獲，且在該處有強(qiáng)烈的上升運(yùn)動，有時(shí)在緊鄰該上升運(yùn)動的上游還有強(qiáng)烈的下沉運(yùn)動（如4.1小節(jié)試驗(yàn)1第6小時(shí)）；這樣的強(qiáng)內(nèi)波波動對潛艇在該處附近的活動影響較大；而在該俘獲波包的下游，則存在一系列的波列。

在山脈地形與定常背景流共同作用下，在演變過程開始（第1小時(shí)），在緊鄰海底山脈下坡處，亦有較強(qiáng)的內(nèi)波波包中心存在，在該處有強(qiáng)下沉運(yùn)動，數(shù)小時(shí)后則該波包中心向下游傳播，而在該地形附近則波動減弱，這點(diǎn)與臺階地形不同?？磥韮?nèi)波波包中心是否被地形俘獲取決于地形陡峭度，山脈地形的陡峭度不如臺階地形，故其不為俘獲。

在有定常背景流時(shí)，在地形的左側(cè)（上游）與右側(cè)（下游），內(nèi)波波包有不同的結(jié)構(gòu)?？偟恼f來，下游內(nèi)波較強(qiáng)，且以正環(huán)流圈占優(yōu)，而上游則反之；這與無背景流的情況有明顯不同（參見圖2、3），定常背景流的影響就體現(xiàn)于此。

當(dāng)背景流隨時(shí)間發(fā)生變化時(shí)，當(dāng)背景流向著臺階地形流動（正向流動）時(shí)，在該地形躍階處易產(chǎn)生較強(qiáng)的內(nèi)波波包中心；當(dāng)背景流改變方向，逆著臺階地形流動（反向流動）時(shí)，在該地形躍變處則無明顯強(qiáng)波動產(chǎn)生。在該臺階地形躍階處，背景流強(qiáng)度的變化與該處內(nèi)波強(qiáng)度的變化不是同步的，內(nèi)波強(qiáng)度的變化比背景流強(qiáng)度的變化要滯后。

5.5 非線性的作用

當(dāng)積分時(shí)間較長時(shí)，非線性的影響將逐步明顯，表現(xiàn)為內(nèi)波波包正、負(fù)環(huán)流圈出現(xiàn)非對稱性，正環(huán)流圈強(qiáng)度要大于負(fù)環(huán)流圈強(qiáng)度；在流體內(nèi)部與該正、負(fù)環(huán)流圈相鄰的交界處，有強(qiáng)上升或下沉運(yùn)動（具體是上升還是下沉要看這兩個(gè)環(huán)流圈的配置），這表明該處有強(qiáng)內(nèi)波波包中心，從試驗(yàn)結(jié)果看，強(qiáng)上升要出現(xiàn)得更多些；這種帶有強(qiáng)上升（下沉）的內(nèi)波波包對潛艇活動有明顯影響。以上情況的發(fā)生是由于非線性追趕效應(yīng)所造成的，因該效應(yīng)使最初的簡諧波波包演變成橢圓余弦波波包（此時(shí)出現(xiàn)明顯的正、負(fù)環(huán)流圈的不對稱性），最終使該內(nèi)波波包演化為孤立內(nèi)波，此時(shí)負(fù)環(huán)流圈將消失，僅剩下正環(huán)流圈。

6 結(jié)語

本文采用Boussinesq近似，研制了一個(gè)準(zhǔn)三維、非線性、非靜力的斜壓數(shù)值模式，利用該模式做了背景流與海底地形對水平尺度較短的海洋內(nèi)波影響的數(shù)值試驗(yàn)。主要結(jié)論有：

（1）在本文數(shù)值試驗(yàn)中，在初始1小時(shí)內(nèi)，有適應(yīng)過程，并最終形成垂直方向?yàn)閱稳Νh(huán)流的海洋內(nèi)波波包，以后進(jìn)入海洋內(nèi)波的演變過程，并在非線性作用下形成橢圓余弦波波包，并向孤立內(nèi)波演變。

（2）當(dāng)無背景流時(shí)，地形會激發(fā)出兩支傳播方向相反的海洋內(nèi)波波包，波包傳播速度與水深有關(guān)，水深大處波包傳播較慢，反之亦然。

（3）當(dāng)存在背景流時(shí)，在其與地形共同作用下，在地形處會激發(fā)出較強(qiáng)的內(nèi)波波包，若地形足夠陡峭，則內(nèi)波波包在該陡峭處會被截陷，造成該處強(qiáng)烈的垂直運(yùn)動。背景流的強(qiáng)度和方向都影響到內(nèi)波的生成和演變。

（4）本文計(jì)算了試驗(yàn)中內(nèi)波的散度渦度比，驗(yàn)證了該內(nèi)波的性質(zhì)屬重力慣性波。

本文給出了以上數(shù)值試驗(yàn)得到的一些結(jié)論，并發(fā)現(xiàn)該模式能較好地模擬背景流和地形對內(nèi)波的影響。因模式范圍所限，本文主要討論了背景流與海底地形對水平尺度較短的海洋內(nèi)波的影響。因模擬是在準(zhǔn)三維模式下進(jìn)行，故模擬的結(jié)果還不完善，對海洋內(nèi)波的討論亦不夠深入，在這方面今后還須做更多的工作。

[1]方欣華,杜濤.海洋內(nèi)波基礎(chǔ)和中國海內(nèi)波[M].青島:中國海洋大學(xué)出版社,2005,1-78.

[2]Thorpe.S A.On wave interaction in a stratified fluid[J].J.Fluid Mech.,1966,24:737-751.

[3]Thorpe.S A.Experiments on instability and turbulence in a stable stratified shear flow[J].J.Fluid Mech,1973,61:731-751.

[4]Thorpe.S A.The excitation dissipation and interaction of internal waves in the deep ocean[J].J Geophys Res.,1975,80(3):328-338.

[5]趙俊生,耿世江,孫洪亮.黃海內(nèi)波場的觀測和研究[R].中國海洋學(xué)文集,北京,海洋出版社,1992,78-99.

[6]陳震遠(yuǎn),許榮中,陳震武.孤立內(nèi)波遇海脊障礙物的能量消散和速度分布[J].西南交通大學(xué)學(xué)報(bào),2006,41(4):433-437.

[7]Xu,Z.H.,Yin,B.S.,and Hou,Y.J.Highly nonlinear internal solitary waves over continental shelf of northwestern South China Sea[J].Chin.J.Oceanol.Limonol,2010,28(5):1049-1054.

[8]Liu A K,Chang Y S,Hsu M K,et al.Evolution of nonlinear internal waves in the east and south China Seas[J].J.Geophys Res.,1998,103(C4):1613-1624.

[9]HOLLOWAY P E,MERRIFIELD M A.Internal tide generation by seamounts,ridges and island[J].Journal of Geophysical Research,1999,104:25937-25951

[10]CUMMINS P F.Stratified flow over topography:time-dependent comparisons between model solutions and observation[J]. Dynam.Atmos.Oceans,2000,33:43-47

[11]CAI Shu-qun,LONG Xiao-min and HUANG Qi-zhou.Preliminary numerical study on the generation condition of induced internal solitons in the northern South China Sea[J].Acta Oceanologica Sinica,2003,25(4):119-124.

[12]GRIMSHAW R.Internal solitary waves in a variable medium[J]. GAMM-Mitteilungen,2007,30(1):96-109.

[13]葉篤正,巢紀(jì)平.論大氣運(yùn)動的多時(shí)態(tài)特征—適應(yīng)、發(fā)展和準(zhǔn)定常演變[J].大氣科學(xué),1998,22(4):385-398

Abstract：Base on Boussinesq approximation,a quasi-three-dimensional,nonlinear,non-static baroclinic numerical model was developed and used to study the effects of background currents and topography on the high-frequency oceanic internal wave.The main conclusions include:the adaptation process lasts for an hour and then the single-circle internal wave packet occurs in the vertical for each numerical experiment;.during the evolution process of ocean internal wave,the internal wave was steepened into oval cosine waveform via the nonlinear effects. Without the constant background currents,topography effects could trigger two wave packages propagating in the opposite direction and wave packet propagate slower in deep water than that in shallow water.With the constant background currents,strong internal wave packages are generated under the influence of topography and background current.When the terrain is steep enough,the internal wave packets are trapped in this region,resulting in strong vertical motions.When background current changes with time,it has different effects on internal wave packet in the step-like topography.

Key words：ocean internal wave;topography;background current

Preliminary study of background current and topography offects on ocean internal wave

ZHANG Xiang1,DENG Bing2,ZHANG Ming3
(1.the Hydrometeorological Center of Navy of P.L.A.Beijing 100161China；2.the Institute of Applied Meteorology in Beijing,Beijing 100029 China；3.Laboratory of Atmospheric Circulation and Short-range Climate Forecast,Meteorological College,P.L.A.University of Science and Technology,Nanjing 211101China)

P731

A

1003-0239（2012）03-0026-09

2011-05-02

國家自然科學(xué)基金資助（90411006）

張翔(1963-)，男，高級工程師，主要從事物理海洋應(yīng)用研究。E-mail:zhxiang0325@163.com