張菀倫, 劉志亮, 王世紅

?

南黃海渦動能及其譜特征的季節(jié)性變化

張菀倫1, 2, 劉志亮1, 王世紅1, 2

(1. 中國科學(xué)院海洋研究所海洋環(huán)流與波動重點開放實驗室, 山東青島 266071; 2. 中國科學(xué)院大學(xué), 北京 100049)

為進(jìn)一步了解南黃海渦動能的季節(jié)變化特點, 本文利用氣候態(tài)月平均數(shù)據(jù)分析了南黃海渦動場的季節(jié)變化規(guī)律, 并通過計算動能通量譜探討了相關(guān)的變化機(jī)制。結(jié)果顯示, 南黃海渦動能在季節(jié)變化尺度上存在雙峰結(jié)構(gòu), 9月和1月各存在一個峰值, 9月峰值較大, 1月峰值較小。而且, 上半年南黃海存在的渦主要為經(jīng)向伸長, 下半年主要為緯向伸長。動能通量譜的計算結(jié)果顯示, 9月存在的渦動能峰值主要是斜壓不穩(wěn)定過程直接作用的結(jié)果, 而1月的峰值主要與渦渦相互作用相關(guān)。本文從一個新的角度(渦動能)對南黃海進(jìn)行了探究, 有利于加深人們對該海域海水動力機(jī)制的認(rèn)識。

南黃海; 渦動能; 動能通量譜

南黃海由山東半島、蘇北平原和朝鮮半島所包圍, 是以山東半島的成山頭與朝鮮半島的長山串連線為北界, 以長江口至濟(jì)州島連線為南界的半封閉海域, 總面積30多萬平方公里, 平均水深46 m。黃海環(huán)流主要受大氣變動產(chǎn)生的風(fēng)應(yīng)力、潮余流、海面熱交換、周邊河流的淡水輸入以及海區(qū)南部開邊界等條件的影響[1-2]。南黃海緊連我國大陸, 并且沿海分布有很多我國重要的漁場、生態(tài)資源豐富, 因此南黃海的水文狀況對我國沿海生態(tài)環(huán)境及經(jīng)濟(jì)發(fā)展有重大影響。南黃海環(huán)流存在明顯的季節(jié)性變化, 冬、夏季環(huán)流的基本形態(tài)有很大差別。冬季環(huán)流基本上是由南黃海中部北上的黃海暖流與其兩側(cè)向南的沿岸流所組成, 黃海暖流及其余脈分別與東西沿岸流之間形成反氣旋和氣旋式環(huán)流[3-4]。而黃海夏季環(huán)流一般認(rèn)為是以熱鹽環(huán)流為主要成分。每當(dāng)夏季來臨, 逐漸產(chǎn)生與黃海底層冷水團(tuán)相適應(yīng)的黃海夏季環(huán)流, 其水平環(huán)流主要表現(xiàn)為上層氣旋式邊緣環(huán)流[5-12]。一般認(rèn)為, 該水平環(huán)流在靠近我國大陸一側(cè)大約沿40~50 m等深線向南流動。而在40 m以淺海域, 劉志亮等[13]認(rèn)為存在風(fēng)驅(qū)的北向流, 并且該觀點已基本被接受[14]。南黃海環(huán)流的季節(jié)變化, 必然與該海域不同運(yùn)動尺度間能量的傳輸過程密切相關(guān)。本文的目的是對南黃海渦動能的季節(jié)變化特征進(jìn)行分析, 并討論相關(guān)的變化機(jī)制。

1 研究區(qū)域及數(shù)據(jù)來源

為減小邊界以及周邊淡水輸入的影響, 并包含沿40~50 m等深線熱鹽環(huán)流及南黃海渦場的主要信號, 本文所選擇的研究區(qū)域為33.5~35.7°N, 121.5°~ 124.1°E的范圍(圖1黑框中所標(biāo)出的區(qū)域)。

由于近海捕撈作業(yè)等因素的影響, 使得近海海流的直接觀測很困難, 長時間的連續(xù)觀測資料更是稀少, 因此對近海環(huán)流的研究常常采用模式數(shù)據(jù)。為研究南黃海渦場的季節(jié)性變化及不同尺度渦之間的能量運(yùn)輸, 本文流場和溫度場數(shù)據(jù)采用的是區(qū)域海洋模式ROMS(Regional Ocean Modeling System)[15]輸出的在東中國海氣候態(tài)模擬的結(jié)果(模式模擬區(qū)域如圖1)。該模式以大氣風(fēng)場(QuikSCAT風(fēng)應(yīng)力)和熱通量(COADS 熱通量)為主要外強(qiáng)迫場, 考慮M2、S2、O1、K1 等10個影響黃海環(huán)流的主要分潮, 采用美國海軍1/12°太平洋海盆模擬結(jié)果為邊界條件進(jìn)行氣候態(tài)模擬。模式水平分辨率為1/12°, 垂直方向分為26層。

圖2a、圖2b分別為夏季垂向平均的黃海環(huán)流模擬結(jié)果和劉志亮等[13]利用夏季ADP實測資料給出的南黃海環(huán)流模式概念圖。兩圖對比可以看出模擬結(jié)果與劉志亮等[13]的結(jié)果比較吻合, 說明模擬結(jié)果較為可信。

2 南黃海渦動能季節(jié)性變化

為研究南黃海渦動能季節(jié)性變化的特點, 我們首先對整個研究區(qū)域的渦動能(Eddy Kinetic Energy, EKE)進(jìn)行平均。渦動能由緯向速度和經(jīng)向速度計算得出:

圖3繪制的是區(qū)域平均的渦動能的時間序列, 從圖中可以看出全年渦動能呈現(xiàn)一個雙峰值的結(jié)構(gòu), 峰值分別出現(xiàn)在冬季和夏季, 最大峰值出現(xiàn)在9月, 約為1.34×10–4(m/s)2, 次峰值出現(xiàn)在1月; 最低值出現(xiàn)在5月, 約為1.34×10–5(m/s)2, 最大值與最小值幾乎相差一個量級, 差值大約為1.2×10–4(m/s)2。

為了更好地研究南黃海渦動能在不同方向上的分布特征, 我們計算了該區(qū)域氣候態(tài)月平均的EKE能量譜密度(圖4), 其中k是緯向波數(shù),k是經(jīng)向波數(shù),是月份。

圖4顯示了1月~6月(除3月外)渦信號在k上比在k上有更大的能量水平, 表明這段期間渦更多的是經(jīng)向伸長。然而, 7月~12月渦信號在k上有較大的能量水平, 表明這段期間渦更多的是緯向伸長。從整個的量值上來看, 在9月和1月渦動能有較高水平, 這與圖3相符。

為觀察不同方向上渦動能所占總渦動能的比例, 本文利用Qiu[16]定義的兩類渦動能, 計算了經(jīng)向伸長的渦動能m(k>k)和緯向伸長的渦動能z(k<k), 如公式(2)和公式(3)所示,cpkm(cycles per kilometer),為網(wǎng)格間距,是每個方向上的格點數(shù)(本文中=32)。水平波數(shù)矢量,。

(3)

圖5繪制了緯向伸長的渦動能和經(jīng)向伸長的渦動能之差與總渦動能的比值, 表示的是緯向伸長或經(jīng)向伸長渦動能所占總渦動能的比例。如果, 則表示該階段渦動能更多的是緯向分布, 反之則更多的是經(jīng)向分布。經(jīng)向伸長的渦主要為剛產(chǎn)生的渦, 緯向伸長的渦主要是經(jīng)過了非線性渦渦相互作用調(diào)制后的渦[16]。

在接下來兩個部分, 本文將探討渦動能的這種特性分別與斜壓不穩(wěn)定和渦渦相互作用(不同尺度渦動之間能量的傳輸)之間的關(guān)系。

3 渦動發(fā)展與斜壓不穩(wěn)定

為了探討南黃海EKE存在雙峰結(jié)構(gòu)的機(jī)制, 本部分分析了南黃海斜壓結(jié)構(gòu)(與斜壓不穩(wěn)定過程相關(guān))與EKE的變化關(guān)系。圖6和圖7分別是利用ROMS模式氣候態(tài)月平均數(shù)據(jù)繪制的本文研究區(qū)域在35°N處全年溫度剖面圖和平均經(jīng)向流的速度剖面圖。二者在夏季與冬季存在著明顯不同的結(jié)構(gòu)特征: 夏季溫度存在顯著的季節(jié)性溫躍層, 深度在10 m左右, 7月、8月溫躍層最強(qiáng), 并且在溫躍層上方存在只有在夏季出現(xiàn)的北向流, 溫躍層下方受南黃海冷水團(tuán)環(huán)流左支的影響出現(xiàn)南向流; 而在冬季, 由于本文研究海區(qū)深度較淺并且受強(qiáng)勁的西北風(fēng)影響, 整個深度上溫度混合均勻, 并且?guī)缀跞繛槟舷蛄鳌４藭r,斜壓結(jié)構(gòu)很弱, 很難為渦動場直接提供能量。

為了分析斜壓不穩(wěn)定過程對EKE的影響, 本文選取2 m處和10 m處的流場數(shù)據(jù)計算了研究區(qū)域的速度垂向切變(圖8實線所示), 圖8中虛線為上文提到的本文研究區(qū)域的氣候態(tài)月平均EKE的時間序列。從圖中可以看出二者變化趨勢幾乎一致。與EKE類似, 垂向切變也存在雙譜峰結(jié)構(gòu), 7月存在較大的譜峰, 1月存在較小的譜峰, 4月和11月分別存在一個極小值, 但垂向切變的季節(jié)變化領(lǐng)先了EKE大約2個月。

由圖8可看出, 南黃海EKE的變化受斜壓不穩(wěn)定過程的影響較為顯著, 特別是在夏季: 當(dāng)3月~7月垂向切變增強(qiáng)時, 渦動能從5月~9月也逐漸增加, 而且這時經(jīng)向伸長的渦比較多(圖5); 當(dāng)8月~11月垂向切變減小時, EKE在2個月后也開始減小, 此時緯向伸長的渦比較多。所以, 我們認(rèn)為, 夏季存在的EKE譜峰主要是由于夏季較強(qiáng)的斜壓不穩(wěn)定過程直接導(dǎo)致的。

4 不同尺度過程的相互作用

4.1 譜演化方程

因為冬季南黃海幾乎沒有斜壓結(jié)構(gòu), 所以1月份存在的譜峰應(yīng)該與斜壓不穩(wěn)定過程無關(guān)。為了探討1月EKE譜峰的產(chǎn)生機(jī)制, 本部分我們通過計算動能通量譜[17-19]分析渦渦相互作用(不同尺度渦之間的能量傳輸過程)對EKE變化的影響。

計算方法如下, 考慮旋轉(zhuǎn)地球上笛卡爾坐標(biāo)系下水平動量方程, 并對其做離散傅里葉變換后整理得到譜密度的發(fā)展方程:

其中譜密度:

(5)

(7)

(8)

則有:

上橫線表示長期平均。

方程(7)和(9)為解釋通量譜提供了重要依據(jù)。根據(jù)方程(7), 正向(負(fù)向)通量表示動能從小(大)波數(shù)向大(小)波數(shù)傳輸, 也就是正向(反向)串級。另外根據(jù)方程(9), 通量譜的斜率為正(負(fù)), 則強(qiáng)迫項強(qiáng)(弱)于耗散項, 表示長期平均在該尺度上的一個能量源(匯)?？梢约僭O(shè)耗散項主要發(fā)生在高波數(shù)處, 在低波數(shù)時渦動能的譜運(yùn)輸主要和強(qiáng)迫項相平衡。

4.2 計算結(jié)果

由以上計算方法可以得出不同尺度的渦動之間的能量傳輸情況。圖9給出了本文計算的南黃海區(qū)域的動能通量譜。圖中可以看出幾乎每個月份都在空間尺度為128.9~154.2 km(即波數(shù)=6.485×10–3~ 7.758×10–3cpkm)范圍內(nèi)存在一個零點, 也就是說在此零點處相對應(yīng)波數(shù)的渦動處于一個相對穩(wěn)定的狀態(tài), 既不向大尺度傳輸能量也不向小尺度傳輸能量。除了10月和11月, 全年在空間尺度為122.1 km (= 8.19×10–3cpkm)處存在一個很大的負(fù)值, 即存在一個反向能量串級。也就是說, 在這些月份中波數(shù)小于零點波數(shù)時, 存在正的能量串級, 即波的能量由大尺度向小尺度傳輸; 波數(shù)大于零點波數(shù)時, 則存在負(fù)的能量串級即波的能量由小尺度向大尺度傳輸, 在零點波數(shù)處達(dá)到一個相對平衡。然而在10月和11月則相反, 小于零點波數(shù)時, 為負(fù)的能量串級, 即波的能量由較大尺度繼續(xù)向更大尺度傳輸; 而波數(shù)大于零點波數(shù)時, 則存在正的能量串級即波的能量由小尺度向更小的尺度傳輸, 我們認(rèn)為這是一個不穩(wěn)定的狀態(tài)。以12月、1月、2月為主的冬季反向能量串級幅度較大, 1月最大值達(dá)到–5.82×10–5m2/s3; 而在7月、8月、9月為主的夏季則與冬季相比反向能量串級幅度較小; 在春季則為過渡階段。

圖9和圖10都顯示, 冬季渦渦之間能量傳輸比較活躍, 而夏季渦渦之間能量傳輸很弱, 這說明冬季渦渦相互作用比較強(qiáng), 而夏季比較弱。結(jié)合圖8我們認(rèn)為, 與9月EKE的譜峰的產(chǎn)生機(jī)制不同, 1月EKE譜峰的存在主要是渦渦相互作用的結(jié)果。

5 總結(jié)與討論

本文采用區(qū)域海洋模式ROMS氣候態(tài)月平均的數(shù)據(jù), 研究了南黃海EKE的季節(jié)性變化并基于動能通量譜探討了相關(guān)變化機(jī)制, 得出以下主要結(jié)論:

1) 南黃海渦動能的季節(jié)性變化呈現(xiàn)雙峰值的特征, 最大峰值出現(xiàn)在夏末的9月, 較小峰值出現(xiàn)在1月, 春季和秋季的量值較低;

2) 南黃海EKE的變化受斜壓不穩(wěn)定過程的影響較為顯著, 特別是在夏季: 當(dāng)3月~7月垂向切變增強(qiáng)時, 渦動能從5月~9月也逐漸增加, 而且這時經(jīng)向伸長的渦比較多; 當(dāng)8月~11月垂向切變減小時, EKE在2個月后也開始減小, 此時緯向伸長的渦比較多。所以, 我們認(rèn)為, 夏季存在的EKE譜峰主要是由于夏季較強(qiáng)的斜壓不穩(wěn)定過程直接導(dǎo)致的;

3) 冬季EKE的變化主要與渦渦相互作用過程相關(guān), 此時不同尺度渦之間的能量傳輸比較活躍。

[1] 郭炳火. 黃海物理海洋學(xué)的主要特征[J]. 海洋科學(xué)進(jìn)展, 1993(3): 7-18. Guo Binghuo. Major features of the physical oceanography in the Yellow Sea[J]. Advances in Marine Science, 1993(3): 7-18.

[2] 徐丹亞, 趙保仁. 青島—石島近海反氣旋中尺度渦旋存在證據(jù)及數(shù)值模擬[J]. 海洋學(xué)報, 1999, 21(2): 18-26. Xu Danya, Zhao Baoren. A numerical study on the mesoscale anticyclonic eddy in Qingdao-Shidao offshore in winter[J]. Acta Oceanologica Sinica(Chinese Version), 1998, 21(2): 64-67.

[3] 湯毓祥, 鄒娥梅, 李興宰, 等. 南黃海環(huán)流的若干特征[J]. 海洋學(xué)報, 2000, 22(1): 1-16. Tang Yuxiang, Zou Emei, Lie Heung-Jae, et al. Some features of circulation in the southern Huanghai Sea[J]. Acta Oceanologica Sinica(Chinese Version), 2000, 22(1): 1-16.

[4] 劉興泉. 黃海冬季環(huán)流的數(shù)值模擬[J]. 海洋與湖沼, 1996, 5: 546-555. Liu Xingquan. Numerical simulation of the Yellow Sea wintertime circulation[J]. Oceanologia et Limnologia Sinica[J].1996, 5: 546-555.

[5] 管秉賢. 黃海冷水團(tuán)的水溫變化以及環(huán)流特征的初步分析[J]. 海洋與湖沼, 1963, 5(4): 225-284. Guan Bingxian. A preliminary study of the temperature variations and the characteristics of the circulation of the cold water mass of the Yellow Sea[J]. Oceanologia et Limnologia Sinica, 1963, 5(4): 225-284.

[6] 袁業(yè)立. 黃海冷水團(tuán)環(huán)流——I.冷水團(tuán)中心部分的熱結(jié)構(gòu)和環(huán)流特征[J].海洋與湖沼, 1979, 10(3): 187-199. Yuan Yehli. A preliminary study on the circulation related of the cold water-mass of the Yellow Sea Ⅰ. The thermal structure and the characteristics of the circulation in the central part of the[J]. Oceanologia et Limnologia Sinica, 1979, 10(3): 187-199.

[7] Hu D, Cui M, Yongxiang L I, et al. On the Yellow Sea cold water mass-related circulation[J]. The Yellow Sea, 1991, 4: 79-88.

[8] 繆經(jīng)榜. 陸架上熱生渦旋的線性模型[J]. 力學(xué)學(xué)報, 1989, 3: 257-265. Miao Jingbang. A linear model for the thermal eddy on the continental shelf[J]. Acta Mechanica Sinica, 1989, 3: 257-265.

[9] 李惠卿, 袁業(yè)立.黃海冷水團(tuán)熱結(jié)構(gòu)及其環(huán)流解析研究[J]. 海洋與湖沼, 1992, 1: 7-13. Li Huiqing, Yuan Yeli. Theoretical study on the thermal structure and circulation pattern related to cold watermass of Yellow Sea[J]. Oceanologia et Limnologia Sinica, 1992, 1: 7-13.

[10] 蘇紀(jì)蘭, 黃大吉. 黃海冷水團(tuán)的環(huán)流結(jié)構(gòu)[J]. 海洋與湖沼, 1995, 26(5): 1-7. Su Jilan, Huang Daji. On the current field associated with the Yellow Sea cold water mass[J]. Oceanologia et Limnologia Sinica, 1995, 26(5): 1-7.

[11] Xu Dongfeng, Yuan Yaochu, Liu Yuan. The baroclinic ciculation structure of Yellow Sea Cold Water Mass[J]. Science in China (Science D: Earth Scieces), 2003, 46(2): 117-126.

[12] 張啟龍, 楊玉玲, 程明華. 南黃海春季溫、鹽結(jié)構(gòu)特征分析[J]. 海洋科學(xué), 1994, 6: 50-55.Zhang Qilong, Yang Yuling, Cheng Minghua.An analysis of characteristics of thermohaline structures in the south Yellow Sea in spring[J].Marine Sciences, 1994, 6: 50-55.

[13] 劉志亮, 胡敦欣. 黃海夏季近岸海區(qū)環(huán)流的初步分析及其與風(fēng)速的關(guān)系[J]. 海洋學(xué)報, 2009, 31(2): 1-7. Liu Zhiliang, Hu Dunxin.Preliminary study on the Huanghai Sea coastal current and its relationship with local wind in summer[J]. Acta Oceanologica Sinica, 2009, 31(2): 1-7.

[14] 韋欽勝, 于志剛, 冉祥濱, 等. 黃海西部沿岸流系特征分析及其對物質(zhì)輸運(yùn)的影響[J]. 地球科學(xué)進(jìn)展, 2011, 26(2): 145-156. Wei Qingsheng, Yu Zhigang, Ran Xiangbin, et al. Characteristics of the western coastal current of the Yellow Sea and its impacts on material transportation[J]. Advances In Earth Science, 2011, 26(2): 145-156.

[15] Shchepetkin A F, Mcwilliams J C. The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model[J]. Ocean Modelling, 2005, 9(4): 347-404.

[16] Qiu B, Scott R B, Chen S, et al. Length scales of eddy generation and nonlinear evolution of the seasonally modulated south Pacific subtropical countercurrent[J]. Journal of Physical Oceanography, 2008, 38(7): 1515-1528.

[17] Frisch U. Turbulence. The legacy of A. N. Kolmogorov.[J]. Physics Today, 1995, 49(11): 521-523.

[18] Scott R B, Wang F. Direct evidence of an oceanic inverse kinetic energy cascade from satellite altimetry[J]. Jphys.oceanogr, 2005, 35(9): 1650-1666.

[19] Wang S, Liu Z, Pang C. Geographical distribution and anisotropy of the inverse kinetic energy cascade, and its role in the eddy equilibrium processes[J]. Journal of Geophysical Research Oceans, 2015, 120(7): 4891-4906.

(本文編輯: 李曉燕)

Seasonal variations of eddy kinetic energy and its spectral characteristics in the southern Yellow Sea

ZHANG Wan-lun1, 2, LIU Zhi-liang1, WANG Shi-hong1, 2

(1. Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. University of Chinese Academy of Science, Beijing 100049, China)

Seasonal variations of the eddy kinetic energy (EKE) and its related mechanisms were analyzed in the southern Yellow Sea based on a climatological model output. We found that EKE had a double-peak structure on seasonal scales in this region, with a stronger peak appearing in September and a weaker peak in January. Furthermore, the eddies were found to be more meridionally elongated in the first half of the year, whereas they were zonally elongated in the second half of the year. The spectral characteristics of EKE showed that the EKE peak in September was mainly due to baroclinic instabilities in summer while that in January was possibly due to eddy?eddy interactions.Thus, we discussed the seasonal variations of EKE in the southern Yellow Sea, which provides an understanding of the dynamic mechanisms in this area.

southern Yellow Sea; eddy kinetic energy; spectral kinetic energy flux

Mar.18, 2016

[National Natural Science Foundation of China, No.41276026; Strategic Priority Research Program of the Chinese Academy of Sciences, No.XDA11020301]

P733.1

A

1000-3096(2017)03-0130-08

10.11759/hykx20160318002

2016-03-18;

2017-02-25

國家自然科學(xué)基金(41276026); 中國科學(xué)院戰(zhàn)略先導(dǎo)專項(XDA11020301)

張菀倫(1990-), 女, 遼寧沈陽人, 碩士研究生, 從事中國近海環(huán)流研究, 電話: 13730920942, E-mail: wanlunzh@126.com; 劉志亮(1977-), 通信作者, 男, 山東濰坊人, 副研究員, 從事中國近海環(huán)流研究, 電話: 0532-82898840, E-mail: zhlliu@qdio.ac.cn