北太平洋中層水模擬的改進(jìn)?

劉奇奇， 范植松??， 胡瑞金， 劉海龍， 楚合濤

(1.中國(guó)海洋大學(xué)海洋與大氣學(xué)院，山東 青島 266100; 2.中國(guó)科學(xué)院大氣物理研究所，北京 100029; 3.萊西氣象局，山東 萊西 266600)

?

北太平洋中層水模擬的改進(jìn)?

劉奇奇1， 范植松1??， 胡瑞金1， 劉海龍2， 楚合濤3

(1.中國(guó)海洋大學(xué)海洋與大氣學(xué)院，山東 青島 266100; 2.中國(guó)科學(xué)院大氣物理研究所，北京 100029; 3.萊西氣象局，山東 萊西 266600)

通過(guò)對(duì)北太平洋中層水模擬和分析的3種結(jié)果的比較：在稀疏網(wǎng)格的LICOM全球大洋環(huán)流模式中使用Canuto等的湍流混合方案(簡(jiǎn)稱(chēng)T方案)并添加范植松等提出的慣性?xún)?nèi)波破碎混合方案(簡(jiǎn)稱(chēng)F方案)，在LICOM模式中使用T方案，以及WOA09資料的分析結(jié)果，證明在添加F方案的情況下，北太平洋中層水的模擬有較明顯的改進(jìn)。慣性?xún)?nèi)波破碎混合可能是維持通風(fēng)過(guò)程的重要機(jī)制之一，慣性頻率風(fēng)場(chǎng)的強(qiáng)迫作用可能是海洋內(nèi)部混合的外部作用源。在WOA09資料的分析結(jié)果與數(shù)值模擬結(jié)果中均顯示出北太平洋中層水的明顯的季節(jié)變化，這意味著除表層水和模態(tài)水之外，北太平洋中層水在氣候變化中扮演了重要的角色。模擬結(jié)果分析表明，為獲得北太平洋中層水模擬的進(jìn)一步改進(jìn)，高分辨率的海洋模式和鋒面與羽流的參數(shù)化方案是必要的。

垂直混合;慣性?xún)?nèi)波;通風(fēng)過(guò)程;北太平洋中層水

引用格式：劉奇奇， 范植松， 胡瑞金, 等. 北太平洋中層水模擬的改進(jìn)[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016, 46(10): 1-9.

LIU Qi-Qi, FAN Zhi-Song, HU Rui-Jin, et al. An improvement of modeling the North Pacific intermediate water[J]. Periodical of Ocean University of China, 2016, 46(10): 1-9.

北太平洋中層水(North Pacific Intermediate Water，簡(jiǎn)稱(chēng)NPIW)的主要特征是具有整個(gè)水柱中鹽度的最小值，其中心密度大約為1026.8kg/m3[1-2]。它在太平洋經(jīng)向翻轉(zhuǎn)環(huán)流中扮演著重要角色，其來(lái)源是Okhotsk海[3]和Alaska灣流[4]。Hasumi等曾提到“盡管北太平洋中層水非常重要，但對(duì)其模擬的研究卻很少。粗網(wǎng)格模式對(duì)于模擬北太平洋中層水是不夠的，因?yàn)槠湫纬缮婕傲撕芏嗑?xì)的過(guò)程，比如Okhotsk海水和北太平洋水的交換，Okhotsk海源頭水的運(yùn)輸以及在混合水區(qū)域中的混合過(guò)程”[5]。Ishizaki和Ishikawa曾用過(guò)高分辨率模式來(lái)模擬北太平洋中層水[6]，并進(jìn)行改進(jìn)[7]。除了上述Hasumi等人指出的精細(xì)過(guò)程，作者認(rèn)為海洋內(nèi)部的垂直混合(跨等密度面混合)過(guò)程的參數(shù)化是模擬北太平洋經(jīng)向翻轉(zhuǎn)環(huán)流的關(guān)鍵，也是模擬北太平洋中層水的關(guān)鍵。眾所周知，海洋內(nèi)部的垂直混合過(guò)程和風(fēng)生南大洋Ekman上升流驅(qū)動(dòng)了全球的經(jīng)向翻轉(zhuǎn)環(huán)流[8-10]。然而，不同于北大西洋，在北太平洋，風(fēng)生南大洋Ekman上升流對(duì)經(jīng)向翻轉(zhuǎn)環(huán)流影響很小。

在F方案中考慮了慣性?xún)?nèi)波破碎產(chǎn)生的細(xì)結(jié)構(gòu)的擾動(dòng)引起的垂直效應(yīng)[13-14]。使用Canuto等的湍流混合方案(簡(jiǎn)稱(chēng)為T(mén)方案)[29]，并添加F方案，對(duì)全球海洋的溫度和鹽度的模擬獲得顯著的改進(jìn)，特別是在中高緯度區(qū)域的中層水和深層水的模擬。如文獻(xiàn)[13]中圖7～10所示，相對(duì)于初始場(chǎng)(WOA05資料)對(duì)太平洋的中層水和深層水模擬的鹽度和溫度的年平均相對(duì)誤差和均方差根誤差的統(tǒng)計(jì)結(jié)果顯示明顯的改進(jìn)。依據(jù)這些模擬結(jié)果，本文進(jìn)行模擬北太平洋中層水的比較研究。

1　研究基礎(chǔ)

1.1 LICOM模式

LICOM是中國(guó)科學(xué)院大氣物理研究所(IAP)大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國(guó)家重點(diǎn)實(shí)驗(yàn)室(LASG)開(kāi)發(fā)的全球海洋環(huán)流模式(OGCM)[30]。它的水平分辨率是1(°)×1(°)，垂直分為30層(300m深度之上為12層)。目前LICOM模式采用的混合方案是T方案，將在1.2節(jié)中介紹。依據(jù)該模式進(jìn)行了氣候平均態(tài)全球大洋數(shù)值試驗(yàn)。模擬時(shí)間為200年，初始無(wú)運(yùn)動(dòng)溫度與鹽度場(chǎng)為WOA05資料。在本次試驗(yàn)中，采用了2種不同混合方案：?jiǎn)为?dú)使用T方案和T方案添加F方案，F(xiàn)方案在1.3節(jié)中給出。模式在運(yùn)行到151a時(shí)達(dá)到基本穩(wěn)定。主要考慮鹽度、溫度和其它參數(shù)在151～200a的平均值。

1.2 湍流混合方案(T方案)

通過(guò)改進(jìn)Mellor和Yamada模式，Canuto等提出OGCM中比較新的湍流混合方案(T方案)[29，31-32]。T方案最顯著的改進(jìn)是將Mellor和Yamada模式的臨界Richardson數(shù)由Ricr=0.2改為Ricr～1。另外，在該新方案中考慮了鹽指的發(fā)生。

1.3 慣性?xún)?nèi)波破碎混合方案(F方案)

不同于湍流混合方案，F(xiàn)方案中考慮了由細(xì)尺度近慣性?xún)?nèi)波(可逆細(xì)結(jié)構(gòu))所產(chǎn)生斑塊的垂直徙動(dòng)。所使用的內(nèi)波能量譜E(β,θ,ω)為[13]：

(1)

(2)

(3)

(4)

(5)

A(β)=

(6)

(7)

細(xì)結(jié)構(gòu)的水平動(dòng)能譜Hsf(ω)為：

(8)

近慣性頻段的水平動(dòng)能譜Ek為：

(9)

細(xì)結(jié)構(gòu)的位移譜Dsf(ω)為：

(10)

近慣性頻段的勢(shì)能Ep為：

(11)

細(xì)結(jié)構(gòu)的剪切譜Φsf(β,ω)為：

(12)

由慣性?xún)?nèi)波破碎混合產(chǎn)生的動(dòng)量的垂直渦旋粘性系數(shù)Av和溫度/鹽度垂直渦旋擴(kuò)散系數(shù)Kv分別由以下兩個(gè)公式計(jì)算：

(13)

(14)

這里：

(15)

βc=0.1cpm=0.1×2π,(rad/m)。

(16)

h=4500m；N0=5.235988×10-3rad/s；b=1300m；π=3.14159；g=9.8m/s2；C=1500m/s。

(17)

因?yàn)樵谥鳒剀S層中，相對(duì)于等鹽度面，等位勢(shì)密度面更接近于等溫度面，所以沿等位勢(shì)密度面產(chǎn)生的破碎產(chǎn)物中溫度的分布比起鹽度來(lái)更為均勻。因此，對(duì)應(yīng)于斑塊垂直徙動(dòng)的溫度的垂直渦旋擴(kuò)散系數(shù)比鹽度的相應(yīng)值更小。本文假定在主溫躍層中的溫度的垂直渦旋擴(kuò)散系數(shù)Kt為：

Kt=Ks/5 。

(18)

這里的Ks是深度小于1500m的水層中鹽度垂直擴(kuò)散系數(shù)Kv，計(jì)算公式為(14)式。本模式的坐標(biāo)原點(diǎn)(z=0)設(shè)立在混合層底部。公式(18)是F方案獨(dú)有的，與湍流混合方案不同。通常在OGCM中只考慮湍流混合，并采用鹽度擴(kuò)散系數(shù)與溫度擴(kuò)散系數(shù)一致。例如，Nash和Moum依據(jù)海洋測(cè)量結(jié)果獲得如下的關(guān)系[33]：

dt=Kts/Ktt,0.6

(19)

這里dt是海洋中湍流混合鹽度擴(kuò)散系數(shù)Kts和溫度擴(kuò)散系數(shù)Ktt的比例。而在層化柵格湍流的室內(nèi)實(shí)驗(yàn)中，Rehmann和Koseff的結(jié)果為[34]：

0.67

(20)

但誤差范圍較大。

在赤道海洋(15°S～15°N)，由于強(qiáng)剪切的背景流和內(nèi)波的水平各向異性[35-36]，內(nèi)波能量譜E(β,θ,ω)的計(jì)算公式(1)～(7)不適用。另外，考慮到風(fēng)生慣性?xún)?nèi)波主要在中高緯度海域盛行[16，28]，本文建議在(15°S～15°N)之外的區(qū)域?qū)方案添加到T方案，在赤道區(qū)域僅用T方案。

2　數(shù)值模擬

2.1 北太平洋中層水的垂直剖面結(jié)構(gòu)

本文中選取北太平洋中層水團(tuán)中密度值為1026.8kg/m3的等密度面附近的水柱中鹽度最小值曲面作為該水團(tuán)的中心表示。為了揭示F方案對(duì)模擬北太平洋中層水及對(duì)其形成和傳輸機(jī)制理解的優(yōu)勢(shì)，本文中采用了2個(gè)方法，即北太平洋中層水的垂直剖面結(jié)構(gòu)的比較研究和其水平分布形態(tài)的比較研究，并分別以1月的模擬結(jié)果代表冬季，7月模擬結(jié)果代表夏季。

圖1為北太平洋中層水的垂直剖面結(jié)構(gòu)的比較，分別為在稀疏網(wǎng)格的LICOM模式中使用T方案添加F方案(表示為“LICOM+F”)，在LICOM模式中僅使用T方案(表示為“LICOM”)，以及WOA09資料的分析結(jié)果(表示為“WOA09”)。對(duì)應(yīng)剖面的經(jīng)度分別為150°E,180°E和210°E，使用1月份1500m以上深度的數(shù)據(jù)。圖2與圖1類(lèi)似，為沿20°N,30°N和40°N的1月份緯度剖面結(jié)構(gòu)的比較。圖1和2清楚地顯示使用F方案顯著改善了對(duì)北太平洋中層水的模擬，而且鹽度模擬的改進(jìn)更為明顯。該成果意味著使用F方案能更好的模擬通風(fēng)過(guò)程。如圖1～6所示，從Okhotsk海和Alaska灣流流入的低溫低鹽水向下俯沖混合并且在低緯度地區(qū)向上涌升，俯沖的最深深度能達(dá)到1000m左右。源頭水通過(guò)副極地環(huán)流圈、混合水區(qū)域的混合過(guò)程以及副熱帶環(huán)流圈在北太平洋散布開(kāi)來(lái)。在這個(gè)通風(fēng)過(guò)程中不能僅考慮湍流混合(T方案)，也應(yīng)考慮慣性?xún)?nèi)波破碎混合(F方案)。而且在F方案中鹽度垂直渦旋擴(kuò)散系數(shù)Ks的特殊性質(zhì)(公式(18))是適宜的。

稀疏網(wǎng)格LICOM模式的局限性在圖1中150°E剖面的高緯度地區(qū)清晰地顯示出來(lái)，而且在圖5中Okhotsk海及其鄰近西北太平洋區(qū)域也顯示出來(lái)。在這些區(qū)域，模擬結(jié)果相當(dāng)不同于觀(guān)測(cè)結(jié)果(WOA09)，甚至沒(méi)有合理的值。引起該缺陷的原因有2個(gè)：粗網(wǎng)格分辨率和海冰模式。在1月份(冬季)，Okhotsk海以及西北太平洋附近有大量的海冰存在，然而粗網(wǎng)格分辨率LICOM模式和海冰模式不能適當(dāng)?shù)啬M出該區(qū)域的海況，以及Okhotsk海和北太平洋的海水交換。另外，圖2中40°N剖面圖中的日本海寬度清楚地小于真實(shí)情況，這也將會(huì)影響日本海和西北太平洋的模擬結(jié)果。

(在稀疏網(wǎng)格的LICOM模式中使用T方案添加F方案(表示為“LICOM+F”)，在LICOM模式中僅使用T方案(表示為“LICOM”)，以及WOA09資料的分析結(jié)果(表示為“WOA09”)。對(duì)應(yīng)剖面的經(jīng)度分別為150°E,180°E和210°E，使用1月份1500m以上深度的數(shù)據(jù)。圖中彩色背景為鹽度(psu)；黑色實(shí)線(xiàn)為等溫線(xiàn)(℃)；藍(lán)色點(diǎn)線(xiàn)是等位密度線(xiàn)(kg/m3)；白色虛線(xiàn)是混合層底部深度(m)；紅色實(shí)線(xiàn)是選取的NPIW中心表示；橫坐標(biāo)是北半球緯度；縱坐標(biāo)為水深(m)。Depicted are the results simulated by performing the coarse-resolution LICOM with T-scheme adding F-scheme (denotation “LICOM+F”), those simulated with T-scheme alone (denotation “LICOM”), and the WOA09 data (denotation “WOA09”), against latitude along longitude 150°E, 180°E and 210°E respectively, in January above the depth of 1500m. In each panel, salinity is denoted with background color (psu)； Black solid curve represents isotherm (℃); Blue dotted curve represents potential isopycnal (kg/m3); White dotted curve represents the mixed layer base (m); Red solid curve is the centre representation of NPIW; Abscissa is latitude in the north hemisphere, and ordinate is the depth of water (m).)

圖1北太平洋中層水的經(jīng)度垂直剖面結(jié)構(gòu)的比較

Fig.1Comparison among simulation results of vertical profile structure along longitude of NPIW

(與圖1相似，但是分別為20°N,30°N和40°N的緯度剖面結(jié)構(gòu)的比較。As in Fig.1, but against longitude along latitude 20°N, 30°N and 40°N respectively.) 圖2　北太平洋中層水的緯度垂直剖面結(jié)構(gòu)的比較Fig.2　Comparison among simulation results of vertical profile structure along latitude of NPIW

(與圖1相似，但是為7月份。As in Fig.1, but in July.) 圖3　北太平洋中層水的經(jīng)度垂直剖面結(jié)構(gòu)的比較Fig.3　Comparison among simulation results of vertical profile structure along longitude of NPIW

(與圖2相似，但是為7月份。As in Fig.2, but in July.) 圖4　北太平洋中層水的緯度垂直剖面結(jié)構(gòu)的比較Fig.4　Comparison among simulation results of vertical profile structure along latitude of NPIW

(在稀疏網(wǎng)格的LICOM模式中使用T方案添加F方案(表示為“LICOM+F”)，在LICOM模式中僅使用T方案(表示為“LICOM”)，以及WOA09資料的分析結(jié)果(表示為“WOA09”)。對(duì)應(yīng)的鹽度(psu)、溫度(℃)和中心表示層深度(m)的水平分布圖所依據(jù)的數(shù)據(jù)均為1月份數(shù)據(jù)。Depicted are the simulation results of horizontal distribution pattern of salinity (psu), temperature (℃) and the depth (m) of the centre representation of NPIW by performing the coarse-resolution LICOM with T-scheme adding F-scheme (denotation “LICOM+F”), those simulated with T-scheme alone (denotation “LICOM”), and the WOA09 data (denotation “WOA09”) respectively in January.)

圖5北太平洋中層水的水平分布形態(tài)的比較

Fig.5Comparison among simulation results of horizontal distribution pattern of NPIW

(與圖5相似，但是為7月份。As in Fig.5, but in July.) 圖6　北太平洋中層水的水平分布形態(tài)的比較Fig.6　Comparison among simulation results of horizontal distribution pattern of NPIW

與圖1和2類(lèi)似，7月份的模擬結(jié)果示于圖3和4。北太平洋中層水的主要特性以及F方案的優(yōu)勢(shì)從1月維持到了7月，但從圖1～4中可觀(guān)察到NPIW明顯的季節(jié)變化。觀(guān)測(cè)資料(WOA09)和模擬結(jié)果都呈現(xiàn)出季節(jié)變化，這是因?yàn)閬?lái)自O(shè)khotsk海和Alaska灣流的源頭水具有強(qiáng)烈的季節(jié)變化特性。該季節(jié)變化暗示除了表層水和模態(tài)水，NPIW在氣候變化中扮演重要的角色[5]。NPIW對(duì)氣候變化的響應(yīng)及其在北太平洋長(zhǎng)期變化中的作用是有待研究的重要課題。

2.2 北太平洋中層水的水平分布形態(tài)

另一個(gè)研究方法是北太平洋中層水的水平分布形態(tài)的比較研究。圖5給出了1月(冬季)北太平洋中層水的水平分布形態(tài)的比較，分別為在稀疏網(wǎng)格的LICOM模式中使用T方案添加F方案(表示為“LICOM+F”)，在LICOM模式中僅使用T方案(表示為“LICOM”)，以及WOA09資料的分析結(jié)果(表示為“WOA09”)。對(duì)應(yīng)的鹽度(psu)、溫度(℃)和中心表示層深度(m)的水平分布圖所依據(jù)的數(shù)據(jù)均為1月份在中心表示曲面上的數(shù)據(jù)。圖6與5相似，是7月(夏季)的比較結(jié)果。與上述的NPIW垂直剖面結(jié)構(gòu)比較相似，使用F方案在NPIW的鹽度、溫度、中心表示層深度的水平分布模擬方面的改進(jìn)是明顯的。另外，在圖5和6中NPIW的季節(jié)變化也清楚地顯示出來(lái)。

值得注意的是，在圖5中Okhotsk海和西北太平洋附近的模擬結(jié)果中存在明顯的空白區(qū)域，此與圖1中150°E剖面圖中高緯度區(qū)域相對(duì)應(yīng)，這是由2.1節(jié)中提到的不足之處所引起。

3　結(jié)論與討論

通過(guò)以上的比較研究可得出結(jié)論如下：

(1)在稀疏網(wǎng)格的LICOM全球大洋環(huán)流模式中使用湍流混合方案(T方案)添加慣性?xún)?nèi)波破碎混合方案(F方案)，模擬北太平洋中層水的成果表明F方案是強(qiáng)有效的。

(2)為獲得北太平洋中層水模擬的進(jìn)一步改進(jìn),高分辨率的海洋模式和鋒面與羽流的參數(shù)化方案是必要的。

[1]Talley L D. An Okhotsk Sea water anomaly: implications for ventilation in the North Pacific [J]. Deep-Sea Res, 1991, 38: S171-S190.

[2]Talley L D. Distribution and formation of North Pacific intermediate water [J]. J Phys Oceanogr, 1993, 13: 517-537.

[3]Yasuda I. The origin of the North Pacific intermediate water [J]. J Geophys Res, 1997, 102: 893-909.

[4]You Y, Suginohara N, Fukasawa M, et al. Roles of the Okhotsk Sea and Gulf of Alaska in forming the North Pacific intermediate water [J]. J Geophys Res, 2000, 105: 3253-3280.

[5]Hasumi H, Tatebe H, Kawasaki T, et al. Progress of North Pacific modeling over the past decade [J]. Deep-Sea Res II, 2010, 57: 1188-1200. doi: 10. 1016/j. dsr2. 2009. 12. 008.

[6]Ishizaki H, Ishikawa I. Simulation of formation and spreading of salinity minimum associated with NPIW using a high-resolution model [J]. J Oceanogr, 2004, 60: 463-485.

[7]Ishikawa I, Ishizaki H. Importance of eddy representation for modeling of the intermediate salinity minimum in the North Pacific in comparison between eddy-resolving and eddy-permitting models [J]. J Oceanogr, 2009, 65: 407-426.

[8]Munk W, Wunsch C. Abyssal recipes II: Energetics of tidal and wind mixing [J]. Deep Sea Res Part I, 1998, 45(12): 1977-2010.

[9]Wunsch C, Ferrari R. Vertical mixing, energy and the general circulation of the oceans [J]. Annu Rev Fluid Mech, 2004, 36(1): 281-314.

[10]Kuhlbrodt T, Griesel A, Montoya M, et al. On the driving processes of the Atlantic meridional overturning circulation [J]. Rev Geophys, 2007, 45: RG2001, doi: 10. 1029/2004RG000166.

[11]Nasmyth P. Oceanic Turbulence [D]. Vancouver, Canada: Institute of Oceanography, University of British Columbia, 1970.

[12]Oakey N S. Determination of the rate of dissipation of turbulent energy from simultaneous temperature and velocity shear microstructure measurements [J]. J Phys Oceanogr, 1982, 12(3): 256-271.

[13]Fan Z S, Shang Z Q, Zhang S W, et al. A parameterization scheme of vertical mixing due to inertial internal wave breaking in the ocean general circulation model [J]. Acta Oceanol Sin, 2015, 34(1): 11-22. doi: 10. 1007/s13131-015-0591-1.

[14]尚真琦, 劉海龍, 范植松, 等. LICOM模式中不同混合方案的比較研究 [J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 45(10): 1-6. doi: 10. 16441/j. cnki. hdxb. 20140288.

Shang Z Q, Liu H L, Fan Z S, et al. A comparative study of different mixing schemes in framework of the LICOM (in Chinese)[J]. Periodical of Ocean University of China, 2015, 45(10): 1-6. doi: 10. 16441/j. cnki. hdxb. 20140288.

[15]Müller P, Olbers D, Willebrand J. The IWEX spectrum [J]. J Geophys Res, 1978, 83(C1): 479-500.

[16]Gregg M C, D’Asaro E A, Shay T J, et al. Observations of persistent mixing and near-inertial internal waves [J]. J Phys Oceanogr, 1986, 16(5): 856-885.

[17]Kunze E, Briscoe M G, Williams III A J. Interpreting shear and strain fine structure from a neutrally buoyant float [J]. J Geophys Res, 1990, 95(C10): 18111-18125.

[18]Kunze E, Williams III A J, Briscoe M G. Observations of shear and vertical stability from a neutrally buoyant float [J]. J Geophys Res, 1990, 95(C10): 18127-18142.

[19]Pinkel R, Sherman J, Smith J, et al. Strain: observations of the vertical gradient of isopycnal vertical displacement [J]. J Phys Oceanogr, 1991, 21: 527-540.

[20]Shcherbina A Y, Sundermeyer M A, Kunze E, et al. The LatMix summer campaign: submesoscale stirring in the upper ocean [J]. Bulletin of the American Meteorological Society, 2014. doi: 10. 1175/BAMS-D-14-00015. 1.

[21]范植松, 方欣華. 旋轉(zhuǎn)向量水平分量對(duì)大洋內(nèi)波方程的影響[J]. 海洋學(xué)報(bào), 1998, 20(3): 129-133.

Fan Z S, Fang X H. Effect of horizontal component of rotation vector on equations for oceanic internal waves[J]. Acta Oceanologica Sinica (in Chinese), 1998, 20(3)： 129-133.

[22]范植松, 方欣華. 考慮旋轉(zhuǎn)向量水平分量的大洋內(nèi)波方程的一個(gè)漸近解 [J]. 海洋學(xué)報(bào), 1998, 20(4): 1-8.

Fan Z S, Fang X H. An asymptotic solution of equations for oceanic internal waves under considering horizontal component of rotation vector[J]. Acta Oceanologica Sinica (in Chinese), 1998, 20(4)： 1-8.

[23]Fan Z S, Fang X H. A possible mechanism of ocean finestructure. Part I: Energy and coherence [J]. Journal of Ocean University of Qingdao, 1999, 29(2): 207-214.

[24]Fan Z S, Fang X H, Xu Q C. A possible mechanism of ocean finestructure. Part II: Shear and strain [J]. Journal of Ocean University of Qingdao, 1999, 29(3): 405-414.

[25]Fan Z S, Fang X H. A possible mechanism of ocean finestructure. Part III: Estimation of the kinetic energy dissipation in mixing [J]. Journal of Ocean University of Qingdao, 2000, 30(1): 7-14.

[26]Fan Z S. Is the small-scale turbulence an exclusive breaking product of oceanic internal waves [J]. Acta Oceanol Sin, 2011, 30(6): 1-11.

[27]Garrett C J R. What is the “near-inertial” band and why is it different from the rest of the internal wave spectrum? [J]. J Phys Oceanogr, 2001, 31: 962-971.

[28]Alford M H, Cronin M F, Klymak J M. Annual cycle and depth penetration of wind-generated near-inertial internal waves at Ocean Station Papa in the Northeast Pacific [J]. J Phys Oceanogr, 2012, 42(6): 889-909. doi: 10. 1175/JPO-D-11-092. 1.

[29]Canuto V M, Howard A, Cheng Y, et al. Ocean turbulence, III: new GISS vertical mixing scheme [J]. Ocean Modelling, 2010, 34(3-4): 70-91. doi: 10. 1016/j. ocemod. 2010. 04. 006.

[30]劉海龍, 俞永強(qiáng), 李薇， 等. LASG/LAP 氣候系統(tǒng)海洋模式(LICOM1. 0)參考手冊(cè)[M]. 北京: 科學(xué)出版社, 2004: 107.

Liu H, Yu Y, Li W, and Zhang X. 2004. Reference Manual for LASG/IAP Climate System Ocean Model (LICOM1. 0). Beijing: Science Press, China, 2004, 107 (in Chinese).

[31]Canuto V M, Howard A, Cheng Y, et al. Ocean turbulence. Part I: one point closure model, momentum and heat vertical diffusivities [J]. J Phys Oceanogr, 2001, 31(6): 1413-1426.

[32]Canuto V M, Howard A, Cheng Y, et al. Ocean turbulence. Part II: vertical diffusivities of momentum, heat, salt, mass and passive scalars[J]. J Phys Oceanogr, 2002, 32: 240-264.

[33]Nash J D, Moum J N. Microstructure estimates of turbulent salinity flux and the dissipation spectrum of salinity [J]. J Phys Oceanogr, 2002, 32: 2312-2333.

[34]Rehmann C R, Koseff J R. Mean potential energy change in stratified grid turbulence [J]. Dynamics of Atmospheres and Oceans, 2004, 37: 271-294.

[35]Eriksen C C. Some characteristics of internal gravity waves in the equatorial Pacific [J]. J Geophys Res, 1985, 90(C4): 7243-7255.

[36]Boyd T J, Luther D S, Knox R A, et al. High-frequency internal waves in the strongly sheared currents of the upper equatorial Pacific: Observations and a simple spectral model [J]. J Geophys Res, 1993, 98(C10): 18089-18107.

[37]Hibiya T, Nagasawa M, Niwa Y. Global mapping of diapycnal diffusivity in the deep ocean based on the results of expendable current profiler (XCP) surveys [J]. Geophys Res Lett, 2006, 33: L03611. doi: 10. 1029/2005GL025218.

責(zé)任編輯龐旻

An Improvement of Modeling the North Pacific Intermediate Water

LIU Qi-Qi1, FAN Zhi-Song1, HU Rui-Jin1, LIU Hai-Long2, CHU He-Tao3

(1.College of Oceanic and Atmospheric Sciences, Ocean University of China, Qingdao 266100, China; 2.State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China; 3.Laixi Meteorological Bureau, Laixi 266600, China)

The comparison among results simulated by performing the coarse-resolution LICOM model over the global ocean with the turbulence mixing scheme of Canuto et al. (T-scheme for short) adding the inertial internal wave breaking mixing scheme (F-scheme for short) put forward by Fan et al. recently, those simulated with T-scheme alone, and the WOA09 data, demonstrates that a notable improvement of modeling the North Pacific Intermediate Water (NPIW) has been attained by using F-scheme, which illustrates the robustness and effectiveness of F-scheme. The inertial internal wave breaking mixing may be one of important mechanisms maintaining the ventilation process, as well as inertial forcing of wind may be an external mechanical source of mixing in the ocean interior. The obvious seasonal variability of NPIW occurs in both the observations (WOA09) and the simulation results, which implies that NPIW plays an important role in climate changes besides the surface water and Mode water. Finally analysis of the simulation results indicates that to achieve further progress of modeling NPIW the high-resolution ocean model and the parameterization scheme of fronts and filaments are necessary.

vertical mixing; inertial internal wave; ventilation process; North Pacific intermediate water

國(guó)家自然科學(xué)基金面上項(xiàng)目(41275084)資助

2016-03-07；

2016-04-24

劉奇奇(1990-)，女，碩士生。E-mail: lqqi1990@163.com

??通訊作者：E-mail：fanzhs@hotmail.com

P722

A

1672-5174(2016)10-001-09

10.16441/j.cnki.hdxb.20160056

Supported by the National Natural Science Foundation of China(41275084)