簡單地球模型的研究進(jìn)展

王勇劉蘇峽邵亞平彭公炳,4

（1 中國科學(xué)院地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點(diǎn)實(shí)驗(yàn)室，北京 100101；2 南京信息工程大學(xué)遙感學(xué)院，南京 210044；3 德國科隆大學(xué)地球物理與氣象研究所，科隆 50923；4 國際歐亞科學(xué)院中國科學(xué)中心，北京 100107）

簡單地球模型的研究進(jìn)展

王勇1,2劉蘇峽1邵亞平3彭公炳1,4

（1 中國科學(xué)院地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點(diǎn)實(shí)驗(yàn)室，北京 100101；2 南京信息工程大學(xué)遙感學(xué)院，南京 210044；3 德國科隆大學(xué)地球物理與氣象研究所，科隆 50923；4 國際歐亞科學(xué)院中國科學(xué)中心，北京 100107）

以簡單概念性氣候模式（SCMs）和中等復(fù)雜程度的地球系統(tǒng)模式（EMICs）為對象，介紹了簡單地球模型的研究進(jìn)展。在簡單地球模型中，簡單概念性氣候模式已經(jīng)成為了解氣候理論和氣候變化規(guī)律的有效工具；中等復(fù)雜程度的地球系統(tǒng)模式涵蓋了地球系統(tǒng)的大多數(shù)組成部分，對各部分之間的反饋和過程描述比較詳細(xì)，同時(shí)計(jì)算要求比較低，使其應(yīng)用幾乎覆蓋了簡單概念性氣候模式以及耦合氣候系統(tǒng)模式的所有研究領(lǐng)域。EMICs必定會有更大的發(fā)展，從而在地球氣候模擬中起到更加重要的作用。

簡單地球模型，簡單概念性氣候模式，地球系統(tǒng)模式（EMICs）

1 引言

當(dāng)今的地球,氣候變化無常,這已為世人矚目。人們一直不懈地從不同角度去探求其原因,尋找合適的應(yīng)對措施。在所有的歸因研究方法中,氣候模式和地球系統(tǒng)模式等(這里統(tǒng)稱為地球模型)是氣候變化歸因研究的重要工具。

地球模型的共性是在地球系統(tǒng)中的動力、物理、化學(xué)和生物過程基礎(chǔ)上建立起來的數(shù)學(xué)方程組(包括動力學(xué)方程組和參數(shù)化方案)來確定其各個(gè)部分(大氣圈、水圈、冰雪圈、巖石圈、生物圈)的性狀,由此構(gòu)成地球系統(tǒng)的數(shù)學(xué)物理模型,然后用數(shù)值的方法進(jìn)行求解,編制成一種大型綜合性計(jì)算程序,并通過計(jì)算機(jī)付諸實(shí)現(xiàn)對地球系統(tǒng)復(fù)雜行為和過程的模擬與預(yù)測的科學(xué)工具[1-2]。

20世紀(jì)60年代末,隨著計(jì)算機(jī)技術(shù)的發(fā)展,復(fù)雜耦合氣候系統(tǒng)模式出現(xiàn)了,這類模式包括了發(fā)展成熟的大氣模式、海洋模式、陸面模式、海冰模式等模塊,用以研究包括海洋狀況、冰雪過程、土壤溫濕等在內(nèi)的氣候系統(tǒng)變化規(guī)律,是目前研究大氣、海洋及陸地之間復(fù)雜相互作用的主要工具[1,3]。這類模式基本出發(fā)點(diǎn)在于完善地描述氣候系統(tǒng)的各種過程,模型極其復(fù)雜,計(jì)算時(shí)間和花費(fèi)比較大。同時(shí),根據(jù)研究的實(shí)際需要,陸續(xù)產(chǎn)生了大量不同種類的模式,主要有概念模型、能量平衡模式(EBMs)、輻射對流模式(RCMs)、統(tǒng)計(jì)動力模式(SDMs)和地球系統(tǒng)模式(EMICs)等。目前,依據(jù)復(fù)雜程度的不同,廣

泛得到應(yīng)用的地球模式可分為3類[4-5]：(1)簡單概念性氣候模式(SCMs)；(2)耦合氣候系統(tǒng)模式(CGCMs)；(3)復(fù)雜程度介于前2類之間的中等復(fù)雜程度的地球系統(tǒng)模式(EMICs)。

相對于CGCMs來說,概念模型、能量平衡模式、輻射對流模式和統(tǒng)計(jì)動力模式屬于SCMs,這些模式對氣候系統(tǒng)高度整合,強(qiáng)調(diào)氣候系統(tǒng)的主要?dú)夂蛱卣骱臀锢磉^程,空間分辨率較低,對計(jì)算機(jī)要求不高,能進(jìn)行長尺度的模擬,故在地球模擬研究中得到廣泛應(yīng)用,在研究氣候變化成因上發(fā)揮了重要的作用。在20世紀(jì)90年代末期,中等復(fù)雜程度的EMICs發(fā)展起來的,涵蓋了CGCMs中的大多數(shù)過程,不過采用相對CGCMs更加簡化的方法來表述,這里把這兩類模式統(tǒng)稱為簡單地球模型。本文以簡單概念性氣候模式、中等復(fù)雜程度的地球系統(tǒng)模式為對象,介紹它們的應(yīng)用以及存在的不足。

2 簡單概念性氣候模式簡介

簡單概念性氣候模式中的能量平衡模式最先由前蘇聯(lián)的Budyko提出,美國的Sellers在1969年發(fā)表了相類似的一維能量平衡模式,計(jì)算出的氣候要素隨緯度的分布與實(shí)測值吻合得很好。20世紀(jì)70年代是動力氣候研究的春天,各種類型的模式相繼出現(xiàn),如輻射對流模式、統(tǒng)計(jì)動力模式、概念模型等。下文將對部分簡單概念性氣候模式進(jìn)行詳細(xì)介紹。

2.1 能量平衡模式

2.1.1 零維能量平衡模式

公式(1)未考慮大氣中的輻射過程,如果考慮大氣成分(CO2和水汽等)不僅吸收輻射能,而且發(fā)射長波輻射,使大氣氣溫增加,地球表面的氣溫Ts應(yīng)該比有效輻射溫度高。即：

如果進(jìn)入大氣的輻射和射出大氣的輻射量不平衡,那么地球氣溫將發(fā)生變化,變化率?Ts/?t可以表示為：

式中,c為地球的比熱容,τa為大氣紅外傳輸率,ε為比例參數(shù)。

2.1.2 一維能量平衡模式

一維模式不再將地球視為一點(diǎn),引入地球氣候系統(tǒng)的主要差異,即緯度變化特征,某一緯度帶(φ)能量平衡公式為：

式中,F(Tφ)是附加項(xiàng),為射出大氣的輻射量,Sφ和Tφ都隨緯度變化。顯然,一維能量平衡模式是零維能量平衡模式的擴(kuò)展。

將模式中附件項(xiàng)和射出大氣的輻射量經(jīng)過近似和簡化,得到如下公式：

再利用逐次迭代法,計(jì)算出各個(gè)緯度帶的氣溫。

除了以上兩種模式外,還可以引入二維(水平或者經(jīng)向/垂直)的能量平衡模式,另外,還有盒式能量平衡模式(Box-EBM),Box-EBM將大氣作為一層,陸面和海洋混合層為一層,海洋中間層和深海為另外兩層,和其他EBM模式相比,Box-EBM強(qiáng)調(diào)深海內(nèi)部的熱量過程以及不同深度層上的能量交換。

2.2 輻射對流模式

輻射對流模式[7](RCM)把大氣簡化為一個(gè)鉛直的大氣柱,詳細(xì)考慮大氣柱內(nèi)的輻射過程,通過對流調(diào)整來獲取大氣的垂直溫度分布,是一種考慮時(shí)間變化的一維模式,此模式建立在以下兩個(gè)原理上：(1)在任何高度上的太陽輻射和長波輻射通量與對流熱通量保持平衡；(2)因輻射差異引起的溫度垂直分布的不穩(wěn)定由對流調(diào)整而達(dá)到平衡。

2.2.1 無對流調(diào)整的輻射平衡模式

假定大氣氣溫的垂直分布是由輻射收支決定的,則溫度變化方程為：

式中,右邊第一項(xiàng)、第二項(xiàng)分別代表長波輻射和短波輻射引起的溫度變化。

2.2.2 有對流調(diào)整的輻射平衡模式

基本假定：(a)在大氣層頂,凈入射短波輻射等于射出長波輻射；(b)大氣的凈輻射冷卻作用等于大氣長波輻射與短波輻射之差；(c)溫度垂直遞減率小于規(guī)定值時(shí),氣層維持局地輻射平衡；(d)當(dāng)氣溫直減率大于規(guī)定值ε時(shí),對氣溫分布進(jìn)行調(diào)整,以達(dá)到規(guī)定值。

利用雙極晶體管在中子輻照下的硅位移損傷特性，成功研制了晶體管直流增益在線測試系統(tǒng)。該系統(tǒng)以LabVIEW為開發(fā)平臺，通過多任務(wù)并行處理的模塊化程序設(shè)計(jì)，實(shí)現(xiàn)了快速響應(yīng)、高效率同步采集、實(shí)時(shí)數(shù)據(jù)處理、存儲及顯示；采用回讀技術(shù)建立了電壓補(bǔ)償方法，解決了遠(yuǎn)程監(jiān)測中因電壓損耗帶來的測量誤差，系統(tǒng)的測試精度高達(dá)0.2%。

模式方程為：

式中,右邊第二項(xiàng)為對流調(diào)整項(xiàng),當(dāng)不考慮對流調(diào)整時(shí)該項(xiàng)為零。在對流層質(zhì)量守恒條件下有：

在地表面應(yīng)滿足熱量平衡方程,即：

同時(shí)可定義各層的短波輻射加熱率和長波輻射冷卻率分別為：

2.3 統(tǒng)計(jì)動力模式

將大氣沿著緯度圈進(jìn)行平均,用緯度和高度組成的網(wǎng)格點(diǎn)來表示大氣,就構(gòu)成了緯向平均動力模式(ZADM),該模式包括基本的動力和物理過程,是介于一維大氣模式和三維大氣模式之間,連接一維和三維模式的橋梁,因此這類模式在氣候模擬的研究中起著重要的作用。緯向平均模式的主要困難是渦旋輸送的參數(shù)化。由于渦旋輸送的處理是建立在統(tǒng)計(jì)近似的基礎(chǔ)上,所以這類模式也常稱為統(tǒng)計(jì)動力模式[6]。

把渦旋動量通量寫成：

式中,KT和Km分別為熱量和動量傳輸系數(shù),可以視為常數(shù)或?yàn)闇囟忍荻鹊暮瘮?shù),緯向平均量用表示。渦旋輸送的重要性實(shí)際上在于通過渦旋熱通量和渦旋動量通量將大氣積蓄的斜壓和正壓不能穩(wěn)定能量釋放出來,這必然對大氣環(huán)流和氣候變化發(fā)生重要的作用。

2.4 各模式之間的關(guān)系

通過圖1來表示簡化氣候模式間及其與CGCMs之間的關(guān)系,圖中三角錐的3條邊代表模式的3個(gè)基本要素,即輻射、地面過程和動力學(xué)。垂直方向表示模式的復(fù)雜程度,越接近椎體頂端表示復(fù)雜程度越大。當(dāng)?shù)竭_(dá)頂端時(shí),模式已有了足夠的分辨率,所有因素都被正確而合理地考慮到模式里去了,達(dá)到一個(gè)逼真的氣候模式。在椎體的底部是最簡單的氣候模式,它僅僅考慮一個(gè)主要的過程。同時(shí),椎體底部被認(rèn)為是空心的,因?yàn)樵谧刁w的各邊之間基本上沒有相互作用。沿著金字塔越往上,各主要過程之間的相互作用越大。

3 簡單概念性氣候模式的應(yīng)用現(xiàn)狀

圖1 氣侯模式的三角錐[8]（EBM：能量平衡模式；RCM：對流輻射模式；SDM：統(tǒng)計(jì)動力模式；CGCMs：耦合氣候系統(tǒng)模式）

簡單概念性氣候模式在氣候模擬研究中應(yīng)用廣泛,在研究氣候變化成因上發(fā)揮了重要的作用。從第一、二次IPCC報(bào)告十幾個(gè)氣候模式參加評估,到第三次報(bào)告34個(gè)模式被引用,再到第四次和第五次評估報(bào)告更多模式被采用,IPCC歷次報(bào)告所引用模式的數(shù)量在不斷增加,其中,簡單氣候模式(SCM)一直未缺席(IPCC第四次和第五次報(bào)告采用EBMs簡單氣候模式的結(jié)果)。利用簡單氣候模式能夠估測由于熱膨脹導(dǎo)致的全球平均溫度的變化和海平面的上升,此類模式以一套全球或者半球框的形式代表海洋—大氣系統(tǒng),進(jìn)而利用能量平衡方程預(yù)測全球地表溫度,能夠與地球生物化學(xué)圈簡化模式耦合,而可以快速估測氣候?qū)Ω鞣N氣候情景的響應(yīng)[9-10]。Gilliland[11]用零維能量平衡模式,考慮在太陽輻射變化、溫室氣體以及火山氣溶膠強(qiáng)迫情況下,得出太陽活動變化、CO2變化可以解釋北半球氣候變化中低頻變化的93%左右。Schlesinger等[12]也用類似的模式模擬得到溫室氣體、太陽活動和平流氣溶膠三者可以解釋1886—1985年全球氣溫變化的60%。畢鳴等[13]用海氣耦合的能量平衡模式研究了130年來全球氣溫的變化,考慮溫室氣體、太陽活動、火山活動及ENSO的情況下,模擬結(jié)

果與觀測序列的相關(guān)系數(shù)達(dá)0.884,可解釋77.8%的溫度變化。Crowley[14]利用EBM開展的研究指出,太陽輻射和火山活動可以解釋工業(yè)化以前北半球氣溫低頻變化的41%～64%,20世紀(jì)變暖主要來自溫室氣體的增加,自然變率的貢獻(xiàn)只占25%。Zhao等[15]用一維能量平衡模式研究下墊面改變對降水和溫度造成的影響,取得一定的成功。王紹武等[16]利用簡化的能量平衡模式,綜合考慮太陽、火山、溫室氣體的輻射效應(yīng),發(fā)現(xiàn)EBM能夠較好模擬近千年的全球平均氣溫序列,能夠成功再現(xiàn)中世紀(jì)暖期、小冰期和20世紀(jì)變暖等特征。Bertrand等[17]利用二維緯向平均模式所開展的研究,同樣表明太陽和火山活動可以解釋小冰期和中世紀(jì)暖期,而20世紀(jì)后期的加速變暖,則只有溫室氣體可以解釋；1850年以后氣溫的合理模擬,還需要考慮人為硫酸鹽氣溶膠強(qiáng)迫和森林采伐的影響。劉玉芝等[18]利用一維輻射對流氣候模式,詳細(xì)研究了云量、云光學(xué)厚度以及云高等要素的變化對大氣頂和地面太陽短波輻射和紅外長波輻射通量以及云的輻射強(qiáng)迫的影響。丁仲禮等[19]就晚更新世東亞古季風(fēng)變化方面建立起一個(gè)概念模型,模型考慮了太陽輻射變化、中高緯、青藏高原、邊緣海以及熱帶海洋等因素,得出東亞古季風(fēng)主要由全球冰量通過極地冰蓋擴(kuò)張、邊緣海變化以及青藏高原冰川作用3個(gè)過程來影響的,太陽輻射變化的直接控制作用相對較小,較合理地解釋了黃土所記錄的古季風(fēng)變化特征。

從上述模擬來看,簡單概念性氣候模式涵蓋了主要?dú)夂蛱卣骱椭饕锢磉^程,如熱力動力過程、輻射對流過程或關(guān)鍵反饋過程等,同時(shí)模式對計(jì)算機(jī)要求很低,能夠進(jìn)行長時(shí)間的模擬試驗(yàn)。不過,簡單概念性氣候模式通常采用緯向平均,最多二維,只能反映大區(qū)域的一個(gè)平均過程,空間分辨率不高。隨著當(dāng)前氣候變化動力機(jī)制逐漸明確、考慮更多的物理過程的情況下,這種模式存在一定的局限性,體現(xiàn)在難以給出氣候變化的高分辨率空間特征,同時(shí)還不包含氣候系統(tǒng)內(nèi)部變率過程。如果僅考慮“簡單”背景下的研究,其研究結(jié)果具有一定的借鑒意義。

4 中等復(fù)雜程度的地球系統(tǒng)模式簡介

中等復(fù)雜程度的地球系統(tǒng)模式(EMICs)于20世紀(jì)90年代末期建立,它一般具有二維(兩個(gè)水平維向或垂直和一個(gè)水平維向)空間分辨率,能夠描述包括人類和自然相互作用的自然地球系統(tǒng),EMICs涵蓋了CGCMs中的大多數(shù)過程,不過采用相對CGCMs更加簡化的方法去表述,與CGCMs模擬的分辨率也有差別。EMICs的復(fù)雜程度介于簡單概念性氣候模式和CGCMs模式之間,整合化程度沒有簡單概念性氣候模式高,但比CGCMs要高,所包含的過程和描述較CGCMs更加簡化,過程比簡單概念性氣候模式要多,細(xì)節(jié)描述也較之更加豐富,兼顧了簡單概念性氣候模式和CGCMs的部分優(yōu)點(diǎn)(圖2)。不難看出,相對于CGCMs模式,EMICs也是一種簡化氣候模式。

EMICs中的大氣模塊主要采用的是能量水汽平衡模式、準(zhǔn)地轉(zhuǎn)模式、統(tǒng)計(jì)動力學(xué)模式和大氣環(huán)流模式。大氣模塊的分辨率通常比較低,大多采用緯向平均和能量水氣平衡模式。由表1可見,EMICs在大氣模塊構(gòu)成和簡化程度方面存在明顯的差異,根據(jù)這種差異,模式可以分為兩類：一類是基于CGCMs進(jìn)行大氣或者海洋過程的簡化,使其總體的復(fù)雜程度較高,對氣候過程及反饋描述比較完善,對氣候過程描述比較細(xì)膩,如德國漢堡大學(xué)Planet Simulator模式[21-22]；而另一類基于低分辨率且大量參數(shù)化的模式,這類EMICs在更長時(shí)間尺度(千年范圍乃至冰期間冰期循環(huán))氣候變化的研究中有著重要作用,如CLIMBER-2的大氣模塊是二維統(tǒng)計(jì)動力學(xué)模式。從模式的復(fù)雜程度上看,前者更像是簡化了的CGCMs,而后者則更接近于簡單概念性氣候模式。綜述所述,從整體上簡化了EMICs的復(fù)雜程度,降低了對計(jì)算能力的要求,有助于在古氣候資料的強(qiáng)迫下模擬出不同于今天的長尺度古氣候狀態(tài),充分顯示出了EMICs模式的優(yōu)點(diǎn)；但EMICs的分辨率比較粗糙,對地球系統(tǒng)細(xì)節(jié)的描述也相當(dāng)簡化,因而模擬無法達(dá)到CGCMs的高度。

5 中等復(fù)雜程度的地球系統(tǒng)模式的應(yīng)用現(xiàn)狀

當(dāng)前EMICs在氣候研究中得到廣泛的應(yīng)用,雖然其對氣候系統(tǒng)的描繪不如CGCMs細(xì)致,但是EMICs具有與簡單概念性模式相當(dāng)?shù)挠?jì)算效率,所以EMICs主要注重長期氣候模擬,以足夠的細(xì)節(jié)來刻畫地球系統(tǒng)的長期演化。

圖2 EMICs 的圖解定義[20]

表1 主要EMICs及相關(guān)介紹[5]

目前,國際古氣候模擬比較計(jì)劃(PMIP)已經(jīng)歷了PMIP1(1991—2001年)和PMIP2(2002—2007年)兩個(gè)階段,并正在開展第三階段PMIP3(2008年至今)的研究。PMIP3階段致力于利用地球系統(tǒng)模式對末次盛冰期和中全新世(距今最近的、最有代表性的寒冷和溫暖的氣候狀態(tài))的氣候進(jìn)行模擬。在古氣候模擬的研究方面,產(chǎn)生了不少研究成果。Gallcee等[23]和Berger等[24]的研究表明,氣候模式LLN-2D在變化的太陽輻射強(qiáng)迫下,當(dāng)指定的常量CO2濃度低于220×10-6后,能夠重現(xiàn)觀測到的冰期—間冰期循環(huán)；而固定軌道參數(shù),讓CO2依據(jù)Vostok重建記錄變化,則LLN-2D無法再現(xiàn)冰期—間冰期循環(huán)[25]。EMICs對冰期的開始和結(jié)束也進(jìn)行了各種試驗(yàn),Yoshimori等[26]利用Uvic模式進(jìn)行了末次冰期結(jié)束對軌道參數(shù)和CO2強(qiáng)迫的敏感性試驗(yàn),Meissner等[27]使用UVic研究了植被對末次冰期開始的作用,Kageyama等[28]使用CLIMBER-2探討了末次冰期的開始特別是各種不同冰川的反饋?zhàn)饔?。Wang等[29]使用MPM-2對冰期氣候突變和D-O震蕩的機(jī)制進(jìn)行了深入探討,發(fā)現(xiàn)北大西洋經(jīng)向翻轉(zhuǎn)流(AMOC)狀態(tài)的快速轉(zhuǎn)換導(dǎo)致冰期氣候突變,這一突變被海冰面積的擴(kuò)張進(jìn)一步強(qiáng)化。

另外,對歷史、當(dāng)今、未來氣候模擬的應(yīng)用研究,在IPCC AR4中,EMIC已較前幾次評估報(bào)告被更多人使用并得到深入的評估,在海平面上升歸因分析上,利用一個(gè)低分辨率CGCM 和數(shù)個(gè)以氣候敏感性最佳估計(jì)值3℃為基礎(chǔ)的EMICs進(jìn)行推估,得出海平面上升的促成因素是海洋熱膨脹,并對未來進(jìn)行了預(yù)估[9]。在2013年出版的IPCC AR5中,EMICs在古氣候模擬、未來氣候預(yù)估等方面發(fā)揮了重要的作用[10]。Loutre等[30]使用LLN-2D對下次冰期的開始也進(jìn)行了探討,研究了未來的氣候及其對不同CO2情景的敏感性,發(fā)現(xiàn)目前的氣候很可能經(jīng)歷較長的間冰期(大約5萬年),而在未來的6萬年和10萬年將可能出現(xiàn)較小和較大的冰期。Bertrand等[31]和Bauer等[32]分別使用MoBiDiC和CLIMBER-2對過去1000年的各種氣候強(qiáng)迫年進(jìn)行了評估,與高分辨率的重建數(shù)據(jù)進(jìn)行比較,發(fā)現(xiàn)火山和太陽輻射作為自然強(qiáng)迫,是工業(yè)革命時(shí)期以前的主要?dú)夂驈?qiáng)迫,能夠很好地解釋中世紀(jì)暖期和小冰期。但自然強(qiáng)迫不足以產(chǎn)生近代以來的全球變暖,人類活動排放的溫室氣體是造成全球變暖的主要強(qiáng)迫,而且研究結(jié)果表明所有的氣候強(qiáng)迫對模擬過去千年氣候變化都是重要的。Petoukhov等[33]對目前氣候和CO2加倍模擬的EMICs進(jìn)行比較,充分顯示EMICs對目前氣候的模擬非常接近觀測結(jié)果,而且EMICs模擬誤差均在CGCMs誤差的范圍之內(nèi)。Wang等[34]和Marsh等[35]分別使用MPM-2和GENIE對溫鹽環(huán)流進(jìn)行研究,發(fā)現(xiàn)溫鹽環(huán)流這個(gè)洋流系統(tǒng)最少有兩種穩(wěn)定的氣候狀態(tài),其中一種是目前所出現(xiàn)的,即環(huán)流非常強(qiáng)勁和向極地方向傳輸大量熱能,從而帶來溫暖氣候的狀態(tài)；另一種則是環(huán)流非常微弱,甚至接近于停止的狀態(tài),從而帶來寒冷的氣候狀態(tài)。數(shù)值模擬試驗(yàn)還顯示,這兩種狀態(tài)之間是可能出現(xiàn)相互轉(zhuǎn)化的。另外,利用EMICs對陸地、海洋、全球碳循環(huán)等多個(gè)方面的問題也進(jìn)行了探討。

結(jié)合當(dāng)前EMICs的應(yīng)用和發(fā)展,可見EMICs應(yīng)用范圍很廣,模擬能力很強(qiáng),未來如在下面幾個(gè)方面取得突破,必定會有更大的發(fā)展,從而在地球系統(tǒng)模擬中起到更加重要的作用：(1)加入人類活動、考慮固體地球過程和太陽活動的影響,使EMICs成為真正意義上的地球模式；(2)注重氣候突變機(jī)理的研究,區(qū)分自然和人為因素在突變中的相對貢獻(xiàn)；(3)改進(jìn)模式算法,提高硬件設(shè)備的運(yùn)算能力,使其耗費(fèi)更少的計(jì)算時(shí)間,能夠進(jìn)行更長時(shí)間尺度上的模擬研究。

6 結(jié)論

本文介紹了以簡單概念性氣候模式(SCMs)和中等復(fù)雜程度的地球系統(tǒng)模式(EMICs)為代表的簡單地球模型,它們在氣候模擬方面都有不凡表現(xiàn)。簡單概念性氣候模式將對氣候變化的研究拓展到了古氣

候領(lǐng)域,已經(jīng)成為了解氣候理論和氣候變化規(guī)律的有效工具。EMICs涵蓋了地球系統(tǒng)的大多數(shù)組成部分,對各部分之間的反饋和過程描述比較詳細(xì),同時(shí)計(jì)算要求比較低,使其應(yīng)用幾乎覆蓋了SCMs以及耦合氣候系統(tǒng)模式(CGCMs)的所有研究領(lǐng)域。EMICs必定會有更大的發(fā)展,從而在地球氣候模擬中起到更加重要的作用。

EMICs也有其自身的缺點(diǎn),分辨率比較低,對地球系統(tǒng)細(xì)節(jié)的描述也相當(dāng)?shù)暮喕?尚不能達(dá)到CGCMs的高度。但是,EMICs能夠刻畫出氣候變化的大體特征,更能進(jìn)行CGCMs在目前計(jì)算條件下無法比擬的長期氣候變化模擬。2005年,德國漢堡大學(xué)完成了Planet Simulator模式,此模式屬于典型的EMICs模式,由全球大氣環(huán)流模式、海洋—海冰模塊和陸地土壤—生物圈模塊組成,以并行計(jì)算機(jī)為硬件支撐,支持地球和類地球行星的氣候動力學(xué)數(shù)值試驗(yàn),同時(shí)還具備結(jié)構(gòu)簡單、可視化界面和開放源代碼便于改造等特點(diǎn),依托此模式使改進(jìn)模式、模擬長尺度的古氣候和開展敏感性試驗(yàn)成為了可能。

EMICs從出現(xiàn)到發(fā)展到今天不過短短的十幾年,尚處于發(fā)展的初始階段,隨著計(jì)算機(jī)硬件水平的提高、模式的不斷改進(jìn)和對氣候系統(tǒng)認(rèn)識的不斷加深,在未來的幾十年或者更長的時(shí)期內(nèi),對氣候時(shí)間尺度(千年甚至更長)的氣候變化的研究,EMICs很可能依然是比較理想的選擇。今后一段時(shí)間,如何利用EMICs就一些重要驅(qū)動因子(如天文和地球因子)進(jìn)行敏感性試驗(yàn),如何采用混沌理論研究其非線性問題[36],如何改進(jìn)對重點(diǎn)區(qū)域和重點(diǎn)過程的模擬,將是氣候數(shù)值模擬領(lǐng)域面臨的重要課題。

[1]王斌, 周天軍, 俞永強(qiáng), 等. 地球系統(tǒng)模式發(fā)展展望. 氣象學(xué)報(bào), 2008, 66(6): 857-869.

[2]《10000個(gè)科學(xué)難題》地球科學(xué)編委會. 10000個(gè)科學(xué)難題·地球科學(xué)卷. 北京: 科學(xué)出版社, 2010.

[3]周天軍, 王在志, 宇如聰, 等. 基于LASG/IAP大氣環(huán)流譜模式的氣候系統(tǒng)模式. 氣象學(xué)報(bào), 2005, 63(5): 702-715.

[4]Houghton J T, Ding Y, Griggs D J, et al. Climate change 2001: the scientific basis. United Kingdom and New York: Cambridge University Press, 2001.

[5]張冉, 李力, 郭慶春, 等. 古氣候研究中氣候模式的發(fā)展與應(yīng)用. 干旱區(qū)研究, 2007, 24(5): 704-711.

[6]李崇銀. 氣候動力學(xué)引論(第二版). 北京: 氣象出版社, 2000.

[7]于革. 古氣候動力模擬. 北京: 高等教育出版社, 2007.

[8]曹鴻興. 氣候模式與氣候模擬. 氣象, 1987, 13(11): 3-6.

[9]IPCC. IPCC fourth assessment report (AR4). Cambridge: Cambridge University Press, 2007.

[10]IPCC. IPCC fifth assessment report (AR5). Cambridge: Cambridge University Press, 2013.

[11]Gilliland R L. Solar, volcanic and CO2forcing of recent climatic changes. Climate Change, 1982, 4: 11-13.

[12]Schlesinger M E, Ramalakutty N. Implication for globalwarming of intercycle solar irradiance variations. Nature, 1992, 360: 330-333.

[13]畢鳴, 王紹武, 李曉東. 130年來全球氣溫變化的模擬. 第四紀(jì)研究, 1997(1): 21-27.

[14]Crowley T J. Causes of climate change over the past 1000 years. Science, 2000, 289: 270-277.

[15]Zhao M, Zeng X. A theoretical analysis on the local climate change induced by the change of landuse. Adv Atm Sci, 2002, 19: 45-63.

[16]王紹武, 謝志輝, 蔡靜寧, 等. 近千年全球平均氣溫變化的研究.自然科學(xué)進(jìn)展, 2002, 12(11): 1145-1149.

[17]Bertrand C, Loutre M F, Crucifix M, et al. Climate of the last millennium: a sensitivity study. Tellus, 2002, 54A: 221-244.

[18]劉玉芝, 石廣玉, 趙劍琦. 一維輻射對流模式對云—輻射強(qiáng)迫的數(shù)值模擬研究.大氣科學(xué), 2007, 21(3): 486-493.

[19]丁仲禮, 劉東生. 晚更新世東亞古季風(fēng)變化動力機(jī)制的概念模型. 科學(xué)通報(bào), 1998, 43(2): 122-132.

[20]Claussen M. Earth system models //Ehlers E, Krafit T, eds. Understanding the Earth system: compartments, processes and interactions. Heidelberg, Berlin, New York: Springer, 2000.

[21]Fraedrich K, Jansen H, Kirk E, et al. The planet simulator: Green planet and desert world. Meteorologische Zeitschrift, 2005, 14(3): 305-314.

[22]Fraedrich K, Jansen H, Kirk E, et al. The planet simulator: towards a user friendly model. Meteorologische Zeitschrift, 2005, 14(3): 299-304.

[23]Gallcee H, Ypersele J P, Fichefet T, et al. Simulation of the last glacial cycle by a coupled, sectorially averaged climate-ice sheet model II. Response to insolation and CO2variation. Journal of Geophysical Research, 1992, 97(15):713-715.

[24]Berger A, Loutre M F, Gallcee H. Sensitivity of the LLN climate model to the astronomical and CO2forcings over the last 200kyr. Climate Dynamics, 1998, 14: 615-629.

[25]Loutre M F, Berger A. No glacial-interglacial cycle in the ice volume simulated under a constant astronomical forcing and avariable CO2. Geophysical Research Letters, 2000, 27(6): 783-786.

[26]Yoshimori M, Weaver A J, MarshallS J, et al. Glacial terminations: sensitivity to orbital and CO2forcing in a coupled climate system model. Climate Dynamics, 2001, 17: 571-588.

[27]Meissner K J, Weaver A J, Matthews H D, et al. The role of land surface dynamics in glacial inception: a study with the Uvic Earth System Model. Climate Dynamics, 2003, 21: 515-537.

[28]Kageyama M, Charbit S, Ritz C,et al. Quantifying ice-sheet feedbacks during the last glacial inception. Geophysical Research Letters, 2004, 31, L24203, do:i 10. 1029/2004GL021339.

[29]Wang Z, Mysak L A. Glacial abrupt climate changes and Dansgaard-Oeschger oscillations in a coupled climatemodel. Paleoceanography, 2006, 21: PA2001.

[30]Loutre M F, Berger A. Future climatic changes: Are we entering an exceptionally long Interglacial. Climatic Change, 2000, 46: 61-90.

[31]Bertrand C, Loutre M F, Crucifix M, et al. Climate of the last millennium: a sensitivity study. Tellus, 2002, 54A: 221-244.

[32]Bauer E, Claussen M, Brovkin V. Assessing climate forcings of the Earth system for the past millennium. Geophysical Research Letters, 2003, 30(6): 1276-1280.

[33]Petoukhov V, Claussen M, Berger A, et al. EMIC intercomparison project (EMIP-CO2): Comparative analysis of EMIC simulations of currentclimate and equilibrium and transient responses to atmospheric CO2doubling. Climate Dynamics, 2005, 25: 363-385.

[34]Wang Z. Two climatic states and feedbacks on thermohaline circulation in an Earth system model of intermediate complexity. Climate Dynamics, 2005, 25: 299-314.

[35]Marsh R, Yool A, Lenton T M, et al. Bistability of the thermohaline circulation identified through comprehensive 2-parameter sweeps of an efficient climate model. Climate Dynamics, 2004, 23: 761-777.

[36]Shao Y, Peng G, Leslie L M. The environmental dynamic system //Peng G, Leslie LM, Shao Y. Environmental modelling and prediction. Berlin: Springer-Verlag, 2002.

Progress in Simple Earth Models

Wang Yong1,2, Liu Suxia1, Shao Yaping3, Peng Gongbing1,4
(1 Key Laboratory of Water Cycle and Related Land Surface Processes, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101; 2 School of Remote Sensing, Nanjing University of Information Science and Technology, Nanijing 210044; 3 Institute for Geophysics and Meteorology, University of Cologne, Cologne 50923; 4 China Science Center, International Eurasian Academy of Sciences, Beijing 100107)

In this paper, by taking Simple conceptual Climate Models (SCMs) and Earth system Models of Intermediate Complexity (EMICs) as research objects, simple climate models and their application in climate modeling are introduced. Simple conceptual climate models have become effective tools to understand the theory of climate and climate-change law. Earth system Models of Intermediate Complexity (EMICs) cover most of the components in the Earth system, and describe the process and feedback among the various components in more details. Because EMICs have low computational requirements, modeling applications of EMICs cover almost all fi elds of simple conceptual climate models and CGCMs simulations. EMICs will surely achieve a greater development and play a more important role in the Earth's climate simulation.

sample earth model, simple conceptual climate models, EMICs

10.3969/j.issn.2095-1973.2014.03.004

2013年12月17日；

2014年4月14日

王勇（1979—），Email：wyong0210@igsnrr.ac.cn

資助信息：國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目（2012CB957802）