一個(gè)較高分辨率的基于通量訂正的海-氣耦合環(huán)流模式的設(shè)計(jì)和性能評估

薛洪斌，楊玉震，楊金波，趙艷玲，宋樹剛

（1、61741部隊(duì)，北京 100094；2、92493部隊(duì)，遼寧 葫蘆島 125000；3、95949部隊(duì)，河北 滄州 061736）

1 引言

當(dāng)前氣候變化研究認(rèn)為氣候變化不僅僅是大氣自身的演變，還包括海洋、冰雪、陸面、生態(tài)植被等的相互作用。其中，海-氣相互作用對氣候變化的影響非常重要，也是目前研究的重點(diǎn)。海-氣耦合環(huán)流模式為這種研究提供了平臺，被各個(gè)研究所和業(yè)務(wù)中心所關(guān)注。世界氣候研究計(jì)劃（WCRP）氣候變率與可預(yù)報(bào)性研究（CLIVAR）中心議題就是耦合模式的發(fā)展。為此，WCRP聯(lián)合科學(xué)委員會（JSC）和CLIVAR科學(xué)指導(dǎo)小組（SSG）聯(lián)合建立了耦合模擬工作組（WGCM），其目的就是評估和支持耦合氣候系統(tǒng)模式的發(fā)展。

從上世紀(jì)80年代末開始，我國自行設(shè)計(jì)了大氣環(huán)流模式、大洋環(huán)流模式及海-氣耦合模式，為我國的氣候變化模擬和預(yù)測研究提供了工具[1-5]。在發(fā)展的初期，受計(jì)算能力所限模式的分辨率相對較低。隨著計(jì)算能力的飛速進(jìn)步，發(fā)展高分辨率環(huán)流模式的需求越來越強(qiáng)烈。高分辨率的環(huán)流模式特別是海-氣耦合模式，可以更細(xì)致地描述海-氣相互作用，對提高模式的性能具有重要意義。

本文基于“日通量距平”耦合方案將兩個(gè)具有較高分辨率的大氣和海洋環(huán)流模式耦合，耦合積分了10年，結(jié)果顯示耦合模式比較成功地抑制了“氣候漂移”現(xiàn)象，對于模式的性能評估這里僅僅給出了平均氣候態(tài)和季節(jié)變化。

2 耦合模式簡介

海-氣耦合模式包括大氣環(huán)流模式和大洋環(huán)流模式兩部分。其中，大氣環(huán)流模式為具有較高分辨率的IAPAGCM-3并行版大氣環(huán)流格點(diǎn)模式。模式框架在水平方向?yàn)榍蛎娼?jīng)緯網(wǎng)格坐標(biāo)，垂直方向取σ坐標(biāo)，模式分辨率在水平方向?yàn)?°×1°，垂直方向按σ坐標(biāo)不等距分為21層，大氣物理過程主要包括：云和輻射過程，大尺度凝結(jié)降水，對流參數(shù)化（含深對流和淺對流），動量、熱量和水汽的水平和垂直擴(kuò)散，地形重力波拖曳以及邊界層等物理過程。大洋環(huán)流模式為目前分辨率較高的Licom模式。模式框架在水平方向?yàn)榍蛎娼?jīng)緯網(wǎng)格坐標(biāo)，垂直方向取η坐標(biāo)，模式分辨率在水平方向?yàn)?.5°×0.5°，垂直方向按η坐標(biāo)不等距分為30層。海洋環(huán)流模式在提高分辨率的同時(shí)對模式地形、水平粘性方案和熱力強(qiáng)迫方案等都進(jìn)行重新設(shè)計(jì)，同時(shí)改進(jìn)對海冰的處理，增加了熱力學(xué)海冰模式。雖然大氣環(huán)流模式和海洋環(huán)流模式單獨(dú)運(yùn)行都能成功模擬和預(yù)測出大尺度大氣和海洋的許多特征，但當(dāng)二者耦合后，耦合界面上各種通量的誤差在耦合過程中不穩(wěn)定增長，將會產(chǎn)生氣候漂移現(xiàn)象。如何把這種“氣候漂移”減至最小，是設(shè)計(jì)一個(gè)成功的海氣耦合模式的主要挑戰(zhàn)之一。為了抑制這種現(xiàn)象的出現(xiàn)，早在80年代后期，Sausen[6]就提出了所謂的“通量訂正”的方法。其后，人們又對這個(gè)方法進(jìn)行改進(jìn)并應(yīng)用到不同的耦合模式中。1992年，張學(xué)洪[7]采用“月平均通量異常”成功地將二層大氣環(huán)流模式和四層海洋模式耦合。其后，1995年中國科學(xué)院大氣物理研究所又發(fā)展了一個(gè)基于“線性統(tǒng)計(jì)修正”同步耦合方案的熱帶太平洋環(huán)流模式和全球兩層大氣環(huán)流模式的耦合模式[8]。

本文采用了所謂的日通量距平耦合方案實(shí)現(xiàn)了海洋和大氣模式的成功耦合。具體步驟如下：

（1）對未耦合的海洋模式長期積分，模式積分五十年，基本能達(dá)到平衡。

（2）對未耦合的大氣模式進(jìn)行積分。模式積分五年達(dá)到平衡。用大氣模式的平衡態(tài)作初值，用未耦合的海洋模式達(dá)到平衡時(shí)的SST代替洋面氣候平均的SST強(qiáng)迫大氣模式積分五年，并計(jì)算其五年平均海氣界面的熱通量和風(fēng)應(yīng)力，這就是本文所取得參考通量（見圖1）。

（3）開始耦合積分，大氣模式和海洋模式每天交換一次通量。以海洋環(huán)流模式計(jì)算出的前一天的海表溫度（SST）作為大氣模式的下邊界強(qiáng)迫條件。

（4）大氣模式計(jì)算出當(dāng)天的海表風(fēng)應(yīng)力和熱通量根據(jù)公式（1）修正后作為海洋模式的外強(qiáng)迫條件，而影響海水鹽度的淡水通量仍采用氣候平均值。

式中Fc為OGCM和AGCM直接耦合時(shí)海洋得到的熱通量，稱為“耦合通量”；Fu為OGCM在觀測的大氣強(qiáng)迫下海洋得到的熱通量，稱為“觀測通量”；為未耦合的AGCM計(jì)算得到的熱通量，稱為“參考通量”，它的選取對耦合模式的成敗至關(guān)重要，本模式中該量選取為用OGCM算出的SST代替氣候平均的SST積分一定時(shí)間，并平均為12個(gè)月的資料作為。

（5）對耦合模式長期積分，本文積分了15年。

3 耦合模式基本氣候態(tài)的模擬

本節(jié)主要對耦合模式積分最后10年的結(jié)果進(jìn)行分析，并與未耦合的大氣環(huán)流模式積分20年和海洋環(huán)流模式積分50年的結(jié)果進(jìn)行了對比。這里僅簡要對比大氣和海洋環(huán)流的基本變量，以檢測耦合模式的性能及對“氣候漂移”的控制。

3.1 長期趨勢

SST對海氣界面上通量的變化極為敏感，全球平均SST的時(shí)間序列是衡量耦合模式長期“氣候漂移”趨勢的通用指標(biāo)。圖1給出了耦合模式積分15年的全球平均SST變化趨勢。15年時(shí)間序列表明全球平均SST沒有明顯的漂移現(xiàn)象，SST做為大氣模式的下邊界條件直接影響模式大氣，全球平均SST沒有明顯漂移趨勢保證了大氣模式?jīng)]有明顯“氣候漂移”。

圖1 耦合模式積分15年的逐月全球平均SST序列（單位/℃）

3.2 大氣環(huán)流部分

我們利用耦合模式積分后10年的結(jié)果，主要分析了耦合模式對全球氣候平均態(tài)的特征模擬。將NOAA-CIRES Climate Diagnostics Center生成的2.5°×2.5°多年平均的NCEP再分析氣候資料作為海平面氣壓場和高度場的實(shí)測資料；將2.5°×2.5°的CMAP（the CPC（Climate Prediction Center,National Oceanic and Atmospheric Administration（NOAA），U.S.A.）Merged Analysis of Precipitation Enhanced）氣候月平均降水資料作為降水的實(shí)測資料。

圖2給出了1月和7月海平面氣壓場的全球分布，可以看出，耦合模式模擬出了1月和7月海平面氣壓場全球分布的基本特征。例如，無論1月或7月對亞索爾高壓的位置、范圍和走向的模擬與再分析資料基本一致，對南緯副熱帶地區(qū)高壓范圍、中心位置等也模擬得比較準(zhǔn)確，對1月亞洲大陸高壓的強(qiáng)度、位置以及7月北太平洋高壓的中心位置、走向等模擬也比較成功。另外與未耦合模式的比較可以看出，兩者模擬能力相當(dāng)，對7月亞索爾高壓的模擬似有所改善。不足之處是無論耦合或未耦合，系統(tǒng)總體模擬偏強(qiáng)。

圖3給出了1月和7月500 hPa位勢高度場的全球分布，可以看出，兩模式（耦合和未耦合，下同）均較好地再現(xiàn)了500 hPa高度場的基本分布結(jié)構(gòu)，模擬出了冬季和夏季主要的槽脊分布特征，并較好地刻畫了北半球的季節(jié)變化情況，尤其是夏季副熱帶高壓的分布特征模擬比較成功。對南半球平直的位勢高度等值線形態(tài)，兩模式都模擬得較好。

圖4給出了1月、7月模擬和觀測的降水全球分布。1月最主要的特征是熱帶存在一條環(huán)球雨帶，在北太平洋和北大西洋分別存在一條雨帶。兩模式模擬的結(jié)果能夠再現(xiàn)冬季雨帶的這兩個(gè)主要特征，兩大洋上的雨帶位置、強(qiáng)度和走向都與觀測相當(dāng)一致，對基本環(huán)繞全球的熱帶雨帶的模擬也較好；模擬的南太平洋輻合帶降水比較清晰，與觀測較為一致；對南非大陸的降水有一定的刻畫能力；對澳大利亞干區(qū)和南半球大陸西邊海上少雨區(qū)、以及非洲北部至西亞的干旱區(qū)的模擬均很準(zhǔn)確。不足之處是兩模式對赤道印度洋北側(cè)的雨帶模擬偏強(qiáng)，南側(cè)雨帶范圍模擬偏小，南半球中緯度地區(qū)降水帶過于明顯；此外，兩模式對南美北部的較強(qiáng)降水模擬都偏弱。

圖2 一月（左）和七月（右）耦合模擬(a)、未耦合模擬(b)和NCEP再分析(c)海平面氣壓場的全球分布(單位/hPa)

圖3 1月（左）和7月（右）耦合模擬(a)、未耦合模擬(b)和NCEP再分析(c)500 hPa位勢高度場的全球分布（單位/gpm）

圖5 耦合模式(a)、海洋環(huán)流模式模擬(b)和觀測（c）的年平均海表溫度

7月，實(shí)際大氣的雨帶普遍北移，由于北半球西風(fēng)減弱，兩大洋北部雨帶也明顯減弱，而東亞出現(xiàn)強(qiáng)大的季風(fēng)降水。從圖上可以看出，兩模式能夠模擬出上述基本特征，其基本形態(tài)與CMAP的觀測比較一致。對東亞雨帶，兩模式模擬結(jié)果類似；能基本模擬出南半球太平洋輻合帶雨帶；模擬的南半球南美大陸、澳大利亞和南非南端的干區(qū)也與實(shí)測基本一致；模擬出了赤道中東太平洋ITCZ雨帶的緯向帶狀特征。此外，夏季南半球西風(fēng)帶加強(qiáng)，相應(yīng)的雨帶也比冬季強(qiáng)，模式能夠模擬出這一現(xiàn)象，但雨區(qū)的范圍有點(diǎn)過大。

3.3 海洋環(huán)流部分

耦合模式中的一個(gè)重要變量是海表溫度，它直接反應(yīng)海氣間的相互作用，做為大氣環(huán)流的下邊界條件，其模擬的好壞直接關(guān)系到耦合的成敗。圖5給出耦合模擬后5年海表溫度年平均模擬結(jié)果與未耦合海洋環(huán)流模式運(yùn)行20年的后5年年平均結(jié)果及觀測的海表溫度（NOAA Extended Reconstructed SST V3b）的2006年1月—2010年12月5年年平均結(jié)果對比分析?？梢钥闯?，兩者模擬的結(jié)果基本相似都能再現(xiàn)觀測SST的基本特征，特別是熱帶太平洋，對赤道西太平洋暖池和東太平洋的冷舌的位置和形狀具有一定的刻畫能力，但兩者對東太平洋冷舌模擬過強(qiáng)，西伸過大。

4 小結(jié)

本文基于日通量距平耦合方案，把中國科學(xué)院大氣物理研究所發(fā)展的具有較高分辨率大氣環(huán)流模式和大洋環(huán)流模式耦合起來。對耦合模式積分10年的結(jié)果分析表明，耦合模式對氣候平均態(tài)和季節(jié)變化特征具有一定的模擬能力，現(xiàn)將其主要性能歸納為以下幾點(diǎn)：

（1）耦合模式可以穩(wěn)定地長期積分，并且沒有出現(xiàn)明顯的“氣候漂移”現(xiàn)象；

（2）耦合模式能模擬出與未耦合的海-氣模式接近的氣候平均態(tài)及其季節(jié)變化；

（3）耦合模式繼承了未耦合模式的特征，既包括未耦合模式的優(yōu)點(diǎn)，同時(shí)也有缺點(diǎn)，如對副高的模擬耦合與未耦合都偏強(qiáng)。

發(fā)展耦合模式的目的是研究氣候變化，本文僅僅對耦合模式的基本性能—?dú)夂驊B(tài)進(jìn)行了初步的評估，對于耦合過程中的各種時(shí)間尺度的氣候變率還要做進(jìn)一步分析。同時(shí)，分量模式本身也需進(jìn)一步完善，如改善副高模擬偏強(qiáng)等等。

[1]梁信忠.IAP GCM的設(shè)計(jì)與氣候模擬[D].北京：中國科學(xué)院大氣物理研究所,1986：25-30.

[2]畢訓(xùn)強(qiáng).IAP九層大氣環(huán)流模式及氣候數(shù)值模擬[D].北京：中國科學(xué)院大氣物理研究所，1993:34-40.

[3]左瑞亭.新一代高分辨率格點(diǎn)大氣環(huán)流模式的研制[D].南京：解放軍理工大學(xué)氣象學(xué)院，2003:7-10.

[4]劉海龍.高分辨率海洋環(huán)流模式和熱帶太平洋上層環(huán)流的模擬研究[D].北京：中國科學(xué)院大氣物理研究所，2002：9-13.

[5]俞永強(qiáng)，張學(xué)洪.一個(gè)修正的海氣通量距平耦合方案[J].科學(xué)通報(bào),1998,43:866-870.

[6]Sausen R,Barthel K,Hasselmann K.Coupled ocean-atmosphere models with flux correction[J].Clim.Dyn,1988,m2:145-163.

[7]Zhang X H,Bao N,Yu R C,et al.Coupling scheme experiments based on an atmospheric and oceanic GCM[J].Chinese J Atmos Sci,1992,16:129-144.

[8]Zhou G Q,Zeng Q C.Prediction of ENSO with a coupled atmosphere-ocean general circulation model[J].Adv Atmos Sci,2001,18(4):587-601.

[9]薛洪斌.一個(gè)耦合環(huán)流模式的研制和應(yīng)用[D]，南京：解放軍理工大學(xué)氣象學(xué)院，2005:64-98，

[10]周廣慶,李崇銀.西太平洋暖池次表層海溫異常與ENSO關(guān)系的CGCM模擬結(jié)果[J].氣候與環(huán)境研究,1994,4:346-352.

[11]YuYQ,GuoYF. The interannual variability of climate in a coupled ocean-atmospheremodel[J].Adv.Atmos.Sci,1995,12:273-288.

[12]Zhang R H,Endoh M.A free surface general circulation model for the tropical Pacific Ocean[J].J Geophys Res,1992,97(C7):11237-11255.