李秀鎮(zhèn), 盛立芳, 劉　騫, 劉玉芝, 王力群

(1.中國(guó)衛(wèi)星海上測(cè)控部，江蘇 江陰 214431； 2.中國(guó)海洋大學(xué)海洋與大氣學(xué)院，山東 青島 266100； 3.中國(guó)海洋大學(xué)物理海洋教育部重點(diǎn)實(shí)驗(yàn)室，海洋-大氣相互作用與氣候?qū)嶒?yàn)室，山東 青島 266100； 4.蘭州大學(xué)大氣科學(xué)學(xué)院，甘肅 蘭州 730000)

?

基于SBDART輻射傳輸模式的晴天地面總輻射模擬誤差分析*

李秀鎮(zhèn)1, 盛立芳2,3**, 劉騫2, 劉玉芝4, 王力群1

(1.中國(guó)衛(wèi)星海上測(cè)控部，江蘇 江陰 214431； 2.中國(guó)海洋大學(xué)海洋與大氣學(xué)院，山東 青島 266100； 3.中國(guó)海洋大學(xué)物理海洋教育部重點(diǎn)實(shí)驗(yàn)室，海洋-大氣相互作用與氣候?qū)嶒?yàn)室，山東 青島 266100； 4.蘭州大學(xué)大氣科學(xué)學(xué)院，甘肅 蘭州 730000)

基于SBDART輻射傳輸模式，利用POM-02型天空輻射計(jì)觀測(cè)和反演獲得的氣溶膠光學(xué)參數(shù)計(jì)算地面總輻射，并與CM21輻射表的觀測(cè)結(jié)果對(duì)比以探討晴天地面總輻射模擬的誤差來(lái)源。結(jié)果表明：地面總輻射對(duì)氣溶膠單次散射比和水汽柱含量的敏感性較強(qiáng)，模擬中對(duì)這2個(gè)因子減少50%會(huì)分別造成地面總輻射量增加56.99W/m2和減少25.45W/m2。在準(zhǔn)確輸入水汽柱含量基礎(chǔ)上變換大氣廓線(xiàn)對(duì)地面總輻射計(jì)算結(jié)果影響甚微。相關(guān)性及多元線(xiàn)性回歸分析表明，天頂角、氣溶膠單次散射比和光學(xué)厚度是地面總輻射相對(duì)誤差的主要來(lái)源。利用上述3個(gè)因子建立回歸方程，回歸結(jié)果與地面總輻射相對(duì)誤差的相關(guān)系數(shù)為0.486(通過(guò)了置信度0.01的顯著性檢驗(yàn))，說(shuō)明回歸方程能夠準(zhǔn)確的計(jì)算相對(duì)誤差，利用該方法可以對(duì)本站模式計(jì)算結(jié)果進(jìn)行修訂。

輻射傳輸模式； 探空資料； 總輻射表； 氣溶膠光學(xué)厚度

引用格式：李秀鎮(zhèn), 盛立芳, 劉騫， 等. 基于SBDART輻射傳輸模式的晴天地面總輻射模擬誤差分析[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016, 46(8): 13-18.

LI Xiu-Zhen, SHENG Li-Fang, LIU Qian, et al. Error in calculation of surface radiation based on SBDART radiative transfer model[J]. Periodical of Ocean University of China, 2016, 46(8): 13-18.

太陽(yáng)輻射是地表面最主要的能量熱源[1]，到達(dá)地面的太陽(yáng)輻射受到大氣中氣體吸收、氣溶膠和云散射等作用的影響，地面太陽(yáng)輻射的變化也會(huì)對(duì)地表溫度、空氣濕度和水循環(huán)等造成影響，研究地面總輻射的特征變化具有重要的氣候意義。中國(guó)的地面輻射測(cè)站尚存在數(shù)量少、時(shí)空分布不均的不足，在某些地區(qū)尚不能滿(mǎn)足相關(guān)需求。國(guó)內(nèi)外快速發(fā)展的輻射傳輸模式為開(kāi)展相關(guān)研究開(kāi)辟了新的途徑，隨之而來(lái)的對(duì)模式性能的檢驗(yàn)受到科學(xué)家的頗多關(guān)注[2-9]。

早期研究發(fā)現(xiàn)，模式計(jì)算的地面直接輻射與觀測(cè)值吻合較好，而模式計(jì)算的散射輻射卻遠(yuǎn)遠(yuǎn)高于觀測(cè)值[2,10-11]。相關(guān)研究將模擬誤差歸因于作為標(biāo)準(zhǔn)的輻射觀測(cè)值[12-15]，但對(duì)觀測(cè)值進(jìn)行修訂后模擬值依然偏高[16-18]。誤差原因進(jìn)而轉(zhuǎn)移到模式輸入?yún)?shù)(水汽、臭氧、氣溶膠等)的數(shù)據(jù)質(zhì)量上來(lái)。Michalsky等[19]使用高質(zhì)量氣溶膠光學(xué)參數(shù)、水汽、臭氧觀測(cè)數(shù)據(jù)分別對(duì)6類(lèi)輻射傳輸模式進(jìn)行檢驗(yàn)，平均誤差減小到小于1%和小于1.9%，其他一些研究也表明，提高輸入?yún)⒘康木葘⑻岣吣Ｊ降挠?jì)算準(zhǔn)確度[7,20]。在眾多的研究中，觀測(cè)手段、地點(diǎn)的不同所造成的觀測(cè)數(shù)據(jù)質(zhì)量的差異導(dǎo)致了模式計(jì)算精度的不同,其中在對(duì)水汽、臭氧等參數(shù)在大氣中的垂直分布——大氣廓線(xiàn)的使用上，有的學(xué)者僅采用模式提供的標(biāo)準(zhǔn)大氣廓線(xiàn)[21-25]，有的站點(diǎn)則采用了實(shí)時(shí)觀測(cè)獲得的探空廓線(xiàn)[2-4,6,10]，其計(jì)算結(jié)果的誤差精度也存在差異。

在實(shí)際工作中，高質(zhì)量數(shù)據(jù)的時(shí)間長(zhǎng)度通常僅限于某次實(shí)驗(yàn)過(guò)程或某個(gè)單一站點(diǎn)，大氣廓線(xiàn)的實(shí)時(shí)獲取也難以長(zhǎng)期保證，提高現(xiàn)有數(shù)據(jù)情況下模式的計(jì)算準(zhǔn)確度顯得尤為重要。本文以長(zhǎng)期觀測(cè)的高分辨率天空輻射計(jì)數(shù)據(jù)等資料作為模式輸入?yún)?shù)，使用SBDART(Santa Barbara DISORT Atmospheric Radiative Transfer)輻射傳輸模式計(jì)算地面總輻射，與CM21輻射表觀測(cè)數(shù)據(jù)進(jìn)行比較，區(qū)分各參數(shù)對(duì)模式計(jì)算結(jié)果的影響程度，確定輸入?yún)?shù)及大氣廓線(xiàn)的最優(yōu)組合方案。利用相關(guān)性及多元回歸分析方法研究造成本站地面總輻射計(jì)算誤差的影響因子，探討提高準(zhǔn)確度的方法，以期為改善現(xiàn)有數(shù)據(jù)情況下地面總輻射計(jì)算質(zhì)量提供依據(jù)。

1　資料與方法

1.1 氣溶膠和地面總輻射資料來(lái)源

氣溶膠光學(xué)參數(shù)(光學(xué)厚度，單次散射比，不對(duì)稱(chēng)因子)和地面總輻射觀測(cè)點(diǎn)位于中國(guó)海洋大學(xué)校園內(nèi)八關(guān)山(120°20.134′E，36°3.809′N(xiāo))氣象臺(tái)樓頂。八關(guān)山海拔高度68m，氣象臺(tái)樓頂距地面15m。氣溶膠觀測(cè)儀器來(lái)自PREDE POM-02型11波段天空輻射計(jì)，最小視場(chǎng)角為3°，11個(gè)波段的中心波長(zhǎng)分別位于340、380、400、500、675、870、940、1020、1225、1600和2200nm，濾光片的半波寬度為10nm。該輻射計(jì)可以對(duì)太陽(yáng)自動(dòng)跟蹤，并按照設(shè)定的程序每一定時(shí)間間隔進(jìn)行太陽(yáng)直接輻射和天空散射輻射觀測(cè)。觀測(cè)時(shí)間2002—2009年，其中2005年缺測(cè)。該輻射計(jì)每年3月份定標(biāo)一次。光學(xué)厚度(AOD: Aerosol Optical Depth)、單次散射比(SSA: Single-Scattering Albedo)、不對(duì)稱(chēng)因子等數(shù)據(jù)由SKYRAD.PACK4.2分析獲得。地面總輻射由Kipp&Zonen公司的CM21輻射表進(jìn)行觀測(cè)。該儀器是一種高性能的總輻射表，用于測(cè)量水平面上的輻射通量，即直接輻射和上半球2 π角弧度天空內(nèi)向下的散射輻射之和。儀器響應(yīng)時(shí)間5s，方向性誤差±5W/m2，傾斜誤差±0.2%，光譜范圍200～3600nm。觀測(cè)時(shí)間2003—2009年，其中2005年缺測(cè)。

晴天數(shù)據(jù)的選取依據(jù)Khatri和Takamura[26]給出的云剔除算法。首先依據(jù)晴空條件下地面總輻射日變化趨勢(shì)與有云條件下地面總輻射日變化趨勢(shì)的差異將有云影響和無(wú)云條件下的觀測(cè)時(shí)段分離開(kāi)來(lái)，再采用改進(jìn)的SVA算法進(jìn)行光譜變率的檢驗(yàn)，最后利用統(tǒng)計(jì)分析的方法剔除極端值,獲得36d晴天數(shù)據(jù)。選取氣溶膠光學(xué)參數(shù)及地面觀測(cè)輻射同時(shí)存在的時(shí)刻，總計(jì)1274個(gè)時(shí)次。

1.2 模式介紹

SBDART輻射傳輸模式是美國(guó)加利福尼亞大學(xué)地球系統(tǒng)計(jì)算學(xué)院的地球空間研究組開(kāi)發(fā)的，用來(lái)計(jì)算晴空和有云狀況下平面平行輻射傳輸[27]。模式采用了5S模型中的6種大氣廓線(xiàn)(熱帶、中緯度夏季、中緯度冬季、副極地夏季、副極地冬季及US62大氣)，LOWTRAN7和MODTRAN3模型中的標(biāo)準(zhǔn)大氣氣溶膠模型，5種基本的地表類(lèi)型(海水、湖水、植被、雪面、沙地)，允許用戶(hù)自定義大氣廓線(xiàn)和氣溶膠特性。程序提供了水汽、臭氧、甲烷、一氧化碳、一氧化二碳的混合比垂直廓線(xiàn)及其他微量氣體的垂直廓線(xiàn)。輻射傳輸方程利用離散縱坐標(biāo)法程序(DISORT)。

本研究中模式輸入的參數(shù)包括時(shí)間、經(jīng)緯度、大氣廓線(xiàn)、氣溶膠光學(xué)參數(shù)(AOD，SSA，不對(duì)稱(chēng)因子)、地表反照率。本次研究主要考慮了可見(jiàn)光波段，設(shè)定計(jì)算波長(zhǎng)范圍為250～3500nm。

基于探空資料和模式自帶的中緯度大氣廓線(xiàn)建立模式所需要的大氣溫度、壓強(qiáng)、水汽密度、臭氧密度廓線(xiàn)(探空廓線(xiàn))。大氣柱高度設(shè)定為100km。探空高度范圍以?xún)?nèi)各層上的溫度、氣壓、水汽密度取探空值，臭氧密度由中緯度大氣廓線(xiàn)中相鄰高度層的臭氧數(shù)據(jù)通過(guò)線(xiàn)性?xún)?nèi)插獲得。探空高度范圍之外的溫度、壓強(qiáng)、臭氧密度、水汽密度由中緯度大氣廓線(xiàn)補(bǔ)充。臭氧柱含量數(shù)據(jù)取自TOMS臭氧數(shù)據(jù)(http://toms.gsfc.nasa.gov)。

文中使用了兩類(lèi)水汽柱含量數(shù)據(jù)，包括探空獲得的大氣可降水量數(shù)據(jù)(Precipitable water for entire sounding)以及通過(guò)積分探空廓線(xiàn)數(shù)據(jù)間接獲得的水汽柱含量[28]。探空觀測(cè)時(shí)間分別為00(UTC)時(shí)和12時(shí)，水汽柱含量數(shù)據(jù)取2個(gè)時(shí)次的平均值。

由于550nm的AOD沒(méi)有直接觀測(cè)，利用400、500、675和870nm的數(shù)值通過(guò)內(nèi)插獲得。計(jì)算公式參考?ngstrom關(guān)系式：

τ=aλ-n。

其中:τ是氣溶膠光學(xué)厚度;a為常數(shù);n是?ngstrom指數(shù)。

1.3 誤差分析方法

利用平均偏差(MBE)、平均絕對(duì)偏差(MABE)、平均絕對(duì)百分比誤差(MAPE)、相對(duì)均方根誤差σ、相對(duì)誤差(RE)進(jìn)行了誤差分析。

2　結(jié)果與分析

2.1 參數(shù)敏感性試驗(yàn)

為了檢驗(yàn)地面總輻射模擬對(duì)輸入?yún)⒘康拿舾行?，通過(guò)單獨(dú)改變某一輸入?yún)⒘康臄?shù)值，比較計(jì)算結(jié)果的變化，進(jìn)而說(shuō)明由于該參量的不確定性引起的模擬誤差。所選晴天數(shù)據(jù)中AOD、SSA、不對(duì)稱(chēng)因子、水汽柱含量、臭氧柱含量取值范圍分別為0.06～1.02，0.36～1，0.47～0.76，0.1～3.0g/cm2，0.272～0.409atm-cm。

首先確定一組基準(zhǔn)參量。大氣溫度、壓強(qiáng)、水汽密度、臭氧密度廓線(xiàn)取自2008年11月19日的探空資料(該天光學(xué)厚度日變化最小)，地表反照率取0.2[29]。AOD、SSA、不對(duì)稱(chēng)因子、水汽柱含量、臭氧柱含量基準(zhǔn)值分別取0.14、0.98、0.63、0.35g/cm2、0.33atm-cm(MBE=-2.29%)。在上述參量取值范圍內(nèi)，保持其他參量數(shù)值不變，單獨(dú)改變某一參量，其中修改步長(zhǎng)為各參量對(duì)應(yīng)基準(zhǔn)值的10%。地面總輻射計(jì)算結(jié)果隨各參量的變化見(jiàn)圖1，橫坐標(biāo)數(shù)值為輸入?yún)⒘恐蹬c對(duì)應(yīng)基準(zhǔn)值的比值(比值大于1時(shí)為正誤差，反之為負(fù)誤差)，紅色橫線(xiàn)表示總輻射基準(zhǔn)值。各參量單獨(dú)改變時(shí)，總輻射計(jì)算值均有明顯改變。AOD、水汽柱和臭氧柱含量輸入值大于基準(zhǔn)值時(shí)將導(dǎo)致地面總輻射出現(xiàn)負(fù)誤差(計(jì)算值偏小)，反之出現(xiàn)正誤差；SSA、不對(duì)稱(chēng)因子、地表反照率存在誤差時(shí)所引起的誤差情況與前三者相反。

圖1　敏感性試驗(yàn)輸入?yún)⒘颗c輻射計(jì)算值變化

取定曲線(xiàn)上任意一點(diǎn)，曲線(xiàn)在該處的斜率表征了輸入?yún)⒘康牟淮_定性所引起的總輻射的模擬誤差，斜率越大誤差越大。除水汽柱含量曲線(xiàn)外各輸入?yún)⒘康挠?jì)算曲線(xiàn)近似呈線(xiàn)性分布，地表反照率、AOD、SSA、不對(duì)稱(chēng)因子、臭氧柱含量的斜率分別為13.88、-22.95、113.98、31.96和-10.92。故各參量存在相同的百分比誤差時(shí)所引起的地面總輻射的計(jì)算誤差大小排序?yàn)镾SA>不對(duì)稱(chēng)因子>AOD>地表反照率>臭氧柱含量。當(dāng)水汽柱含量存在負(fù)誤差時(shí)的斜率為-50.91，其引起的計(jì)算誤差小于SSA而大于不對(duì)稱(chēng)因子；正誤差情況下斜率為-6.45，計(jì)算誤差小于其余各參量。通過(guò)對(duì)個(gè)參量斜率的分析可以看出，地面總輻射對(duì)SSA和水汽主含量的變化最為敏感，當(dāng)模式中這2個(gè)因子減少50%時(shí)將會(huì)分別造成地面總輻射量增加56.99W/m2和減少25.45W/m2。

2.2 大氣廓線(xiàn)對(duì)模擬結(jié)果的影響

水汽柱含量和臭氧柱含量對(duì)模式的影響在敏感性試驗(yàn)中已得到證實(shí)，在二者總量確定的情況下，其在大氣中的不同分布對(duì)計(jì)算結(jié)果有無(wú)影響，大氣溫度、氣壓等參數(shù)的垂直分布對(duì)模式的影響如何？分別采用探空廓線(xiàn)以及標(biāo)準(zhǔn)大氣廓線(xiàn)(中緯度)進(jìn)行了地面總輻射的計(jì)算，其中水汽柱含量分別采用積分計(jì)算水汽柱含量與可降水量數(shù)據(jù)。共設(shè)計(jì)4類(lèi)計(jì)算方案(見(jiàn)表1)，對(duì)1274個(gè)時(shí)刻的計(jì)算結(jié)果進(jìn)行比較。

相關(guān)分析表明，各方案的計(jì)算結(jié)果與觀測(cè)值間的相關(guān)系數(shù)均大于0.997，通過(guò)了置信度0.01的顯著性檢驗(yàn)。就各類(lèi)偏差(MBE、MABE、MAPE、σ)而言，方案2<方案4<方案1<方案3。分別采用標(biāo)準(zhǔn)大氣廓線(xiàn)與探空廓線(xiàn)的情況下，可降水量數(shù)據(jù)的計(jì)算結(jié)果均優(yōu)于通過(guò)積分獲得的水汽柱含量數(shù)據(jù)，即可降水量數(shù)據(jù)比水汽柱含量數(shù)據(jù)更具代表性。使用相同的水汽柱含量的情況下，方案1比方案3的MBE低0.7，方案2的MBE比方案4低0.2；標(biāo)準(zhǔn)大氣廓線(xiàn)計(jì)算誤差略低于探空廓線(xiàn)，但二者的偏差值不大，這與余予等人[7]的結(jié)果相一致。分析標(biāo)準(zhǔn)大氣廓線(xiàn)計(jì)算結(jié)果優(yōu)于探空廓線(xiàn)的可能原因，探空廓線(xiàn)由08時(shí)和20時(shí)的探空數(shù)據(jù)計(jì)算獲得，該廓線(xiàn)不能準(zhǔn)確代表全天的大氣狀態(tài)。

表1　總輻射計(jì)算值與觀測(cè)值比較

2.3 影響地面總輻射計(jì)算誤差的因子分析

選取最優(yōu)方案2對(duì)誤差來(lái)源進(jìn)行了初步探討，為排除輻射值本身的影響選取相對(duì)誤差進(jìn)行分析。各輸入?yún)⒘颗c相對(duì)誤差的時(shí)間序列見(jiàn)圖2。氣溶膠AOD、SSA與不對(duì)稱(chēng)因子呈現(xiàn)倒“V”型日變化(見(jiàn)圖2a)。與之相類(lèi)似，相對(duì)誤差在一定程度上也呈現(xiàn)倒“V”型日變化趨勢(shì)(見(jiàn)圖2c)。相對(duì)誤差為負(fù)值(總輻射計(jì)算值<觀測(cè)值)的情況占總樣本數(shù)83.28%；相對(duì)誤差為正值的情況較少，集中在2004年1、6月和2009年5月。天頂角(模式通過(guò)時(shí)間、經(jīng)緯度計(jì)算獲得)存在與相對(duì)誤差相反的日變化趨勢(shì)，即呈“V”狀分布(見(jiàn)圖2b)。

圖2　輸入?yún)⒘?氣溶膠AOD、SSA、不對(duì)稱(chēng)因子及天頂角)與相對(duì)誤差的時(shí)間序列

相關(guān)分析表明，相對(duì)誤差與SSA、天頂角的相關(guān)系數(shù)分別為0.37、-0.45，通過(guò)了置信度0.01的顯著性檢驗(yàn)，分別呈現(xiàn)正相關(guān)和負(fù)相關(guān)關(guān)系。相對(duì)誤差與AOD、不對(duì)稱(chēng)因子間的相關(guān)系數(shù)分別為-0.11、-0.13，均未通過(guò)顯著性檢驗(yàn)。在觀測(cè)時(shí)間段內(nèi)相對(duì)誤差主要來(lái)源于SSA和天頂角?？赡艿恼`差原因有二：地面總輻射對(duì)SSA的敏感性較強(qiáng)，SSA的觀測(cè)誤差導(dǎo)致了一定的計(jì)算誤差；在天頂角較小即入射輻射值較大時(shí)SBDART輻射傳輸模式的模擬性能較差。由于地面總輻射對(duì)AOD、不對(duì)稱(chēng)因子的敏感性不如SSA，二者的觀測(cè)誤差對(duì)模式誤差的貢獻(xiàn)小于SSA。

通過(guò)回歸分析的方法對(duì)模式的計(jì)算結(jié)果進(jìn)行修訂，選取回歸模型為

y=β0+β1x1+β2x2+…+βixi+ε,ε～N(0,σ2)。

其中：y為地面總輻射相對(duì)誤差；xi為預(yù)報(bào)因子(影響誤差的輸入?yún)⒘?。

回歸分析設(shè)定4組，預(yù)報(bào)因子組合如下：天頂角(Ⅰ)；天頂角和SSA(Ⅱ)；天頂角、SSA與AOD(Ⅲ)；天頂角、SSA、AOD和不對(duì)稱(chēng)因子(Ⅳ)。選擇前950個(gè)時(shí)刻的樣本進(jìn)行多元回歸，利用回歸方程計(jì)算剩余時(shí)刻的相對(duì)誤差與模式計(jì)算誤差比較檢驗(yàn)回歸效果。

4組回歸分析結(jié)果與原誤差間的相關(guān)系數(shù)分別為0.0936(Ⅰ)、0.2896(Ⅱ)、0.4763(Ⅲ)、0.4766(Ⅳ)。僅采用天頂角的情況下計(jì)算效果比較差，誤差來(lái)源不能僅考慮天頂角的影響。第Ⅱ組添加了SSA的影響，相比第一組相關(guān)系數(shù)有了較大變化，說(shuō)明SSA對(duì)相對(duì)誤差有著非常重要的影響。第Ⅲ組又增加了AOD的影響，相比第Ⅱ組相關(guān)系數(shù)變化也比較大，AOD與地面總輻射相對(duì)誤差間的相關(guān)性雖然不顯著，但回歸分析的結(jié)果表明其對(duì)相對(duì)誤差的影響是不可忽略的，而且影響是顯著的。第Ⅳ組的結(jié)果與第Ⅲ組相差甚微，表明不對(duì)稱(chēng)因子對(duì)相對(duì)誤差的影響是可以忽略不計(jì)的。

通過(guò)以上分析，確定回歸方程為:

y=-0.09964-0.001938×x1+0.2333×

x2-0.1552×x3。

其中x1、x2、x3分別是天頂角、SSA、AOD。利用上式對(duì)剩余374個(gè)時(shí)刻的相對(duì)誤差進(jìn)行計(jì)算，與模式計(jì)算誤差進(jìn)行比較(見(jiàn)圖3)，結(jié)果表明，回歸方程計(jì)算結(jié)果與模式誤差存在較好的一致性，二者相關(guān)系數(shù)為0.485，通過(guò)了置信度0.01的顯著性檢驗(yàn)，利用回歸方程計(jì)算相對(duì)誤差是可行的。

3　結(jié)論與討論

本文利用SBDART輻射傳輸模式計(jì)算地面總輻射，通過(guò)與觀測(cè)值的比較分析得到如下幾點(diǎn)認(rèn)識(shí)：

圖3　輸入?yún)⒘颗c相對(duì)誤差

(1)SBDART輻射傳輸模式對(duì)氣溶膠SSA及水汽柱含量最為敏感，提高其數(shù)據(jù)質(zhì)量是保證計(jì)算準(zhǔn)確度的有效途徑。敏感性試驗(yàn)結(jié)果主要體現(xiàn)了模式本身的特點(diǎn)，因此在使用SBDART輻射傳輸模式對(duì)不同時(shí)間和地點(diǎn)進(jìn)行模擬時(shí)，本研究的結(jié)果依然具有參考價(jià)值。

(2)單位氣柱中水汽和臭氧的總量對(duì)模式的影響大于其垂直分布(大氣廓線(xiàn))帶來(lái)的影響。當(dāng)探空數(shù)據(jù)不能表征全天大氣狀態(tài)時(shí)，使用標(biāo)準(zhǔn)大氣廓線(xiàn)可以更準(zhǔn)確地計(jì)算地面總輻射。以中國(guó)為例，氣象業(yè)務(wù)中探空觀測(cè)為北京時(shí)08和20時(shí)，研究顯示利用這一探空資料得到的大氣廓線(xiàn)進(jìn)行模擬，其結(jié)果準(zhǔn)確度并不高于使用模式中標(biāo)準(zhǔn)大氣廓線(xiàn)的結(jié)果。

(3)通過(guò)相關(guān)性和回歸分析對(duì)地面總輻射誤差來(lái)源進(jìn)行分析表明，天頂角、SSA、AOD是影響誤差來(lái)源的主要因子?；貧w方程計(jì)算結(jié)果與模式計(jì)算誤差間相關(guān)系數(shù)達(dá)0.485，利用回歸方程對(duì)地面總輻射計(jì)算結(jié)果進(jìn)行修訂是可行的。但是由于相關(guān)性和多元回歸分析對(duì)數(shù)據(jù)存在一定的“依賴(lài)性”，該結(jié)論目前僅適用于與青島類(lèi)似的沿海地區(qū)。

由于僅采用了SBDART輻射傳輸模式,文中確定的影響模式計(jì)算效果的因子以及提出的誤差修訂方法對(duì)于其他模式是否適用尚待進(jìn)一步研究。

[1]左大康, 王懿賢, 陳建綏. 中國(guó)地區(qū)太陽(yáng)總輻射的空間分布特征[J]. 氣象學(xué)報(bào), 1963, 33(1): 78-96.

Tso T K, Wang Y X, Chen C X. Characteristics of the distribution of total radiation in China [J]. Acta Meteorologica Sinica, 1963, 33(1): 78-96.

[2]Halthore R N, Stephen E S, Joseph J M, et al. Comparison of model estimated and measured direct-normal solar irradiance [J]. J Geophys Res, 1997, 102(D25): 29991-30002.

[3]Halthore R N, Stephen E S. Comparison of model-estimated and measured diffuse downward irradiance at surface in cloud-free skies [J]. J Geophys Res, 2000, 105(D15): 20165-20177.

[4]Kato S, Ackerman T P, Dutton E G, et al. A comparison of modeled and measured surface shortwave irradiance for a molecular atmosphere [J]. J Quant Spectrosc Radiat Transfer, 1999, 61(4): 493-502.

[5]Barnard J C, Powell D M. A comparison between modeled and measured clear-sky radiative shortwave fluxes in Arctic environments, with special emphasis on diffuse radiation [J]. J Geophys Res, 2002, 107(D19): 4383.

[6]Xia X, Li Z, Wang P, et al. Estimation of aerosol effects on surface irradiance based on measurements and radiative transfer model simulations in northern China [J]. J Geophys Res, 2007, 112(D22S10): 1-11.

[7]余予, 夏祥鰲, 陳洪濱, 等. 晴空大氣太陽(yáng)短波輻射觀測(cè)與模式比較[J]. 太陽(yáng)能學(xué)報(bào), 2007, 28(3): 233-240.

Yu Y, Xia X A, Chen H B, et al. A comparison between measured and modeled clear-sky surface solar irradiance [J]. Acta Energiae Solaris Sinica, 2007, 28(3): 233-240.

[8]權(quán)維俊, 陳洪濱, 郭文利, 等. 兩種反演地表凈太陽(yáng)輻射的參數(shù)化模式在北京地區(qū)的比較[J]. 氣象學(xué)報(bào), 2009, 67(4): 623-630.

Quan W J, Chen H B, Guo W L, et al. Net surface solar radiation over Beijing area derived by two parameterization models: A comparison [J]. Acta Meteorologica Sinica, 2009, 67(4): 623-630.

[9]王明歡, 賴(lài)安偉, 陳正洪, 等. WRF模式模擬的地表短波輻射與實(shí)況對(duì)比分析[J]. 氣象, 2012, 38(5): 585-592.

Wang M H, Lai A W, Chen Z H, et al. Comparison of WRF forecast downward shortwave radiation with observations-A pilot study [J]. Meteorological Monthly, 2012, 38(5): 585-592.

[10]Kato S, Ackerman T P, Clothiaux E E, et al. Uncertainties in modeled and measured clear-sky surface shortwave irradiances [J]. J Geophys Res, 1997, 102(25): 881-25, 898.

[11]Halthore R N, Nemesure S, Schwartz S E, et al. Models overestimate diffuse clear-sky irradiance: A case for excess atmospheric absorption [J]. Geophys Res Lett, 1998, 25: 3591-3594.

[12]Bush B C, Valero F P J, Simpson A S, et al. Characterization of thermal effects in pyranometers: A data correction algorithm for improved measurement of surface insolation [J]. J Atmos Ocean Technol, 2000, 17: 165-175.

[13]Haeffelin M, Kato S, Smith A M, et al. Determination of the thermal offset of the Eppley precision spectral pyranometer [J]. Appl Opt, 2001, 40: 472-484.

[14]Dutton E G, Michalsky J J, Stoffel T, et al. Measurement of broadband diffuse solar irradiance using current commercial instrumentation with a correction for thermal offset errors [J]. J Atmos Ocean Technol, 2001, 18: 297-314.

[15]Michalsky J J, Dolce R, Dutton E G, et al. Results from the first ARM diffuse horizontal shortwave irradiance comparison [J]. J Geophys Res， 2003, 108(D3)： 4108. doi: 10. 1029/2002JD002825.

[16]Halthore R N, Schwartz S E. Comparison of model estimated and measured diffuse downward irradiance at surface in cloud-free skies [J]. J Geophys Res, 2000, 105(20): 165-2177.

[17]Halthore R N, Miller M A, Ogren J A, et al. Further developments in closure experiments for diffuse irradiance under cloud-free skies at a continental site [J]. Geophys Res Lett, 2004, 31, L07111, doi: 10. 1029/2003GL019102.

[18]Henzing J S, Knap W H, Stammes P, et al. Effect of aerosols on the downward shortwave irradiances at the surface: Measurements versus calculations with MODTRAN4. 1 [J]. J Geophys Res, 2004, 109, D14204. Doi: 10. 1029/2003JD004142.

[19]Michalsky J J, Anderson G P, Barnard J, et al. Shortwave radiative closure studies for clear skies during the atmospheric radiation measurements 2003 aerosol intensive observation period [J]. J Geophys Res, 2006, 111： D14S90.

[20]傅炳珊, 陳渭民, 張鳳英. 利用TOVS 資料計(jì)算我國(guó)東南地區(qū)的太陽(yáng)直接輻射和散射輻射[J]. 南京氣象學(xué)院學(xué)報(bào), 2002, 25(6): 807-815.

Fu B S, Chen W M, Zhang F Y. Evaluation of solar direct and scattered radiation over Southeast China using Tiros Operational Vertical Sounder data [J]. Journal of Nanjing Institute of Meteorology, 2002, 25(6): 807-815.

[21]Yu S, Zender C S, Saxena V K. Direct radiative forcing and atmospheric absorption by boundary layer aerosols in the southeastern US: Model estimates on the basis of new observations [J]. Atmos Environ, 2001(35): 3967-3977.

[22]劉勇洪, 權(quán)維俊, 夏祥鰲, 等. 基于MODTRAN模式與衛(wèi)星資料的晴空凈太陽(yáng)輻射模擬[J]. 高原氣象, 2008, 27(6): 1410-1415.

Liu Y H, Quan W J, Xia X A, et al. Net surface solar radiation in the clear sky based on MODTRAN model and satellite data [J]. Plateau Meteorology, 2008, 27(6): 1410-1415.

[23]沈鐘平, 張華. 影響地面太陽(yáng)輻射及其譜分布的因子分析[J]. 太陽(yáng)能學(xué)報(bào), 2009, 30(10): 1209-1215.

Sheng Z P, Zhang H. Analysis on the factors affecting surface solar radiation and its spectral distribution [J]. Acta Energiae Solaris Sinica, 2009, 30(10): 1209-1215.

[24]劉香娥, 王廣河. FY-2C靜止氣象衛(wèi)星紅外通道輻射特性研究[J]. 氣象, 2010, 36(2): 39-47.

Liu X E, Wang G H. Analyses of the radiation characteristics from FY-2C mid-infrared channel[J]. Meteorological Monthly, 2010, 36(2): 39-47.

[25]李浩, 孫學(xué)金, 陳峰, 等. 雙波段大氣向下紅外輻射云遙感數(shù)值模擬[J]. 氣象科技, 2010, 38(2): 222-225.

Li H, Sun X J, Chen F, et al. Numerical simulation of double-band cloud remote sensing based on atmospheric downward infrared radiation [J]. Meteorological Science and Technology, 2010, 38(2): 222-225.

[26]Khatri P, Takamura T. An algorithm to screen cloud-affected data for sky radiometer data analysis [J]. J Meteor Soc Japan, 2009, 87(1): 189-204.

[27]Ricchiazzi P, Yang S, Gautier C, et al. SBDART: A research and teaching software tool for plane-parallel radiative transfer in the Earth’s atmosphere [J]. Bull Amer Meteor Soc, 1998, 79(10): 2101-2114.

[28]Zhai P M, Eskridge R E. Atmospheric water vapor over China [J]. J Climate, 1997, 10: 3643-3652.

[29]張瓊, 錢(qián)永甫. 用NCEP/NCAR再分析輻射資料估算月平均地表反照率[J]. 地理學(xué)報(bào), 1999, 54(4): 309-317.

Zhang Q, Qian Y F. Monthly mean surface albedo estimated from NCEP /NCAR reanalysis radiation data [J]. Acta Geographica Sinica, 1999, 54(4): 309-317.

責(zé)任編輯龐旻

Error in Calculation of Surface Radiation Based on SBDART Radiative Transfer Model

LI Xiu-Zhen1, SHENG Li-Fang2, 3, LIU Qian2, LIU Yu-Zhi4, WANG Li-Qun1

(1.China Satellite Maritime Tracking and Control Department, Jiangyin 214431, China； 2.College of Oceanic and Atmospheric Sciences, Ocean University of China, Qingdao 266100, China； 3.The Key Laboratory of Physical Oceanography, Ministry of Education, Ocean-Atmosphere Interaction and Climate Laboratory, Ocean University of China, Qingdao 266100, China； 4.College of Atmospheric Sciences of Lanzhou University, Lanzhou 730020, China)

Solar irradiance was calculated by SBDART radiative transfer model using ground-based aerosol optical data from POM-02 pyranometer, and the source of error in calculation of surface radiation was analyzed by comparing with observational solar irradiance from CM21 radiometer. It is shown that the sensitivity of SSA and water vapor to surface radiation is strong, when these two factors have a 50% decrease, the surface radiation will increase 56.99 W/m2and decrease 25.45 W/m2respectively. Changing atmospheric profile with accurate water vapor, the difference can be very small. According to correlation and multiple linear regression analysis, solar zenith angle, SSA and AOD are main sources of errors in solar irradiance calculation, the results of regression equation based on these 3 factors has a positive correlation (0.486) with the relative errors of surface radiation. With regression equation the results can be revised better.

radiative transfer model； radiosonde data； pyranometer； aerosol optical depth

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

2015-08-05；

2015-12-17

李秀鎮(zhèn)(1985-)，男，工程師。E-mail: zqlixiuzhen@163.com

**通訊作者：E-mail：shenglf@ouc.edu.cn

P422

A

1672-5174(2016)08-013-06

10.16441/j.cnki.hdxb.20150279

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