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

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

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

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

引用格式：李秀鎮(zhèn), 盛立芳, 劉騫， 等. 基于SBDART輻射傳輸模式的晴天地面總輻射模擬誤差分析[J]. 中國海洋大學(xué)學(xué)報(自然科學(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.

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

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

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

1　資料與方法

1.1 氣溶膠和地面總輻射資料來源

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

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

1.2 模式介紹

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

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

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

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

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

τ=aλ-n。

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

1.3 誤差分析方法

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

2　結(jié)果與分析

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

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

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

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

2.2 大氣廓線對模擬結(jié)果的影響

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

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

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

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

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

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

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

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

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

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

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

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

通過以上分析，確定回歸方程為:

y=-0.09964-0.001938×x1+0.2333×

x2-0.1552×x3。

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

3　結(jié)論與討論

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

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

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

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

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

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

[1]左大康, 王懿賢, 陳建綏. 中國地區(qū)太陽總輻射的空間分布特征[J]. 氣象學(xué)報, 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]余予, 夏祥鰲, 陳洪濱, 等. 晴空大氣太陽短波輻射觀測與模式比較[J]. 太陽能學(xué)報, 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)維俊, 陳洪濱, 郭文利, 等. 兩種反演地表凈太陽輻射的參數(shù)化模式在北京地區(qū)的比較[J]. 氣象學(xué)報, 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]王明歡, 賴安偉, 陳正洪, 等. WRF模式模擬的地表短波輻射與實(shí)況對比分析[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ì)算我國東南地區(qū)的太陽直接輻射和散射輻射[J]. 南京氣象學(xué)院學(xué)報, 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)星資料的晴空凈太陽輻射模擬[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]沈鐘平, 張華. 影響地面太陽輻射及其譜分布的因子分析[J]. 太陽能學(xué)報, 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]張瓊, 錢永甫. 用NCEP/NCAR再分析輻射資料估算月平均地表反照率[J]. 地理學(xué)報, 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

國家自然科學(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)