劉波 李陽(yáng)春 徐永福 范廣洲

1成都信息工程學(xué)院大氣科學(xué)學(xué)院，成都610225

2中國(guó)科學(xué)院大氣物理研究所大氣邊界層物理和大氣化學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100029

3成都信息工程學(xué)院大氣科學(xué)學(xué)院高原大氣與環(huán)境四川省重點(diǎn)實(shí)驗(yàn)室，成都610225

1 引言

中尺度渦旋包含了海洋中大部分的動(dòng)能，對(duì)海洋環(huán)流和氣候的影響意義重大，雖然目前海洋環(huán)流模式的研發(fā)和改進(jìn)取得了很大進(jìn)展，但在分辨率不足以分辨中尺度渦旋的模式中，對(duì)中尺度渦旋的描述仍然存在很多的不確定性，所以對(duì)渦旋的參數(shù)化的研究依然非常重要。Gent and McWilliams（1990）和 Gent et al.（1995）提出的等密度傳輸參數(shù)化（GM90）方案是中尺度渦旋對(duì)大尺度溫鹽以及其他示蹤物傳輸?shù)拿枋觥Ｅc其他次網(wǎng)格混合參數(shù)化方案不同的是，GM90方案在平流速度內(nèi)加入了渦旋誘導(dǎo)速度，該方案使模擬的溫躍層變陡、深海變冷、經(jīng)向熱輸送增加（Danabasoglu et al., 1994）。

GM90方案中的兩個(gè)渦旋擴(kuò)散系數(shù)，等密度擴(kuò)散系數(shù) Ai和厚度擴(kuò)散系數(shù) Aith，是兩個(gè)重要的參數(shù)。早期模式對(duì)這兩參數(shù)通常采用常系數(shù)的方式處理，而從觀測(cè)資料可推導(dǎo)出渦動(dòng)能隨著海洋深度的增加是遞減的（Ledwell et al., 1998；Marshall et al.,2006），這些觀測(cè)事實(shí)也告訴我們渦旋擴(kuò)散系數(shù)不應(yīng)該是一個(gè)一成不變的常數(shù)，而是隨時(shí)間空間變化的變量。很多學(xué)者利用模式對(duì)這兩個(gè)參數(shù)做了大量的研究，Danabasoglu and McWilliams（1995）通過(guò)全球海洋環(huán)流模式對(duì)中尺度示蹤物傳輸參數(shù)化方案做了一系列敏感性試驗(yàn)，包括測(cè)試等密度和厚度擴(kuò)散系數(shù)隨深度指數(shù)遞減的方案，方案中系數(shù)在表層取2.0×103m2s-1，深海取0.5×103m2s-1，認(rèn)為它們不同的取值方式影響著海洋的溫度和流場(chǎng)輸送。另外，不同的海洋模式不同的系數(shù)取值在國(guó)內(nèi)也有很多研究（Xu et al., 2006；李陽(yáng)春等, 2007）。以上這些試驗(yàn)取值往往比較依賴經(jīng)驗(yàn)，Danabasoglu and Marshall（2007）利用 Ferreira et al.（2005）提出的在非絕熱層以下渦度擴(kuò)散系數(shù)隨浮力頻率變化的方案后，發(fā)現(xiàn)這種新方案使得渦旋擴(kuò)散系數(shù)比其他方案更符合客觀事實(shí)，同時(shí)新方案使得渦旋誘導(dǎo)速度的垂直結(jié)構(gòu)及其累積輸送量比使用常系數(shù)方案得到的結(jié)果更符合觀測(cè)事實(shí)。此外，他們的模擬結(jié)果（Danabasoglu and Marshall, 2007）顯示，采用新方案后，模式的模擬能力有很大的提高，比如對(duì)北太平洋上層海洋的垂直結(jié)構(gòu)模擬得更好了，渦旋速度的傳輸也有了提高，對(duì)上層海洋模擬偏暖的情況有了改善，等等。隨后，Danabasoglu et al.（2012）對(duì)該方案的系數(shù)進(jìn)行了進(jìn)一步的修改，并在NCAR（National Center for Atmospheric Research）通用氣候系統(tǒng)模式 CCSM4（Community Climate System Model 4.0）的海洋分量模式中應(yīng)用了該新方案。

海洋模式模擬能力的好壞不能僅憑模式對(duì)物理場(chǎng)模擬能力來(lái)判斷，在國(guó)際上海洋被動(dòng)示蹤物通常被用來(lái)考察海洋模式的模擬能力（Matsumoto et al., 2004），這是因?yàn)楹Ｑ髮?duì)被動(dòng)示蹤物輸送能力的描述影響到海洋模式是否適合用來(lái)研究碳循環(huán)以及氣候變化。目前用于海洋模式的被動(dòng)示蹤物主要有氟氯烴（CFCs）、14C、氚等。England et al.（1997）利用等密度面擴(kuò)散系數(shù)隨深度指數(shù)遞減的方式，即從表層的 5.0×103m2s-1變化到深海 1.0×103m2s-1，并以CFCs作為被動(dòng)示蹤物，研究了次網(wǎng)格混合機(jī)制對(duì)海洋模式通風(fēng)過(guò)程的作用，得出GM90方案相對(duì)Laplace混合等方案在對(duì) CFCs的模擬上有很大改進(jìn)的結(jié)論。Xu et al.（2006）研究了中尺度傳輸參數(shù)化對(duì)被動(dòng)示蹤物CFCs和14C在北太平洋海盆模式中的分布和儲(chǔ)藏的敏感性，指出增大等密度面擴(kuò)散系數(shù)能增大在西北太平洋的 CFCs儲(chǔ)藏量。李陽(yáng)春等（2007）使用全球海洋模式L30T63，研究了不同次網(wǎng)格混合系數(shù)對(duì)CFCs模擬的影響，指出擴(kuò)大等密度擴(kuò)散系數(shù)能改善CFC-11（一氟三氯甲烷）模擬穿透深度不足的問(wèn)題。到目前為止，新的渦旋擴(kuò)散系數(shù)隨浮力頻率變化的方案（Danabasoglu and Marshall, 2007）并沒(méi)有利用海洋被動(dòng)示蹤物來(lái)檢驗(yàn)過(guò)它的性能。

本文的目的是利用中國(guó)科學(xué)院大氣物理研究所（IAP）大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國(guó)家重點(diǎn)實(shí)驗(yàn)室（LASG）發(fā)展的海洋環(huán)流模式LICOM（張學(xué)洪等，2003；劉海龍等, 2004；Liu et al., 2004,2012），來(lái)探討渦旋擴(kuò)散系數(shù)隨浮力頻率變化（浮力方案）的次網(wǎng)格尺度混合參數(shù)化方案對(duì) LICOM模式的影響。通過(guò)模擬結(jié)果與示蹤物（包括主動(dòng)以及被動(dòng)）觀測(cè)資料比較，考察新的參數(shù)化浮力方案對(duì)海洋物理場(chǎng)模擬的影響，并檢驗(yàn)其能否對(duì)LICOM的模擬結(jié)果有所改進(jìn)。

2 模式介紹與試驗(yàn)設(shè)計(jì)

本文選用的海洋模式版本為L(zhǎng)ICOM1.0（Liu et al., 2004），其經(jīng)緯度范圍分別為 0°～360°和 78°S～90°N。由于碳循環(huán)需要較長(zhǎng)時(shí)間的積分，模式選為粗網(wǎng)格2°×2°的形式。模式垂直方向采用η坐標(biāo)，共30層，上面15層每隔10 m均勻變化，模式最深達(dá)5600 m，示蹤物層最深深度為5243 m。模式采用自由海表面，水平網(wǎng)格采用Arakawa B網(wǎng)格。本文動(dòng)量方程仍采用Laplace常系數(shù)粘性方案，水平粘性系數(shù)為 8.0×104m2s-1，垂直粘性系數(shù)為1.0×10-3m2s-1，并在 30°S～30°N 范圍內(nèi)使用依賴Richardson數(shù)形式的湍流粘性項(xiàng)。模式對(duì)中尺度次網(wǎng)格的描述采用GM90方案，在海表以及等密度面相對(duì)水平面坡度（isopycnal slope）較大時(shí)使用LDD97方案對(duì)其削弱（Large et al., 1997）。模式更多詳情見(jiàn)LICOM手冊(cè)（劉海龍等, 2004）。

模式使用的初始場(chǎng)資料和用于海表恢復(fù)的溫鹽資料均來(lái)自美國(guó)國(guó)家海洋資料中心（NODC）提供的 WOA09資料，模式強(qiáng)迫場(chǎng)來(lái)自 MPI-OMIP（Max-Planck-Institute—Ocean Model Intercomparison Project）的氣候態(tài)月平均資料，包括風(fēng)應(yīng)力、輻射通量等。模式物理場(chǎng)在以上強(qiáng)迫場(chǎng)的強(qiáng)迫下從靜止態(tài)開(kāi)始積分達(dá)到準(zhǔn)穩(wěn)態(tài)后，將CFC-11作為被動(dòng)示蹤物放入模式中同步積分，從1944年至2007年，共64年。大氣中CFC-11濃度是由Atmospheric Oxygen Research Group（AORG）提供的 1944～2007年的資料（https://bluemoon.ucsd.edu/pub/cfchist/[2011-6-27]）。CFC-11海氣通量公式為

其中，F(xiàn)為海氣交換通量；k為傳輸速度；Ceqm代表大氣CFC-11在海表的平衡濃度，Ceqm= αCa，Ca為大氣濃度，α為 CFC-11在海水中的溶解度，是海表溫鹽的函數(shù)（Warner and Weiss, 1985）；Cocn代表模擬的海表CFC-11濃度；ε為模式網(wǎng)格內(nèi)海冰所占的面積比例，資料來(lái)源于Hadley中心氣候態(tài)月平均結(jié)果（Rayner et al., 2003）。k的計(jì)算，本工作采用趙琦等（2012）認(rèn)為最適合 LICOM 模式模擬CFC-11的Wanninkhof（1992）傳輸速度公式：

其中，u為海表10 m高度處的風(fēng)速；Sc代表Schimit數(shù)。

本文針對(duì)GM90方案設(shè)計(jì)了兩個(gè)不同的試驗(yàn)：對(duì)照試驗(yàn)和浮力試驗(yàn)。對(duì)照試驗(yàn)采用常系數(shù)方案，等密度擴(kuò)散系數(shù) Ai和厚度擴(kuò)散系數(shù) Aith均取值為1.0×103m2s-1。浮力方案是參考 Ferreira et al.（2005）和Danabasoglu and Marshall（2007）提出的方案，使渦旋擴(kuò)散系數(shù)Aρ隨浮力頻率變化：

式中，Aρ代表等密度擴(kuò)散系數(shù)Ai和厚度擴(kuò)散系數(shù)Aith；Aref取3.0×103m2s-1；γ為浮力系數(shù)，浮力頻率（N）的平方為

其中，g為重力加速度，取g=9.806 m s-2；ρ為海水位密度；z取向上為正。當(dāng)N2＞0時(shí)，海洋是層結(jié)穩(wěn)定的，當(dāng)N2＜0時(shí)，海洋層結(jié)不穩(wěn)定。為參考層的N2，取非絕熱層底層DLD（diabatic layer depth）為參考層，非絕熱層一般是指與外界存在熱量傳輸?shù)乃畬?，Danabasoglu et al.（2008）形式上將非絕熱層定義為海洋邊界層（boundary layer depth）和海洋過(guò)渡層（transition layer thickness）之和，并建議使用海洋混合層深度MLD（mixed layer depth）能較好的代替非絕熱層。因此，本文采用混合層底層并保證N2≥0的情況下作為浮力頻率的參考層，若該層 N2＜0，則取混合層以下最先穩(wěn)定的那層作為參考層。在模擬試驗(yàn)中，淺層區(qū)域即參考層以上的浮力系數(shù)γ取為0.33，表示非絕熱層中渦旋擴(kuò)散系數(shù)沒(méi)有垂直變化且保證其值約為1×103m2s-1。在參考層與海底之間，我們確保 γmin≤γ≤1.0，γmin是指定的系數(shù)下限，按照 Danabasoglu and Marshall（2007）的建議，取 γmin= 0 .1，并規(guī)定參考層= 0的地方，γ= 0.1。本文每一模式時(shí)間步長(zhǎng)均計(jì)算γ，也就是說(shuō)渦旋擴(kuò)散系數(shù)是隨時(shí)間空間變化的變量。

除了對(duì)GM90方案中的等密度擴(kuò)散系數(shù)和厚度擴(kuò)散系數(shù)處理之外，本文試驗(yàn)中其他參數(shù)均保持一致。由于參數(shù)取值不同，模式物理場(chǎng)達(dá)到穩(wěn)態(tài)的時(shí)間也有所差異，對(duì)照試驗(yàn)和浮力試驗(yàn)分別運(yùn)行了3200年和4100年后達(dá)到準(zhǔn)穩(wěn)態(tài)（模式最后100年結(jié)果每層溫度變化均小于 0.009°C，平均鹽度變化均小于0.001 psu）。本文用兩個(gè)試驗(yàn)各自穩(wěn)態(tài)后的十年平均結(jié)果作比較分析（下文若沒(méi)單獨(dú)指出，年平均結(jié)果均為最后十年平均結(jié)果）。

3 結(jié)果與討論

3.1 浮力參數(shù)

從圖 1d可以看出浮力系數(shù)γ緯向平均呈現(xiàn)隨深度增加而遞減的垂直分布狀態(tài)，熱帶海域垂直變化最快，南北兩副熱帶海區(qū)擴(kuò)散系數(shù)變化較熱帶海區(qū)較慢，在熱帶海區(qū)400 m深度內(nèi)系數(shù)很快從3.0變成0.3，而副熱帶海區(qū)0.3等值線一直延伸到1800 m 左右，由公式（3）可知，浮力系數(shù)γ的分布特征即渦旋擴(kuò)散系數(shù)Aρ的分布特征。上述這些參數(shù)的特征與Danabasoglu and Marshall（2007）得到的結(jié)果基本一致。

對(duì)于被動(dòng)示蹤物的吸收，混合層往往起了很重要的作用，但在GM90方案中，因?yàn)榛旌蠈訛榉墙^熱層，不滿足使用條件，因而模式采用了 LDD97方案（Large et al., 1997）對(duì)GM90方案所產(chǎn)生的擴(kuò)散過(guò)程進(jìn)行了削弱。此外，本文兩個(gè)試驗(yàn)在混合層所采用的系數(shù)是一致的?？紤]到大洋內(nèi)部溫鹽場(chǎng)的改變?cè)谝欢ǔ潭壬峡梢杂绊懟旌蠈拥椎暮Ｋ芏冗M(jìn)而可能導(dǎo)致兩個(gè)試驗(yàn)混合層深度產(chǎn)生的差異，本文首先對(duì)這一影響海氣物質(zhì)能量交換的水層進(jìn)行對(duì)比分析。圖2代表觀測(cè)、對(duì)照試驗(yàn)1月份、7月平均混合層深度以及浮力試驗(yàn)與對(duì)照試驗(yàn)相應(yīng)月份混合層深度差值圖，其中觀測(cè)結(jié)果是使用WOA09溫鹽資料得出（本文使用與海表位密度差Δσθ=0.125 kg m-3判斷混合層的最底層深度）。對(duì)照試驗(yàn)1月份北半球MLD模擬形態(tài)基本與觀測(cè)分布一致（圖 2b），主要差異出現(xiàn)在南大洋 40°S～60°S之間混合帶沒(méi)有模擬出來(lái)，在北大西洋MLD模擬偏深，最大深度達(dá)2000 m。浮力試驗(yàn)一月份的MLD模擬與對(duì)照試驗(yàn)相似（圖 2e），在北半球副熱帶海區(qū)有少部分海區(qū)有差異，相對(duì)差異在 5%以內(nèi)，同時(shí)對(duì)對(duì)照試驗(yàn)在南大洋MLD模擬偏弱的情況稍有改善。對(duì)照試驗(yàn)七月份的MLD分布整體與觀測(cè)相似（圖2f），差異同樣主要在南大洋。在40°S～50°S南太平洋區(qū)域，對(duì)照試驗(yàn)?zāi)M的MLD偏弱，而對(duì)60°S以南靠近南極大陸的海區(qū)，對(duì)照試驗(yàn)MLD模擬偏強(qiáng)。與對(duì)照試驗(yàn)相比，浮力試驗(yàn)在7月份MLD模擬上與之的主要差異也出現(xiàn)在南大洋（圖 2f）。浮力試驗(yàn)對(duì)南太平洋40°S以南130°E～130°W海區(qū)MLD模擬增強(qiáng)，其余大部分海區(qū)模擬減弱?？傮w來(lái)說(shuō)，兩個(gè)試驗(yàn)MLD模擬在南大洋的差異較大，這是由于浮力試驗(yàn)在混合層以下應(yīng)用渦旋擴(kuò)散系數(shù)隨浮力頻率變化的方案，使得混合層下溫鹽發(fā)生改變，溫鹽的改變會(huì)影響海洋層結(jié)穩(wěn)定結(jié)構(gòu)，從而使得模式垂直混合強(qiáng)度發(fā)生變化。比如在南太平洋高緯西側(cè)海區(qū)，渦旋擴(kuò)散系數(shù)垂直變化的浮力試驗(yàn)對(duì)對(duì)照試驗(yàn)?zāi)蠘O中層水模擬過(guò)強(qiáng)的情況有所改善，使得在該區(qū)域中層水海溫相對(duì)對(duì)照試驗(yàn)增大，從而海洋層結(jié)穩(wěn)定度降低，導(dǎo)致該區(qū)域垂直混合加強(qiáng)，最大加深達(dá)75 m以上。

圖1 （a）觀測(cè)及（b）浮力試驗(yàn)得到的參考層浮力頻率平方的分布（單位：10-4s-2）；（c）上層海洋（0～876 m）平均渦旋擴(kuò)散系數(shù)分布（單位：m2 s-1）；（d）年平均γ的緯向平均分布Fig. 1 Global distribution of buoyancy frequency ( N r2e f)in the reference layer from (a)observations and (b)BUOYANCY (the test in which the eddy diffusion coefficients varied in the vertical direction along with the N2 below the diabatic layer depth); (c)eddy diffusion coefficient averaged in the upper-layer ocean (0–876 m); (d)annual mean zonally averagedγ. The units are 10-4 s-2 in (a)and (b), and m2 s-1 in (c)

3.2 環(huán)流特征

圖3是各試驗(yàn)渦旋誘導(dǎo)經(jīng)向流函數(shù)和歐拉流速加上渦旋誘導(dǎo)流速經(jīng)向流函數(shù)（文中簡(jiǎn)稱全經(jīng)向流函數(shù)）。在所有的試驗(yàn)中，可看出由渦旋誘導(dǎo)速度得到的經(jīng)向流函數(shù)表現(xiàn)為南半球 30°S～65°S一個(gè)從海表到底部的反時(shí)針?lè)D(zhuǎn)和北半球副熱帶附近一個(gè)較弱且淺的順時(shí)針?lè)D(zhuǎn)。從圖3a，b對(duì)比可以發(fā)現(xiàn)，渦旋擴(kuò)散系數(shù)垂直變化的浮力試驗(yàn)導(dǎo)致南半球200 m深度30°S～60°S附近流函數(shù)強(qiáng)度變?nèi)酰?5°S附近浮力試驗(yàn)輸送強(qiáng)度出現(xiàn)低值約為14 Sv（1 Sv=106m3s-1），而對(duì)照試驗(yàn)為 18 Sv。在南半球600～2000 m深度處，浮力試驗(yàn)渦旋誘導(dǎo)流函數(shù)很明顯較對(duì)照試驗(yàn)弱，浮力試驗(yàn)在200 m以下的大值（絕對(duì)值）中心深度明顯淺于對(duì)照試驗(yàn)，且最大值僅為16 Sv，小于對(duì)照試驗(yàn)的18 Sv。總的來(lái)說(shuō)，通過(guò)與對(duì)照試驗(yàn)對(duì)比可以看出，渦旋擴(kuò)散系數(shù)垂直變化的浮力試驗(yàn)使得渦旋誘導(dǎo)經(jīng)向流函數(shù)南半球的反時(shí)針?lè)D(zhuǎn)上層以及深海減弱。Danabasoglu and Marshall（2007）指出厚度擴(kuò)散系數(shù)垂直變化對(duì)渦旋誘導(dǎo)速度有很大的影響，且上層海洋300 m以上厚度擴(kuò)散系數(shù)垂直梯度很大，因而在使用了存在空間變化的厚度擴(kuò)散系數(shù)后模式在溫躍層的模擬結(jié)果有了很大的提升，此外等密度擴(kuò)散系數(shù)直接和間接的影響等密度斜坡（isopycnal slopes）的分布，從而對(duì)渦旋誘導(dǎo)經(jīng)向流函數(shù)的模擬產(chǎn)生影響。與本文上述模擬結(jié)果稍有不同，Danabasoglu and Marshall（2007）得到渦旋擴(kuò)散系數(shù)隨浮力頻率變化使得渦旋誘導(dǎo)經(jīng)向流函數(shù)南半球的反時(shí)針?lè)D(zhuǎn)上層加強(qiáng)。這也說(shuō)明參數(shù)化方案對(duì)不同模式的影響存在一定差異。

圖2 （a，b）1月和（c，d）7月平均混合層深度，（a，c）為觀測(cè)，（c，d）為對(duì)照試驗(yàn)；浮力試驗(yàn)與對(duì)照試驗(yàn)的（e）1月份、（f）7月份差值。單位：mFig. 2 The (a, b)January and (c, d)July mean mixed layer depth from (a, c)observations and (b, d)CONTROL (the test in which the eddy diffusion coefficients were spatially constant); the difference between BUOYANCY and CONTROL in (e)January and (f)July. Units: m

圖3c、d分別為對(duì)照試驗(yàn)和浮力試驗(yàn)的全經(jīng)向流函數(shù)分布圖，從圖 3c可以看出，全經(jīng)向流函數(shù)的主要特征為：在赤道表層兩側(cè)存在由Ekman輻散造成的兩個(gè)淺的翻轉(zhuǎn)，均由赤道流向兩極；在南半球主要存在一個(gè)較大迪肯流環(huán)（Deacon Cell），從海表伸向2000 m深度；在北半球30°～60°N之間的北大西洋深水（NADW），最大輸送強(qiáng)度為12 Sv；象征南極底層水的反時(shí)針?lè)D(zhuǎn)環(huán)流等等。上述全經(jīng)向流函數(shù)特征，浮力試驗(yàn)與對(duì)照試驗(yàn)基本一致，兩個(gè)試驗(yàn)的主要差異在南半球Deacon Cell上。由于浮力試驗(yàn)渦旋誘導(dǎo)經(jīng)向流函數(shù)對(duì)南半球的反時(shí)針?lè)D(zhuǎn)上層以及深海減弱，使得浮力試驗(yàn)全經(jīng)向流函數(shù)的 Deacon Cell環(huán)流在上層以及深海加強(qiáng)，在400～800 m深度有所減弱。比如，在南半球45°S上層海區(qū)，浮力試驗(yàn)在200 m深度中心最大值為22 Sv，該處對(duì)照試驗(yàn)僅為14 Sv；在250～750 m區(qū)域，對(duì)照試驗(yàn)出現(xiàn)Deacon Cell環(huán)流大值中心，最大約為16 Sv，該處比浮力試驗(yàn)大2 Sv；在深海，以等值線4 Sv為例，浮力試驗(yàn)延伸最大深度達(dá)3000 m，對(duì)照試驗(yàn)約為1700 m左右。由于Deacon Cell深水環(huán)流的增強(qiáng)，導(dǎo)致浮力試驗(yàn)?zāi)蠘O底層水模擬變?nèi)?，從圖 3d可以明顯看出浮力試驗(yàn)?zāi)蠘O底層水環(huán)流輸送比對(duì)照試驗(yàn)弱得多。Danabasoglu and Marshall（2007）也得到類似的結(jié)論，并指出，南極底層水的減弱會(huì)加強(qiáng)北大西洋深水向下的輸送。這在我們的浮力試驗(yàn)中也有所體現(xiàn)，浮力試驗(yàn)中象征北大西洋深水的經(jīng)向流函數(shù)大值中心12 Sv等值線的范圍比對(duì)照試驗(yàn)的有一定的增大，環(huán)流向下向南擴(kuò)散增強(qiáng)。

圖3 （a）對(duì)照試驗(yàn)和（b）浮力試驗(yàn)?zāi)昶骄鶞u旋誘導(dǎo)經(jīng)向流函數(shù)，（c）對(duì)照試驗(yàn)和（d）浮力試驗(yàn)?zāi)昶骄?jīng)向流函數(shù)（渦旋誘導(dǎo)流速＋歐拉流速）。單位: Sv，等值線間隔均為2Fig. 3 Annual mean eddy-induced meridional stream function from (a)CONTROL and (b)BUOYANCY, and the full meridional stream function(eddy-induced plus Euler’s velocity)from (c)CONTROL and (d)BUOYANCY. Units: Sv; contour interval: 2

圖4是各試驗(yàn)正壓流函數(shù)分布圖，正壓流函數(shù)能描述海洋整層輸送的無(wú)輻散部分。從圖4a，b可看出兩個(gè)試驗(yàn)均能很好地模擬出大洋環(huán)流的西向強(qiáng)化特點(diǎn)，比如北太平洋副熱帶西邊界流（黑潮）強(qiáng)度達(dá) 40 Sv，北大西洋也存在一條強(qiáng)度稍弱的西邊界流，同樣，在南半球也存在較強(qiáng)的西邊界流，其中印度洋副熱帶西邊界流強(qiáng)度可達(dá)60 Sv。兩個(gè)試驗(yàn)的主要差異在南大洋。由于南極繞流的存在，南大洋海區(qū)的等值線非常密集，流函數(shù)梯度很大。Whitworth and Peterson（1985）通過(guò)觀測(cè)指出在南大洋輸送強(qiáng)度約為135±13 Sv，對(duì)照試驗(yàn)在南大洋高緯輸送強(qiáng)度模擬偏弱，最大僅為100 Sv，浮力試驗(yàn)整體模擬增強(qiáng)，南大洋大部分海區(qū)增大約為對(duì)照試驗(yàn)的20%～30%，在南大西洋部分最大值約為130 Sv，比對(duì)照試驗(yàn)增大約30%，更加接近觀測(cè)事實(shí)。渦旋擴(kuò)散系數(shù)垂直變化的浮力試驗(yàn)在這方面相對(duì)對(duì)照試驗(yàn)有顯著的改進(jìn)。

3.3 位密度與溫度分布

圖4 年平均正壓流函數(shù)（單位：Sv）：（a）對(duì)照試驗(yàn)；（b）浮力試驗(yàn)。等值線間隔均為20Fig. 4 Annual mean barotropic stream function from (a)CONTROL and (b)BUOYANCY. Units: Sv; contour interval: 20

圖5 （a）觀測(cè)位密度緯向平均分布以及（b）對(duì)照試驗(yàn)－觀測(cè)、（c）浮力試驗(yàn)－觀測(cè)的差值分布。單位：kg m-3；（a）等值線間隔為0.2,（b）和（c）等值線間隔為0.04Fig. 5 Zonally averaged potential density: (a)Observations, (b)difference between CONTROL and observations, (c)difference between BUOYANCY and observations. Units: kg m-3; contour interval: 0.2 in (a); 0.04 in (b, c)

位密度的變化能反映海水溫度和鹽度變化的綜合效果，圖5是基于WOA09溫鹽觀測(cè)資料計(jì)算得到的海水位密度（異常值）緯向平均圖以及各試驗(yàn)與其差值圖。由于密度最劇烈的變化出現(xiàn)在上層海洋，所以我們?nèi)『１碜鳛閰⒖紝佑?jì)算海水位密度。位密度分布從上往下逐漸增大，在1000 m以上呈“W”分布型，越往下等值線越平坦（圖5a）。對(duì)照試驗(yàn)與浮力試驗(yàn)位密度整體分布形狀與觀測(cè)類似（圖略）。相對(duì)觀測(cè)資料，除北半球中低緯表層水、60°N附近600 m以上海水，以及30°N附近600～1000 m深度和南半球小部分海區(qū)位密度模擬偏重外，對(duì)照試驗(yàn)大部分海區(qū)海水位密度模擬偏輕（圖 5b）。浮力試驗(yàn)與觀測(cè)的偏差和對(duì)照試驗(yàn)的類似（圖 5c）。兩個(gè)試驗(yàn)之間差異主要體現(xiàn)在南半球15°S以南海區(qū)。由于對(duì)照試驗(yàn)里過(guò)強(qiáng)的南極中層水輸送在浮力試驗(yàn)中有所遏制，使得南極中層水覆蓋區(qū)海水位密度有所增加，大部分區(qū)域與觀測(cè)更為接近，但對(duì)45°S附近600～1000 m深度位密度增的稍大，這是由于Deacon Cell環(huán)流在該處的減弱導(dǎo)致暖水輸送不夠所致。由于浮力試驗(yàn)對(duì) 40°S以南2000～3000 m深度溫度模擬升高，以及鹽度模擬降低導(dǎo)致在該深度區(qū)位密度模擬有所降低。

海表溫度 SST是檢驗(yàn)物理海洋模式非常重要的指標(biāo)（張學(xué)洪等，2013），同時(shí)SST對(duì)海洋吸收被動(dòng)示蹤物也有直接影響。圖6是WOA09年平均SST、對(duì)照試驗(yàn)結(jié)果與觀測(cè)差值以及浮力試驗(yàn)與對(duì)照試驗(yàn)?zāi)M結(jié)果差值分布圖。從圖 6a可以看出，年平均SST顯示赤道海區(qū)溫度較高，在赤道印度洋東側(cè)和西太平洋附近為暖池，最高溫度超過(guò)28°C，在東赤道太平洋有一個(gè)明顯的冷舌伸向西側(cè)。溫度由赤道向周圍遞減，南北極SST達(dá)到最低可達(dá)-1°C以下。北半球海盆由于邊界流的影響出現(xiàn)西側(cè)高于東側(cè)的現(xiàn)象，而在南半球 30°S以南，溫度經(jīng)向梯度較大，由于南極繞流的存在等溫線緯向分布平直。兩個(gè)試驗(yàn)的模擬結(jié)果與觀測(cè)在分布形態(tài)以及數(shù)值上均相似（圖略）。從圖6b可以看出，對(duì)照試驗(yàn)與觀測(cè)的相對(duì)差異主要分布在海洋垂直運(yùn)動(dòng)較強(qiáng)的海區(qū)，比如南大洋混合區(qū)、西北大西洋、西北太平洋以及赤道上升流區(qū)，其中在西北大西洋和西北太平洋邊界差值較大，分別達(dá)4°C和3°C。為了驗(yàn)證兩個(gè)試驗(yàn)差異是否為試驗(yàn)方案差異引起，我們對(duì)兩個(gè)試驗(yàn)十年月平均SST差值圖做了顯著性t檢驗(yàn)（120個(gè)樣本），通過(guò)95%顯著性檢驗(yàn)的區(qū)域主要分布在南大洋部分海區(qū)（圖 6c）。從上文對(duì)環(huán)流場(chǎng)的分析可知，浮力試驗(yàn)的Deacon Cell在南大洋上層有所加強(qiáng)，這有利于南大洋Deacon Cell下沉支所處海表海水溫度降低，而上升支海表海水溫度增加。除此以外，從表1可以看出，兩個(gè)試驗(yàn)間年平均 SST模擬的改變跟其年平均混合層深度的變化存在較大的正相關(guān)，混合層深度增加的海區(qū)浮力試驗(yàn)海表的SST就會(huì)有所增加。通過(guò)南大洋150～300 m（年平均混合層深度以下）不同深度層次的年平均渦旋擴(kuò)散系數(shù)水平分布來(lái)看（圖略），渦旋擴(kuò)散系數(shù)的變化與對(duì)應(yīng)層的溫度的變化有一定的正相關(guān)關(guān)系，即與對(duì)照試驗(yàn)相比，浮力試驗(yàn)渦旋擴(kuò)散系數(shù)增大的海區(qū)，有利于對(duì)應(yīng)深度海區(qū)溫度的增加，反之亦然。比如在通過(guò)顯著性檢驗(yàn)的南太平洋高緯160°E～160°W附近（圖6c），浮力試驗(yàn)相對(duì)對(duì)照試驗(yàn)在中層水區(qū)域渦旋擴(kuò)散系數(shù)有所增大，模擬水溫升高，而垂直運(yùn)動(dòng)加強(qiáng)使得升溫的水與海表混合，從而導(dǎo)致SST增大。

圖6 （a）WOA09觀測(cè)年平均海表溫度以及（b）對(duì)照試驗(yàn)－觀測(cè)、（c）浮力試驗(yàn)－對(duì)照試驗(yàn)的差值分布。單位：°C；（a）等值線序列為 –1，0，2，4，6，8，12，16，18，20，22，24，26，28；（b）等值線序列為-2，-1，-0.5，0，0.5，1，2，4；（c）等值線間隔：0.25Fig. 6 Global distribution of annual mean sea surface temperature: (a)WOA09 observation; (b)difference between CONTROL and the observation;(c)difference between BUOYANCY and the observation. Units: °C; contours are drawn at –1, 0, 2, 4, 8, 12, 16, 18, 20, 22, 24, 26, and 28 in (a), and -2, -1, -0.5,0, 0.5, 1, 2, and 4 in (b); contour interval is 0.25 in (c)

表1 浮力試驗(yàn)與對(duì)照試驗(yàn)?zāi)昶骄旌蠈由疃萂LD（Mixed Layer Depth）差值與兩個(gè)試驗(yàn)其他變量差值的空間相關(guān)系數(shù)（所有結(jié)果均通過(guò)99%信度檢驗(yàn)）Table 1 Pattern correlation coefficients between the difference in annual mean mixed layer depth (MLD)and the difference in other variables from the two runs(BUOYANCY minus CONTROL)

表2是南大洋（40°S～78°S）不同深度各試驗(yàn)與 WOA09年平均海溫的均值以及各試驗(yàn)與WOA09之間的空間均方根誤差。從表2可以看出，兩個(gè)試驗(yàn)?zāi)昶骄１頊囟萐ST的均值與WOA09觀測(cè)資料的相差不大，且各試驗(yàn)?zāi)昶骄?SST與WOA09之間的空間均方根誤差也較小，說(shuō)明兩個(gè)試驗(yàn)對(duì)南大洋年平均 SST的水平空間分布模擬均較好。對(duì)南大洋其他深度而言，對(duì)比兩個(gè)試驗(yàn)與WOA09的均值可以發(fā)現(xiàn)，模擬結(jié)果與觀測(cè)是比較接近的，但是也需要指出的是，空間均方根誤差顯示出在空間分布上，南大洋的模擬結(jié)果與觀測(cè)資料還是有一定偏差的，比如在0～150 m深度，對(duì)照試驗(yàn)?zāi)M結(jié)果與WOA09的空間均方誤差達(dá)1.26，浮力試驗(yàn)的為 1.24。對(duì)比兩個(gè)試驗(yàn)，從兩個(gè)試驗(yàn)上層（0～150 m）深度以及在150～3500 m深度的空間均方根誤差來(lái)看，浮力試驗(yàn)相對(duì)對(duì)照試驗(yàn)有一定改進(jìn)；在3500 m以下深度，由于浮力試驗(yàn)?zāi)蠘O底層水溫度模擬偏低，導(dǎo)致與觀測(cè)資料的偏差稍微加大。

表2 南大洋不同層次各試驗(yàn)與WOA09年平均海溫場(chǎng)的均值以及各試驗(yàn)與WOA09資料之間的空間均方根誤差Table 2 Averages of annual mean temperature from the runs and the observations and spatial mean square error between the runs and the observations at the different layers of the Southern Ocean

3.4 CFC-11的吸收、儲(chǔ)藏與輸送

本文選取的 CFC-11觀測(cè)資料為 GLODAP（Global Ocean Data Analysis Project）提供（Key et al.,2004）。GLODAP資料綜合了 WOCE（World Ocean Circulation Experience）、JGOFS（Joint Global Ocean Flux Study）、OACES（Ocean Atmosphere Carbon Exchange Study）等大洋環(huán)流計(jì)劃資料，其中CFC-11數(shù)據(jù)中0～1200 m深度的CFC-11數(shù)據(jù)是來(lái)自于 WOCE航測(cè)數(shù)據(jù)，該計(jì)劃航測(cè)時(shí)間集中在1991～1997年，故我們?nèi)≡囼?yàn)1991～1997年的平均結(jié)果表示1990年代與之比較。

圖7是各試驗(yàn)以及觀測(cè)的1990年代CFC-11海表濃度分布以及試驗(yàn)間的差值圖。從圖 7a觀測(cè)資料可以看出，CFC-11海表濃度從低緯向高緯增加，這是由于CFC-11從大氣通過(guò)海氣交換進(jìn)入海洋，而海洋吸收CFC-11主要受海溫影響，海溫低則CFC-11溶解度高的緣故。在北半球30°N附近由于邊界流場(chǎng)的影響，西邊界海水從赤道帶來(lái)低濃度的CFC-11，而東邊界從高緯帶來(lái)高濃度的CFC-11，使得北半球海盆同一緯度西邊濃度較東邊低。上述特征各試驗(yàn)?zāi)M結(jié)果與觀測(cè)相似。從圖7b可看出，對(duì)照試驗(yàn)與觀測(cè)的偏差主要是赤道以及熱帶大部分海區(qū)模擬偏低，尤其在赤道太平洋西海岸附近水域，最大偏差達(dá)-1 pmol kg-1，模擬值約為觀測(cè)的50%，這是由于模式模擬暖池溫度偏高的緣故。在其他部分海區(qū)對(duì)照試驗(yàn)?zāi)M偏高，在西北太平洋、西北大西洋以及南大西洋高緯60°W附近存在最大偏差為1 pmol kg-1，偏高觀測(cè)約20%～25%。浮力試驗(yàn)與對(duì)照試驗(yàn)?zāi)M基本相似，兩者的差異能通過(guò)95%顯著性t檢驗(yàn)（84個(gè)樣本）的區(qū)域在南大洋附近和北冰洋（沒(méi)有觀測(cè)資料），浮力試驗(yàn) 60°S以南海區(qū)CFC-11濃度模擬較對(duì)照試驗(yàn)高（圖7c）。一方面，SST的差異會(huì)引起海表CFC-11濃度的分布的差異，SST的增大會(huì)減小海表CFC-11的溶解度。另一方面，從表 1可以看出，MLD的變化與海表 CFC-11在南大洋 40°S以南的相關(guān)系數(shù)達(dá)-0.22（通過(guò)99%置信度檢驗(yàn)），這是因?yàn)镸LD的改變會(huì)影響CFC-11垂直輸送的能力。比如在南太平洋高緯混合帶，與混合層深度加深海區(qū)對(duì)應(yīng)，浮力試驗(yàn)與對(duì)照試驗(yàn)CFC-11差值在南印度洋東側(cè)、南太平洋西側(cè)以及南大西洋 50°S緯度帶附近存在負(fù)值區(qū)。兩個(gè)試驗(yàn)?zāi)M與觀測(cè)均很接近，對(duì)照試驗(yàn)海表CFC-11濃度與觀測(cè)值之間的空間均方根誤差為0.29 pmol kg-1，浮力試驗(yàn)的為0.31 pmol kg-1，對(duì)照試驗(yàn)在海表濃度的模擬上稍好于浮力試驗(yàn)。

圖7 1990年代（a）觀測(cè)CFC-11海表濃度以及（b）對(duì)照試驗(yàn)減去觀測(cè)、（c）浮力試驗(yàn)減去對(duì)照試驗(yàn)的差值分布。單位：pmol kg-1；（a）等值線序列為 0.5，1，2，3，4，5，6，6.5；（b）等值線序列為-1，-0.5，-0.2，0，0.2，0.5，1；（c）為間隔 0.4Fig. 7 Global distribution of CFC-11 sea surface concentration in the 1990s: (a)Observations; (b)difference between CONTROL and the observation;(c)difference between BUOYANCY and CONTROL. Units: pmol kg-1; contours are drawn at 0.5, 1, 2, 3, 4, 5, 6, and 6.5 in (a), and -1, 0.5, -0.2, 0, 0.2, 0.5, and 1 in (b); contour interval is 0.4 in (c)

將1944年1月至1994年12月海洋對(duì)CFC-11吸收通量進(jìn)行時(shí)間積分，可求出51 a海洋對(duì)CFC-11的累積吸收量。圖8是兩個(gè)試驗(yàn)截止至1994年12月的累計(jì)吸收量及其差值分布圖。從圖 8a對(duì)照試驗(yàn)累積吸收量可以看出，西北太平洋、西北大西洋、南半球副熱帶以南是海洋吸收 CFC-11的主要區(qū)域。對(duì)照試驗(yàn)在西北太平洋最大累積吸收量達(dá) 6 μmol m-2，在西北大西洋的最大值為12 μmol m-2，這種空間差異跟大洋的垂直運(yùn)動(dòng)強(qiáng)度有關(guān)，在西北大西洋北半球冬季存在較強(qiáng)的垂直混合，且北大西洋深水的下沉也會(huì)將海表吸收的 CFC-11帶至深海，從而使得海表處于較低濃度狀態(tài)，因此西北大西洋累積海氣通量比較大。在55°S的20°E和100°E附近，對(duì)照試驗(yàn)存在高值中心，強(qiáng)度與西北大西洋相當(dāng)，但南大洋廣闊的海域使得其成為CFC-11最大的吸收海區(qū)。浮力試驗(yàn)的累積吸收量大體跟對(duì)照試驗(yàn)類似（圖 8b），但兩者之間的差異也較明顯，圖 8c是兩個(gè)試驗(yàn)累積吸收量的差值分布圖，可以看出，南大西洋高緯和南太平洋高緯度西側(cè)大部分海區(qū)以及西北大西洋部分海區(qū)，浮力試驗(yàn) CFC-11累積吸收量相對(duì)對(duì)照試驗(yàn)要高，而在南印度洋副熱帶以及南太平洋副熱帶東側(cè)要低。海洋累計(jì)吸收量跟海洋垂直混合強(qiáng)度有著很大的關(guān)系，從表1可以看出，混合層深度與CFC-11累計(jì)吸收量呈現(xiàn)一定的正相關(guān)。這是因?yàn)榇怪被旌蠌?qiáng)度的加大會(huì)促使海表CFC-11向下運(yùn)動(dòng)，從而加強(qiáng)海洋吸收量，反之亦然。比如，在南太平洋副極地，其東側(cè)浮力試驗(yàn)MLD的減小直接影響到累積吸收量的減小，最大達(dá)-1 μmol m-2，而西側(cè)則相反。在65°S南大西洋50°W 附近，由于浮力試驗(yàn)在該處出現(xiàn)了一個(gè)相對(duì)對(duì)照試驗(yàn)深的多的混合層（比對(duì)照試驗(yàn)深約 1000 m），使得該處吸收了較多的CFC-11。

圖8 （a）對(duì)照試驗(yàn)、（b）浮力試驗(yàn)截止至1994年12月CFC-11累計(jì)吸收量以及（c）浮力試驗(yàn)減去對(duì)照試驗(yàn)差值分布（單位：μmol m-2）。（a）和（b）等值線序列為 1，2，4，6，8，10；（c）等值線序列為-2，-1，-0.5，0，0.5，1，2，4，6，8，10Fig. 8 Global distribution of CFC-11 cumulative uptake up to December 1994: (a)CONTROL; (b)BUOYANCY; (c)the difference between BUOYANCY and CONTROL. Units: μmol m-2; contours are drawn at 1, 2, 4, 8, and 10 in (a, b), and -2, -1, -0.5, 0, 0.5, 1, 2, 4, 6, 8, and 10 in (c)

對(duì)CFC-11濃度進(jìn)行整層垂直積分得到其海洋水柱總量的分布（圖9）。圖9a為觀測(cè)資料1990年代CFC-11水柱總量分布，可以看出北太平洋副熱帶、西北大西洋、南半球副熱帶以南大部分海區(qū)為CFC-11主儲(chǔ)藏區(qū)，在赤道海區(qū)附近儲(chǔ)藏較少。從圖9b可以看出，對(duì)照試驗(yàn)1990年代的水柱總量模擬結(jié)果與觀測(cè)基本一致，主要差異是對(duì)照試驗(yàn)在北太平洋副熱帶以及南大洋40°S～50°S模擬偏低，西北大西洋模擬偏高。在北太平洋副熱帶，觀測(cè)資料存在最大值2 μmol m-2，對(duì)照試驗(yàn)?zāi)M較低，最大值僅約為1.2 μmol m-2，這是由于對(duì)照試驗(yàn)在太平洋 30°N附近中層水向下輸送偏弱，導(dǎo)致在太平洋副熱帶儲(chǔ)藏量偏低。在南大洋高緯，海洋吸收CFC-11后沿等密度面向下向北輸送，在40°S～50°S達(dá)到最大深度，使得該附近海區(qū)水柱總量較大，觀測(cè)資料最大值達(dá) 4 μmol m-2，高于對(duì)照試驗(yàn)的3 μmol m-2。在西北大西洋，由于該處冬季存在較強(qiáng)的垂直混合，且北大西洋深水的下沉可將CFC-11輸送至2000 m以下深度，導(dǎo)致該處成為全球海洋水柱總量最大的海區(qū)，觀測(cè)資料為10 μmol m-2，由前文分析可知對(duì)照試驗(yàn)在該處MLD模擬偏深，導(dǎo)致對(duì)照試驗(yàn)該處最大儲(chǔ)藏量模擬偏大，達(dá)12 μmol m-2，同時(shí)可以看出該處對(duì)照試驗(yàn)高值區(qū)范圍較觀測(cè)的?。ū热绲戎稻€5 μmol m-2），這也說(shuō)明試驗(yàn)在該處南向輸送強(qiáng)度偏弱的問(wèn)題（對(duì)照試驗(yàn)北大西洋深水經(jīng)向最大輸送僅為12 Sv，相對(duì)觀測(cè)偏弱）。在赤道海區(qū)，尤其在暖池附近，由于對(duì)照試驗(yàn)暖池海表溫度模擬偏高，使得CFC-11海表濃度模擬偏低，導(dǎo)致赤道附近水柱總量模擬稍有偏低。

浮力試驗(yàn)?zāi)M 1990年代的水柱總量與對(duì)照試驗(yàn)有著類似的結(jié)果（圖9b，c），但對(duì)南大洋部分海域偏低的情況改善明顯，這跟浮力試驗(yàn)改變了南大洋累積吸收量有關(guān)。兩個(gè)試驗(yàn)差異通過(guò)95%顯著性檢驗(yàn)的海區(qū)主要在南大洋 50°S以南高緯海區(qū)（圖9d），該區(qū)域?qū)φ赵囼?yàn)大部分海區(qū)偏低于觀測(cè)資料，其中對(duì)照試驗(yàn)在60°S以南0～70°E附近為1.5 μmol m-2低值區(qū)，比觀測(cè)值偏低40%左右，而浮力試驗(yàn)在該處儲(chǔ)藏CFC-11有明顯的增強(qiáng)，比對(duì)照試驗(yàn)大部分區(qū)域增強(qiáng) 0.5 μmol m-2，從該區(qū)域斷面CFC-11濃度資料分析可知（圖略），這是由于浮力試驗(yàn)加強(qiáng)了該剖面60°S以南CFC-11向下的滲透以及高的海表濃度導(dǎo)致水柱總量增大。同樣在南大西洋 40°S～60°S，浮力試驗(yàn)的大值中心 3.5 μmol m-2明顯大于對(duì)照試驗(yàn)，與觀測(cè)資料更接近，而在70°S南大西洋60°W～30°W之間，由于累計(jì)吸收通量的過(guò)大，浮力試驗(yàn)存在較小區(qū)域模擬偏高，最大偏差達(dá) 5 μmol m-2。

利用模式截止至1994年12月的CFC-11海洋累積吸收量減去1994年12月的水柱總量就可以反映出其在海洋內(nèi)部的輸送情況。給海洋粗略劃分成13個(gè)區(qū)域（圖10），分別計(jì)算出各區(qū)域的總累積吸收量以及各緯度帶上的輸送量（包括平流以及擴(kuò)散作用輸送的量）。從圖10的輸送量可以看出，南北半球高緯海區(qū)輸出量較大。在南半球 34°S以南至整個(gè)南大洋海區(qū)，對(duì)照試驗(yàn)中該處共吸收了29.73×107mol CFC-11，同時(shí)該區(qū)域存在較大的輸出量，輸出量達(dá)4.54×107mol，而在兩半球低緯至副熱帶海區(qū)輸入量較大，尤其在南半球10°S～34°S海區(qū)，對(duì)照試驗(yàn)輸入量達(dá)4.65×107mol。北大西洋高緯略有不同，可以看出北大西洋高緯為吸收大值區(qū)，對(duì)照試驗(yàn)在40°N～75°N之間吸收CFC-11量達(dá)5.44×107mol，同時(shí)從輸送量來(lái)看，該處還為較弱的輸入?yún)^(qū)，說(shuō)明該處同時(shí)儲(chǔ)藏了較多的 CFC-11。對(duì)北大西洋進(jìn)一步細(xì)分，在40°N～50°N海區(qū)可以算出存在較弱的輸出區(qū)，輸出量達(dá)0.04×107mol。浮力試驗(yàn)與對(duì)照試驗(yàn)整體模擬基本相似，如前文所述，由于渦旋擴(kuò)散系數(shù)的空間變化，使得浮力試驗(yàn)的混合強(qiáng)度與對(duì)照試驗(yàn)產(chǎn)生差異，進(jìn)而對(duì) CFC-11的吸收和輸送產(chǎn)生影響。兩個(gè)試驗(yàn)的差異主要體現(xiàn)在南大洋海區(qū)。在南大洋 60°S以南，浮力試驗(yàn)總累積吸收比對(duì)照試驗(yàn)大1.39×107mol，輸出量增大了 11%；在 34°S～60°S緯度帶，該區(qū)域?yàn)楹Ｑ驝FC-11主儲(chǔ)藏區(qū)，其儲(chǔ)藏量大于該區(qū)域的吸收量，對(duì)比兩個(gè)試驗(yàn)，盡管浮力試驗(yàn)在該區(qū)域的吸收量比對(duì)照試驗(yàn)要少，但由于浮力試驗(yàn)從 60°S以南輸送過(guò)來(lái)的CFC-11量比對(duì)照試驗(yàn)要大，從34°S向南輸出的量要小，因此該區(qū)域的凈輸入量浮力試驗(yàn)比對(duì)照試驗(yàn)約增大 116%，使得兩個(gè)試驗(yàn)在該區(qū)域CFC-11的儲(chǔ)藏量是相當(dāng)?shù)摹?/p>

圖9 （a）觀測(cè)、（b）對(duì)照試驗(yàn)、（c）浮力試驗(yàn)1990年代CFC-11水柱總量，（d）浮力試驗(yàn)減去對(duì)照試驗(yàn)差值分布。單位：μmol m-2；（a）、（b）和（c）等值線序列為：0，0.25，0.5，1，1.5，2，2.5，3.5，5，8，12；（d）等值線序列為-0.5，-0.1，0，0.5，1，1.5，2，2.5，3，4（深灰陰影為通過(guò)95%信度檢驗(yàn)區(qū)域）Fig. 9 Global distributions of CFC-11 column inventory in the 1990s from (a)observations, (b)CONTROL, and (c)BUOYANCY, and (d)the difference between BUOYANCY and CONTROL. Units: μmol m-2; contours are drawn at 0, 0.25, 0.5, 1, 1.5, 2, 2.5, 3.5, 5, 8, and 12 in (a–c), and -0.5, -0.1, 0, 0.5, 1,1.5, 2, 2.5, 3, and 4 in (d). Dark gray shading indicates the 95% confidence levels based on a two-sided Student’s t-test

圖10 截止至1994年12月對(duì)照試驗(yàn)和浮力試驗(yàn)中海洋CFC-11的累計(jì)吸收量與輸送，單位均為107 mol。箭頭代表該緯度帶的輸送量（箭頭大小不代表數(shù)值大小），箭頭方向代表輸送方向；黑色粗體數(shù)值代表該區(qū)域的總吸收量，黑色細(xì)體數(shù)值代表該緯度帶的輸送量；括號(hào)內(nèi)的為浮力試驗(yàn)對(duì)應(yīng)的值Fig. 10 Cumulative uptake and transport of CFC-11 in the ocean by December 1994 in units of 107 mol. An upward (downward)arrow indicates the northward (southward)movement of CFC-11; black bold type numbers represent cumulative uptake, and thin type numbers near arrows indicate transport. The values in brackets are from BUOYANCY

由上述分析可知試驗(yàn)在南大洋的差異最大，故選取由WOCE 航測(cè)提供的南太平洋170°W斷面資料，該斷面從15°S～75°S，時(shí)間為1996年1～3月。各試驗(yàn)均取1996年2月平均值結(jié)果與其對(duì)比。圖11分別為觀測(cè)、對(duì)照試驗(yàn)以及浮力試驗(yàn)在該斷面的分布情況，背景表示對(duì)應(yīng)的位密度分布。由于高緯海表溫度較低，吸收大量的 CFC-11，海表濃度較高，隨著緯度的減小，海表溫度越高，海表CFC-11濃度逐漸降低。由于被動(dòng)示蹤物CFC-11主要是沿等密度面輸送，從觀測(cè)資料可以看出，CFC-11的分布與等密度面分布吻合，例如在 62°S附近的兩側(cè)，150 m以下等密度面有很明顯的向下彎曲，這勢(shì)必使沿著等密度面輸送的被動(dòng)示蹤物向下輸送。等密度面在 40°S～50°S存在最大深度，也使得CFC-11的穿透深度達(dá)到了最大。在分布形態(tài)上，對(duì)照試驗(yàn)與浮力試驗(yàn)位密度以及CFC-11剖面的整體分布基本相似，兩個(gè)試驗(yàn)相比，浮力試驗(yàn)稍好于對(duì)照試驗(yàn)，例如在海表 65°S附近，觀測(cè)資料與浮力試驗(yàn)均出現(xiàn)明顯的高值中心，為6.5 pmol kg-1，并有明顯的向北向下輸送，而對(duì)照試驗(yàn)海表的高值中心處于 70°S以南海表，且向下輸送不明顯。同時(shí)在高值中心下200 m以下海區(qū)，對(duì)照試驗(yàn)向下穿透過(guò)強(qiáng)，比如0.5 pmol kg-1等值線，觀測(cè)資料在63°S附近的最淺深度在 250 m左右，浮力試驗(yàn)在 66°S附近最淺深度能控制在 500 m左右，對(duì)照試驗(yàn)在65°S附近達(dá)到最淺深度為800 m。基于以上原因，對(duì)照試驗(yàn)相對(duì)浮力試驗(yàn)在南半球高緯吸收的CFC-11在60°S以南的海區(qū)儲(chǔ)藏得更多，而浮力試驗(yàn)向 60°S以北的海區(qū)輸送的較多。浮力試驗(yàn)對(duì)南半球高緯上層海水的位密度模擬偏低情況有所改善使得該剖面相對(duì)對(duì)照試驗(yàn)更符合觀測(cè)事實(shí)。此外需要指出的是，兩個(gè)試驗(yàn)CFC-11分布與各自位密度的吻合程度要弱于觀測(cè)資料，說(shuō)明模式在南半球高緯度低估了沿等密度面擴(kuò)散的作用，這提醒我們，模式其他物理過(guò)程與等密度擴(kuò)散之間的關(guān)系仍需進(jìn)一步研究分析。

4 結(jié)論

本文使用中國(guó)科學(xué)院大氣物理研究所發(fā)展的LICOM海洋環(huán)流模式，考察了GM90參數(shù)化方案中不同渦旋擴(kuò)散系數(shù)取值對(duì)模式模擬結(jié)果的影響。做了兩個(gè)試驗(yàn)，一個(gè)為常系數(shù)方案，即渦旋擴(kuò)散系數(shù)均取1.0×103m2s-1的對(duì)照試驗(yàn)；另一個(gè)為在非絕熱層下渦旋擴(kuò)散系數(shù)隨浮力頻率變化的浮力試驗(yàn)。通過(guò)對(duì)比兩個(gè)試驗(yàn)物理場(chǎng)以及對(duì)被動(dòng)示蹤物CFC-11的分布、吸收、儲(chǔ)藏以及輸送的情況，具體得出以下結(jié)論：

（1）次網(wǎng)格參數(shù)取值的不同對(duì)環(huán)流場(chǎng)的影響是顯著的。兩個(gè)試驗(yàn)對(duì)全經(jīng)向流函數(shù)的模擬主要差異出現(xiàn)在南半球，浮力試驗(yàn)通過(guò)采用渦旋擴(kuò)散系數(shù)的空間變化方案，使得模擬的渦旋誘導(dǎo)流場(chǎng)發(fā)生變化，導(dǎo)致南半球Deacon Cell上層以及深海的環(huán)流加強(qiáng)，對(duì)過(guò)強(qiáng)的中層水有一定的抑制，對(duì)南極底層水的模擬有所減弱。同時(shí)浮力試驗(yàn)在南大洋正壓流函數(shù)的模擬上明顯好于對(duì)照試驗(yàn)，比對(duì)照試驗(yàn)增大約20%～30%，更加接近觀測(cè)事實(shí)。

圖11 1996年2月南太平洋170°W斷面CFC-11濃度的分布（單位：pmol kg-1）：(a)觀測(cè)資料；(b)對(duì)照試驗(yàn)；(c)浮力試驗(yàn)。填色圖為該剖面對(duì)應(yīng)資料位密度分布（單位：kg m-3）Fig. 11 The vertical distribution of CFC-11 concentrations along the 170°W section in the South Pacific in Feb 1996 from (a)observations, (b)CONTROL,and (c)BUOYANCY (contour, units: pmol kg-1); shading indicates the distribution of potential density (shaded , units: kg m-3)

（2）通過(guò)與觀測(cè)資料位密度的對(duì)比，由于浮力試驗(yàn)對(duì)對(duì)照試驗(yàn)?zāi)M過(guò)強(qiáng)的南極中層水輸送有所削弱，使得南極上層海水位密度有所增加，大部分區(qū)域與觀測(cè)更為接近。在對(duì)海溫的模擬上，兩個(gè)試驗(yàn)同樣在南大洋差異較為顯著，兩者海表模擬結(jié)果接近，在海表以下至3500 m深度浮力試驗(yàn)?zāi)M結(jié)果更接近觀測(cè)資料。

（3）通過(guò)與CFC-11海表分布、水柱總量觀測(cè)資料的對(duì)比可知，兩個(gè)試驗(yàn)均能很好的模擬CFC-11在海洋中的分布以及儲(chǔ)藏情況。模式與觀測(cè)資料的差異主要表現(xiàn)在：兩個(gè)試驗(yàn)均對(duì)西北大西洋、西北太平洋以及南大洋大部分海表CFC-11濃度模擬偏高，對(duì)赤道以及熱帶海區(qū)模擬偏低；對(duì)南大洋40°S～50°S主儲(chǔ)藏區(qū)水柱總量模擬偏低。兩個(gè)試驗(yàn)之間在屬于主吸收區(qū)和主儲(chǔ)藏區(qū)的南大洋的相對(duì)差異較大。在南大洋，海表CFC-11濃度與SST和混合層深度的變化均呈現(xiàn)較強(qiáng)的反相關(guān)關(guān)系，如浮力試驗(yàn)在高緯度的南印度洋東側(cè)和南太平洋西側(cè)以及南大西洋50°S緯度帶附近SST較對(duì)照試驗(yàn)偏高，模擬的海表CFC-11濃度也有所偏低，而在靠近南極大陸的海區(qū)模擬偏高。

（4）從對(duì)CFC-11的輸送量可以看出，浮力試驗(yàn)在靠近南極大陸海區(qū)向南半球主儲(chǔ)藏區(qū)的輸送量增大，增大約為對(duì)照試驗(yàn)的11%，使得浮力試驗(yàn)對(duì)水柱總量的模擬稍好于對(duì)照試驗(yàn)。從CFC-11剖面圖可以看出，由于浮力試驗(yàn)對(duì)南大洋高緯上層位密度結(jié)構(gòu)的模擬比對(duì)照試驗(yàn)要好，如位密度的穿透深度以及位密度向北向下延伸結(jié)構(gòu)等的模擬上，浮力試驗(yàn)相對(duì)對(duì)照試驗(yàn)與觀測(cè)資料更接近，使得浮力試驗(yàn)從高緯向北輸送的 CFC-11比對(duì)照試驗(yàn)的要強(qiáng)，因而在南大洋對(duì)CFC-11的模擬上更接近觀測(cè)資料。總的來(lái)說(shuō)，通過(guò)對(duì)比分析物理場(chǎng)和 CFC-11結(jié)果后，可認(rèn)為Danabasoglu et al.（2007）提出的浮力方案是一個(gè)很有前景的中尺度參數(shù)化方案，包括方案中更符合實(shí)際情況的渦旋擴(kuò)散系數(shù)分布以及對(duì)模式模擬能力的提升等等，但也存在不足的地方，對(duì)該方案的進(jìn)一步調(diào)試和檢驗(yàn)仍有很大的必要。

（

）

Danabasoglu G, McWilliams J C. 1995. Sensitivity of the global ocean circulation to parameterizations of mesoscale tracer transports [J]. J.Climate, 8 (12): 2967-2987.

Danabasoglu G, Marshall J. 2007. Effects of vertical variations of thickness diffusivity in an ocean general circulation model [J]. Ocean Modelling, 18(2): 122-141.

Danabasoglu G, Ferrari R, McWilliams J C. 2008. Sensitivity of an ocean general circulation model to a parameterization of near-surface eddy fluxes [J]. J. Climate, 21 (6): 1192-1208.

Danabasoglu G, McWilliams J C, Gent P R. 1994. The role of mesoscale tracer transports in the global ocean circulation [J]. Science, 264 (5162):1123-1126.

Danabasoglu G, Bates S C, Briegleb B P, et al. 2012. The CCSM4 ocean component [J]. J. Climate, 25 (5): 1361-1389.

England M H, Hirst A C. 1997. Chlorofluorocarbon uptake in a world ocean model: 2. Sensitivity to surface thermohaline forcing and subsurface mixing parameterizations [J]. J. Geophys. Res.: Oceans, 102 (C7):15709-15731.

Ferreira D, Marshall J, Heimbach P. 2005. Estimating eddy stresses by fitting dynamics to observations using a residual-mean ocean circulation model and its adjoint [J]. J. Phys. Oceanogr., 35 (10): 1891-1910.

Gent P R, McWilliams J C. 1990. Isopycnal mixing in ocean circulation models [J]. J. Phys. Oceanogr., 20 (1): 150-155.

Gent P R, Willebrand J, Mcdougall T J, et al. 1995. Parameterizing eddy-induced tracer transports in ocean circulation models [J]. J. Phys.Oceanogr., 25 (4): 463-474.

Key R M, Kozyr A, Sabine C L, et al. 2004. A global ocean carbon climatology: Results from Global Data Analysis Project (GLODAP)[J].Global Biogeochem. Cycles, 18 (4), doi:10.1029/2004GB002247.

Large W G, Danabasoglu G, Doney S C, et al. 1997. Sensitivity to surface forcing and boundary layer mixing in a global ocean model: Annual-mean climatology [J]. J. Phys. Oceanogr., 27 (11): 2418-2447.

Ledwell J R, Watson A J, Law C S. 1998. Mixing of a tracer in the pycnocline [J]. J. Geophys. Res.: Oceans, 103 (C10): 21499-21529.

李陽(yáng)春, 徐永福, 趙亮, 等. 2007. 全球海洋模式中CFC-11吸收對(duì)次網(wǎng)格尺度混合參數(shù)化的敏感性 [J]. 海洋學(xué)報(bào), 29 (3): 31-38. Li Yangchun, Xu Yongfu, Zhao Liang, et al. 2007. Sensitivity of CFC-11 uptake in a global ocean model to subgrid-scale mixing parameterizations[J]. Acta Oceanologica Sinica (in Chinese), 29 (3): 31-38.

劉海龍, 俞永強(qiáng), 李薇, 等. 2004. LASG/IAP 氣候系統(tǒng)海洋模式 (LICOM1. 0)參考手冊(cè) [M]. 北京: 科學(xué)出版社, 107pp. Liu Hailong, Yu Yongqiang,Li Wei, et al. 2004. LASG/IAP Climate System Ocean Model (LICOM 1.0):User Manual (in Chinese)[M]. Beijing: Science Press, 107pp.

Liu Hailong, Zhang Xuehong, Li Wei, et al. 2004. An eddy-permitting oceanic general circulation model and its preliminary evaluations [J]. Adv.Atmos. Sci., 21 (5): 675-690.

Liu Hailong, Lin Pengfei, Yu Yongqiang, et al. 2012. The baseline evaluation of LASG/IAP Climate system Ocean Model (LICOM)version 2 [J]. Acta Meteor. Sinica, 26 (3): 318-329.

Marshall J, Shuckburgh E, Jones H, et al. 2006. Estimates and implications of surface eddy diffusivity in the Southern Ocean derived from tracer transport [J]. J. Phys. Oceanogr., 36 (9): 1806-1821.

Matsumoto K, Sarmiento J L, Key R M, et al. 2004. Evaluation of ocean carbon cycle models with data-based metrics [J]. Geophys. Res. Lett., 31(7), doi:10.1029/2003GL018970.

Rayner N A, Parker D E, Horton E B, et al. 2003. Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century [J]. J. Geophys. Res.: Atmos, 108 (D14),doi:10.1029/2002JD002670.

Wanninkhof R. 1992. Relationship between wind-speed and gas-exchange over the ocean [J]. J. Geophys. Res.: Oceans, 97 (C5): 7373-7382.

Warner M J, Weiss R F. 1985. Solubilities of chlorofluorocarbon-11 and chlorofluorocarbon-12 in water and seawater [J]. Deep Sea Res., Part I,32 (12): 1485-1497.

Whitworth T, Peterson R G. 1985. Volume transport of the Antarctic circumpolar current from bottom pressure measurements [J]. J. Phys.Oceanogr., 15 (6): 810-816.

Xu Y F, Aoki S, Harada K. 2006. Sensitivity of the simulated distributions of water masses, CFCs, and bomb C-14 to parameterizations of mesoscale tracer transports in a model of the North Pacific [J]. J. Phys. Oceanogr.,36 (3): 273-285.

張學(xué)洪, 俞永強(qiáng), 劉海龍. 2003. 海洋環(huán)流模式的發(fā)展和應(yīng)用 I. 全球海洋環(huán)流模式 [J]. 大氣科學(xué), 27 (4): 607-617. Zhang Xuehong, Yu Yongqiang, Liu Hailong. 2003. The development and application of the oceanic general circulation models. Part I. The global oceanic general circulation models [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 27 (4): 607-617.

張學(xué)洪, 俞永強(qiáng), 周天軍, 等. 2013. 大洋環(huán)流和海氣相互作用的數(shù)值模擬講義 [M]. 北京: 氣象出版社, 298pp. Zhang Xuehong, Yu Yongqiang, Zhou Tianjun, et al. 2013. Lecture Notes on Numerical Simulations of Ocean Circulation and Air–Sea Interaction (in Chinese)[M]. Beijing: China Meteorological Press, 298pp.

趙琦, 陳中笑, 徐永福, 等. 2012. 全球海洋 CFC-11吸收對(duì)傳輸速度的敏感性 [J]. 大氣科學(xué), 36 (6): 1253-1268. Zhao Qi, Chen Zhongxiao,Xu Yongfu, et al. 2012. Sensitivity of CFC-11 uptake in a global ocean model to air–sea gas transfer velocity [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 36 (6): 1253-1268.