二流—四流球諧函數(shù)譜展開累加輻射傳輸方案在全球氣候模式中的應(yīng)用

張華 盧鵬 荊現(xiàn)文

?

二流—四流球諧函數(shù)譜展開累加輻射傳輸方案在全球氣候模式中的應(yīng)用

張華1, 2盧鵬3, 1荊現(xiàn)文1, 2

1中國氣象局氣候研究開放實(shí)驗(yàn)室/國家氣候中心，北京100081;2南京信息工程大學(xué)氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心，南京210044;3江蘇省氣候中心，南京210008

本文首先構(gòu)建了二流—四流球諧函數(shù)譜展開累加輻射傳輸?shù)男路桨?，然后將其?yīng)用于國家氣候中心第二代大氣環(huán)流模式BCC_AGCM2.0.1的新版本中，并與模式中原有的Eddington累加方案進(jìn)行了比較。由于新方案本質(zhì)上是單層Eddington近似方案在四流上的推廣。因此新方案在計算精度上要優(yōu)于原方案。通過在全球氣候模式中的應(yīng)用與比較，本文發(fā)現(xiàn)新方案對氣候模擬會產(chǎn)生比較大的影響。在晴空條件下，新方案計算的在南緯30°到60°區(qū)間、北大西洋東北部以及非洲北部的撒哈拉沙漠區(qū)域的地表向下年平均短波輻射通量要小于原方案結(jié)果，最大差別可以達(dá)到3.5 W/m2；同時，新方案計算的在南緯30°到60°區(qū)間和北大西洋東北部的大氣頂向上年平均短波輻射通量要大于原方案結(jié)果，最大差別達(dá)到3 W/m2。在有云大氣情況下，新方案計算的地表向下年平均短波輻射通量要小于原方案結(jié)果，并隨著緯度的增加，新舊兩種方案的差別逐漸變大，在南北極時達(dá)到最大5.5 W/m2；同時，新方案計算的在赤道區(qū)域的大氣頂?shù)哪昶骄滩ㄏ蛏陷椛渫恳∮谠桨附Y(jié)果，最大差別為2.5 W/m2，而在南北緯30°到60°區(qū)間，新方案計算的在大氣頂?shù)哪昶骄滩ㄏ蛏陷椛渫縿t要大于原方案結(jié)果，最大差別為1.5 W/m2。新方案計算的年平均短波加熱率普遍高于原方案結(jié)果，特別是在800 hPa到地表之間的低層大氣以及50 hPa到100 hPa的高層大氣，最大差別可達(dá)0.03 K/d。因此，新方案有助于改善全球氣候模式中普遍存在的赤道平流層中下層的溫度冷偏差現(xiàn)象。

四流球諧函數(shù) 輻射傳輸方案 輻射通量 加熱率 冷偏差

1 引言

輻射過程是大氣中最重要的物理過程之一。一方面，地氣系統(tǒng)的輻射收支決定了長期的全球平均氣候狀況；另一方面，太陽輻射能和熱輻射能在大氣中的輻散、輻合，形成非絕熱冷、熱源，造成大氣水平和垂直層結(jié)的不穩(wěn)定（石廣玉，2007）。在具體的大氣輻射模式中，前者通過大氣頂和地表的輻射通量表征；后者通過大氣加熱率反映。

氣候模式中的輻射模式比較計劃（簡稱：ICRCCM，the Intercomparison of Radiation Codes in Climate Models）第一階段第二部分（Fouquart et al., 1991）比較了26個短波輻射模式。其中在晴空條件下共有21個大氣輻射模式參加了比較（不包括逐線積分模式），結(jié)果表明：各模式間的均方根差 別在4%左右。在有云大氣條件下，均方根的差 別隨著云層光學(xué)厚度的變化而變化，從4%到10%。ICRCCM第三階段（Barker et al., 2003）比較了25個短波輻射模式，其中包括2個逐線積分模式、4個蒙特卡洛光子輸運(yùn)算法。在熱帶大氣廓線條件下，在太陽直射時，無論有云與否，大多數(shù)短波輻射模式都會低估15～25 W/m2的大氣吸收。Halthore et al.（2005）評估了16個短波輻射模式，結(jié)果表明，在晴空不考慮氣溶膠大氣情況下，各種短波輻射模式計算的地表直接輻射強(qiáng)度之間的均方根方差在1%以內(nèi)，但在濕大氣情況下，均方根方差達(dá)到5%左右。從上述比較可以看出，各種短波輻射算法很難達(dá)到一致，短波輻射模式的精度還需要進(jìn)一步的提高。

目前大多數(shù)氣候模式仍然采用二流近似輻射傳輸方案。該方案的優(yōu)點(diǎn)在于計算速度快，但是與四流近似輻射傳輸方案相比，該方案在有云大氣的情況下的計算誤差比較大。在美國標(biāo)準(zhǔn)大氣廓線情況下，二流近似輻射傳輸方案與48流DISORT離散縱坐標(biāo)的參考結(jié)果相比，云頂加熱率有可能被低估約6%（Zhang et al., 2013；Zhang and Li, 2013）。

Shibata and Uchiyama（1992）評估了四流矩陣算子法，當(dāng)太陽天頂角小于72.5°時，相對誤差小于5%；當(dāng)太陽天頂角大于72.5°時，相對誤差小于10%。Chou（1992）利用單層四流離散坐標(biāo)法結(jié)合累加法計算輻射通量和加熱率，與傳統(tǒng)的通過求解線性方程組的四流離散坐標(biāo)法相比，明顯地節(jié)省了計算時間，因此可以應(yīng)用于氣候模式中。Kay et al. （2001）分別評估了二流和四流離散縱坐標(biāo)法以及矩陣求逆法，結(jié)果表明在光化通量的計算中，四流離散縱坐標(biāo)法略優(yōu)于矩陣求逆法，而兩者都要優(yōu)于二流離散縱坐標(biāo)法。而在輻射通量的比較中，四流離散縱坐標(biāo)法要優(yōu)于矩陣求逆法和二流離散縱坐標(biāo)法。Ayash et al.（2008）在第三代加拿大氣候中心大氣環(huán)流模式（Scinocca et al., 2008）（簡稱：CCC GCM3，the third-generation Canadian Climate Center Atmosphere GCM）中評估了二流—四流混合離散縱坐標(biāo)法（Liou et al., 1988）與二流近似輻射傳輸算法對輻射通量的影響。Liu et al.（2009）將Fu-Liou四流離散縱坐標(biāo)輻射傳輸算法加入美國海軍海氣耦合中尺度預(yù)報系統(tǒng)（簡稱：COAMPS，Coupled Ocean–Atmosphere Mesoscale Prediction System）中，發(fā)現(xiàn)對三天短期天氣預(yù)報的溫度場有很大的改進(jìn)。

Li and Ramaswamy（1996）通過比較表明，對于求解單層輻射傳輸方程而言，四流球諧函數(shù)輻射傳輸方法的透過率和反射率的精度都要略高于四流離散縱坐標(biāo)方法的結(jié)果。本文利用單層四流球諧函數(shù)譜展開代替原有的二流近似解來計算均勻的單層介質(zhì)的的透射率和反射率，其采用的累加過程還是原有的二流累加法（Coakley et al., 1983）。這種單層四流球諧函數(shù)譜展開和二流近似累加法結(jié)合的方法，即為二流—四流球諧函數(shù)譜展開累加法（下文中簡稱：新方案）。

本文將以上新方案應(yīng)用到國家氣候中心氣候模式BCC_AGCM2.0.1新的版本（荊現(xiàn)文和張華，2012；Zhang et al., 2014）中進(jìn)行檢驗(yàn)和評估，以改進(jìn)該模式對大氣輻射的計算精度。

本文的第二部分給出了在非均勻多層大氣中四流球諧函數(shù)輻射傳輸方案的構(gòu)建；第三部分簡要介紹了所用的全球氣候模式和數(shù)值檢驗(yàn)方法；第四部分比較分析了新方案與Eddington近似方案在氣候模式中的表現(xiàn)及其模擬差異；最后在第五部分給出了本文的主要結(jié)論。

2 四流球諧函數(shù)輻射傳輸方案構(gòu)建

輻射傳輸方程可以寫成以下形式：

其中，為局地天頂角，是輻射強(qiáng)度，是光學(xué)厚度，是多次散射項，0是單次散射項。四流球諧函數(shù)輻射傳輸算法將輻射強(qiáng)度展開成球諧函數(shù)如下：

將（2），（3），（4）式代入（1）整理，并將展開到第4項，可得常系數(shù)線性微分方程組。通過求解可得輻射強(qiáng)度的各展開項的數(shù)值。利用公式（2）可得到輻射強(qiáng)度，同時也可以得到透過率和反射率。詳細(xì)的推導(dǎo)過程也可見Li and Ramaswamy（1996）及Zhang and Li（2013）。

在求解方程時如果不考慮直接入射項，則可以通過求解齊次常系數(shù)線性微分方程組得到漫射的透過率和反射率。對于兩層總的透過率和反射率則可通過Coakley et al.（1983）的公式求得：

通過上述公式（5）到（8），可以得到任意2層的反射率和透過率，從而得到非均勻大氣每層的向上通量和向下通量。

本文在Li and Ramaswamy（1996）基礎(chǔ)上，在數(shù)值上實(shí)現(xiàn)了二流—四流球諧函數(shù)譜展開累加計算，并將之應(yīng)用于全球氣候模式中進(jìn)行比較和評估。

3 全球氣候模式介紹和數(shù)值檢驗(yàn)方法

本文采用國家氣候中心第二代大氣環(huán)流模式BCC_AGCM2.0.1（Wu et al., 2008，2010）的新版本。其水平分辨率為T42（近似于2.8度乘以2.8度），垂直方向采用混合坐標(biāo)，共26層。新版本的重要特征是更新了BCC_AGCM2.0.1的輻射過程和云的垂直重疊處理，新的輻射過程采用Zhang et al.（2003，2006a，2006b）發(fā)展的BCC-RAD輻射模塊（Randles et al., 2013；Zhang et al., 2014）。云的垂直重疊采用McICA蒙特卡洛隨機(jī)次網(wǎng)格柱方法處理（荊現(xiàn)文和張華，2012；Jing and Zhang，2013）。氣溶膠過程采用Gonget al.（2002，2003）開發(fā)的CUACE_Aero氣溶膠理化數(shù)值模式（Zhang et al., 2012）。Wuet al.（2010）對BCC_AGCM2.0.1模式的模擬性能進(jìn)行了系統(tǒng)的評估。

BCC-RAD輻射模塊對于氣體吸收采用相關(guān)分布算法（Shi，1981；張華，1999）；對于輻射傳輸方程的求解，本文采用了Eddington二流累加法和二流—四流球諧函數(shù)譜展開累加兩種輻射傳輸算法來計算大氣輻射通量和加熱/冷卻率。

本文通過離線診斷和在線模擬兩種方式比較Eddington近似和本文提出的新算法在氣候模式中的模擬效果。離線診斷是在運(yùn)行Eddington二流輻射傳輸算法的同時，加入新方案，但新方案的輸出量僅作為診斷量，不產(chǎn)生氣候反饋。離線診斷主要用于檢驗(yàn)在相同大氣條件下輻射傳輸算法本身的差別。在線模擬是在相同的初始場條件下分別運(yùn)行Eddington近似和新方案，以檢驗(yàn)兩種方法各自對氣候模擬的影響。首先離線診斷了兩種方法的輻射場差異，重點(diǎn)給出了地表向下的短波輻射通量、大氣頂向上短波輻射通量和短波輻射加熱率的結(jié)果，它們分別反映了以上不同輻射傳輸算法下大氣的透射率、反射率和大氣加熱情況。然后給出了兩種方法對氣候模擬的影響，重點(diǎn)分析了兩種輻射傳輸算法對云輻射強(qiáng)迫的影響。

離線診斷和在線運(yùn)行都是從1949年9月1日開始，積分52個月。其中輻射方案每小時調(diào)用一次（模式的積分步長為20分鐘，也就是3個模式積分步長調(diào)用一次）。前16個月為Spin-up時間，取后三年（1951～1953年）的結(jié)果進(jìn)行分析。所用海溫資料為多年平均的月平均氣候態(tài)數(shù)據(jù)（Hurrell and Trenberth, 1999）。

4 結(jié)果分析

4.1 晴空輻射通量

圖1a、b分別給出了兩種輻射傳輸算法得到 的晴空條件下地表向下和大氣頂向上年平均短波輻射通量差異。從圖1a可以看出，兩種方案在晴空大氣條件下地表向下短波輻射通量的差別較大的區(qū)域主要發(fā)生在南緯30°到60°區(qū)間的海洋表面以及非洲北部的撒哈拉沙漠，分別處于海鹽氣溶膠和沙塵氣溶膠含量較高的區(qū)域，這些區(qū)域Eddington近似方法都高估了地表向下的短波輻射通量。Li and Ramaswamy（1996）的結(jié)果表明二流近似與新方案相比，不論是粒子單次散射比為1還是0.9都高估了大氣的透射率，因此，在氣溶膠濃度高的區(qū)域，Eddington近似會高估對晴空地表向下的輻射通量。

從圖1b可以看出，Eddington近似與新方案相比，對大氣頂向上短波輻射通量的最大低估區(qū)主要位于南北緯30°到60°的海洋區(qū)域，但在非洲北部的撒哈拉沙漠卻僅出現(xiàn)微弱的低估和高估。Li and Ramaswamy（1996）表明：與Eddington近似相比，不論是單次散射比為1還是0.9，新方案的反射率在太陽天頂角余弦比較大時存在高估，而在太陽天頂角比較小時存在低估。而在南北緯30°到60°的海洋區(qū)域正好對應(yīng)太陽天頂角余弦值比較小的時候，因此大氣頂?shù)南蛏贤繒坏凸?；而對于撒哈拉沙漠區(qū)域，當(dāng)太陽在北半球時則對應(yīng)太陽天頂角余弦值比較大的情況，當(dāng)太陽在南半球時則對應(yīng)太陽天頂角余弦值比較小的情況，因此年平均的通量值正好被這種高估和低估相互抵消，僅出現(xiàn)很微弱的低估和高估。

4.2 有云大氣輻射通量

圖2與圖1相似，但是為有云條件的相應(yīng)結(jié)果。從圖2a中可以看出，兩種方案在有云大氣條件下地表向下短波輻射通量的差值從赤道向兩極逐漸變大，南北兩極的差別最大。由于Eddington近似與新方案相比，整體都是高估單層透過率（Li and Ramaswamy，1996），同時在相同光學(xué)厚度條件下，兩種算法計算的透過率差別隨著太陽天頂角的增大而增大。而太陽天頂角隨著緯度的增加而增加，正好與有云情況下地表向下短波輻射通量差值的分布相符。Ayash et al.（2008）的結(jié)果也表明，在有云大氣情況下，二流輻射傳輸算法與二流—四流離散縱坐標(biāo)累加算法相比到達(dá)地面的輻射通量偏多，特別是在高緯地區(qū)。

圖2b表明兩種方案在有云大氣條件下大氣頂向上短波輻射通量在熱帶區(qū)域存在正的差別，在南北半球中高緯度的海洋區(qū)域存在負(fù)的差別。這兩個區(qū)域分別對應(yīng)對流云量和層云量出現(xiàn)比較大的區(qū)域，云的光學(xué)厚度較大。由于當(dāng)太陽天頂角余弦較大的時候，Eddington二流近似算法在反射率上存在高估，而在太陽天頂角余弦值較小的時候，Eddington二流近似算法在反射率上存在低估（Li and Ramaswamy，1996），因此在高緯度區(qū)域，對應(yīng)太陽天頂角余弦值較小的情況，二流近似算法相對新方案低估了有云大氣的大氣頂向上短波輻射；在熱帶區(qū)域，對應(yīng)太陽天頂角余弦值較大的情況，二流近似算法相對新方案高估了有云大氣的大氣頂向上短波輻射。Ayash et al.（2008）比較了二流近似與二流—四流離散縱坐標(biāo)法也得到了相同的結(jié)果。

圖1 用Eddington近似與新方案計算的晴空年平均輻射通量（W/m2）差值的全球分布（Eddington近似減去新方案）：（a）地表向下短波輻射通量；（b）大氣頂向上短波輻射通量

圖2 同圖1，但為有云大氣條件

4.3 短波加熱率

圖3給出了Eddington方法與新方案計算的年平均短波加熱率的差別的緯度—高度分布?？梢钥闯觯瑑煞N方法造成的短波加熱率差異主要發(fā)生在800 hPa到地表的低層大氣以及50 hPa到100 hPa的高層大氣。對800 hPa到地表的區(qū)域，偏差最大超過0.2 K/d，這主要是受低云的影響。對于100 hPa以上的區(qū)域，Eddington近似方案相對新方案，加熱率也有一個負(fù)偏差。在100 hPa以上的區(qū)域，主要處于平流層，目前全球氣候模式對溫度場的模擬在熱帶上空的平流層中下層區(qū)域幾乎都存在一個冷偏差，這是全世界氣候模擬普遍存在的問題（Forster et al., 2011），而Eddington近似和新方案計算的加熱率差值表明，在該區(qū)域如果采用新方案計算，將會比Eddington近似計算的加熱率強(qiáng)，有助于改善氣候模式的冷偏差問題。

圖3 用Eddington方法與新方案計算的年平均短波加熱率差值（K/d）的緯度—高度分布（前者減去后者）

圖4給出Eddington近似方法、新方案與ECMWF資料溫度差值的緯度—高度分布圖。從圖4可以發(fā)現(xiàn)，采用新方案計算的溫度場結(jié)果在50～100 hPa赤道地區(qū)的的確要優(yōu)于Eddington近似方案，對原有模式中熱帶上空平流層的冷偏差有所改進(jìn)。而該區(qū)域大氣成分的分布及變化對于認(rèn)識氣候長期變化極為重要，因?yàn)樵搮^(qū)域的臭氧、水汽、卷云和氣溶膠對太陽短波輻射和地球長波輻射有很強(qiáng)的調(diào)節(jié)作用（陳洪濱等，2006）。

4.4 云輻射強(qiáng)迫

圖5給出Eddington近似方法和新方案模擬的短波區(qū)間云輻射強(qiáng)迫與CERES資料的差值圖，可以看出Eddington近似與CERES資料的全球平均值相比低估了―1.32 W/m2, 新方案與CERES資料的全球平均值相比低估了―0.33 W/m2。原因可能是由于新方案計算的加熱率要大于Eddington近似方案，尤其是在加熱率最大的云頂處，因此云頂?shù)募訜崧首兇?，將抑制云的發(fā)展從而減少云量，從而減少反射到大氣頂?shù)妮椛渫?，使得短波區(qū)間云的負(fù)輻射強(qiáng)迫變小、云輻射強(qiáng)迫與CERES資料的結(jié)果更加接近。

圖4 Eddington近似和新方案與ECMWF資料溫度場差值的緯度—高度分布：（a）Eddington近似方法與ECMWF資料的差值；（b）新方案與ECMWF資料的差值

圖5 Eddington近似方法和新方案與CERES資料大氣頂短波云輻射強(qiáng)迫的差值場：（a）Eddington近似方法與CERES資料的差值（a）新方案與CERES資料的差值

從全球分布來看，新方案與Eddington近似方法相比，在北太平洋中部和南緯45°到60°海域，對模式模擬的短波云頂輻射強(qiáng)迫的誤差有所減少。丁守國等（2005）利用ISCCP月平均云氣候資料研究了全球云量分布，從年平均的經(jīng)向平均圖來看，云量的分布有三個峰值帶，分別位于北緯10°，和南北緯60°附近，云量都在70%以上。而新方案與Eddington方法相比，改善最大的地方也是在北緯10°和南緯60°附近，表明，本文提出的新方案在云量較多的地方對模式模擬的短波大氣頂云輻射強(qiáng)迫有較大改進(jìn)。

5 主要結(jié)論

本文利用單層四流球諧函數(shù)譜展開算法結(jié)合二流累加法，構(gòu)造了一種精度介于二流輻射方案和四流矩陣求解法之間的二流—四流球諧函數(shù)譜展開累加輻射傳輸算法。并將其應(yīng)用于中國氣象局國家氣候中心的全球氣候模式BCC_AGCM2.0.1的新版本中（荊現(xiàn)文和張華，2012；Zhang et al., 2014），檢驗(yàn)了本文構(gòu)建的新方案的模擬效果并與原方案進(jìn)行了比較。得到的主要結(jié)論如下：

（1）在晴空情況下，Eddington近似方案和本文提出的新方案在大氣頂和地表的輻射通量差異主要是由氣溶膠引起的。其中對于晴空短波地表向下輻射通量，海鹽氣溶膠和沙塵氣溶膠分別是造成Eddington近似方案和新方案差異的主要來源。對于晴空短波大氣頂向上輻射通量，海鹽氣溶膠是造成這種差異的主要來源。

（2）在有云大氣情況下，Eddington近似方案和新方案在大氣頂和地表的輻射通量差異主要是由于太陽天頂角決定的。對于有云大氣短波地表向下輻射通量而言，Eddington近似方案和新方案的差值從赤道向兩極逐漸變大，南北兩極的差別最大。對于有云大氣大氣頂短波向上輻射通量，Eddington近似方案與新方案相比，在高緯度地區(qū)低估了大氣頂短波向上輻射通量，而在熱帶地區(qū)高估了大氣頂短波向上輻射通量。

（3）對于短波加熱率而言，新方案的加熱率普遍要強(qiáng)于Eddington近似方案，尤其在800 hPa到地面區(qū)域和50 hPa到100 hPa區(qū)域。其中800 hPa下方的加熱率，差異在云頂最大，因此會抑制云的生成，從而增加云的短波輻射通量；而100 hPa上方的正加熱偏差則有助于改善在大多數(shù)模式中普遍存在的赤道平流層中下層的冷偏差現(xiàn)象。

此外，本文提出的新方案與二流累加輻射傳輸方案相比的最大優(yōu)勢是顯著提高了對單層反射率和透過率的模擬精度。當(dāng)然本文構(gòu)造的新方案也有一定的局限性，單層的四流球諧函數(shù)方法可以用來計算各個方向的輻射強(qiáng)度，但受限于二流累加法，必須將各個方向的輻射強(qiáng)度積分到向上和向下兩個方向的輻射通量。在這個過程中將損失部分精度，因此下一步工作我們期望將新研制的四流球諧函數(shù)輻射傳輸累加法（Zhang and Li，2013）應(yīng)用 于本文的全球氣候模式中并對模擬精度和計算效率進(jìn)行對比研究。

（References：）

Ayash T, Gong S L, Jia C Q. 2008. Implementing the delta-four-stream approximation for solar radiation computations in an atmosphere general circulation model [J]. J. Atmos. Sci., 65 (7): 2448–2457.

Barker H W, Stephens G L, Partain P T, et al. 2003. Assessing 1D atmospheric solar radiative transfer models: Interpretation and handling of unresolved clouds [J]. J. Climate, 16 (16): 2676–2699.

陳洪濱, 卞建春, 呂達(dá)仁. 2006. 上對流層—下平流層交換過程研究的進(jìn)展與展望 [J]. 大氣科學(xué), 30 (5): 813–820. Chen H B, Bian J C, Lü D R. 2006. Advances and prospects in the study of stratosphere exchange [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 30(5): 813–820.

Chou M D. 1992. A solar radiation model for use in climate studies [J]. J. Atmos. Sci., 49 (9): 762–772.

Coakley J A Jr, Cess R D, Yurevich F B. 1983. The effect of tropospheric aerosols on the Earth’s radiation budget: A parameterization for climate models [J]. J. Atmos. Sci., 40 (1): 116–138.

丁守國, 趙春生, 石廣玉, 等. 2005. 近20年全球總云量變化趨勢分析 [J]. 應(yīng)用氣象學(xué)報, 16 (5): 670–677. Ding S G, Zhao C S, Shi G Y, et al. 2005. Analysis of global total cloud amount variation over the past 20 years [J]. Journal of Applied Meteorological Science (in Chinese), 16(5): 670–677.

Forster P M, Fomichev V I, Rozanov E, et al. 2011. Evaluation of radiation scheme performance within chemistry climate models [J]. J. Geophys. Res., 116 (D10): 10302, doi:10.1029/2010JD015361.

Fouquart Y, Bonnel B, Ramaswamy V. 1991. Intercomparing shortwave radiation codes for climate studies [J]. J. Geophys. Res., 96 (D5): 8955– 8968.

Gong S L, Barrie L A, Lazare M. 2002. Canadian Aerosol Module (CAM): A size-segregated simulation of atmospheric aerosol processes for climate and air quality models 2. Global sea-salt aerosol and its budgets [J]. J. Geophys. Res., 107 (D24): AAC 13-1–AAC 13-14, doi:10.1029/ 2001JD002004.

Gong S L, Barrie L A, Blanchet J P, et al. 2003. Canadian Aerosol Module: A size-segregated simulation of atmospheric aerosol processes for climate and air quality models 1. Module development [J]. J. Geophys. Res., 108 (D1): AAC 3-1–AAC 3-16, doi:10.1029/2001JD002002.

Halthore R N, Crisp D, Schwartz S E, et al. 2005. Intercomparison of shortwave radiative transfer codes and measurements [J]. J. Geophys. Res., 110(D11): D11206, doi:10.1029/2004JD005293.

Hurrell J W, and Trenberth K E. 1999. Global sea surface temperature analyses: multiple problemsand their implications for climate analysis, modeling, and reanalysis [J] Bull. Amer. Met. Soc., 80 (12):2661?2678

荊現(xiàn)文, 張華. 2012. McICA云—輻射方案在國家氣候中心全球氣候模式中的應(yīng)用與評估 [J]. 大氣科學(xué), 36 (5): 945–958. Jing X W, Zhang H. 2012. Application and evaluation of McICA Cloud-Radiation framework in the AGCM of the National Climate Center [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 36(5): 945–958.

Jing X W, and Zhang H. 2013. Application and evaluation of McICA scheme in BCC_AGCM2.0.1 [C]. AIP Conf. Proc. 1531, 756, doi:10.1063/1.4804880.

Kay M J, Box M A, Trautmann T, et al. 2001. Actinic flux and net flux calculations in radiative transfer—A comparative study of computational efficiency [J]. J. Atmos. Sci., 58 (24): 3752–3761.

Li J, Ramaswamy V. 1996. Four-stream spherical harmonic expansion approximation for solar radiative transfer [J]. J. Atmos. Sci., 53 (8): 1174–1186.

Liou K N, Fu Q, Ackerman T P. 1988. A simple formulation of the delta-four-stream approximation for radiative transfer paramterization [J]. J. Atmos. Sci., 45 (13): 1940–1947.

Liu M, Nachamkin J E, Westphal D L. 2009. On the improvement of COAMPS weather forecasts using an advanced radiative transfer model [J]. Wea. Forecasting, 24 (1): 286–306.

Randles C A, Kinne S, Myhre G, et al. 2013. Intercomparison of shortwave radiative transfer schemes in global aerosol modeling: Results from the AeroCom Radiative Transfer Experiment [J]. Atmos. Chem. Phys., 13: 2347–2379.

Scinocca J F, McFarlane N A, Lazare M, et al. 2008. Technical Note: The CCCma third generation AGCM and its extension into the middle atmosphere [J]. Atmos. Chem. Phys., 8, 7055–7074, doi:10.5194/ acp-8–7055-2008.

Shibata K, Uchiyama A. 1992. Accuracy of the delta-four-stream approximation in inhomogeneous scattering atmospheres [J]. J. Meteor. Soc. Japan, 70 (6): 1097–1109.

Shi G Y, 1981. An accurate calculation and representation of the infrared transmission function of the atmospheric constituents [D]. Ph.D. dissertation, Dept. of Science, Tohoku University of Japan, 71pp.

石廣玉. 2007. 大氣輻射學(xué)[M]. 北京: 科學(xué)出版社, 1. Shi Guangyu. 2007. Atmospheric Radiation (in Chinese) [M]. Beijing: Science Press, 1.

Wu T W, Yu R C, Zhang F. 2008. A modified dynamic framework for the atmospheric spectral model and its application [J]. J. Atmos. Sci., 65 (7): 2235–2253.

Wu T W, Yu R C, Zhang F, et al. 2010. The Beijing Climate Center atmospheric general circulation model: Description and its performance for the present-day [J]. Climate Dyn., 34 (1): 123–147.

Zhang F, Li J. 2013. Doubling-adding method for delta-four-tream spherical harmonic expansion approximation in radiative transfer parameterization [J]. J. Atmos. Sci., 70: 3084–3101.

Zhang F, Shen Z P, Li J N, et al. 2013. Analytical delta-four-stream doubling–adding method for radiative transfer parameterizations [J]. J. Atmos. Sci., 70 (3): 794–808.

張華. 1999. 非均勻路徑相關(guān)-分布方法的研究 [D]. 中國科學(xué)院大氣物理研究所博士學(xué)位論文, 169pp. Zhang Hua. 1999. On the study of a new correlated-distribution method for nongray gaseous absorption in the inhomogeneous scattering atmosphere [D]. Ph. D. dissertation (in Chinese), Institute of Atmospheric Physics, Chinese Academy of Sciences, 169pp.

Zhang H, Nakajima T, Shi G Y, et al. 2003. An optimal approach to overlapping bands with correlateddistribution method and its application to radiative calculations [J]. J. Geophys. Res., 108(D20): 4641, doi:10.1029/2002JD003358.

Zhang H, Shi G Y, Nakajima T, et al. 2006a. The effects of the choice of the-interval number on radiative calculations [J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 98 (1): 31–43.

Zhang H, Suzuki T, Nakajima T, et al. 2006b. Effects of band division on radiative calculations [J]. Optical Engineering, 45 (1): 016002.

Zhang H, Wang Z L, Wang Z Z, et al. 2012. Simulation of direct radiative forcing of aerosols and their effects on East Asia climate using an interactive GCM–aerosol coupled system [J]. Climate Dyn., 38 (7): 1675–1693.

Zhang H, Jing X W, Li J N. 2014. Application and evaluation of a new radiation code under McICA scheme in BCC_AGCM2.0.1 [J]. Geosci. Model Dev., 7: 737–754.

張華, 盧鵬, 荊現(xiàn)文. 2015. 二流—四流球諧函數(shù)譜展開累加輻射傳輸方案在全球氣候模式中的應(yīng)用[J]. 大氣科學(xué), 39 (1): 137?144, doi:10.3878/j.issn. 1006-9895.1404.13316. Zhang Hua, Lu Peng, Jing Xianwen. 2015. Application of two-four stream spherical harmonic expansion approximation in a global climate model [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 39 (1): 137?144.

Application of Two?Four Stream Spherical Harmonic Expansion Approximation in a Global Climate Model

Zhang Hua1, 2, LU Peng3, 1, and Jing Xianwen1, 2

1,,,100081;2,,210044;3,210008

In this study, a new scheme for the radiative transfer algorithm, calledthe two?four stream spherical harmonic expansion approximation, is built and applied in the new version of the Beijing Climate Center atmospheric general circulation model (BCC_AGCM2.0.1). It is then compared with the original Eddington approximation scheme. Because this new scheme expands the Eddington approximation to solve radiative transfer through the atmosphere, it has better accuracy. We found the new scheme to have a great effect on climatic simulation. In a clear sky, the new scheme reduces the shortwave downward radiative flux in the surface in the 30°–60°S regions, in the Northeast Atlantic, and in the Sahara desert, with the largest reduction being 3.5 W/m2. Meanwhile, it increases the shortwave upward radiative flux at the top of the atmosphere (TOA) in the 30–60°S regions and in the Northeast Atlantic, with the largest increase being 3 W/m2. For all-sky cases, the new scheme reduces the shortwave downward radiative flux, and the difference between the two schemes becomes larger with increasing latitude. The largest difference reaches 5.5 W/m2in the two polar regions. The new scheme also reduces the shortwave upward radiative flux at the TOA in the tropics, with the largest difference being 2.5 W/m2, but increases this flux in the 30°–60°S regions, with the largest difference being 1.5 W/m2. Moreover, the new scheme increases the shortwave heating rate within the atmosphere generally, especially for the levels between 800 hPa and the surface and between 50 and 100 hPa where the largest difference reaches 0.03 K/d. Therefore, the new scheme is useful in global climate modeling for improving the so-called temperature cold bias phenomena generally existing in the lower parts of the stratosphere above the tropics.

Four-stream spherical harmonic expansion, Radiative transfer, Radiative flux, Heating rate, Temperature cold bias

1006?9895(2015)01?0137?08

P422

A

10.3878/j.issn.1006-9895.1404.13316

2013?03?16；網(wǎng)絡(luò)預(yù)出版日期2014?05?07

國家自然科學(xué)基金項目41375080，科技部公益性行業(yè)（氣象）科研專項項目GYHY201406023，國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃項目2011CB403405

張華，女，1965年出生，博士，研究員，主要從事氣溶膠—云—輻射相互作用及其氣候效應(yīng)、大氣輻射傳輸模式等方面的研究。E-mail: huazhang@cma.gov.cn