曠芳芳，潘愛軍*，張俊鵬，黃獎(jiǎng)，蔡尚湛

( 1.自然資源部第三海洋研究所，福建 廈門 361005)

1 引言

由于深海潛標(biāo)能提供長時(shí)間序列的海流的直接觀測(cè)資料，對(duì)深海海流的研究具有重要作用。于克哲[1]使用磁錄海流計(jì)在千米潛標(biāo)系統(tǒng)南海試驗(yàn)期間獲取的深海海流資料，對(duì)該海域深層海流進(jìn)行了初步分析，提出了潮流和余流的運(yùn)動(dòng)規(guī)律；Zhai等[2]使用菲律賓以東的深海潛標(biāo)數(shù)據(jù)分析了近底層海流和潮流的季節(jié)變化特征；在南海，多名學(xué)者通過定點(diǎn)觀測(cè)海流資料分析了觀測(cè)海區(qū)的正壓潮和內(nèi)潮等特征[3–9]。

深層海洋與上層海洋存在著不同的動(dòng)力環(huán)境和機(jī)制，近年來深海在全球海洋氣候變化中的作用得到越來越多海洋學(xué)家的重視，深海對(duì)海洋熱量的再分配引起了全球變暖的停滯[10]，深刻影響著全球的氣候變化。菲律賓海位于副熱帶西太平洋，該海域上層有北赤道流、副熱帶逆流等多種流系交匯，還是黑潮的源地，在西太平洋海氣相互作用以及氣候變化等的研究中具有重要作用；針對(duì)該海域表層和次表層已經(jīng)有非常多的研究成果，但對(duì)深海的綜合環(huán)境信息仍知之甚少。為了解呂宋海峽東側(cè)深層海水的流動(dòng)特征及其對(duì)深海生態(tài)的影響，我們?cè)谠摵Ｓ虿挤帕?套深海潛標(biāo)進(jìn)行長期觀測(cè)，獲得了長達(dá)1年的連續(xù)觀測(cè)資料。以下將對(duì)潛標(biāo)觀測(cè)數(shù)據(jù)進(jìn)行分析，重點(diǎn)討論該海域潮流的垂向分布以及時(shí)間變化特征。

2 數(shù)據(jù)和方法

2.1 數(shù)據(jù)介紹

潛標(biāo)位于西太平洋呂宋海峽以東約500 km的菲律賓海（圖1a），位置為 19.75°N，126.75°E。全球地形數(shù)據(jù)GEBCO顯示，潛標(biāo)所在位置水深約5 300 m，周邊50 km的范圍內(nèi)水深為4 800～5 810 m（圖1b），局地地形西北?東南走向。在潛標(biāo)上層水深160 m配置了1個(gè)ADCP海流剖面儀向上觀測(cè)，層間距為4 m，有效觀測(cè)深度為60～160 m，本文中選取100 m和160 m兩個(gè)層次代表次表層；在810 m、1 550 m、2 560 m和4 040 m掛載了單點(diǎn)海流計(jì)觀測(cè)深海和近底層的海流。海流的觀測(cè)時(shí)間為2015年6月至2016年6月，時(shí)間頻率為1 h?1。另外，使用WOA2013提供的季節(jié)平均的溫鹽數(shù)據(jù)計(jì)算觀測(cè)點(diǎn)的浮力頻率（圖2），其計(jì)算公式為

圖1 潛標(biāo)位置（a）和水深分布（b）Fig.1 Location of the mooring (a) and water depth around the mooring (b)

式中，N為浮力頻率；g為重力加速度；ρ為根據(jù)水溫和鹽度數(shù)據(jù)計(jì)算的位勢(shì)密度。

可看到在300 m以淺冬季的層結(jié)相對(duì)較弱，躍層位于50～200 m的深度，300 m以深季節(jié)變化不明顯。

圖2 通過WOA2013計(jì)算出的浮力頻率剖面Fig.2 Vertical structure of buoyancy frequency calculated from data in WOA2013

2.2 潮流調(diào)和分析

使用T_tide程序[11]對(duì)逐時(shí)的觀測(cè)海流數(shù)據(jù)進(jìn)行潮流的調(diào)和分析，計(jì)算主要分潮的潮流橢圓要素，包括長軸、短軸、傾角、相位和橢圓率。其中，長軸代表分潮的最大流速；短軸代表分潮的最小流速，當(dāng)短軸為正值時(shí)表示逆時(shí)針旋轉(zhuǎn)，短軸為負(fù)值時(shí)表示順時(shí)針旋轉(zhuǎn)；傾角代表最大流速與東方向的夾角；橢圓率為短軸和長軸的比值。

2.3 動(dòng)能計(jì)算公式

潮動(dòng)能（Ketide）和總動(dòng)能（Keall）的計(jì)算公式為

式中，u和v為觀測(cè)的緯向流速和經(jīng)向流速；utide和vtide為經(jīng)T_tide程序提取的潮流流速。

2.4 潮流性質(zhì)判斷方法

通過計(jì)算潮流的形態(tài)數(shù)來判斷潮流的性質(zhì)，其計(jì)算公式為

式中，a為潮流形態(tài)數(shù)；WK1、WO1和WM2分別為K1、O1和M2分潮流的平均最大流速（即為潮流橢圓的長軸）。若a≤0.5，為正規(guī)半日潮流；若 0.5

3 結(jié)果與討論

3.1 海流高頻變化頻譜特征分析

圖3是各層次海流在高頻段（3 d以內(nèi)）的旋轉(zhuǎn)譜分析結(jié)果。如圖所示，各層次海流的高頻波動(dòng)以順時(shí)針（CW）旋轉(zhuǎn)為主，譜密度峰值分別對(duì)應(yīng)S2(12 h)、M2(12.4 h)、K1(23.9 h)、O1(25.8 h)4個(gè)分潮的潮周期，其中M2和K1分潮相對(duì)S2和O1分潮更為顯著。另外各層次海流還表現(xiàn)出31.8～35.4 h的近慣性周期，經(jīng)計(jì)算本地理論慣性周期為35.5 h，可以看到上層（100 m和160 m）的近慣性周期與理論慣性周期非常接近，而中層和深層（810～4 040 m）為32 h左右，略小于理論慣性周期。Zhai等[2]對(duì)西太平洋潛標(biāo)（8°N）的分析中發(fā)現(xiàn)該海域近底層海流的慣性振蕩周期略大于理論慣性周期，與本文的分析結(jié)果有所不同。其原因還有待進(jìn)一步研究。

圖3 各層次海流旋轉(zhuǎn)譜分析Fig.3 Rotary spectrum analysis of currents at different depth

3.2 潮流橢圓要素及潮流能量

使用全年的海流觀測(cè)資料用T_tide程序計(jì)算各層次 4 個(gè)主要分潮（M2、S2、K1和 O1）的潮流橢圓要素，計(jì)算后的結(jié)果見圖4。可見除1 550 m處的O1分潮呈逆時(shí)針旋轉(zhuǎn)外，其他層次的4個(gè)分潮都為順時(shí)針旋轉(zhuǎn)。在次表層，160 m層各分潮的潮流橢圓與100 m層相似，但流速相對(duì)較??；4個(gè)主要分潮中，M2分潮最大（100 m和160 m的最大流速分別為8.58 cm/s和7.14 cm/s）、K1分潮次之（100 m和160 m的最大流速分別為6.01 cm/s和5.51 cm/s）、S2分潮最小；各分潮的傾角為148°～174°，即各分潮的最大流速方向均為西北?東南向；橢圓率為?0.75～?0.5，表現(xiàn)出旋轉(zhuǎn)流的特征。在中層（810 m），4個(gè)主要分潮中，K1分潮最大（最大流速為 2.86 cm/s）、M2次之（1.85 cm/s）、S2分潮最小，K1和S2分潮的最大流速方向?yàn)闁|北?西南向，M2和 O1為西北?東南向。在深層（1 550～4 040 m），K1分潮的最大流速最大，M2次之，S2分潮最??；1 550 m層K1和S2分潮的最大流速方向?yàn)闁|北?西南向，M2為東?西向，O1為西北?東南向；在2 560 m和4 040 m層，K1分潮的最大流速為西北?東南向，M2為東北?西南向，O1為南?北向，S2為東-西向。

圖5為各層次的潮動(dòng)能、總動(dòng)能以及潮動(dòng)能占總動(dòng)能之比。由圖可見次表層（100～160 m）的潮動(dòng)能比中層（810 m）和深層（1 550～4 040 m）大 1～2個(gè)數(shù)量級(jí)，在810～2 560 m潮動(dòng)能隨深度減小，近底層（4 040 m層）的潮動(dòng)能大于深層（1 550 m和2 560 m）。次表層（100 m、160 m）和近底層（4 040 m）的潮動(dòng)能占總動(dòng)能的20%以上，在中層和深層（810～2 560 m）的潮動(dòng)能占總動(dòng)能的10%左右。

另外，通過計(jì)算各層次潮流的形態(tài)數(shù)，可判斷次表層（100～160 m）為不正規(guī)半日潮流，中層和深層（810～4 040 m）為不正規(guī)全日潮流。

3.3 潮流的季節(jié)變化

為了探討潮流的季節(jié)變化，在每個(gè)時(shí)間點(diǎn)使用30 d的數(shù)據(jù)（前后各15 d）用T_tide程序進(jìn)行潮流的調(diào)和分析，獲得了主要分潮M2和K1的長軸的時(shí)間序列（圖6），以及計(jì)算了潮動(dòng)能的時(shí)間序列及其小波譜（圖7）。

由圖6可見，次表層M2分潮的變化范圍較大，最大流速在100 m層可接近20 cm/s，在160 m層可接近13 cm/s，在冬季（2–4 月）和夏季（8–9 月）增強(qiáng)，冬季增強(qiáng)尤為明顯（圖6a）；次表層K1分潮在秋季（9–10月）明顯減弱，可由最大時(shí)的10 cm/s減弱到3 cm/s（圖6b）；在中層的 810 m，M2分潮均在夏季（7–8月）和冬季（2–4月）增強(qiáng)（圖6a），K1分潮也在夏季（8–9 月）和冬季（1–3 月）增強(qiáng)（圖6b）；深層（1 550～4 040 m）各分潮的變化范圍相對(duì)較小，1 550 m處M2分潮在夏季（7–8 月）和冬季（2–3月）有所增強(qiáng)，K1分潮在秋季（10–12月）有所減弱；2 560 m處K1和M2分潮的變化范圍均很??；4 040 m處K1分潮在夏季（7–8月）略有增強(qiáng)。對(duì)于M2分潮的時(shí)間變化，可能與月球繞地的橢圓形軌道引起的地月距離的變化有關(guān)；經(jīng)查詢，觀測(cè)時(shí)段（2015年6月至2016年5月）月球的近地點(diǎn)時(shí)間為 2015年的 8–10月和 2016年的 4–5月[12]，文中M2分潮動(dòng)能在夏季的增強(qiáng)可以用該時(shí)段地月距離的變短來解釋，但2016年2–3月M2分潮動(dòng)能的增強(qiáng)應(yīng)該是其他的原因?qū)е隆?/p>

圖4 各層次各分潮的潮流橢圓Fig.4 Tidal ellipse at different levels

圖7表明各層次潮動(dòng)能表現(xiàn)出較顯著的季節(jié)變化特征。從圖中可以看到，各層次潮動(dòng)能均在2–3月份增強(qiáng)，在其他時(shí)間的變化則不太一致。在次表層和中層，潮動(dòng)能在夏季（8–9月）和冬季（2–3月）顯著增強(qiáng)，可達(dá)到年均值的2倍以上（圖7a，圖7b），這主要與M2分潮和K1分潮在夏季和冬季的增強(qiáng)相對(duì)應(yīng)；在深層，潮動(dòng)能均在夏季（6–8月）和冬季（2–3月）增強(qiáng)，除此之外，在2 560 m處潮動(dòng)能在1月出現(xiàn)另一個(gè)峰值，而在4 040 m處潮動(dòng)能在11月也出現(xiàn)峰值。小波分析結(jié)果顯示，次表層的潮動(dòng)能全年表現(xiàn)出了近100 d的振蕩周期；中層的潮動(dòng)能表現(xiàn)出87 d左右的振蕩周期，在夏半年明顯；深層和近底層則為47 d和114 d左右的振蕩周期。從季節(jié)平均的浮力頻率剖面來看（圖2），在300 m以淺的上層，冬季躍層較深且較弱，其他季節(jié)躍層較淺且相對(duì)較強(qiáng)，而在深層季節(jié)變化不明顯，潮動(dòng)能的季節(jié)變化特征似乎與躍層的變化并不完全一致，其背后的機(jī)制還有待將來的進(jìn)一步分析。

圖6 各層次M2和K1分潮流長軸的時(shí)間序列Fig.6 Time series of M2 and K1 constituent major axes at different layers

圖5 各層次潮動(dòng)能及其占比Fig.5 Tidal kinetic energy and its portion to total kinetic energy at different levels

圖7 各層次潮動(dòng)能的時(shí)間序列及其小波分析Fig.7 Time series and wavelet analysis of tidal kinetic energy at different levels

另外值得注意的是，在做潮流調(diào)和分析時(shí)由于所使用的海流時(shí)間序列的長短不同，所計(jì)算的結(jié)果也有所不同。使用30 d的數(shù)據(jù)比使用全年的數(shù)據(jù)計(jì)算的分潮的潮流橢圓相對(duì)較大、潮動(dòng)能也相對(duì)較強(qiáng)。上節(jié)使用全年的數(shù)據(jù)來分析分潮及潮動(dòng)能的整體特征，而這節(jié)則重點(diǎn)關(guān)注其時(shí)間變化。

4 結(jié)論

通過分析一套深海潛標(biāo)觀測(cè)海流的高頻變化特征，本文得到以下結(jié)論：

（1）各層次海流的高頻波動(dòng)以順時(shí)針旋轉(zhuǎn)為主，次表層海流的近慣性周期與當(dāng)?shù)乩碚搼T性周期非常接近，而中層和深層（810～4 040 m）略小于當(dāng)?shù)乩碚搼T性周期。

（2）次表層（100～160 m）的潮動(dòng)能比中層（810 m）和深層（1 550～4 040 m）大 1～2個(gè)數(shù)量級(jí)，在 810～2 560 m潮動(dòng)能隨深度減小，4 040 m層的潮動(dòng)能大于1 550 m和2 560 m層。表層（100～160 m）為不正規(guī)半日潮流，中層和深層（810～4 040 m）為不正規(guī)全日潮流。

（3）各層次潮動(dòng)能均在夏季（6–9月）和冬季（2–3月）增強(qiáng)，這主要與M2分潮和K1分潮在夏季和冬季的增強(qiáng)相對(duì)應(yīng)；除此之外，在2 560 m處潮動(dòng)能在1月增強(qiáng)，而在4 040 m處潮動(dòng)能在11月也出現(xiàn)峰值。小波分析結(jié)果顯示，次表層的潮動(dòng)能全年表現(xiàn)出了近100 d的振蕩周期；中層的潮動(dòng)能表現(xiàn)出87 d左右的振蕩周期，在夏半年明顯；深層和近底層則為47 d和114 d左右。

以上研究表明，深海潮流在不同的深度層次具有不同特征，這些特征背后的動(dòng)力機(jī)制如何，還需要在將來通過更多的深海觀測(cè)以及數(shù)值模擬進(jìn)一步地研究。