丁睿彬，陳大可，姜良紅

（衛(wèi)星海洋環(huán)境動力學(xué)國家重點(diǎn)實(shí)驗(yàn)室，國家海洋局 第二海洋研究所，浙江 杭州 310012）

0 引言

起源于北赤道流，在菲律賓東岸15°N左右分叉后形成北分支的黑潮，是北太平洋的西邊界強(qiáng)流，它攜帶大量的高溫、高鹽水北上，給東亞區(qū)域的海洋環(huán)境與氣候帶來重要影響［1－2］。黑潮在流經(jīng)呂宋海峽時失去西邊界的支撐（圖1），因而發(fā)生重大的流態(tài)改變，這種改變直接關(guān)系到太平洋與南海的水交換以及黑潮下游的狀況，因此黑潮在呂宋海峽附近的變化問題一直是區(qū)域性海洋學(xué)關(guān)注和爭議的焦點(diǎn)。這一問題的主要爭議在于黑潮是否入侵南海，或者是以何種方式進(jìn)入南海的。第一種觀點(diǎn)認(rèn)為存在南?！傲魈住保春诔痹趨嗡魏{南邊進(jìn)入南海后從北邊出來，在南海東北部形成一個“流套”［3］。第二種觀點(diǎn)認(rèn)為黑潮有一個“南海分支”，在夏季表現(xiàn)為一支向西流動的海流，直接進(jìn)入南海［4］。第三種觀點(diǎn)是“流套”和“分支”觀點(diǎn)的組合，認(rèn)為黑潮在呂宋海峽中部以南進(jìn)入南海，一部分在臺灣島的西南側(cè)形成“流套”，另一部分則在臺灣島的南邊形成一個氣旋性的“分支”［5］。近年來也有學(xué)者認(rèn)為黑潮水能以中尺度渦或次中尺度渦的形式進(jìn)入南海。

黑潮主流軸（即跨越黑潮速度梯度最大處）的位置及其隨時間的變化在很大程度上能代表黑潮的時空變化。過去已有一些針對黑潮主軸的研究，如日本學(xué)者在日本島附近海域做了不少相關(guān)工作［6］，但這些早期研究主要基于有限的現(xiàn)場觀測資料，其空間覆蓋率不足以獲得黑潮流軸變化的完整圖像；后來有學(xué)者結(jié)合衛(wèi)星高度計(jì)和浮標(biāo)數(shù)據(jù)，較為準(zhǔn)確地確定了局地的黑潮的主軸［7］；最近 LIU et al［8］用衛(wèi)星高度計(jì)資料進(jìn)一步研究了東海黑潮的路徑變化。然而，對于黑潮流軸在呂宋海峽附近的大幅度變化還缺乏系統(tǒng)研究。

本文以（1／12）°分辨率的區(qū)域海洋環(huán)流模式（ROMS）在北太平洋試驗(yàn)9a（2000—2008年）的輸出結(jié)果為基本數(shù)據(jù)，利用自行編制的基于速度場的流軸自動檢測方法，分析了黑潮流軸在呂宋海峽附近的位置及其在年際、季節(jié)和季節(jié)內(nèi)時間尺度上的變化特征，并對引起這些變化的物理過程和機(jī)制進(jìn)行了初步的探討。

圖1 研究區(qū)域地形圖以及基于觀測表層流場（AVISO，黑線）和ROMS模式10m（藍(lán)線）、200m（紅線）流場的黑潮流軸平均位置Fig.1 Map of research area and mean Kuroshio axes based on surface velocities from observation（AVISO，black）and 10m （blue）and 200m （red）velocities from ROMS model

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

1.1 數(shù)據(jù)

本文所用的主要數(shù)據(jù)來自ROMS模式對1個北太平洋區(qū)域（9°S～41°N、99°～146°E）的試驗(yàn)結(jié)果，該模式使用C網(wǎng)格和垂向S坐標(biāo)，水平網(wǎng)格分辨率為（1／12）°，垂向分32層，其中沒有加入同化；NCEP／NCAR再分析風(fēng)場（1958—2000年）和混合的表面風(fēng)場（2000—2008年）被用于驅(qū)動模式；Simple Oceanic Data Assimilation（SODA）資料被用于確定模式初始場和邊界條件。雖然一共有50a的模式輸出結(jié)果，但為了與觀測資料對比，本文僅選取2000—2008年這9a的數(shù)據(jù)進(jìn)行分析。

用來驗(yàn)證模式的表面流場數(shù)據(jù)由法國的Archiving，Validation，and Interpretation of Satellite Oceanographic Data（AVISO）提供，其時間分辨率為7d，空間分辨率為（1／3）°；AVHRR 溫度資料由 NOAA提供，時間分辨率為1d，空間分辨率為（1／4）°。這些數(shù)據(jù)的時間范圍都為2000—2008年。

圖1比較了分別由AVISO表層流場和模式表層及次表層流場計(jì)算得到的9a平均的黑潮主流軸。從圖1中可見，模式中的黑潮次表層流軸與觀測的相當(dāng)一致，但表層流軸系統(tǒng)性地偏西。這可能是因?yàn)槟Ｊ降谋韺恿鲌龀剞D(zhuǎn)流外還包含風(fēng)生艾克曼輸運(yùn)，而后者的平均流向在此區(qū)域?yàn)槲飨颍ㄒ驗(yàn)槎撅L(fēng)強(qiáng)于夏季風(fēng)）；但AVISO的表層流只有地轉(zhuǎn)流，因而與模式的次表層流更加一致?？偟膩碚f，模式結(jié)果基本合理，可以被進(jìn)一步用來研究黑潮的變化。

1.2 方法

黑潮主流軸在某一特定時刻的位置通常是通過尋找一定空間范圍內(nèi)的速度最大值來確定［7－8］。然而，由于黑潮在呂宋海峽附近失去海岸約束，再加上地形和渦旋的作用，流態(tài)變化劇烈，僅用流速最大值難以準(zhǔn)確找出黑潮的真正主軸。因此，本文根據(jù)黑潮的物理特性開發(fā)了一種自動檢測黑潮主軸的方法（圖2），檢測區(qū)域范圍為16.5°（呂宋島以東）～24.8°N（臺灣島以東）、118°～130°E。檢測黑潮主軸的具體步驟如下：

圖2 檢測黑潮主軸的方法示意圖Fig.2 Schematic map of a method to detect Kuroshio axis

（2）在每一個緯向網(wǎng)格線上將搜索半徑限定在黑潮平均主軸（圖1）的左右300km之內(nèi)，進(jìn)而找出搜索半徑內(nèi)V最大值的位置；

（3）為保證這一位置為黑潮主軸位置，還須滿足3個條件（圖2）：（a）與前一個緯向網(wǎng)格線上（上游）黑潮主軸位置的緯向距離不能超過n1；（b）2個連續(xù)時刻找出的黑潮主軸位置間的距離不能超過n2；（c）此處的溫度（T2）和鹽度（S2）與上游黑潮主軸的溫度（T1）和鹽度（S1）必須滿足

的條件，即保持水團(tuán)的相似性；否則將重新計(jì)算。其中參數(shù)n1、n2和n3通過敏感性試驗(yàn)獲得，在本文中分別取為12，12和1.5。敏感性試驗(yàn)的過程通過不停地修改3個參數(shù)的值來比對已經(jīng)得出的黑潮主軸和背景流場的一致性，最終通過觀察得出在研究時間尺度內(nèi)時空上合理的黑潮主軸，然后確定最終參數(shù)。

圖3顯示出采用上述方法尋找得到的2000年1月1日的黑潮流軸。由圖3可見，該方法很好地抓住了黑潮的空間形態(tài)，黑潮流軸（圖3中的黑線）不僅與背景流場對應(yīng)，還與表層溫度和鹽度的分布高度相關(guān)。

2 黑潮流軸變化特征

圖3 由自動檢測方法確定的2000年1月1日的黑潮表層流軸（黑線）Fig.3 Surface axis of Kuroshio（black curve）on January 1，2000captured by the automatic detection method

黑潮流軸在呂宋海峽附近的變化可以從不同層次和不同時間尺度上進(jìn)行研究。本文分別選擇10m和200m這2個代表水深以及年際、季節(jié)和季節(jié)內(nèi)3個時間尺度來分析黑潮流軸的變化特征。為簡單起見，季節(jié)變化定義為9a平均的氣候態(tài)季節(jié)變化，年際變化為減去季節(jié)變化并作3個月滑動平均后的變化，季節(jié)內(nèi)變化則為去除上述季節(jié)和年際變化后剩余的變化。

2.1 年際變化

如圖4所示，黑潮流軸在表層（水深為10m）和次表層（水深為200m）都有明顯的年際變化，但僅限于呂宋海峽（右圖），且表層的變化大于次表層。在表層，2000年和2001年黑潮在呂宋海峽北部向西入侵南海，2003、2006和2008年黑潮在呂宋海峽中部和北部向東漂移。在次表層，黑潮僅2000年和2007年在呂宋海峽北部出現(xiàn)，明顯向東漂移的跡象則只在2008年出現(xiàn)。與Ni~no3指數(shù)（左圖中的黑線）對比顯示，黑潮流軸的年際變化與ENSO存在一定的聯(lián)系，在Ni~no3指數(shù)變化大的年份（2002、2006和2008年）黑潮流軸偏移變化都較其它時期大，有科學(xué)家指出在2006—2008年黑潮流量變化明顯，但由于數(shù)據(jù)長度不夠難以量化表征。

2.2 季節(jié)變化

如圖5所示，黑潮流軸的季節(jié)變化在表層與次表層比較一致，其幅度略大于年際變化，但與年際變化一樣主要出現(xiàn)在呂宋海峽。在表層，黑潮流軸夏季在呂宋海峽北部向東漂移，秋季在呂宋海峽南部入侵，冬季在呂宋海峽北部入侵。在次表層，黑潮流軸夏季東移且較弱，冬季西進(jìn)則更強(qiáng)。季節(jié)變化似乎與局地經(jīng)向風(fēng)（左圖中黑線）有關(guān)：春、夏季西南季風(fēng)盛行時黑潮向東偏移，秋、冬季東北季風(fēng)盛行時黑潮向西入侵，這與 LIANG et al［9］研究的結(jié)果一致。

圖4 黑潮流軸在呂宋海峽（18.5°～22°N）表層（上圖）和次表層（下圖）的年際偏移變化Fig.4 Interannual variations of the Kuroshio axis near the Luzon Strait（18.5°N～22°N）at 10m （upper panels）and 200m （lower panels）depths

2.3 季節(jié)內(nèi)變化

如圖6所示，黑潮流軸的季節(jié)內(nèi)變化比年際和季節(jié)變化強(qiáng)烈得多，且不限于呂宋海峽（圖6右）。黑潮流軸在呂宋海峽的季節(jié)內(nèi)變化（圖6左）似乎存在季節(jié)和年際尺度上的調(diào)制，例如在2003年以后表層與次表層的變化相當(dāng)一致，但在此之前則相差較大?？偟膩碚f，無論在表層還是次表層，黑潮在季節(jié)內(nèi)尺度上向西入侵的幅度大于向東的移動。在呂宋海峽，黑潮流軸的季節(jié)內(nèi)變化占其總變化方差的80%，表明 季節(jié)內(nèi)變化是此處黑潮變化最為重要的分量。

3 機(jī)制討論

3.1 風(fēng)和黑潮慣性的影響

為了探討局地風(fēng)場驅(qū)動的艾克曼輸運(yùn)和黑潮本身慣性對呂宋海峽黑潮流軸變化的可能影響，本文沿120.75°E對呂宋海峽（18.5°～22°N）的黑潮流軸變化和經(jīng)向風(fēng)分別做平均，察看不同水深的流軸與風(fēng)場之間的關(guān)系；同時在呂宋島東側(cè)黑潮進(jìn)入?yún)嗡魏{之前的位置（18°N）取1個流速斷面，觀察該斷面不同水深的黑潮經(jīng)向流速對呂宋海峽黑潮流軸變化的影響。

在年際尺度上，呂宋海峽黑潮流軸的變化與局地的經(jīng)向風(fēng)在表層有較好的正相關(guān)，即北（南）風(fēng)大時黑潮流軸向西（東）偏，相關(guān)系數(shù)隨著深度的增加而降低；黑潮流軸的變化與上游黑潮流速也有較好的正相關(guān)，即向北的黑潮流速大（?。r黑潮流軸偏東（西），相關(guān)系數(shù)卻隨著深度的增加而增大（圖7a）。在季節(jié)尺度上也存在類似的關(guān)系，即黑潮流軸在表層與局地風(fēng)場高度正相關(guān)而在次表層與上游黑潮流速高度正相關(guān)，只是風(fēng)的影響似乎限制在一個更薄的表層內(nèi)（圖7b）。在季節(jié)內(nèi)尺度上，黑潮流軸的變化與局地風(fēng)場和上游黑潮流速的相關(guān)性都很低，說明黑潮在呂宋海峽的高頻變化（也是黑潮的主要變化）基本不受表層艾克曼輸運(yùn)或黑潮本身慣性的控制（圖7c）。

這些結(jié)果可簡單解釋如下：對于黑潮的低頻（季節(jié)到年際）變化來說，在準(zhǔn)地轉(zhuǎn)理論的框架內(nèi)（如SHEREMET［10］所論述），西邊界流（黑潮）在流經(jīng)一個較大的豁口（如呂宋海峽）時，其路徑將在β效應(yīng)驅(qū)使下向西偏，而西偏的程度取決于流速（慣性）的大?。毫魉俸艽髸r可直接躍過豁口，流速小時則進(jìn)入豁口形成“流套”，這就是呂宋海峽黑潮流軸變化在年際和季節(jié)尺度上與上游流速呈正相關(guān)的原因。另一方面，由風(fēng)驅(qū)動的表層艾克曼輸運(yùn)則是呂宋海峽黑潮流軸變化與局地風(fēng)呈正相關(guān)的原因，這也說明為什么與風(fēng)的相關(guān)隨水深遞減而與上游流速的相關(guān)卻隨水深遞增。對于黑潮的高頻（季節(jié)內(nèi)）變化來說，最大的影響可能來自渦旋，而進(jìn)入?yún)嗡魏{的渦旋主要是由環(huán)流的不穩(wěn)定產(chǎn)生的，因而與局地風(fēng)場和黑潮慣性的關(guān)系不大。

圖7 黑潮流軸在呂宋海峽（18.5°～22°N）不同水深中的年際（a）、季節(jié)（b）和季節(jié)內(nèi)（c）變化與當(dāng)?shù)亟?jīng)向（V向）風(fēng)速和上游黑潮（V向）流速的相關(guān)系數(shù)Fig.7 Correlations of interannual，seasonal and intra－seasonal variations of Kuroshio axis at the Luzon Strait（18.5°N～22°N）with local meridional wind and up－stream Kuroshio velocity

3.2 渦旋的分布和影響

海洋中的渦旋以接近于羅斯貝波的速度向西移動［11］，但目前對渦旋到達(dá)西邊界之后的命運(yùn)特別是與西邊界流的相互作用機(jī)制仍然不清楚。前人的研究結(jié)果表明，渦旋在呂宋海峽呈現(xiàn)不同的形態(tài)［12－14］，且與黑潮發(fā)生強(qiáng)烈的相互作用。鑒于渦旋對黑潮季節(jié)內(nèi)變化的重要影響，本文采用一個自動探測渦旋的方法［15－16］察看呂宋海峽周邊（16.5°～30°N、120.5°～130°E）渦旋的分布情況，并探討渦旋運(yùn)動與呂宋海峽黑潮流軸變化的關(guān)系。

圖8 2000年1月1日呂宋海峽周邊海域的渦旋分布Fig.8 Eddy distribution in the vicinity of the Luzon Strait on January 1，2000

作為一個例子，圖8是用自動探測渦旋方法找出的渦旋在2000年1月1日的分布，表明該方法能有效地檢測聚集在呂宋海峽周邊海域的大量渦旋。圖9進(jìn)一步顯示鄰近呂宋海峽的1個方形區(qū)域（18°～23°N、120°～126°E）2000—2008年渦旋的收支情況。從圖9中可見：（1）由于黑潮的平流作用，北邊界出去的渦多于進(jìn)來的渦，南邊界進(jìn)來的渦多于出去的渦；（2）由于渦旋的西傳特性，東邊界進(jìn)來的渦多于出去的渦，西邊界出去的渦多于進(jìn)來的渦；（3）在區(qū)域內(nèi)部生成的渦遠(yuǎn)多于邊界上進(jìn)來的渦，且這些渦的絕大部分都消亡在區(qū)域里，說明這里大部分渦旋的生命周期和移動距離有限；（4）由于區(qū)域內(nèi)生成的渦多于消失的渦，所以從北邊界和西邊界出去的渦分別多于從南邊界和東邊界進(jìn)來的渦。氣旋和反氣旋渦都具有這些收支特征，只是后者在區(qū)域內(nèi)的凈生成量大于前者，因而從邊界出去的也更多。呂宋海峽鄰近區(qū)域既是渦旋的“溫床”，也是渦旋的“墓地”［12］；黑潮里的渦旋有可能進(jìn)入南海［13］，也有可能被黑潮帶走或消亡在黑潮里［14］。

為了進(jìn)一步察看黑潮與渦旋的相互作用，圖10比較了穿過圖9中120°E和123°E兩條邊界的渦旋數(shù)量與相對應(yīng)的黑潮流軸變化次數(shù)。當(dāng)渦旋從西向東運(yùn)動時，無論在120°E或123°E，無論是反氣旋渦還是氣旋渦，黑潮流軸在大多數(shù)（70%～80%）情況下向 東漂移；反之，當(dāng)渦旋從東向西運(yùn)動時，黑潮流軸在大多數(shù)（75%～85%）情況下向西漂移。由于在同一時刻穿過黑潮兩邊的渦旋可能不止1個，呂宋海峽黑潮流軸的變化不一定與單個渦旋的進(jìn)出對應(yīng)，但圖10中渦旋的進(jìn)出與黑潮的變化基本一致，表明渦旋是影響呂宋海峽黑潮流軸擺動的主要原因。又因?yàn)闇u旋以西向運(yùn)動為主，所以黑潮流軸向西擺動的次數(shù)明顯多于向東擺動的次數(shù)。

4 結(jié)論

本文以高分辨率的區(qū)域海洋環(huán)流模式（ROMS）的9a（2000—2008年）輸出結(jié)果為基本數(shù)據(jù)，利用自行編制的基于速度場的流軸自動檢測方法，分析了黑潮流軸在呂宋海峽附近的位置及其在年際、季節(jié)和季節(jié)內(nèi)時間尺度上的變化特征，并初步探討了引起這些變化的物理過程和機(jī)制。得出的主要結(jié)論如下：

（1）黑潮流軸在呂宋海峽附近具有明顯的年際、季節(jié)和季節(jié)內(nèi)變化，其中季節(jié)內(nèi)變化最為強(qiáng)烈；總的來說黑潮流軸在表層的變化大于次表層的變化，向西擺動的次數(shù)多于向東擺動的次數(shù)。

（2）在年際和季節(jié)時間尺度上，黑潮流軸在表層主要受局地風(fēng)場驅(qū)動的艾克曼漂流影響，北風(fēng)強(qiáng)時西偏，南風(fēng)強(qiáng)時東偏；而次表層黑潮流軸的變化則主要由黑潮本身的慣性決定，上游流速小時西偏拐入南海，上游流速大時采取更為直接的路徑越過呂宋海峽北上。

（3）在季節(jié)內(nèi)時間尺度上，黑潮流軸的變化主要由渦旋與黑潮的相互作用控制，當(dāng)渦旋（包括氣旋和反氣旋渦）向東（向西）運(yùn)動時，黑潮流軸在大多數(shù)情況下向東（向西）漂移；由于渦旋以西向運(yùn)動為主，所以黑潮流軸的變化以向西的移動居多。

（4）呂宋海峽附近海域的渦旋收支分析表明，這里既是產(chǎn)生渦的“溫床”，也是渦消失的“墓地”；由于產(chǎn)生的渦多于消失的渦，所以從北邊和西邊出去的渦多于從南邊和東邊進(jìn)來的渦；隨黑潮北上的渦可以對下游產(chǎn)生影響，而向西移動的渦可以將黑潮水帶入南海。

（References）：

［1］GAN J P，LI H，CURCHITSER E N，et al.Modeling South China Sea circulation：Response to seasonal forcing regimes［J］.J Geophys Res，2006，111：C06034.

［2］BARKLEY R A.The Kuroshio current［J］.Science Journal，1970，7302：54－55.

［3］NITANI H.Oceanographic conditions in the sea east of the Philippines and Luzon Strait in the summer of 1965and 1966［C］∥MARR J.The Kuroshio—A Symposium on the Japan Current.Honolulu：East West Center Press，1970：213－232.

［4］ QIU De－zhong，YANG Tian－h(huán)ong，GUO Zhong－xin.A west flowing current in the northern part of the South China Sea in summer［J］.Tropic Oceanology，1984，3（4）：65－73.

仇德忠，楊天鴻，郭忠信.夏季南海北部一支向西流動的海流［J］.熱帶海洋，1984，3（4）：65－73.

［5］HUANG Qi－zhou，ZHENG You－ren.Bashi Channel currents in northeastern of South China Sea in March 1992［C］∥Northeast of South China Sea Circulation Research Group.China oceanography corpus：edition 6.Beijing：Ocean Press，1996：42－52.

黃企洲，鄭有任.1992年3月南海東北部巴士海峽的海流［C］∥南海東北部海區(qū)環(huán)流課題組.中國海洋學(xué)文集：第六集.北京：海洋出版社，1996：42－52.

［6］EBUCHI N，HANAWA K.Influence of mesoscale eddies on variations of the Kuroshio path south of Japan［J］.Oceanogr，2003，59：25－36.

［7］DAISUKE A，SHIRO I，HIROSHI U，et al.Estimating the Kuroshio Axis South of Japan using combination of satellite altimetry and drifting buoys［J］.Journal of Oceanography，2004，60：375－382.

［8］LIU Z Q，GAN J P.Variability of the Kuroshio in the East China Sea derived from satellite altimetry data［J］.Deep－Sea Research I，2012，59：25－36.

［9］LIANG W D，YANG Y J，TANG T Y，et al.Kuroshio in the Luzon Strait［J］.Journal of Geophysical Research，2008，113：17－18，C08048，doi：10.1029／2007JC004609.

［10］SHEREMET V A.Hysteresis of a Western Boundary Current Leaping across a Gap［J］.Phydical Oceanography，2000：1 247－1 255.

［11］CHELTON D B，SCHLAX M G，SAMELSON R M，et al.Global observations of large oceanic eddies［J］.Geophys Res Lett，2007，34：L15606，doi：10.1029／2007GL030812.

［12］ZHAI X，JOHNSON H L，MARSHALL D P.Significant sink of ocean－eddy energy near western boundaries［J］.Nat Geosci，2010，3：608－612.

［13］SHEU W J，WU C R，OEY L Y.Blocking and westward passage of eddies in the Luzon Strait［J］.Deep Sea Res：Part II，2010，57：1 783－1 791.

［14］HE Y，CAI S，WANG S.The correlation of the surface circulation between the Western Pacific and the South China Sea from satellite altimetry data［J］.Int J Remote Sens，2010，31：4 757－4 778.

［15］NENCIOLI F，DONG C，DICKEY T，et al.A vector geometry based eddy detection algorithm and its application to high－resolution numerical model products and high－frequency radar surface velocities in the Southern California Bight［J］.Atmos Ocean Tech， 2010，27 （3）：564－579，doi：10.11.1175／2009JTECHO725.1.

［16］LIU Yu，DONG Chang－ming，GUAN Yu－ping，et al.Eddy analysis in a zonal band of the North Pacific Ocean［J］.Deep－Sea Research I，2012，68：54－67.