譚 偉 左軍成 李 娟 陳美香 楊逸秋

(1. 河海大學(xué) 南京 210098; 2. 國家海洋局南海預(yù)報中心 廣州 510310; 3. 國家海洋環(huán)境預(yù)報中心 北京 100091)

潮流是潮波運動的一種表現(xiàn)形式。20世紀80年代以來, 對東中國海潮流的研究從傳統(tǒng)的觀測分析向數(shù)值模擬發(fā)展(萬振文等, 1998; 沙文鈺, 2000)。采用數(shù)值模式計算東中國海的全日及半日潮流, 證實M2分潮流在北黃海山東北部近海及南黃海北部海域各存在一個圓流點(趙保仁等, 1994)。隨著高性能計算機的發(fā)展, 為了能更精確地反映潮流的運動規(guī)律,潮波數(shù)值模擬已經(jīng)從二維模型發(fā)展到三維模型(萬振文等, 1998; 王凱等, 1999; 沙文鈺, 2000; 朱學(xué)明等,2012)。利用POM三維數(shù)值模式模擬東中國海潮波系統(tǒng), 發(fā)現(xiàn) M2分潮在舟山群島附近存在一個圓流點,M1分潮在黃海北部和臺灣海峽各存在一個圓流點(萬振文等, 1998)。采用高分辨率三維模式進行模擬, 可以解釋遼東灣中M2分潮無潮點的出現(xiàn)以及黃海中出現(xiàn)的流速上下一致的現(xiàn)象(Guo et al, 1998)。

在近岸海域, 由于水深、地形、緯度不同和天體引潮力的驅(qū)動使潮流能量存在顯著的區(qū)域差異(Bruno et al, 2000; Tsimplis, 2000), 例如, 在直布羅陀海峽、海底門檻以半日分潮M2和S2為主, 而從德克薩斯州到美國路易斯安那州的大陸架則以全日分潮K1為主(DiMarco et al, 1998)。北海潮流垂向剖面觀測結(jié)果表明, 潮流振幅隨著水深的增大表現(xiàn)出減小的趨勢, 潮流位相隨著水深的增大而增大, 由于對潮流垂向結(jié)構(gòu)的分辨率較低, 因此不能很好地分析研究潮流的垂向結(jié)構(gòu)及影響因素(Maas et al, 1987)。影響潮流結(jié)構(gòu)的主要因素有摩擦力、層化和低頻海流的波動(Souza et al, 1996), 其中, 底摩擦力顯著影響潮流的垂向結(jié)構(gòu)和流速大小(Fan et al, 2009), 在黃海,底摩擦效應(yīng)會使 M2分潮流的流速迅速減小 20%—40%以及海洋底部的 M2分潮流的位相增加 15—20min。在淺水海域, 地形是引起潮流垂向結(jié)構(gòu)變化的主導(dǎo)因素, 層化主要影響潮流的垂向切變和季節(jié)變化(Lee et al, 2011), 在墨西哥瓜伊馬斯灣, 層化會導(dǎo)致全日潮流在海洋上層逆時針旋轉(zhuǎn), 而在海洋下層則是順時針旋轉(zhuǎn)(Gómez-Valdés et al, 2012)。

在潮流垂直結(jié)構(gòu)中, 分潮流最大流速發(fā)生時刻隨深度的加深而不斷提前, 不僅與摩擦有關(guān), 而且還受到水深以及科氏參量的影響(葉安樂, 1984a)。水深較淺的海域, 上層最大潮流時刻隨深度的變化率絕對值較大, 而水深較深的海域, 上層最大潮流時刻變化率的絕對值較小, 兩者在底層的變化率都比較接近(方國洪, 1984; 葉安樂, 1984a)。當(dāng)科氏參量f大于分潮角速率σ時, 潮流最大流速方向隨著深度的增加向左偏(Fang et al, 1983; 葉安樂, 1984b)。外力強迫下,層結(jié)的海洋會激發(fā)出各種形式的內(nèi)波, 自由傳播的內(nèi)波其低頻部分的構(gòu)成必須與浮力頻率所決定的一組本征模態(tài)相符合(方國洪, 1984)。海洋內(nèi)部波動模態(tài)可以改變海流的垂直結(jié)構(gòu), 使海流表現(xiàn)出很強的斜壓性(葉安樂, 1984c), 在直布羅陀海峽, 垂向剖面中的內(nèi)波模態(tài)可以解釋潮流半長軸的位相隨深度的變化(Tsimplis, 2000)。Zuo等(2007)利用一個月的海流觀測數(shù)據(jù), 分析得到葡萄牙北部大陸架海域的潮流特征, 并使用 PWP模式探討風(fēng)應(yīng)力對流速的影響。在東海大陸架邊緣海域潮流橢圓隨深度旋轉(zhuǎn), 且潮流以順時針旋轉(zhuǎn)為主, 反映了潮流是正壓流的特點(Li et al, 2012)。

這些研究成果加深了人們對東中國海潮波的認識。但已有的東海實測潮流資料主要集中于近岸海域和東北部海域, 東海東南部大陸架邊緣海域的實測潮流資料還比較少。本文利用 973計劃“近海及鄰近大洋動力環(huán)境低頻變異的可預(yù)測性研究”項目投放在東海的海床基觀測系統(tǒng)得到的 6個多月的實測潮流數(shù)據(jù), 研究該海域的潮流垂向結(jié)構(gòu), 探討層化對潮流垂向結(jié)構(gòu)的影響。

圖1 觀測站點的位置和東中國海水深(單位: m)Fig.1 The observation site and bathymetric map of the Bohai Sea, Yellow Sea, and East China Sea

1 數(shù)據(jù)

海流資料來源于東海海床基觀測系統(tǒng)的 ADCP(聲學(xué)多普勒流速剖面儀) 觀測結(jié)果。觀測點位于東海大陸架邊緣海域(26°30.052′N, 122°35.998′E), 水深約為113m(圖1)。采樣時間從2008年9月20日5:30到2009年4月22日11:00, 實際有效觀測時間是從2008年9月20日6:00到09年4月6日7:00, 采樣時間間隔為10min, 有效數(shù)據(jù)共28519組。垂向空間采樣間隔為8m, 共13層, 實際可供分析的連續(xù)資料有 11層, 分別位于 24、32、40、48、56、64、72、80、88、96m和104m層上。由于天氣和儀器原因, 海流數(shù)據(jù)在不同的時間段和不同層次皆有缺測, 缺測部分用鄰近時刻觀測值線性內(nèi)插加以彌補。水深資料來源于海軍航海保證部提供的海圖水深。本文采用的Argo溫鹽數(shù)據(jù)是基于Argo浮標資料、Triangle Trans-Ocean Buoy Network(TRITON)和CTD數(shù)據(jù)進行融合插值后形成逐月的格點數(shù)據(jù), 空間分辨率為 1°×1°,垂向 25層分別為10、20、30、50、75、100、125、150、200、250、300、400、500、600、700、800、900、1000、1100、1200、1300、1400、1500、1750m和 2000m, 時間從 2001年 1月開始, 下載網(wǎng)址為http://www.argo.ucsd.edu/Gridded_ fields.html。本文選取的是觀測點附近的一個格點(26.5°N, 122.5°E)上125m層的溫鹽數(shù)據(jù), 時間段是 2008年 9月至 2009年4月。

2 分析結(jié)果

2.1 密度的垂向結(jié)構(gòu)

夏季, 充分的潮汐能和太陽熱能傳輸?shù)胶Ｑ蟊韺? 導(dǎo)致大陸架海洋層化加強(Howarth, 1998)。圖 2給出了利用觀測點附近的Argo溫鹽數(shù)據(jù)計算得到的位勢密度的垂向分布(2008年9月到次年4月)。除9月份外, 上混合層密度的垂向分布相對比較均勻, 而在密度躍層, 海水密度隨深度的增加顯著增大, 密躍層以深密度增大的趨勢減緩, 到底層密度達到最大,而9月份海水上層溫度仍比較高, 密度較低, 層化加強, 密躍層強度大, 幾乎沒有上混合層的存在, 海水密度在整個深度上遞增。隨著時間推移, 海水上層溫度降低, 密度逐月增大, 出現(xiàn)上混合層, 并且隨著對流混合的加強, 上混合層有逐漸加深的趨勢, 到次年2月份密度達到最大, 密躍層幾乎消失, 海水密度在整個深度上變化較小。3月份密度垂向結(jié)構(gòu)與2月份相近, 但上層海水密度開始減小, 到4月份上層海水密度已明顯減小, 躍層重新出現(xiàn)。

圖2 觀潮點附近的位勢密度垂向分布Fig.2 Profiles of the potential density nearby the observation site

圖3 觀測點的浮力頻率垂向分布Fig.3 Profiles of the buoyancy frequency nearby the observation site

2.2 平均流

本文分別對海流的x方向流速U分量和y方向流速V分量在觀測期間躍層強度最大的9月份和最弱的2月份做時間平均(圖4), 在各觀測層上U分量明顯大于V分量, 海流基本為偏東的東北向流。海洋密度層結(jié)對流速結(jié)構(gòu)的影響非常明顯, 9月份幾乎不存在上混合層, 在56m以淺U分量隨深度的增加而增大, 在躍層最強處達到最大, 躍層以深流速變化較小,在底層, 流速因受底摩擦的影響而迅速減小; 2月份,位勢密度躍層最弱(圖1), U分量只在近底層迅速減小,在其它深度上變化不大(圖4)。與U分量類似, V分量的變化同樣以56m為界。56m以淺V隨深度略微減小; 在躍層以深, 9月份V分量隨深度的增大而增大,而2月份變化很小; 9月份和2月份V分量在底層都迅速減小。以上結(jié)果表明該海域海流的動量受到密度層結(jié)的約束, 層化加強減小了動量和質(zhì)量的垂向湍流交換, 導(dǎo)致密度垂向微分在層化處加強(Howarth,1998)。圖 3中浮力頻率在密躍層深度處最大也可以由此得到解釋。

2.3 觀測點的潮流調(diào)和分析

2.3.1 垂向平均潮流 原始觀測的流速中包含潮流成分, 潮流的調(diào)和分析可以對實測潮流按各天文分潮的頻率進行分離。為了除去高頻噪音, 取 60min滑動平均, 采用Pawlowicz等(2002) 提供的潮流調(diào)和分析, 在 95%的置信度下得到垂向平均的 O1、K1、M2和S2這4個主要分潮的潮流橢圓要素(表1), 并繪制潮流橢圓(圖 4)。潮流橢圓短軸為正值, 表示潮流矢量隨時間逆時針旋轉(zhuǎn), 短軸為負表示潮流矢量隨時間順時針旋轉(zhuǎn)。從表1和圖5中可以明顯看出, 觀測站點以半日潮流為主, M2分潮最大, 其次為S2分潮,K1和 O1分潮則相對較小。潮汐性質(zhì)比為0.12, 小于0.5, 說明該觀測點所在海域為正規(guī)半日潮海域。除了 O1分潮外, 其余 3個分潮的潮流橢圓均做順時針旋轉(zhuǎn)。這4個主要分潮的長軸指向基本一致, 均是西北—東南方向。

圖4 觀測點海流的U、V分量在9月份和2月份的垂向分布Fig.4 Profiles of the time-averaged U and V components of observed currents in the September and February nearby the observation site

表1 主要分潮的垂向平均潮流橢圓參數(shù)Tab.1 The vertical average tidal current ellipse parameters of the main constituent

2.3.2 潮流的垂向結(jié)構(gòu) 在淺海海域?qū)踊茱@著影響潮汐混合和潮流的垂向結(jié)構(gòu), 層化的影響是減小密度躍層內(nèi)的擾動/湍流, 導(dǎo)致潮流振幅、位相和傾角在通過密度躍層時發(fā)生較大改變(Howarth, 1998)。

本文對9月、12月、次年2月和4月這4個典型月份的潮流進行調(diào)和分析, 得到主要分潮M2和K1的潮流橢圓要素, 并給出橢圓要素的垂向分布。

(1)橢圓長軸的垂向結(jié)構(gòu)

圖6給出了M2和K1分潮的潮流橢圓長軸的垂向分布。半日分潮M2潮流的最大流速(30—60cm/s) 遠遠大于全日分潮K1(2—12cm/s), 9月份, M2分潮流的長軸在80m以淺的深度上變化不大, 在43cm/s 附近振蕩, 從表層開始先隨深度的增加略微增大, 在 32m處存在一個極大值, 然后隨深度的增大而減小, 在密躍層深度(48m)處有一極小值, 再然后隨深度顯著增大, 在 88m處達到最大值為 47.44cm/s, 由于受底摩擦的影響, 在近底層迅速減到最小為37.07cm/s(圖6),這與方國洪得到的理論結(jié)果吻合(方國洪, 1984)。12月份, M2分潮流的長軸的垂向結(jié)構(gòu)與9月份相近, 但隨著冬季的到來, 對流混合加強, 密度躍層強度減弱,相比9月份, M2分潮流的長軸在整個深度上減小, 且在80m以淺的深度上垂向變化與9月份類似, 但變化幅度略小(圖6)。2月份, M2分潮流的長軸在垂向上幾乎沒有變化, 大小和12月份相當(dāng); 2月過后, 整個深度上的M2長軸開始增大, 4月份的長軸已明顯大于9月, 80 m以淺的垂向變化率也在增大, 并且伴隨著躍層的形成及加深, 極大、極小值出現(xiàn)的深度也隨之下移(圖 6); 另外, M2分潮流的長軸由于受底摩擦的影響在近底層迅速減小(Souza et al, 1996)。

圖5 垂向平均潮流橢圓Fig.5 Vertical average current ellipse nearby the observation site

圖6 潮流橢圓半長軸的垂向分布Fig.6 Profiles of major axes of current ellipse nearby the observation site

圖7 潮流橢圓率的垂向分布Fig.7 Profiles of tide ellipticity nearby the observation site

全日分潮流 K1的長軸大體表現(xiàn)為隨深度增加而遞減的趨勢, 各個月份的K1長軸在64m或72m層處有一個明顯的極小值, 然后隨深度增加而增大, 接近海底的深度上迅速減小, 最小可達2—3cm/s。全日潮流的強流會在上混合層出現(xiàn)的部分原因是日周期的冷卻、加熱以及風(fēng)混合的共同作用(Price et al, 1987) ,這在一定程度上也解釋了圖6中上混合層的K1長軸在2月最小, 之后開始增大這一現(xiàn)象。

(2)橢圓率的垂向結(jié)構(gòu)

橢圓長半軸與短半軸的比值被定義為橢圓率。橢圓率為正表示潮流按逆時針方向旋轉(zhuǎn), 為負表示潮流按順時針方向旋轉(zhuǎn), 橢圓率為0表示往復(fù)流, 此時潮流橢圓退化成線段, 為1或–1時表示橢圓變成圓周(DiMarco et al, 1998)。潮流橢圓率k隨深度增加而增大, 當(dāng)上層潮流左旋(k＞0) 時, 到下層橢圓將變的較圓(k≈1); 當(dāng)上層潮流右旋(k＜0) 時, 到下層橢圓將變的較扁(k～0) , 甚至有可能變成左旋(k＞0) ( Fang et al,1983)。M2半日分潮的橢圓率始終為負值, 并且變化很小, 均在–0.2-0之間變化, 說明 M2潮流橢圓在整個深度上順時針旋轉(zhuǎn), 并且橢圓較扁, 9月份的M2潮流在密度躍層(56m) 處幾乎退化為往復(fù)流(圖7)。

K1全日分潮的橢圓率隨接近海底而呈現(xiàn)出波動增加的趨勢, 上層橢圓率 k為負, 潮流順時針旋轉(zhuǎn),在底層 9月份的橢圓率接近 0, 其余月份均為正, 即底層潮流橢圓變的較扁, 甚至是逆時針旋轉(zhuǎn)(圖 7),與方國洪(1984) 得到的理論結(jié)果一致。在 9月和 4月層化較強的月份, 橢圓率在密躍層以深到 88m處在0附近出現(xiàn)震蕩。

(3)橢圓傾角的垂向結(jié)構(gòu)

對于固定觀測點, 由于密度水平差異產(chǎn)生的斜壓海流導(dǎo)致潮流橢圓傾角θ隨水深發(fā)生旋轉(zhuǎn), 且傾角的垂向結(jié)構(gòu)與科氏力也有關(guān)。當(dāng)相應(yīng)的天文潮頻率σ大于局地慣性頻率 f時, 橢圓長軸方向變化不大, 如水深足夠大, 從表層到中層潮流橢圓略向右偏轉(zhuǎn), 到底層又有向左偏回來的趨勢; 當(dāng) σ小于或等于 f時,隨著水深的不斷加深, 最大流速方向逐漸向左偏(方國洪, 1984; 葉安樂, 1984)。該觀測點( 26°30.052′N,122°35.998′E) 的局地慣性頻率 f=0.0373/h, 而 M2天文分潮頻率天文潮頻率大于局地頻率f, M2傾角先向右略偏, 密躍層以深至底層左偏; K1天文分潮頻率與局地慣性頻率f相近,除4月份外, K1傾角隨著接近海底明顯左偏, 與上述理論結(jié)果吻合, 4月份相對異常分布的原因還需進一步的探討(圖8)。

圖8 潮流橢圓傾角的垂向分布Fig.8 Profiles of tide ellipse inclination nearby the observation site

3 結(jié)論

觀測點位于東海大陸架邊緣海域, 它的潮流特征能夠反映該海域的潮流特點, 本文利用 Argo數(shù)據(jù)分析觀測點附近位勢密度的垂向結(jié)構(gòu); 基于觀測點 6個多月的ADCP海流數(shù)據(jù)分析研究了觀測點潮流的垂向結(jié)構(gòu)及層化對其垂向結(jié)構(gòu)的影響, 得到以下結(jié)論

(1)海水密躍層在9月份強度最大, 此時, 位勢密度在表層最大, 并隨著深度的增加而逐漸遞減; 隨著冬季的到來, 對流混合逐漸增強, 最終在整個深度上混合均勻, 此時密度躍層消失, 而位勢密度逐漸增加,到次年2月份達到最大; 從4月份開始, 海水位勢密度明顯減小, 密度躍層重新出現(xiàn)。相應(yīng)的浮力頻率在夏季達到最大值冬季最小, 且在各個月份的密度躍層深度處最大。

(2)觀測點海域的潮汐性質(zhì)比為 0.12, 屬于正規(guī)半日潮海域, 其中M2分潮最大, S2次之, K1和O1相對較小; 4個分潮的最大流速方向基本一致, 均為西北－東南方向; 分別對4個典型月份的海流數(shù)據(jù)進行潮流調(diào)和分析, 得到 M2和 K1分潮流的潮流橢圓要素,發(fā)現(xiàn)本文結(jié)果與前人的理論結(jié)果吻合較好。

a. M2分潮流橢圓長軸的垂向結(jié)構(gòu)在各個月份相近, 潮流橢圓長軸在 88m 以深直到底層由于底摩擦的影響迅速減到最小。冬季隨著對流混合增強, 海水層化減弱, 潮流橢圓長軸在整個深度上同步減小; 從3月份開始, 潮流橢圓長軸開始逐漸增大, 到 4月份達到最大。K1分潮的潮流橢圓長軸在上混合層較大,可能與日周期的冷卻、加熱和風(fēng)混合的共同作用有關(guān),與M2相似, K1分潮的潮流橢圓長軸在底層達到最小。

b. M2分潮的潮流橢圓率在整個深度上變化不大,且在各個月份變化相似; 在海水層化效應(yīng)較強的月份, K1分潮的潮流橢圓率在整個深度上在0值附近振蕩隨接近海底而增加, 這與方國洪的理論結(jié)果吻合。

c. 潮流橢圓傾角的垂向結(jié)構(gòu)與相應(yīng)的天文潮頻率 σ和局地慣性頻率f的相對大小有關(guān)。觀測點M2分潮的頻率大于f, M2分潮流的橢圓傾角先向右略偏, 在密躍層以深海域向左偏; K1分潮的頻率與局地慣性頻率f相近, 潮流橢圓傾角接近海底時明顯左偏。

萬振文, 喬方利, 袁業(yè)立, 1998. 渤、黃、東海三維潮波運動數(shù)值模擬. 海洋與湖沼, 29(6): 611—616

王 凱, 方國洪, 馮士笮, 1999. 渤海、黃海、東海M2潮汐潮流的三維數(shù)值模擬. 海洋學(xué)報, 21(4): 1—13

方國洪, 1984. 潮流垂直結(jié)構(gòu)的基本特征——理論和觀測的比較. 海洋科學(xué), (3): 1—11

葉安樂, 1984a. 分潮流最大流速發(fā)生時刻隨深度變化的特征.海洋通報, 3(2): 1—8

葉安樂, 1984b. 頻率低于慣性頻率的潮流之垂直結(jié)構(gòu). 海洋湖沼通報, (3): 1—7

葉安樂, 1984c. 潮流橢圓長軸方向隨深度變化的特征. 海洋湖沼通報, (2): 1—6

朱學(xué)明, 鮑獻文, 宋德海等, 2012. 渤、黃、東海潮汐、潮流的數(shù)值模擬與研究. 海洋與湖沼, 43(6): 1103—1113

沙文鈺, 2000. 中國海潮波研究的進展. 海洋預(yù)報, 17(2):73—77

趙保仁, 方國洪, 曹德明, 1994. 渤、黃、東海潮汐潮流的數(shù)值模擬. 海洋學(xué)報, 16(5): 1—10

Bruno M, Ma?anes R, Alonso J J et al, 2000. Vertical structure of the semidiurnal tidal currents at Camarinal Sill, the strait of Gibraltar. Oceanologica Acta, 23(1): 15—24

DiMarco S F, Reid R O, 1998. Characterization of the principal tidal current constituents on the Texas-Louisiana shelf.Journal of Geophysical Research, 103(C2): 3093—3109

Fan X, Wei H, Yuan Y et al, 2009. Vertical structure of tidal current in a typically coastal raft-culture area. Continental Shelf Research, 29(20): 2345—2357

Fang G H, Ichiy T, 1983. On the vertical structure of tidal currents in a homogeneous sea. Geophysical Journal International, 73(1): 65—82

Gómez-Valdés J, Dworak J A, Vázquez H J et al, 2012. Vertical structure of tidal flows at the entrance to Guaymas Bay,Mexico. Geofísica Internacional, 51-3: 271—279

Guo X Y, Yanagi T, 1998. Three-dimensional structure of Tidal current in the East China Sea and the Yellow Sea. Journal of Oceanography, 54(6): 651—668

Howarth M J, 1998. The effect of stratification on tidal current profiles. Continental Shelf Research, 18(11): 1235—1254

Lee S, Lie H J, Cho C H et al, 2011. Vertical Structure of the M2Tidal Current in the Yellow Sea. Ocean Science Journal,46(2): 73—84

Li L, Jiang W W, Li P L et al, 2012. Vertical Structure of the Tidal Currents on the Continental Shelf of the East China Sea. Journal of Ocean University of China, 11(3): 347—353

Maas L R M, Van Haren J J M, 1987. Observations on the vertical structure of tidal and inertial currents in the central North Sea. Journal of Marine Research, 45(2): 293—318

Pawlowicz R, Beardsley B, Lentz S, 2002. Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE. Computers & Geosciences, 29(8): 929—937 Price J F, Weller R A, Schudlich R R, 1987. Wind-driven ocean currents and Ekman transport. Science, 238(4833): 1534—1538 Souza A J, Simpson J H, 1996. The modification of tidal ellipses by stratification in the Rhine ROFI. Journal of Continental Shelf Research, 16(8): 997—1007

Tsimplis M N, 2000. Vertical structure of tidal currents over the Camarinal Sill at the Strait of Gibraltar. Journal of Geophysical Research, 105(C8): 19709—19728

Zuo J C, Du L, Peliz A et al, 2007. The characteristics of Near-surface Velocity During the upwelling Season on the Northern Portugal Shelf. Journal of Ocean University of China, 6(3): 213—225