朱聿超， 趙小龍， 趙　瑋

(中國海洋大學(xué)物理海洋實驗室， 山東 青島 266100)

?

南海簡正模態(tài)型季節(jié)內(nèi)振蕩的數(shù)值模擬?

朱聿超， 趙小龍， 趙瑋??

(中國海洋大學(xué)物理海洋實驗室， 山東 青島 266100)

摘要：本文基于大洋環(huán)流數(shù)值模式HYCOM對南海季節(jié)內(nèi)振蕩進行模擬，并利用CEOF方法對季節(jié)內(nèi)振蕩的物理特性進行了分析。結(jié)果表明，季節(jié)內(nèi)振蕩起源于南海東邊界并以5.87cm·s-1向西傳播，在南海海盆內(nèi)存在3個完整波型，波長為360km，周期為71d。季節(jié)內(nèi)振蕩引起的流速振幅在表層最大，隨著深度的增加迅速減小。Rossby波理論的分析結(jié)果表明,該振蕩是南海海盆內(nèi)的簡正模態(tài)，屬于Rossby波的第一斜壓模態(tài)。

關(guān)鍵詞：南海； 第一斜壓Rossby波； 復(fù)經(jīng)驗正交函數(shù)

ZHU Yu-Chao, ZHAO Xiao-Long, ZHAO Wei. The intraseasonal basin normal mode in South China Sea envisioned from HYCOM[J]. Periodical of Ocean University of China, 2016, 46(2): 1-5.

南海是西北太平洋最大的邊緣海，平均水深大于1 000m，最大水深超過5000m，主要通過呂宋海峽、臺灣海峽、民都洛海峽等與周圍海區(qū)進行水體交換。其中呂宋海峽水深超過2000m，是連接南海與太平洋的唯一深水通道。

南海地處東亞季風(fēng)區(qū)，季風(fēng)是海洋上層環(huán)流的主要驅(qū)動力。利用早期水文觀測、海表面高度記錄，Wyrtki[1]首先給出南海表層環(huán)流結(jié)構(gòu)：在冬季，南海海盆呈現(xiàn)氣旋型環(huán)流；在夏季，南海南部出現(xiàn)反氣旋型環(huán)流。對于深層，Tian等[2]在呂宋海峽的觀測表明，南海1000m以下的垂向混合率高于太平洋一側(cè)2個數(shù)量級。南海深層的強混合導(dǎo)致太平洋水密度更大，由此產(chǎn)生的水平壓強差驅(qū)動了從太平洋到南海的深層溢流。太平洋深層水通過巴士海峽進入南海后，在深層形成水平氣旋型的環(huán)流[3-4]。在南海內(nèi)區(qū)強混合的作用下，深層水上升并在垂直方向上形成南海深層徑向翻轉(zhuǎn)流。由于缺乏直接觀測，溫鹽診斷與數(shù)值模擬是研究南海深層環(huán)流的主要手段[4-8]。

季節(jié)內(nèi)振蕩對流場的變異有非常重要的影響：WOCE計劃對馬斯克林海盆進行了2年的海流觀測，發(fā)現(xiàn)流場存在顯著的季節(jié)內(nèi)振蕩，周期為60d，流速振幅達到10cm·s-1。振蕩信號以7.0cm·s-1的相速度向西傳播，符合正壓Rossby波的特征[9]；北大西洋深層西邊界流的觀測結(jié)果表明流場存在以40d為周期的季節(jié)內(nèi)振蕩，地形Rossby波是造成振蕩的原因[10]。由于觀測資料的匱乏，對于南海季節(jié)內(nèi)振蕩的研究主要集中在中上層。Zhuang等[11]結(jié)合衛(wèi)星高度計與OFES模式資料，對南海海表面高度的季節(jié)內(nèi)振蕩進行了研究：在南海北部陸架區(qū)及越南東側(cè)海域分別存2個強季節(jié)內(nèi)振蕩區(qū)域，振蕩周期約為10周，西傳Rossby波是產(chǎn)生季節(jié)內(nèi)振蕩的原因之一。Zheng等[12]利用2007年與2009年夏季航次獲得的CTD資料，發(fā)現(xiàn)南海內(nèi)部存在3個沿緯向的完整波型。Zhang等[13]基于1年的潛標(biāo)觀測數(shù)據(jù)，發(fā)現(xiàn)在呂宋海峽中層存在周期為10～30d和60d的季節(jié)內(nèi)振蕩，并且與上層的季節(jié)內(nèi)振蕩呈現(xiàn)負相關(guān)關(guān)系。

地形與緯度、層結(jié)、幾何形狀共同決定了封閉海盆內(nèi)的簡正模態(tài)[14]，而南海海盆內(nèi)的簡正模態(tài)具有季節(jié)內(nèi)振蕩的特征[12]。本文基于HYCOM(Hybrid Coordinate Ocean Model)模式并結(jié)合現(xiàn)有的觀測結(jié)果，對南海簡正模態(tài)引起的季節(jié)內(nèi)振蕩進行研究。

1模式設(shè)置及分析方法

本文對季節(jié)內(nèi)振蕩進行模擬時采用了與Zhao等[5]相同的模式設(shè)置。模式計算區(qū)域包括整個南海及部分西北太平洋區(qū)域(4°N～25°N，105°E～125°E，見圖1)。模式水平分辨率為1/12(°)，垂向分為32層，采用2000m處參考壓強，目標(biāo)位勢位密依次為28.10，28.90，29.70，30.50， 30.95，31.50，32.05，32.60，33.15，33.70，34.25，34.75，35.15，35.50，35.80，36.04，36.20，36.34，36.46，36.56，36.64，36.70，36.74，36.78，36.82，36.84，36.86，36.88，36.92，36.96，37.01和37.06kg·m-3。模式采用GDEM3溫鹽數(shù)據(jù)作為初始場，鑒于Tian等[2]在呂宋海峽的觀測，本文將南海1000m以下垂向混合率設(shè)為10-3m2·s-1。外強迫不是影響海盆內(nèi)簡正模態(tài)的因素，因此忽略海表面強迫，東邊界采用溫鹽松弛邊界條件，其它邊界均為閉邊界。模式在5年之后達到穩(wěn)定狀態(tài)，因此本文采用5～10a的模式結(jié)果進行分析。

通過經(jīng)驗正交函數(shù)(EOF)方法可將隨時間變化的變量場分解為不隨時間變化的空間函數(shù)部分與只依賴時間變化的時間函數(shù)部分。傳統(tǒng)EOF分解獲得的是空間駐波結(jié)構(gòu)，不能刻畫西傳Rossby波。因此，本文采用復(fù)EOF(CEOF)方法對模式流場進行分析。CEOF方法是把實數(shù)隨機向量，利用Hilbert變換后作為虛部構(gòu)成復(fù)值隨機向量。對復(fù)值隨機向量進行EOF分解，得到的結(jié)果不僅反映出狀態(tài)變率，還會反映出在相同時間上狀態(tài)的空間變化。具體原理可參見Barnett[15]的工作。

圖1　模式計算區(qū)域地形

2模式結(jié)果

根據(jù)Barnett[15]的經(jīng)驗，應(yīng)用CEOF方法分析波動時，需要先對原始變量場進行適當(dāng)?shù)臑V波。為了選取合適的帶通頻段，作者先對衛(wèi)星高度計獲得的海表面高度異常(Sea level anomaly, SLA)資料進行分析。文中所用的SLA數(shù)據(jù)來于法國AVISO(Archiving, Validation and Interpretation of Satellite Oceanographic Data)項目,該數(shù)據(jù)融合了TOPEX/Poseidon,ERS-1/2等多種衛(wèi)星測高數(shù)據(jù)，空間分辨率為25km，時間分辨率為7d。本文利用1992—2013年的SLA數(shù)據(jù)，分別計算振蕩周期在 8～120d(見圖2(a))，40～80d(見圖2(b))，14～40d(由于采樣間隔為7d，可分辨最小周期為14d，(見圖2(c))之間的能量比例。

(選擇的頻段為2π/120 ～2π/80d(a)、2π/80 ～2π/40d(b)and 2π/40 ～2π/14d(c)。The selected frequency bands are 2π/120 ～2π/80 d (a), 2π/80 ～2π/40 d (b) and 2π/40 ～2π/14 d (c).)

圖2不同頻段內(nèi)譜能量百分比

Fig.2The spectrum percentage between a selected frequency band

從圖2可以看出80～120d周期振蕩在南海北部的能量比例可以達到30%，在其余海域的能量比例在20%以下；40 ～80d周期振蕩的能量比例最大，在中央深水區(qū)的能量比例可以達到70%；14 ～40d周期振蕩在南海北部陸架區(qū)、泰國灣內(nèi)的能量比例可以達到60%以上，但在南海內(nèi)區(qū)的能量比例很低。雖然SLA反應(yīng)的季節(jié)內(nèi)振蕩可能與大氣強迫有密切的關(guān)系，但大氣強迫如果與簡正模態(tài)具有相同的頻率，會通過共振機制將此頻率的信號增強[9]。因此先假設(shè)40～80d是本文所關(guān)心的振蕩周期，把帶通濾波后流速數(shù)據(jù)進行CEOF分解。選取14°N斷面的結(jié)果進行CEOF分解，展示其垂向結(jié)構(gòu)與傳播特性。CEOF分解后可獲得第一模態(tài)的時間相位函數(shù)與時間振幅函數(shù)(見圖3(a))、空間相位函數(shù)(見圖4(a))與空間振幅函數(shù)(見圖4(b))。第一模態(tài)方差貢獻率為41.6%。

((a)藍色實線為時間振幅函數(shù)，紅色虛線為時間相位函數(shù)；(b)為時間函數(shù)的功率譜分析，其中紅色虛線為95%置信度的顯著性檢驗。(a) The solid blue line denotes the temporal amplitude value and doted red line denotes the temporal phase value between ±π in the left panel. (b) Represents the power spectral analysis for temporal phase function, and doted red line denotes the 95% significance test.)

圖314°N斷面徑向速度的第一模態(tài)時間振幅函數(shù)與時間相位函數(shù)(a)與時間相位函數(shù)的功率譜密度(b)

Fig.3Temporal amplitude and temporal phase for first empirical mode of the meridional

velocities in Lat 14°N (a) and power spectral analysis for the temporal phase function (b)

通過時間振幅函數(shù)與時間相位函數(shù)可得到空間結(jié)構(gòu)隨時間的變化，并可以通過功率譜分析找到振蕩的主要周期。本文對時間相位函數(shù)進行功率譜分析，發(fā)現(xiàn)其主要周期是71d(見圖3(b))。

圖4　14°N斷面徑向速度的第一模態(tài)空間相位

通過空間相位函數(shù)可以獲得波動傳播方向、波長等信息。從圖4(a)可清晰看出，在南海深海盆內(nèi)，波自東向西傳播，緯向上存在3個完整波形，與Zheng[12]的觀測結(jié)果相一致。同相位線的形態(tài)受地形影響較大。例如在114°E～116°E之間地形變化劇烈，1000m以下同相位線發(fā)生扭曲；在116°E～118°E之間地形平坦，同相位線呈現(xiàn)垂直狀。根據(jù)空間振幅函數(shù)可以分析流速振幅的空間分布。從圖4(b)可以看出，上層流速振幅明顯大于下層，流速振幅的最大值出現(xiàn)在南海西邊界附近，1 000m以下振幅變得很小，變化較微弱。上述結(jié)果表明，波動存在垂向結(jié)構(gòu)，與Rossby波的第一斜壓模態(tài)相似。將時間振幅函數(shù)與空間振幅函數(shù)相乘可以計算出振蕩引起的流速變化。在表層，流速振幅可達到1.5cm·s-1，在1000m以下的深層，流速振幅只有mm·s-1量級。

14°N斷面CEOF分解的結(jié)果表明，這種振蕩是一種西傳波動。為了描述波動的傳播特性，本文選取500m深度處經(jīng)向速度做CEOF分解(此處只描述其空間相位函數(shù)，見圖5)，發(fā)現(xiàn)波動起源于南海東邊界，同相位線基本與經(jīng)線平行，表明波動主要是向西傳播，只有在12°N以南以及中沙島附近的海域，同相位線的形態(tài)受地形影響較大，發(fā)生了一定的扭曲。

通過上述分析得到了振蕩的基本性質(zhì)：振蕩的生成源地是南海的東邊界，以行波的形式向西傳播；在南海海盆內(nèi)存在3個完整波型；波動的周期為71d；波動引起的流速振幅在表層最大，隨著深度的增加迅速減小。

圖5　500 m深度徑向速度第一模態(tài)空間相位函數(shù)

3討論

南?？煽醋魇欠忾]的海盆，在其內(nèi)部激發(fā)的Rossby波存在水平結(jié)構(gòu)。正如圖5所示，在南海深海盆內(nèi)存在3個完整的波形，封閉海區(qū)內(nèi)存在準(zhǔn)地轉(zhuǎn)簡正模態(tài)s，需要滿足海盆幾何尺度是Rossby波長整數(shù)倍的條件[14]。在14°N斷面，南海海盆寬度約為1080km，因此可得到緯向波長約為360km。

層結(jié)流體中激發(fā)的Rossby波會存在垂向結(jié)構(gòu)，不同斜壓模態(tài)對應(yīng)著不同的Rossby變形半徑，此處先驗證n=1時第一斜壓模的結(jié)果。無論層結(jié)的垂向結(jié)構(gòu)如何，正壓模態(tài)都可能存在，在這種情況下，變形半徑會足夠大以至于在頻散關(guān)系中，與k2相比，1/R2較小而可以忽略。緯向波長約為360km，則k約為1.75×10-5m-1(見圖5)。對于l的取值，因圖5表明波動主要呈現(xiàn)向西傳播的特性，所以本文在頻散關(guān)系中忽略了與l有關(guān)的項。β在14°N取2.22×10-11s-1，計算獲得正壓周期為57d，計算結(jié)果明顯小于71d的周期，與譜分析的結(jié)果并不一致，所以排除了正壓Rossby波的可能性。對于R1的計算，本文借助WOA 09氣候態(tài)溫鹽數(shù)據(jù)，采用WKB方法計算出南海第一斜壓Rossby變形半徑[16]。南海第一斜壓Rossby變形半徑空間變化較大，對于深海區(qū)域，變形半徑在80～120km之間[17]，如果Rossby變形半徑取100km，將這些參量帶入頻散關(guān)系后發(fā)現(xiàn)計算獲得的周期為75d，這與譜分析獲得的結(jié)果是非常接近的。所以可以確定，模式中獲得的季節(jié)內(nèi)振蕩是Rossby波的第一斜壓模態(tài)。

4結(jié)語與展望

本文利用HYCOM模式對南海環(huán)流進行模擬，模式輸出流場中出現(xiàn)以71d為周期的季節(jié)內(nèi)振蕩。利用CEOF分析方法，證明西傳第一斜壓Rossby波是振蕩產(chǎn)生的原因。波動在緯向上存在3個完整波型，水平波長為360km，對應(yīng)的西傳波速為5.87cm·s-1。通過對南海上升流與下沉流量級的分析，本文推測太平洋深層水通過呂宋海峽入侵南海后下沉，下沉流導(dǎo)致的位渦變化是波動產(chǎn)生的主要機制。以上的機制解釋只是模式結(jié)果的初步探討，在下一步工作中，需要對上升流與下沉流的空間分布進行進一步的研究。

參考文獻：

[1]Wyrtki K. Physical oceanography of the Southeast Asian waters [J]. Scripps Institution of Oceanography, 1961.

[2]Tian J, Yang Q, Zhao W. Observation of enhanced diapycnal mixing in the South China Sea [J]. J Phys Oceanogr, 2009, 39: 3191-3203.

[3]Qu T, Girton J B, Whitehead J A. Deepwater overflow through Luzon strait [J]. J Geophys Res, 2006, 111: C01002. doi: 10. 1029/2005JC003139.

[4]Wang G, Xie S P, Qu T, et al. Deep South China Sea circulation [J]. Geophysical Research Letters, 2011, 38(5). 05601

[5]Zhao W, Zhou C, Tian J, et al. Deep water circulation in the Luzon Strait [J]. J Geophys Res, 2014, 119(2): 790-804.

[6]Dongliang Y. A numerical study of the South China Sea deep circulation and its relation to the Luzon Strait transport [J]. Actaoceanol Ogica Sinica, 2002, 21(2): 187-202.

[7]Lan J, Zhang N, Wang Y. On the dynamics of the South China Sea deep circulation [J]. J Geophys Res, 2013, 118(3): 1206-1210.

[8]Shu Y, Xue H, Wang D, et al. Meridional overturning circulation in the South China Sea envisioned from the high-resolution global reanalysis data GLBa0.08 [J]. Journal of Geophysical Research: Oceans, 2014, 119(5): 3012-3028.

[9]Warren B A, Whitworth T, La Casce J H. Forced resonant undulation in the deep Mascarene Basin [J]. Deep Sea Research, 2002, 49(7): 1513-1526.

[10]Pickart R S, Watts D R. Deep western boundary current variability at Cape Hatteras [J]. Journal of Marine Research, 1990, 48(4): 765-791.

[11]Zhuang W, Xie S P, Wang D, et al. Intraseasonal variability in sea surface height over the South China Sea [J]. J Geophys Res, 2010, 115: C04010. doi: 10. 1029/2009JC005647.

[12]Zheng Q, Hu J, Zhu B, et al. Standing wave modes observed in the South China Sea deep basin [J]. J Geophys Res, 2014, 119(7): 4185-4199.

[13]Zhang Z , Zhao W, Tian J, et al. Spatial structure andtemporal variability of the zonal flow in the Luzon Strait[J]. J Geophys Res Oceans, 2015, 120, doi: 10. 1002/2014JC010308.

[14]Pedlosky J. Geophysical Fluid Dynamics [M]. New York: Springer-Verlag, 1986：144-153.

[15]Barnett T P. Interaction of the monsoon and Pacific trade wind system at interannual time scales Part I: the equatorial zone [J]. Monthly Weather Review, 1983, 111(4): 756-773.

[16]Chelton D B, Deszoeke R A, Schlax M G, et al. Geographical variability of the first baroclinic Rossby radius of deformation [J]. J Phys Oceanogr, 1998, 28(3): 433-460.

[17]Cai S, Long X, Wu R, et al. Geographical and monthly variability of the first baroclinicRossby radius of deformation in the South China Sea [J]. Journal of Marine Systems, 2008, 74(1): 711-720.

[18]Chelton D B, Schlax M G. Global observations of oceanic Rossby waves [J]. Science, 1996, 276: 234-238.

責(zé)任編輯龐旻

The Intraseasonal Basin Normal Mode in South China Sea Envisioned from HYCOM

ZHU Yu-Chao, ZHAO Xiao-Long, ZHAO Wei

(The Lab of Physical Oceanography, Ocean University of China, Qingdao 266100, China)

Abstract:A general ocean circulation model HYCOM isdescribedto usein the modeling of the intraseasonal Basin Normal Mode in the South China Sea (SCS) and a notable intraseasonal oscillation derived from the velocity field is investigated. The Complex Empirical Orthogonal Function(CEOF) analysis is used to characterize the intraseasonal oscillation. The oscillation at a period of 71 days originates along the east boundary of SCS and propagates westward witha wavelength of 360 km, about one third of the basin width. Amplitude of the oscillation reaches its maximum at upper layer and decreases with increasing depth. By utilizing the WOA climatological dataset and Rossbywaves theories, we identified that the wave is first baroclinicRossby wave.

Key words:South China Sea; first baroclinic Rossby wave; complex empirical orthogonal function

DOI：10.6441/j.cnki.hdxb.20150109

中圖法分類號：P731.21

文獻標(biāo)志碼：A

文章編號：1672-5174(2016)02-001-06

作者簡介：朱聿超 (1989-)，男，碩士生。E-mail：zhuyuchao-2008@163.com??通訊作者：E-mail：weizhao@ouc.edu.cn

收稿日期：2015-04-02；

修訂日期：2015-05-21

基金項目：? 國家重點基礎(chǔ)研究發(fā)展計劃項目(2014CB745003)；南海海洋環(huán)流形成變異機理及其氣候效應(yīng)項目(GASI-03-01-01-03)；國家自然科學(xué)基金重大研究計劃項目(91028008)；國家自然科學(xué)基金項目(41176010；41176008)；中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(XDA11010202)；三亞深?？茖W(xué)與工程研究所項目(SIDSSE-201207)資助

引用格式：朱聿超， 趙小龍， 趙瑋.南海簡正模態(tài)型季節(jié)內(nèi)振蕩的數(shù)值模擬[J].中國海洋大學(xué)學(xué)報(自然科學(xué)版)， 2016, 46(2)： 1-5.

Supported by National Basic Reasearch Porogram of China(2014CB745003)；Formation，variations and climatic impacts of the circulation in the South China Sea(GASI-03-01-01-03)；National Natural Science Foundation of China(91028008)；National Natural Science Foundation of China(41176010；41176008)；Strategic Priority Research Program of the Chinese Academy of Sciences(XDA11010202)；Sanya Institute of Deep-sea Science and Engineering(SIDSSE-201207)