改進(jìn)的單次散射相函數(shù)解析表達(dá)式?

程晨 史澤林 崔生成 徐青山

1)(中國(guó)科學(xué)院安徽光學(xué)精密機(jī)械研究所,中國(guó)科學(xué)院大氣成分與光學(xué)重點(diǎn)實(shí)驗(yàn)室,合肥 230031)

2)(中國(guó)科學(xué)技術(shù)大學(xué)研究生院科學(xué)島分院,合肥 230031)

3)(中國(guó)科學(xué)院沈陽(yáng)自動(dòng)化研究所,沈陽(yáng) 110016)

改進(jìn)的單次散射相函數(shù)解析表達(dá)式?

程晨1)2)史澤林3)崔生成1)徐青山1)?

1)(中國(guó)科學(xué)院安徽光學(xué)精密機(jī)械研究所,中國(guó)科學(xué)院大氣成分與光學(xué)重點(diǎn)實(shí)驗(yàn)室,合肥 230031)

2)(中國(guó)科學(xué)技術(shù)大學(xué)研究生院科學(xué)島分院,合肥 230031)

3)(中國(guó)科學(xué)院沈陽(yáng)自動(dòng)化研究所,沈陽(yáng) 110016)

散射,相函數(shù),Henyey-Greenstein相函數(shù),新的Henyey-Greenstein相函數(shù)

1 引 言

大氣輻射傳輸模型是大氣科學(xué)研究的主要課題之一,它是聯(lián)系大氣物理屬性和輻射屬性的重要途徑.大氣輻射傳輸模型的主要研究方法有解析數(shù)值法和統(tǒng)計(jì)模擬法,研究過(guò)程中需要考慮大氣光學(xué)厚度、氣體成分的單次散射反照率以及粒子的單次散射相函數(shù),其中散射相函數(shù)的準(zhǔn)確性和復(fù)雜程度直接影響到輻射傳輸模型的計(jì)算結(jié)果和計(jì)算效率[1?4].因此,對(duì)相函數(shù)的研究主要集中在縮短相函數(shù)的計(jì)算時(shí)間以及提高相函數(shù)的準(zhǔn)確性.

目前常用于輻射傳輸模擬過(guò)程中的相函數(shù)計(jì)算方法有:Legendre多項(xiàng)式展開(kāi)法、Mie理論遞推公式計(jì)算方法、相函數(shù)調(diào)整方法以及Henyey-Greenstein(H-G)相函數(shù)近似法[5,6].Legendre多項(xiàng)式展開(kāi)法主要用于離散縱坐標(biāo)法求解輻射傳輸問(wèn)題中,對(duì)于球形小粒子,相函數(shù)只需要很少的展開(kāi)項(xiàng)就可以求解出它的散射問(wèn)題;對(duì)于非球形粒子及大粒子,它們的前向散射峰異常尖銳,需要數(shù)千項(xiàng)的展開(kāi)式才能收斂,且計(jì)算復(fù)雜、計(jì)算量巨大[7,8].Mie散射計(jì)算過(guò)程中將粒子假設(shè)為球體,對(duì)Mie理論直接編程時(shí)涉及到復(fù)變量Bessel函數(shù)的計(jì)算,因?yàn)榱Ｗ映叨认拗?即使粒子為球形,散射角度也不能由幾何光學(xué)或者格林公式[9]獲得的漸近表達(dá)式足夠準(zhǔn)確的近似[10].相函數(shù)調(diào)整方法的基本思想是將相函數(shù)近似為狄拉克Delta函數(shù)和截?cái)嘞嗪瘮?shù)的線性組合,并通過(guò)相似性原理[11,12]得到新的大氣光學(xué)參數(shù)[13].以上3種相函數(shù)計(jì)算方法都存在計(jì)算時(shí)間長(zhǎng)、效率低的問(wèn)題.為減小計(jì)算時(shí)間、提高計(jì)算效率,研究粒子的輻射傳輸特性時(shí)通常使用相函數(shù)的近似解析表達(dá)式.H-G相函數(shù)是最普遍適用的相函數(shù)近似解析表達(dá)式[14],該相函數(shù)公式簡(jiǎn)單、數(shù)值計(jì)算方便,便于求解出關(guān)于散射角隨機(jī)抽樣的逆函數(shù),可方便地應(yīng)用于蒙特卡羅方法關(guān)于散射方向抽樣的問(wèn)題.但該近似方法仍然存在諸多局限,例如不能很好地描述后向散射問(wèn)題[15].為了更好地描述不同的物理模型及物理問(wèn)題提出了多種改進(jìn)的H-G相函數(shù)近似方法,其中為了能很好地模擬后向散射問(wèn)題,Cornette和Shanks[16]利用Rayleigh相函數(shù)與H-G相函數(shù)推導(dǎo)出Henyey-Greenstein*(H-G*)相函數(shù).盡管如此,有些情況下H-G*相函數(shù)仍不能很好地模擬真實(shí)的散射相函數(shù).為了能夠得到最大程度符合粒子散射特性的散射相函數(shù),結(jié)合Rayleigh相函數(shù)和H-G相函數(shù),本文提出了一種新的散射相函數(shù)解析式(RH-G*相函數(shù)).為了驗(yàn)證RH-G*相函數(shù)在描述單個(gè)粒子及多分散系粒子散射特性的適用性,綜合比較了單個(gè)粒子的RH-G*相函數(shù)與H-G相函數(shù)、H-G*相函數(shù)之間的偏離程度,并分析了多種氣溶膠模式下RH-G*相函數(shù)與多分散系Mie散射相函數(shù)、H-G*相函數(shù)的吻合度,為更加準(zhǔn)確、高效地計(jì)算輻射傳輸特性提供理論支持.

2 散射相函數(shù)

散射相函數(shù)描述的是入射能量經(jīng)粒子散射后的能量角分布.入射方向與散射方向之間的夾角θ稱為散射角,為了更好地描述散射能量隨散射角的變化,有必要引入一個(gè)稱之為相函數(shù)的無(wú)量綱參數(shù)p(cosθ),

按此定義,相函數(shù)是歸一化的.按照定義計(jì)算粒子的散射相函數(shù),可以保證計(jì)算結(jié)果的精度,但計(jì)算效率不夠.為此,研究過(guò)程中常采用相函數(shù)的解析表達(dá)式來(lái)近似計(jì)算散射相函數(shù).相函數(shù)解析表達(dá)式的優(yōu)點(diǎn)是大大提高了計(jì)算效率,且在一定程度上符合實(shí)際的散射規(guī)律,但總體來(lái)說(shuō)計(jì)算精度稍有不足[6,16].這就要求發(fā)展更加適合實(shí)際散射規(guī)律的相函數(shù)解析表達(dá)式.

2.1 Mie散射相函數(shù)

借助于振幅散射矩陣,球形粒子的散射相函數(shù)p表示為[17]

其中,k是波數(shù);是散射截面,X是矢量散射振幅,?是立體角;an和bn為復(fù)函數(shù);S1和S2為復(fù)振幅.S1和S2的計(jì)算公式如下:

其中πn和τn為角度依賴函數(shù),可以通過(guò)遞推關(guān)系求解.對(duì)于滿足f(r)分布的多分散粒子系,散射相函數(shù)為[18]

其中p(λ,r,m,θ)是已知粒子尺寸參數(shù)及復(fù)折射率的單個(gè)粒子散射相函數(shù),λ是波長(zhǎng),r是粒子半徑,m是復(fù)折射率.

2.2 H-G相函數(shù)

Mie散射相函數(shù)的計(jì)算過(guò)程比較復(fù)雜,為了計(jì)算方便,常用表達(dá)形式簡(jiǎn)單、數(shù)值計(jì)算方便的H-G相函數(shù)來(lái)替代[10]:

其中θ是散射角,g是不對(duì)稱因子.H-G相函數(shù)可以較好地模擬前向散射峰,但是不能正確模擬后向散射,且當(dāng)g→0時(shí),該相函數(shù)不能還原成瑞利散射相函數(shù):

2.3 雙H-G相函數(shù)

為了克服H-G相函數(shù)不能正確模擬后向散射這個(gè)缺陷,Cornette和Shanks[16]提出使用雙H-G相函數(shù)(DH-G相函數(shù))來(lái)近似:

DHG相函數(shù)雖然可以較好地體現(xiàn)粒子的散射特征,但是解析式是由3個(gè)參數(shù)確定的:f是一個(gè)正的參數(shù);g1是正值,表示前向散射部分;g2是負(fù)值,表示后向散射部分,使用起來(lái)比較復(fù)雜,不適合實(shí)際應(yīng)用.當(dāng)g→0時(shí),為了將雙H-G相函數(shù)還原成瑞利散射相函數(shù),前向散射峰和后向散射峰需要相等,即g1=?g2,且f=1/2,但此時(shí)兩者之間仍不能很好地符合,最小均方根誤差為0.15[16].且H-G相函數(shù)、DH-G相函數(shù)只用于粒子的尺寸參數(shù)不是太小的情況.

2.4 H-G*相函數(shù)

為了能很好地模擬后向散射問(wèn)題,Cornette和Shanks[16]利用Rayleigh相函數(shù)與H-G相函數(shù)推導(dǎo)出了修正的H-G相函數(shù),即H-G*相函數(shù):

該相函數(shù)在一定程度上彌補(bǔ)了H-G相函數(shù)和DHG相函數(shù)的不足,且只有一個(gè)自由參數(shù)g,在g→0時(shí),與瑞利散射相函數(shù)符合;當(dāng)g→1時(shí),近似于H-G相函數(shù).

3 本文提出的新的散射相函數(shù)

比較H-G*和H-G相函數(shù)會(huì)發(fā)現(xiàn),H-G*在前向散射部分修正較多,而在后向散射部分修正較少,不能很好地模擬后向散射峰值比較大的粒子的散射特征,例如海洋型氣溶膠[18].為了更好地模擬后向散射峰值比較大的粒子的散射特征,對(duì)H-G相函數(shù)解析式與瑞利散射相函數(shù)解析式按照不對(duì)稱因子進(jìn)行修正.將新的相函數(shù)表示成H-G相函數(shù)和瑞利相函數(shù)以及不對(duì)稱因子的多項(xiàng)式的和,其中H-G相函數(shù)和瑞利相函數(shù)的系數(shù)是不對(duì)稱因子g的多項(xiàng)式.要求在角度大于90?時(shí),增大后向散射峰值,且在g→0時(shí),新的相函數(shù)與瑞利散射相函數(shù)符合;當(dāng)g→1時(shí),新的相函數(shù)近似于H-G相函數(shù).基于這些考慮,改進(jìn)的散射相函數(shù)(RH-G*)為

該相函數(shù)在一定程度上增大了后向散射峰,解析式形式簡(jiǎn)單,只有一個(gè)自由參數(shù)g.

圖1給出了不對(duì)稱因子從0到0.9(間隔為0.1)時(shí),H-G相函數(shù)、H-G*相函數(shù)和RH-G*相函數(shù)隨角度變化的數(shù)值結(jié)果.

圖1 (a)—(j)是RH-G*相函數(shù)與H-G*相函數(shù)、H-G相函數(shù)的角分布對(duì)比Fig.1.(a)–(j)Comparison of RH-G*phase function,H-G*phase function and H-G phase function.

由圖1可知:1)三種相函數(shù)角分布變化趨勢(shì)相同,表現(xiàn)在小角度的前向散射峰值較強(qiáng),后向散射峰值較弱,且隨著不對(duì)稱因子g的增大,都有前向散射峰增大、后向散射峰減小的趨勢(shì);2)相對(duì)于H-G相函數(shù)和H-G*相函數(shù),RH-G*相函數(shù)增大了后向散射的峰值;3)當(dāng)g=0時(shí),RH-G*相函數(shù)和H-G*相函數(shù)完全重合,但與H-G相函數(shù)偏離程度較大;當(dāng)g6 0.5時(shí),三種相函數(shù)的前向散射峰值顯現(xiàn)出很大的差異,但RH-G*相函數(shù)與H-G*相函數(shù)后向散射部分符合較好;當(dāng)g>0.5時(shí),則相反.總體來(lái)說(shuō),RH-G*相函數(shù)提高了后向散射峰值.

白璐等[19]提出了由于海洋型氣溶膠的后向散射峰值比較大,采用H-G相函數(shù)及H-G*相函數(shù)難以給出較好的擬合結(jié)果;朱孟真等[20]給出了常見(jiàn)或霾天氣條件下紫外光散射相函數(shù)的比較,結(jié)果顯示H-G相函數(shù)及H-G*相函數(shù)的數(shù)值在120?以后與Mie散射的計(jì)算結(jié)果差別比較大,Mie散射相函數(shù)有比較明顯的后向散射峰,說(shuō)明提高相函數(shù)解析式的后向散射峰值是很有必要的,RH-G*相函數(shù)能夠合理地描述單個(gè)粒子散射輻射能量隨角度的變化.

4 多種氣溶膠模式下Mie散射相函數(shù)與RH-G*相函數(shù)的對(duì)比

顯然,RH-G*相函數(shù)可以很好地描述單個(gè)粒子散射輻射能量的角分布特征.相函數(shù)還必須能夠用于輻射傳輸問(wèn)題的求解,因此需要能夠模擬多分散系粒子散射能量的角分布.

Derimendjian[21,22]利用Mie散射理論描述了三種氣溶膠類型(近似的大陸型霾C(Haze C)、海洋型霾M(Haze M)和積云C.1(Cloud C.1))的散射輻射能量的分布特征.海洋型霾M最初是為海洋線和沿海地區(qū)的氣溶膠設(shè)計(jì)的;積云C.1類似薄的積云.近似大陸型霾C是不連續(xù)分布的;海洋型霾M和積云C.1是按照指數(shù)型公式分布(廣義Γ分布):

其中,a表示粒子濃度;r是粒子半徑;α,b,γ是依氣溶膠種類不同而定的常數(shù).表1列舉了這三種氣溶膠類型尺寸分布的典型參數(shù).

表1 粒子尺度分布參數(shù)Table 1.Particle-size distributions.

非偏振多分散系Mie散射相函數(shù)計(jì)算過(guò)程選取的參數(shù)(波長(zhǎng)、粒子尺度、復(fù)折射率)如表2所列:λ是波長(zhǎng),Re是折射率實(shí)部,Im是折射率虛部,參數(shù)g是散射不對(duì)稱因子.

表2 粒子光學(xué)參數(shù)Table 2.Particle optical properties.

圖2—圖4給出了積云C.1、近似大陸型霾C和海洋型霾M三種模式下多分散系Mie散射相函數(shù)、RH-G*相函數(shù)和H-G*相函數(shù)的角分布.虛線表示Mie散射相函數(shù);虛線加點(diǎn)表示RH-G*相函數(shù);實(shí)線加星號(hào)表示H-G*相函數(shù).

圖2積云C.1模式中,RH-G*相函數(shù)和H-G*相函數(shù)的前向散射峰值很接近;波長(zhǎng)小于5.3μm時(shí),RH-G*相函數(shù)與多分散系Mie散射相函數(shù)在后向散射部分的符合程度高.圖3近似大陸型霾C模式中,總體而言RH-G*相函數(shù),對(duì)多分散系Mie散射相函數(shù)的擬合效果要優(yōu)于H-G*相函數(shù),尤其是在大角度后向散射上.圖4中波長(zhǎng)大于10μm時(shí),海洋型霾M模式的后向散射峰值采用RH-G*相函數(shù)能得到較好的擬合結(jié)果.

總之,1)RH-G*相函數(shù)彌補(bǔ)了H-G*相函數(shù)描述后向散射峰尖銳程度的不足,特別是短波部分;2)RH-G*相函數(shù)提高了與多分散系Mie散射相函數(shù)的符合程度,特別是在前向散射部分;3)由于粒子尺度通常小于λ=100μm很多,散射相函數(shù)近似于瑞利散射相函數(shù),此時(shí)RH-G*相函數(shù)精確地描述了這一現(xiàn)象,其隨角度變化的曲線幾乎與多分散系Mie散射相函數(shù)的曲線重合;4)對(duì)于后向散射特征比較明顯的氣溶膠模式,RH-G*相函數(shù)更接近于Mie散射相函數(shù).

圖2 (a)—(j)積云C.1的Mie散射相函數(shù)、H-G*相函數(shù)與RH-G*相函數(shù)的角度分布對(duì)比Fig.2.(a)–(j)Comparison of Mie phase function for the Cloud C.1 model,H-G*phase function and RH-G*phase function.

圖3 (a)—(j)近似大陸型霾C的Mie散射相函數(shù)、H-G*相函數(shù)與RH-G*相函數(shù)的角度分布對(duì)比Fig.3.(a)–(j)Comparison of Mie phase function for the Haze C model,H-G*phase function and RH-G*phase function.

圖4 (a)—(j)海洋型霾M的Mie散射相函數(shù)、H-G*相函數(shù)與RH-G*相函數(shù)的角度分布對(duì)比Fig.4.(a)–(j)Comparison of Mie phase function for the Haze M model,H-G*phase function and RH-G*phase function.

出現(xiàn)上述結(jié)果的原因如下:1)相同波段不同氣溶膠類型的相函數(shù)計(jì)算結(jié)果不同,這主要是因?yàn)楦鞣N氣溶膠類型的不對(duì)稱因子g存在差異,不對(duì)稱因子主要受粒子尺度因子a(a=πd/λ,d是粒子直徑,λ是入射波長(zhǎng))和光學(xué)常數(shù)折射率的影響,本文中三種氣溶膠類型的粒子尺度和折射率不同,所以計(jì)算結(jié)果存在差異,尺度因子越小,相函數(shù)隨角度變化曲線的對(duì)稱性越強(qiáng),所以波長(zhǎng)為100μm時(shí),各個(gè)曲線的對(duì)稱性是最強(qiáng)的;2)相函數(shù)解析式與Mie散射相函數(shù)的計(jì)算結(jié)果存在差異的原因是Mie散射相函數(shù)是由Mie散射參數(shù)和n階Legendre多項(xiàng)式表示的,Mie散射參數(shù)反映了電場(chǎng)振蕩和磁場(chǎng)振蕩對(duì)散射的影響,而相函數(shù)解析式中將Mie散射參數(shù)近似為與不對(duì)稱因子有關(guān)的多項(xiàng)式,所以RH-G*相函數(shù)和H-G*相函數(shù)隨角度的變化曲線相對(duì)比較平滑,且與Mie散射相函數(shù)的計(jì)算結(jié)果存在差異;3)H-G*相函數(shù)和RH-G*相函數(shù)在相同波段下計(jì)算結(jié)果存在差異時(shí)主要原因在于RH-G相函數(shù)是通過(guò)H-G相函數(shù)與不對(duì)稱因子和散射角有關(guān)的多項(xiàng)式的和來(lái)描述散射過(guò)程的,不對(duì)稱因子增大時(shí),RH-G*相函數(shù)在前向的修正較多,后向的修正較少,但是前向和后向的值都是增大;H-G*相函數(shù)是以不對(duì)稱因子和散射角有關(guān)的多項(xiàng)式作為H-G相函數(shù)的系數(shù)來(lái)描述散射過(guò)程的,不對(duì)稱因子越大,修正的越少,所以RH-G*相函數(shù)與H-G*相函數(shù)之間存在差異.

為了進(jìn)一步說(shuō)明RH-G*相函數(shù)提高大角度后向散射的合理性,表3比較了H-G*相函數(shù)、RH-G*相函數(shù)與多分散系Mie散射相函數(shù)在角度大于90?時(shí)的均方根差值.對(duì)比時(shí),角度間隔為5?,均方根誤差表示為是散射角θ處不同相函數(shù)解析式計(jì)算的值,pi(θ)是對(duì)應(yīng)散射角θ下Mie散射相函數(shù)的計(jì)算值.表中列舉出了三種氣溶膠模式中,RH-G*相函數(shù)與Mie散射相函數(shù)的均方根差值小于H-G*相函數(shù)與Mie散射相函數(shù)的有22個(gè),占總數(shù)的73.3%,說(shuō)明大部分情況下,RH-G*相函數(shù)在大角度后向散射上與Mie散射相函數(shù)的符合程度較好.因此,RH-G*相函數(shù)提高大角度后向散射是合理的、必要的.

表3 多分散系氣溶膠粒子H-G*相函數(shù)和RH-G*相函數(shù)與Mie散射相函數(shù)的方均根誤差比較Table 3.Comparisons of the root-mean-square error between the RH-G*phase function,the H-G*phase function and the Mie-sacttering phase function.

5 結(jié) 論

本文提出了一種新的單次散射相函數(shù)(RH-G*相函數(shù))的解析表達(dá)式,合理地提高了后向散射峰值,公式簡(jiǎn)單,計(jì)算方便.多種氣溶膠模式下,RH-G*相函數(shù)與多分散系Mie散射相函數(shù)的擬合效果優(yōu)于H-G*相函數(shù),尤其是后向大角度方向上.RH-G*相函數(shù)模擬后向散射峰值比較大的粒子的散射特征有明顯的效果,對(duì)準(zhǔn)確模擬輻射傳輸過(guò)程具有重要意義.但新的相函數(shù)不能很好地求解出散射角隨機(jī)抽樣的逆函數(shù),所以不能很好地用于蒙特卡羅方法中.

[1]Si F Q,Liu J G,Xie P H,Zhang Y J,Dou K,Liu W Q 2006Acta Phys.Sin.55 3165(in Chinese)[司福祺,劉建國(guó),謝品華,張玉鈞,竇科,劉文清2006物理學(xué)報(bào)55 3165]

[2]Xie Z L,Zhang R,Xiu X Q,Han P,Liu B,Chen L,Yu H Q,Jiang R L,Shi Y,Zhang Y D 2007Acta Phys.Sin.56 6717(in Chinese)[謝自力,張榮,修向前,韓平,劉斌,陳琳,俞慧強(qiáng),江若鏈,施毅,鄭有炓 2007物理學(xué)報(bào) 56 6717]

[3]Li Y,Song J,Zhang Z M 2003Chin.Phys.Lett.20 1189

[4]Li J N,Barker H,Yang P,Yi B Q 2015J.Geophys.Res.120 128

[5]Liu K,Zhang F,Wu K,Zhou X J,Zhang H,Liu R Q 2016Acta Meteorol.Sin.74 784(in Chinese)[劉琨, 張峰,吳琨,周秀驥,張華,劉仁強(qiáng)2016氣象學(xué)報(bào)74 784]

[6]Toublance D 1996Appl.Opt.35 3270

[7]Chen X H,Liu Q,Wei H L 2007J.Light Scatt.3 283(in Chinese)[陳秀紅,劉強(qiáng),魏合理 2007光散射學(xué)報(bào) 3 283]

[8]Ding S,Xie Y,Yang P,Weng F,Liu Q,Baum B,Hu Y 2009J.Quant.Spectrosc.Radiat.Transfer110 1190

[9]Cai H,Huang N N 2003Chin.Phys.Lett.20 496

[10]Xu L Q,Li H,Xiao Z Y 2008Acta Phys.Sin.57 6030(in Chinese)[徐蘭青,李暉,肖鄭穎 2008物理學(xué)報(bào) 57 6030]

[11]Liou K N(translated by Guo C L,Zhou S J)2004An Introduction to Atmospheric Radiation(2nd Ed.)(Beijing:China Meteorogical Press)pp108–109(in Chinese)[廖國(guó)男 著(郭彩麗,周詩(shī)健 譯)2004大氣輻射導(dǎo)論(第二版)(北京:氣象出版社)第108—109頁(yè)]

[12]Joseph J,Wiscombe W,Weinman J 1976J.Atmos.Sci.33 2452

[13]Wiscombe W 1977J.Atmos.Sci.34 1408

[14]Huang C J,Liu Y F,Wu Z S,Sun Y Q,Long S M 2009Acta Phys Sin.58 2397(in Chinese)[黃朝軍,劉亞峰,吳振森,孫彥清,龍姝明2009物理學(xué)報(bào)58 2397]

[15]Kattawar G W 1975Quant.Spectrosc.Radiat.Transfer15 839

[16]Cornette W M,Shanks J G 1992Appl.Opt.31 3152

[17]Bohre C F,Hu ff man D R 1983Absorption and Scattering of Light by Small Particles(New York:John Wiley)p100

[18]McCartney E J 1988Optics of the Atmosphere(Beijing:Science Press)p276(in Chinese)[麥克特尼 E J著 (潘乃先譯)1988大氣光學(xué)(北京:科學(xué)出版社)第276頁(yè)]

[19]Bai L,Tang S Q,Wu Z S,Xie P H,Wang S M 2010Acta Phys.Sin.59 1749(in Chinese)[白璐,湯雙慶,吳振森,謝品華,汪世美2010物理學(xué)報(bào)59 1749]

[20]Zhu M Z,Zhang H L,Jia H L,Yang J K,Cheng S L 2007J.Light Scatt.3 225(in Chinese)[朱孟真,張海良,賈紅輝,楊建坤,常勝利2007光散射學(xué)報(bào)3 225]

[21]Derimendjian D 1963Scattering and Polarization Properties of Polydispersed Suspensions with Partial Absorption(New York:The MacMillan Company,Pursuant to a Special Arrangement with Pergamon Press)pp4–9

[22]Diermendjian D 1969Electromagnetic Scattering on Spherical Polydispersions(New York:American Elsevier)pp31–78

Modi fied analytic expression for the single-scattering phase function?

Cheng Chen1)2)Shi Ze-Lin3)Cui Sheng-Cheng1)Xu Qing-Shan1)?

1)(Key Laboratory of Atmospheric Composition and Optical Radiation,Anhui Institute of Optics and Fine Mechanics,Chinese Academy of Sciences,Hefei 230031,China)

2)(Science Island Branch of Graduate School,University of Science and Technology of China,Hefei 230031,China)

3)(Shenyang Institute of Automation Chinese Academy of Sciences,Shenyang 110016,China)

6 April 2017;revised manuscript

11 May 2017)

In electromagnetic radiative transfer calculation,the accuracy and the computation time are usually determined by the representation of single-scattering phase function.Accurate calculation is time consuming even for spherical particle,thus,an analytic representation is commonly adopted to approximate the exact phase function and then accelerate the calculation.Most widely used single-scattering phase functions are the Henyey-Greenstein phase function and modi fied Henyey-Greenstein phase function(Henyey-Greenstein*).Although the Henyey-Greenstein phase function and the Henyey-Greenstein*phase function can represent the forward-scattering peak of Mie-scattering phase function well,they fail to reproduce the backscattering behavior,limiting the accuracy of the calculation.In order to better fit exact calculations and simulate the backward-scattering peak,we develop a new analytic expression based on the fundamental theory of electromagnetic scattering and radiation transmission.This phase function is an algebraic expression with one single free parameter(asymmetry factor),and can be expanded into Legendre polynomials.The new phase function converges to the Rayleigh phase function when the asymmetry factor approximates to 0,and it can approach to the Henyey-Greenstein phase function as the asymmetry factor is about 1.We compare the Henyey-Greenstein phase function,the Henyey-Greenstein*phase function,and the new phase function for different asymmetry factors,and f i nd that the new phase function provides a more realistic description for the unpolarized light scattering from small particles.Furthermore,the calculated value for the ratio of the scattering intensity at 90 degree to that in the backward direction is more reasonable.We also investigate the effectiveness by approximating the scattering from polydispersed particles through comparing the new phase function,the Henyey-Greenstein*phase function,and the Mie-scattering phase function for three types of Derimendjian’s polydispersions.Results show that the new phase function fits the Mie-scattering phase function much better than the Henyey-Greenstein*phase function.For the new phase function,the root-mean-square error is small for 73.3%data.By contrast,only 26.7%data fit the Mie-scattering phase function well for the Henyey-Greenstein*phase function.Similarly,the effectiveness of new function is most significant when calculating the ratio of the scattering intensity at 90 degree to that in the backward direction.In summary,the new Henyey-Greenstein*phase function provides a more accurate calculation for the scattering intensity in the backward direction,and is conducive to electromagnetic radiative transfer calculation.Furthermore,because the proposed phase function has the same basic form as the Heny-Greenstein phase function,reformatting radiative transfer model in terms of the new phase function should require relatively little effort.

scattering,phase function,Henyey-Greenstein phase function,new Henyey-Greenstein phase function

PACS:02.30.–f,03.65.Nk,07.60.Rd,42.25.BsDOI:10.7498/aps.66.180201

*Projectsupported by the ChinaSpecial Fund for Meteorological Research in the Public Interest(Grant No.GYHY201106002-03)and the National Natural Science Foundation of China(Grant No.41576185).

?Corresponding author.E-mail:qshxu@aiofm.ac.cn

(2017年4月6日收到;2017年5月11日收到修改稿)

單次散射相函數(shù)對(duì)電磁輻射傳輸模擬過(guò)程的準(zhǔn)確性和計(jì)算效率有重要的影響.基于電磁散射與輻射傳輸中的基本理論,對(duì)單次散射相函數(shù)的解析表達(dá)式進(jìn)行了研究,提出了一種新的單次散射相函數(shù)解析表達(dá)式.比較了單個(gè)粒子的Henyey-Greenstein相函數(shù)、Henyey-Greenstein*相函數(shù)與新的相函數(shù)隨角度的分布,發(fā)現(xiàn)新的散射相函數(shù)提高了后向散射峰值,可以更合理地描述單個(gè)粒子的散射特性.按三種氣溶膠粒子譜分布模式計(jì)算了Henyey-Greenstein*相函數(shù)和新的相函數(shù)對(duì)應(yīng)的數(shù)值結(jié)果,并與多分散系Mie散射相函數(shù)進(jìn)行對(duì)比,發(fā)現(xiàn)新的相函數(shù)提高了與多分散系Mie散射相函數(shù)的符合程度.研究表明,對(duì)于大角度(大于90?)后向散射,新的相函數(shù)與Mie散射相函數(shù)均方根差較小的占73.3%,高于Henyey-Greenstein*相函數(shù)的26.7%,證明了新的相函數(shù)可以顯著提高后向散射峰值.新的相函數(shù)對(duì)準(zhǔn)確模擬輻射傳輸過(guò)程具有重要意義.

10.7498/aps.66.180201

?氣象專項(xiàng)課題(批準(zhǔn)號(hào):GYHY201106002-03)和國(guó)家自然科學(xué)基金(批準(zhǔn)號(hào):41576185)資助的課題.

?通信作者.E-mail:qshxu@aiofm.ac.cn