冀承振， 葉瑞杰， 董濟(jì)海， 張志偉， 田紀(jì)偉

(中國(guó)海洋大學(xué)物理海洋實(shí)驗(yàn)室，山東 青島 266100)

南海中尺度渦邊緣亞中尺度過程模式研究?

冀承振， 葉瑞杰， 董濟(jì)海??， 張志偉， 田紀(jì)偉

(中國(guó)海洋大學(xué)物理海洋實(shí)驗(yàn)室，山東 青島 266100)

基于高分辨率ROMS模式數(shù)據(jù)，對(duì)南海中尺度渦邊緣的亞中尺度過程進(jìn)行了研究。結(jié)果表明：在中尺度渦邊緣處存在明顯的溫度細(xì)絲結(jié)構(gòu)，渦度值要明顯大于1，O～W參數(shù)趨于正值，即在渦旋邊緣存在強(qiáng)速度剪切，加強(qiáng)了渦絲的形成。對(duì)亞中尺度過程動(dòng)能空間分布分析發(fā)現(xiàn)，中尺度渦邊緣區(qū)域動(dòng)能所占的比例大約為88.4%，比中心區(qū)域動(dòng)能所占的比例11.6%要大許多。從動(dòng)能譜的分析來看，在中尺度渦邊緣處動(dòng)能譜的斜率趨近于k-2，與表層準(zhǔn)地轉(zhuǎn)理論(SQG)相吻合，而在中心處呈現(xiàn)k-3，與準(zhǔn)地轉(zhuǎn)理論(QG)相一致。

亞中尺度過程；中尺度渦；動(dòng)能譜；ROMS模式

在現(xiàn)在海洋研究中，人們對(duì)中尺度過程(～100km)已有了深入的認(rèn)識(shí)，中尺度過程主要是以中尺度渦旋的形式存在，其水平尺度有幾百公里，攜帶了海洋中大部分能量，通過衛(wèi)星高度計(jì)等觀測(cè)手段，人們已經(jīng)對(duì)其動(dòng)力機(jī)制及其對(duì)水平熱量、動(dòng)量和物質(zhì)輸運(yùn)的貢獻(xiàn)有了很好的理解；同樣，人們對(duì)尺度小于1km(0.1～100m)的三維湍流過程，及其對(duì)海洋混合和能量耗散的作用也做了大量細(xì)致的研究[1]。然而，處于中間狀態(tài)的亞中尺度過程(時(shí)間尺度為1 d，空間尺度為1～10km)，由于觀測(cè)資料和數(shù)值模擬的限制，往往被人們忽視。實(shí)際上，亞中尺度過程往往以亞中尺度渦旋，渦絲，鋒面等形式廣泛存在上層海洋中[1-8]。在亞中尺度過程發(fā)生區(qū)域，相對(duì)渦度ζ與行星渦度f相當(dāng)，甚至有時(shí)候大于行星渦度f，人們基于羅斯貝數(shù)Ro=ζ/f～O(1)來定義亞中尺度過程[1]。

相對(duì)于中尺度渦和湍流混合觀測(cè)而言，亞中尺度過程這一關(guān)系到從中尺度過程到小尺度混合之間的過渡階段的觀測(cè)工作甚少，其主要原因是：(1)亞中尺度過程水平特征尺度小，目前高度計(jì)分辨率過低難以有效監(jiān)測(cè)到；(2)亞中尺度過程時(shí)間壽命短，難以有效捕捉到它的生成、傳播與消亡過程；(3)亞中尺度過程的生成源具有一定的“隨機(jī)性”，難以有效開展現(xiàn)場(chǎng)斷面觀測(cè)；(4)盡管亞中尺度過程垂向運(yùn)動(dòng)較強(qiáng)，但垂向速度仍在O(mm/s)量級(jí)內(nèi)，對(duì)測(cè)流儀器的精度要求非常高，目前常規(guī)海流觀測(cè)儀器的精度難以滿足。

亞中尺度過程的研究是目前物理海洋研究的前沿之一，其對(duì)海洋環(huán)流的整體運(yùn)行機(jī)制以及海洋能量的收支平衡有很重要的意義。本文將基于高分辨率(空間分辨率為1km)ROMS模式數(shù)據(jù)對(duì)中尺度渦邊緣亞中尺度過程進(jìn)行研究。

1 模式數(shù)據(jù)介紹

本文中的ROMS模式使用了單向離線嵌套的方法，對(duì)南海北部的亞中尺度過程進(jìn)行模擬[11]。大區(qū)海域范圍為15°N～25°N，110°E～126°E，水平分辨率為5 km，該區(qū)域模型命名為R5(見圖1，黑色實(shí)框)；嵌套的小區(qū)區(qū)域范圍為18°N～23°N，115°E～121°E，水平分辨率為1 km，該區(qū)域模型命名為R1(見圖2，藍(lán)色實(shí)框)。地形數(shù)據(jù)使用的是由NOAA提供的分辨率為2′的ETOPO2數(shù)據(jù)，兩個(gè)模型在垂向上分層一致，均分為30層，在表層和底部都進(jìn)行了適當(dāng)?shù)募用芴幚?。湍流封閉模型本文選擇了Generic Length Scale(GLS)模型[12]，而邊界條件則使用的是放射性邊界條件[13]。對(duì)于R5模型，其強(qiáng)迫場(chǎng)數(shù)據(jù)來自于歐洲中期天氣預(yù)報(bào)中心(ECMWF)月平均資料，水平分辨率為0.125°，而開邊界條件數(shù)據(jù)則是來自于HYCOM的月平均資料。為消除不同模式之間的誤差，同時(shí)獲得R5穩(wěn)定的初始場(chǎng)，先以HYCOM以及ECMWF的氣候態(tài)數(shù)據(jù)作為初始、邊界以及強(qiáng)迫條件，對(duì)R5模型模擬了30年。此時(shí)，根據(jù)動(dòng)能和勢(shì)能隨時(shí)間的變化曲線可以發(fā)現(xiàn)，模式結(jié)果已基本穩(wěn)定。然后，再以得到的穩(wěn)定結(jié)果為初始條件，以HYCOM以及ECMWF的2004—2012年的月平均資料為邊界以及強(qiáng)迫條件，對(duì)R5模型進(jìn)行了模擬。在得到R5模型結(jié)果之后，再以其作為初始及邊界

(顏色代表水深。The contour represents water depth.)

條件，以ECMWF的2008—2012年的月平均資料為強(qiáng)迫條件，對(duì)R1模型進(jìn)行了模擬，最終得到R1模型高分辨率的模擬結(jié)果。在分析的過程中，我們主要選取了模式穩(wěn)定后有中尺度渦存在的流場(chǎng)進(jìn)行分析，以來研究中尺度渦邊緣亞中尺度的特征。

2 數(shù)值模擬結(jié)果分析

2.1 溫鹽場(chǎng)、渦度場(chǎng)分布

本文選取了2009年12月第19天的模擬結(jié)果，從海表面高度(SSH，見圖2(a))的結(jié)果來看，在南海北部存在一個(gè)較強(qiáng)的暖渦，渦中心大致位于118°E，21°N，渦旋半徑大約為60km；從海表面溫度(SST，見圖2(b))可以看出明顯的中尺度渦結(jié)構(gòu)，并且在中尺度渦邊緣處有明顯的溫度細(xì)絲結(jié)構(gòu)。本文取圖2(a)中的黑色實(shí)線位置的斷面(“118斷面”)來分析穿越中尺度渦中心溫度和密度的分布。圖3給出了沿?cái)嗝娴臏囟群臀粍?shì)密度的分布，從斷面溫度的分布來看，最高溫度可以達(dá)26℃，出現(xiàn)在渦中心的位置，邊緣處的溫度要低于中心的溫度，由于中尺度渦中心存在下凹的緣故，存在明顯的等溫線下降，在中尺度渦邊緣時(shí)，混合層底的深度為50m左右，由于中心下凹，混合層底被下壓到100m左右的深度上，從位勢(shì)密度的水平分布上可以看到存在明顯的等密度線露頭，上混合層中保持了一定的水平浮力梯度。以上的分析表明，中尺度渦自身的特性支持了海洋上混合層中水平浮力梯度的存在，儲(chǔ)存了有效位能，正是由于上混合層存在的這種不穩(wěn)定，為亞中尺度過程的形成創(chuàng)造了潛在的條件[14]。

((a)中的黑色實(shí)線為選取的“118斷面”位置。The black line in (a) is “118 section” that we choose.)

圖3 沿118斷面的溫度分布(單位：℃) (a)及沿118斷面的位勢(shì)密度分布(b)(單位：kg/m3)

倘若能將我們的設(shè)想與現(xiàn)實(shí)結(jié)合，將會(huì)吸引來大批高校學(xué)生的向往，在互聯(lián)網(wǎng)時(shí)代，什么都很便利了，人們的生活質(zhì)量提高，對(duì)于服務(wù)的追求越來越高，只有以客戶為中心，才能為我們引來更多的人流量，進(jìn)而才能實(shí)現(xiàn)共贏，客戶得到了美的體驗(yàn)，商家得到了利益。然而在這過程中，應(yīng)當(dāng)對(duì)商家進(jìn)行嚴(yán)格的要求，切勿追捧蠅頭小利，而送走了客戶。

圖4 相對(duì)渦度的分布(單位: 10-4 s-1) (a)和Okubo-Weiss 參數(shù)計(jì)算結(jié)果及分布 (Unit:10-8s-1) (b)

2.2 動(dòng)能空間分布

圖5(a)表示水平流速的u的斷面分布，向東為正，由于是一個(gè)暖渦，表現(xiàn)為順時(shí)針旋轉(zhuǎn)，因此在21°N以南，u<0，在21°N以北，u>0。最大流速約0.8m·s-1，發(fā)生在中尺度渦邊緣的表層20m處。相對(duì)于水平流速u來說，流速v要小的多，最大流速僅0.15m·s-1。

圖6給出了歸一化的動(dòng)能空間分布，本文利用Butter濾波的方法將原始流速分解為中尺度波段(MP，尺度大于20km)和亞中尺度波段(SP，尺度小于20km)，對(duì)于中尺度波段來說，動(dòng)能的兩個(gè)峰值分別出現(xiàn)在20.4°N以及21.4°N附近，處于中尺度渦邊緣向中尺度渦內(nèi)部轉(zhuǎn)化的區(qū)域，因?yàn)檫@個(gè)區(qū)域內(nèi)流速出現(xiàn)最大值，因此動(dòng)能較大；相對(duì)于中尺度波段的動(dòng)能來說，亞中尺度波段動(dòng)能分布大大不同，動(dòng)能的峰值主要集中于中尺度渦的兩個(gè)邊緣，通過計(jì)算，兩個(gè)邊緣區(qū)域動(dòng)能所占的比例分別為26.5%和61.9%，比中心區(qū)域動(dòng)能所占的比例11.6%要大許多。亞中尺度波段動(dòng)能的分布表明亞中尺度過程在中尺度渦邊緣處較為活躍，而在中尺度渦內(nèi)部則較為“安靜”。

圖5 水平流速u的斷面分布(a)和水平流速v的斷面分布(b)(單位：m·s-1)

(紅線代表亞中尺度波段的動(dòng)能分布，灰線代表中尺度波段的動(dòng)能分布。縱坐標(biāo)經(jīng)過了各波段最大值歸一化處理。The red line and the gray line represent the kinetic energy distribution of SP and MP, respectively and the kinetic energy is normalized by the maximum value of SP and MP, respectively.)

圖6 從表層到200m深度動(dòng)能垂向積分分布圖

Fig.6 Depth mean distribution of kinetic energy from top layer to 200m

2.3 動(dòng)能譜分析

從上面的分析可以得出初步的結(jié)論，亞中尺度過程在中尺度渦邊緣時(shí)要比中尺度渦內(nèi)部更強(qiáng)一些。作為一種有效的手段，許多高分辨率的數(shù)值模式利用動(dòng)能譜來描述亞中尺度過程[16-20]。在海洋中，處于亞中尺度波段的強(qiáng)非地轉(zhuǎn)流可以從平衡態(tài)中獲取能量，然后傳輸給小尺度，準(zhǔn)地轉(zhuǎn)二維中尺度過程動(dòng)能譜特征是斜率為k-3，而高分辨率的能夠分辨亞中尺度過程的數(shù)值模式發(fā)現(xiàn)，譜的斜率變的平緩，變?yōu)閗-2，受這一啟發(fā)，本文這一小節(jié)從動(dòng)能譜的角度來分析。

本文在模式中選取了一條穿過中尺度渦中心的斷面(118斷面，見圖2(a))，通過對(duì)斷面流速的動(dòng)能譜分析，來探究中尺度渦中心和邊緣的差異。根據(jù)圖6(a)水平流速的分布，作者把穿越中尺度渦中心的118斷面分為三個(gè)部分，RR1，RR2和RR3，其中RR1的區(qū)域范圍為20°N～20.7°N，RR2區(qū)域范圍為20.7°N～21.4°N，RR3區(qū)域范圍為21.4°N～22.1°N，RR1和RR3代表中尺度渦的邊緣區(qū)域，RR2則代表中尺度渦的中心區(qū)域。

圖7給出了RR1，RR2和RR3 3個(gè)區(qū)域不同深度上的水平流速u的動(dòng)能譜分析。從圖7(a)中可以看出，在RR1區(qū)域50、100和150m 3個(gè)不同深度上，水平流速波數(shù)譜的斜率都接近于k-2，通過線性擬合，可以得到在尺度5～20km的亞中尺度波段，50m處的波數(shù)譜斜率為-2.17，100m處為-2.30，150m處為-2.34，說明亞中尺度過程到150m處還比較強(qiáng)；在RR3處也有相類似的結(jié)果(見圖7c)，在5～20km的亞中尺度波段，50，100和150m處的波數(shù)譜斜率分別為-2.15，-2.27，-2.24，斜率隨深度的變化更接近于k-2；但在RR2處，情況有所不同(見圖7(b))，50，100和150m的波數(shù)譜斜率分別為-2.81，-2.82，-2.91，比RR1和RR3區(qū)域的波數(shù)譜斜率要陡一些，斜率隨深度的變化更接近于k-3。

從以上的分析可以得出，中尺度渦邊緣動(dòng)能譜的斜率更接近于k-2，比較符合表層準(zhǔn)地轉(zhuǎn)理論(SQG理論)，而在中尺度渦內(nèi)部，動(dòng)能譜的斜率更接近于k-3，與二維準(zhǔn)地轉(zhuǎn)理論的預(yù)測(cè)相一致。Callies[19]及Rocha[21]等在文章中綜述了內(nèi)潮等也有可能使動(dòng)能譜的斜率變緩，接近于k-2。由于在本文的模式模擬中沒有加入潮的作用，因此，本文分析所得到的在中尺度渦邊緣處斜率平緩接近于k-2的結(jié)果，排除了內(nèi)潮因素的干擾，即中尺度渦邊緣動(dòng)能譜斜率呈現(xiàn)k-2變化是由于存在豐富的亞中尺度過程導(dǎo)致的。

(其中k-2和k-3通過兩條灰色直線標(biāo)出。k2andk-3are warkod by the two grey line.)

圖7 中尺度渦邊緣RR1處50，100和150m處水平流速的波數(shù)譜分析(a)；同子圖(a)，但為中尺度渦中心RR2處的結(jié)果(b)及同子圖(a),但為中尺度渦邊緣RR3的結(jié)果(c)

Fig.7 At the edge of the mesoscale eddy (RR1), wavenumber spectra for horizontal velocity at 50 (red line),100 (green line) and 150m (blue line)(a)，similar to (a) but at RR2 (b) and similar to (a) but at RR3 (c)

3 結(jié)語(yǔ)

本文利用高分辨率ROMS模式數(shù)據(jù)對(duì)南海中尺度渦邊緣亞中尺度過程進(jìn)行研究，溫度場(chǎng)的結(jié)構(gòu)顯示中尺度渦邊緣處有明顯的溫度細(xì)絲結(jié)構(gòu)；上混合層存在水平浮力梯度，導(dǎo)致混合層不穩(wěn)定，為亞中尺度過程的產(chǎn)生奠定了基礎(chǔ)；在中尺度渦邊緣處，正負(fù)渦度值交錯(cuò)，且其值要明顯大于1，O-W參數(shù)在渦邊緣處趨于正值，表明渦邊緣區(qū)域的水平流場(chǎng)切變很強(qiáng)，強(qiáng)速度剪切加強(qiáng)了渦絲的形成；動(dòng)能的空間分布表明中尺度渦邊緣區(qū)域動(dòng)能所占的比例大約為88.4%，比中心區(qū)域動(dòng)能所占的比例11.6%要大許多；從動(dòng)能譜的分析來看，在中尺度渦邊緣處動(dòng)能譜的斜率趨近于k-2，與表層準(zhǔn)地轉(zhuǎn)理論(SQG)相吻合，而在中心處呈現(xiàn)k-3，與準(zhǔn)地轉(zhuǎn)理論(QG)相一致，由于在模式模擬中沒有加入潮的作用，因此排除了內(nèi)潮因素的干擾，即中尺度渦邊緣動(dòng)能譜斜率呈現(xiàn)k-2變化是由于存在豐富的亞中尺度過程導(dǎo)致的。以上只是模式結(jié)果的初步探討，在下一步工作中，需要對(duì)亞中尺度過程產(chǎn)生及消亡機(jī)制進(jìn)一步的研究。

[1] Thomas L N, Tandon A, Mahadevan A. Submesoscale processes and dynamics[J]. Ocean modeling in an Eddying Regime, 2008: 17-38.

[2] Mahadevan A, Tandon A. An analysis of mechanisms for submesoscale vertical motion at ocean fronts[J]. Ocean Modelling, 2006, 14(3): 241-256.

[3] Capet X, McWilliams J C, Molemaker M J, et al. Mesoscale to submesoscale transition in the California Current System. Part I: Flow structure, eddy flux, and observational tests[J]. Journal of Physical Oceanography, 2008, 38(1): 29-43.

[4] Capet X, McWilliams J C, Molemaker M J, et al. Mesoscale to submesoscale transition in the California Current System. Part II: Frontal processes[J]. Journal of Physical Oceanography, 2008, 38(1): 44-64.

[5] Capet X, McWilliams J C, Molemaker M J, et al. Mesoscale to submesoscale transition in the California current system. Part III: Energy balance and flux[J]. Journal of Physical Oceanography, 2008, 38(10): 2256-2269.

[6] Gildor H, Fredj E, Steinbuck J, et al. Evidence for submesoscale barriers to horizontal mixing in the ocean from current measurements and aerial photographs[J]. Journal of Physical Oceanography, 2009, 39(8): 1975-1983.

[7] Klein P, Lapeyre G. The oceanic vertical pump induced by mesoscale and submesoscaleturbulence[J]. Annual Review of Marine Science, 2009, 1: 351-375.

[8] Zhong Y. Submesoscale dynamics and transport properties in the Gulf of Mexico [J].Georgia Insititute of Technology, 2013:121-128.

[9] Lévy M, Iovino D, Resplandy L, et al. Large-scale impacts of submesoscale dynamics on phytoplankton: Local and remote effects[J]. Ocean Modelling, 2012, 43: 77-93.

[10] McGillicuddy D J, Anderson L A, Bates N R, et al. Eddy/wind interactions stimulate extraordinary mid-ocean plankton blooms[J]. Science, 2007, 316(5827): 1021-1026.

[11] Mason E, Molemaker J, Shchepetkin A F, et al. Procedures for offline grid nesting in regional ocean models[J]. Ocean Modelling, 2010, 35(1): 1-15.

[12] Umlauf L, Burchard H. A generic length-scale equation for geophysical turbulence models[J]. Journal of Marine Research, 2003, 61(2): 235-265.

[13] Orlanski I. A simple boundary condition for unbounded hyperbolic flows[J]. Journal of Computational Physics, 1976, 21(3): 251-269.

[14] Boccaletti G, Ferrari R, Fox-Kemper B. Mixed layer instabilities and restratification[J]. Journal of Physical Oceanography, 2007, 37(9): 2228-2250.

[15] Isern-Fontanet J, Font J, García-Ladona E, et al. Spatial structure of anticyclonic eddies in the Algerian basin (Mediterranean Sea) analyzed using the Okubo-Weiss parameter[J]. Deep Sea Research Part II: Topical Studies in Oceanography, 2004, 51(25): 3009-3028.

[16] Capet X, Campos E J, Paiva A M. Submesoscale activity over the Argentinian shelf[J]. Geophysical Research Letters, 2008, 35(15)： 2015-2024.

[17] Wang D P, Flagg C N, Donohue K, et al. Wavenumber spectrum in the Gulf Stream from shipboard ADCP observations and comparison with altimetry measurements[J]. Journal of Physical Oceanography, 2010, 40(4): 840-844.

[18] Qiu B, Chen S, Klein P, et al. Seasonal mesoscale and submesoscale eddy variability along the North Pacific Subtropical Countercurrent[J]. Journal of Physical Oceanography, 2014, 44(12): 3079-3098.

[19] Callies J, Ferrari R. Interpreting energy and tracer spectra of upper-ocean turbulence in the submesoscale range (1-200km)[J]. Journal of Physical Oceanography, 2013, 43(11): 2456-2474.

[20] Callies J, Ferrari R, Klymak J M, et al. Seasonality in submesoscaleturbulence[J]. Nature Communications, 2015, 6.

[21] Rocha C B, Chereskin T K, Gille S T, et al. Mesoscale to submesoscale wavenumber spectra in Drake Passage[J]. Journal of Physical Oceanography, 2015, 46: 601-620.

責(zé)任編輯 龐 旻

The Simulation of Submesoscale Processat the Periphery of a Mesoscale Eddy in the South China Sea

JI Cheng-Zhen, YE Rui-Jie, DONG Ji-Hai, ZHANG Zhi-Wei, Tian Ji-Wei

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

The study investigates thesubmesoscale process at a mesoscale eddy in South China Sea (SCS) based on high-resolution ROMS model. The results show that, at the periphery of themesoscale eddy,there exists obvious temperature filaments, the vorticity value is significantly greater than 1, andO-Wparameter tends to be positive, all of these evidences meanthat the shear of velocity is much stronger and strengthen the formation of vortex filaments. From the view ofkinetic energy’s distribution, for submesoscale process, the proportion at the periphery of the eddy (RR1 and RR3) is 88.4%, which is larger than that at the center (RR2,only 11.6%).The results of kinetic energy spectra indicate that the slope of the spectra at the periphery of the eddy holdsk-2regime, which is more likely with surface quasi-geostrophic (SQG), while the slope exhibitsk-3law at the eddy’s center, consistent with quasi-geostrophic (QG) theory.

submesoscale process; mesoscale eddy; kinetic energy spectra; ROMS model

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(2014CB745003) 資助 Supported by the National Key Basic Research and Development Program of China(2014CB745003)

2016-03-13；

2016-05-06

冀承振(1990-)，男，碩士生。E-mail：jichengzhen@126.com

?? 通訊作者：E-mail：jihaihappy@163.com

P722.7

A

1672-5174(2017)01-001-06

10.16441/j.cnki.hdxb.20160068

冀承振， 葉瑞杰， 董濟(jì)海， 等. 南海中尺度渦邊緣亞中尺度過程模式研究[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 47(1): 1-6.

JI Cheng-Zhen, YE Rui-Jie, DONG Ji-Hai, et al. The simulation of submesoscale processat the periphery of a mesoscale eddy in the South China Sea[J]. Periodical of Ocean University of China, 2017, 47(1): 1-6.