電離層LBH日輝輻射大視場(chǎng)計(jì)算方法

張永超，何飛，張效信，陳波

1中國(guó)科學(xué)院長(zhǎng)春光學(xué)精密機(jī)械與物理研究所，長(zhǎng)春 130033

2中國(guó)科學(xué)院大學(xué)，北京 100049

3中國(guó)氣象局國(guó)家空間天氣監(jiān)測(cè)與預(yù)警中心，北京 100081

1 引言

電離層LBH日輝輻射是由光電子與氮?dú)夥肿优鲎布ぐl(fā)而產(chǎn)生的，是電離層在遠(yuǎn)紫外輻射波段中最強(qiáng)的分子輻射信號(hào).從空間對(duì)此信號(hào)進(jìn)行觀測(cè)，不僅可以監(jiān)測(cè)高層大氣中氮?dú)猓ɑ蜓鯕猓┓肿拥拿芏龋∕eier，1991），還可以獲得電離層中光電子通量的信息（Oran and Strickland，1978），從而為高層大氣狀態(tài)的監(jiān)測(cè)提供了一種強(qiáng)有力的方法.因此，LBH輻射被廣泛用作星載遙感儀器的觀測(cè)目標(biāo)，比如DMSP衛(wèi) 星 上 的SSUSI（Paxton et al.，1992），TIMED衛(wèi)星上的GUVI（Paxton et al.，1999），以及IMAGE衛(wèi)星上的WIC（Mende et al.，2000）等.中國(guó)風(fēng)云三號(hào)（FY－3）衛(wèi)星上即將搭載廣角極光成像儀（視場(chǎng)角為130°），首次對(duì)極光LBH輻射和電離層LBH日輝輻射進(jìn)行探測(cè)，從而對(duì)極光活動(dòng)及電離層的空間天氣狀況進(jìn)行監(jiān)測(cè)預(yù)報(bào).

準(zhǔn)確理解LBH輻射產(chǎn)生的物理過(guò)程，從而建立合理的計(jì)算模型，這是獲得可靠的地球物理參數(shù)信息的重要保證.Mcewen等人（1966）率先討論了LBH波帶的輻射強(qiáng)度分布，但只給出了初步的理論計(jì)算公式.Conway（1992）研究了氮?dú)夥肿拥淖晕招?yīng)對(duì)LBH波帶輻射傳輸過(guò)程的影響，指出氧氣分子的吸收作用遠(yuǎn)大于氮?dú)夥肿拥淖晕兆饔?Dashkevich等人（1993），Eastes（2000）討論了氮?dú)夥肿蛹ぐl(fā)態(tài)之間的耦合效應(yīng)對(duì)LBH波帶輻射的影響.Strickland等（1995），Evans等（1995）在前人工作的基礎(chǔ)上，給出了比較通用的LBH日輝輻射的計(jì)算模型.

當(dāng)前，紫外波段的輻射傳輸模型及相應(yīng)計(jì)算代碼的研制仍然是國(guó)際上研究的熱點(diǎn)（Huffman，1992）.然而，能夠用于遠(yuǎn)紫外波段輻射的計(jì)算代碼很少，比較通用的只有AURIC（Atmospheric Ultraviolet Radiance Integrated Code），即大氣紫外輻射積分代碼，它是由美國(guó)計(jì)算物理公司CPI（Computational Physics，Inc.）與空軍Phillips實(shí)驗(yàn)室聯(lián)合開(kāi)發(fā)的

（Link et al.，1992；Strickland et al.，1992）.AURIC在計(jì)算LBH日輝柱輻射率時(shí)將沿觀測(cè)視線LOS方向的太陽(yáng)天頂角作為常數(shù)，Strickland等人（1999）指出AURIC模式不適用于太陽(yáng)天頂角大于90°的觀測(cè)情況.AURIC僅僅適用于SSUSI等瞬時(shí)掃描型小視場(chǎng)的LBH日輝輻射計(jì)算.以DMSP／SSUSI和TIMED／GUVI的11.8°瞬時(shí)掃描視場(chǎng)為例，設(shè)觀測(cè)點(diǎn)高度為830km，觀測(cè)視場(chǎng)投影至155km海拔處的高層大氣時(shí)，觀測(cè)視場(chǎng)內(nèi)太陽(yáng)天頂角的最大差值約是1.2°（計(jì)算過(guò)程見(jiàn)附錄A）.然而，當(dāng)瞬時(shí)掃描視場(chǎng)比較大時(shí)，對(duì)于FY－3上的廣角極光成像儀，視場(chǎng)角是130°，在同樣的觀測(cè)位置，觀測(cè)視場(chǎng)內(nèi)太陽(yáng)天頂角的最大差值約達(dá)50°（見(jiàn)附錄A），就需要考慮太陽(yáng)天頂角變化帶來(lái)的影響，而不能再將其作為常數(shù)來(lái)處理.

為了滿足大視場(chǎng)LBH日輝輻射特性研究的需要，本文給出了一種改進(jìn)的LBH日輝柱輻射率計(jì)算方法RAURIC（Revised algorithm of AURIC）.該方法考慮了觀測(cè)視線方向上太陽(yáng)天頂角的變化，克服了AURIC對(duì)LBH日輝柱輻射率計(jì)算的局限性.本文第2節(jié)首先對(duì)LBH輻射的譜帶結(jié)構(gòu)作了簡(jiǎn)要分析，然后詳細(xì)描述了LBH日輝柱輻射率的計(jì)算方法.第3節(jié)計(jì)算了140～180nm波段的LBH日輝柱輻射率，同時(shí)將RAURIC計(jì)算結(jié)果與AURIC進(jìn)行了比較，然后對(duì)比較結(jié)果作了詳細(xì)討論.最后，第4節(jié)對(duì)本文工作進(jìn)行了簡(jiǎn)要總結(jié).

2 LBH輻射計(jì)算

2.1 LBH譜帶分析

LBH帶（Lyman－Birge－Hopfield Bands）輻射，如（1）式所示，是氮?dú)夥肿佑杉ぐl(fā)態(tài)（α1Πg）至基態(tài)的躍遷而產(chǎn)生的，此為電偶極禁止躍遷.由選擇定則可知，在這兩個(gè)態(tài)之間可以發(fā)生的是電四極矩躍遷和磁偶極矩躍遷（Conway，1992）.每個(gè)LBH帶由（O，P，Q，R，S）五個(gè)分支組成，其中（P，Q，R）三個(gè)分支是由磁偶極矩躍遷與電四極矩躍遷混合組成的，對(duì)應(yīng)于ΔJ＝－1，0，1；其余兩個(gè)分支（O，S）則為純電四極矩躍遷，對(duì)應(yīng)于ΔJ＝－2，2.

式中，fλ表示某個(gè)LBH帶占總的LBH帶系的輻射組分，j（z）表示海拔高度z處于α1Πg態(tài)氮?dú)夥肿拥漠a(chǎn)生速率.fλ采用Ajello和Shemansky（1985）的數(shù)據(jù)值.

處于α1Πg態(tài)的氮?dú)夥肿拥漠a(chǎn)生速率j（z），可采用電子碰撞直接激發(fā)理論來(lái)計(jì)算，其值等于g因子與氮?dú)夥肿用芏鹊某朔e，即

式中，φ（E，z）代表海拔高度z處的電子通量，φ（E，z）采用Strickland等人（1999），Oran和Strickland（1978）的光電子通量模型來(lái)計(jì)算.σex（E）表示α1Πg態(tài)的激發(fā)截面，本文所采用的α1Πg態(tài)的激發(fā)截面數(shù)據(jù)（Ajello and Shemansky，1985）如圖1所示；E0表示α1Πg態(tài)的激發(fā)勢(shì)能閾值（E0＝8.7eV）（Ajello and Shemansky，1985），Emax表示所考慮的光電子能量的上限（Emax＝800eV）（Strickland et al.，1999）.

圖1 α1Πg態(tài)的激發(fā)截面Fig.1 Excitation cross sections of theα1Πgstate

2.2 LBH輻射計(jì)算方法

從空間對(duì)LBH日輝輻射進(jìn)行觀測(cè)，觀測(cè)量為L(zhǎng)BH日輝柱輻射率.某一觀測(cè)視線LOS方向的柱輻射率Iλ，可由（5）式表達(dá)（單位為Rayleigh）：

式中，Zl，Zu分別代表輻射產(chǎn)生區(qū)域的下上邊界；μ表示觀測(cè)角的余弦值，即μ＝cosθ，觀測(cè)角θ的定義如圖2所示，即觀測(cè)視線LOS方向與當(dāng)?shù)靥斓追较蛑g的夾角（注：此處與AURIC中定義的觀測(cè)角恰好互為補(bǔ)角）；jλ表示某個(gè)波長(zhǎng)的體輻射率；tλ表示氧氣分子對(duì)某個(gè)LBH波長(zhǎng)吸收的光學(xué)厚度；表示氧氣分子對(duì)某個(gè)LBH波長(zhǎng)的吸收截面.

圖2 觀測(cè)幾何Fig.2 Observation geometry

LBH體輻射率的計(jì)算在2.1節(jié)中已經(jīng)討論過(guò)，下一步需要計(jì)算的是LBH帶光子由發(fā)射點(diǎn)傳輸至觀測(cè)點(diǎn)處而不被散射的幾率.LBH帶光子的輻射傳輸過(guò)程，主要考慮氧氣分子的吸收，而忽略氮?dú)夥肿幼陨淼奈眨吹獨(dú)夥肿訉?duì)于LBH帶光子是“光學(xué)薄”的（Conway，1992）.因此LBH帶光子的散射幾率，可通過(guò)計(jì)算氧氣分子的吸收光學(xué)厚度來(lái)解決.本文使用的氧氣分子吸收截面的數(shù)據(jù)（Strickland et al.，1999）如圖3所示.

電離層LBH日輝柱輻射率的具體計(jì)算步驟，如圖4所示.

圖3 氧氣分子的吸收截面Fig.3 Absorption cross sections of molecular oxygen

第一步，確定輸入?yún)?shù).輸入?yún)?shù)包括：日期（Year、Day）、時(shí)間（UT）、地磁活動(dòng)參數(shù)（Ap）、太陽(yáng)10.7cm射電通量值（F10.7）及其81天平均值（〈F10.7〉）.Year、Day、UT、Ap、F10.7和〈F10.7〉用于計(jì)算高層大氣密度分布.為了研究不同太陽(yáng)活動(dòng)下LBH輻射的特征，我們分別選取了太陽(yáng)活動(dòng)極大年和太陽(yáng)活動(dòng)極小年兩種情況，如表1所示.

圖4 柱輻射率計(jì)算流程圖Fig.4 Flowchart of calculation of column emission rates

表1 大氣模型參數(shù)Table 1 Parameters of atmosphere model

第二步，確定LBH輻射產(chǎn)生區(qū)域的邊界.電離層LBH日輝輻射產(chǎn)生于向陽(yáng)側(cè)的高層大氣區(qū)域.因?yàn)楹０胃叨仍?0km以下的低層大氣對(duì)遠(yuǎn)紫外波段的LBH輻射具有強(qiáng)烈的吸收作用，同時(shí)考慮到500～600km的海拔高度時(shí)氮?dú)夥肿右呀?jīng)非常稀薄，從而可以忽略對(duì)LBH日輝輻射的貢獻(xiàn)，故選取Zl＝90km，Zu＝600km作為輻射產(chǎn)生區(qū)域的下、上邊界.在該區(qū)域內(nèi)，采用球面幾何描述的大氣模型，即使用球面對(duì)該區(qū)域進(jìn)行劃分.各個(gè)球面之間間隔的確定，根據(jù)高層大氣密度結(jié)構(gòu)的變化情況，采用“由密到疏”的原則給定（Bush and Chakrabarti，1995）.本文采用AURIC的劃分方法將90～600km范圍內(nèi)的高層大氣區(qū)域分為43層.

第三步，確定觀測(cè)視線LOS的方向.首先需要給定觀測(cè)點(diǎn)的位置坐標(biāo)（H，φ，λ），即觀測(cè)點(diǎn)的海拔高度，緯度及經(jīng)度坐標(biāo).然后需要給定觀測(cè)角的大小，但這時(shí)還不能確定觀測(cè)視線LOS的方向.因?yàn)椴捎玫氖乔驇缀未髿饽Ｐ?，還需要給定具體的觀測(cè)方位，才能確定某一觀測(cè)視線LOS的方向.因此，我們定義一個(gè)觀測(cè)方位角β：以觀測(cè)點(diǎn)的東向觀測(cè)方位為起始線，逆時(shí)針?lè)较蜣D(zhuǎn)向觀測(cè)視線LOS的方向，轉(zhuǎn)過(guò)的角度即觀測(cè)方位角β.

接下來(lái)，需要計(jì)算每個(gè)小柱體上下表面處的經(jīng)度與緯度坐標(biāo).本文將輻射區(qū)域分為43層，那么沿觀測(cè)視線LOS方向就劃分為43個(gè)小柱體.因?yàn)檩椛鋮^(qū)域的分層已經(jīng)完成，所以每個(gè)小柱體上下表面處的海拔高度就確定了.

對(duì)于每個(gè)小柱體上下表面處的經(jīng)度與緯度坐標(biāo)的確定，本文采用的是垂直方位投影技術(shù).此處，以小柱體上表面處為例說(shuō)明經(jīng)度與緯度的確定過(guò)程.由觀測(cè)點(diǎn)向小柱體上表面對(duì)應(yīng)的海拔高度處作一垂直投影平面，根據(jù)觀測(cè)角的大小，以及觀測(cè)視線LOS的方向，可以確定小柱體上表面處在投影平面中所對(duì)應(yīng)的水平坐標(biāo)與垂直坐標(biāo)（x，y）.這時(shí)根據(jù)垂直投影坐標(biāo)的轉(zhuǎn)換公式（詳見(jiàn)附錄B），可以將上表面處的水平坐標(biāo)與垂直坐標(biāo)（x，y）轉(zhuǎn)換為緯度與經(jīng)度坐標(biāo)（Glat，Glon）.對(duì)于小柱體下表面處的經(jīng)度與緯度坐標(biāo)的確定，同樣遵循這個(gè)過(guò)程.

當(dāng)小柱體上下表面處的位置坐標(biāo)確定之后，還需要確定的是每個(gè)小柱體沿觀測(cè)視線LOS方向的觀測(cè)角的大小.本文采用的是球幾何大氣模型，根據(jù)觀測(cè)角的定義，沿觀測(cè)視線LOS方向的各個(gè)小柱體處當(dāng)?shù)氐奶斓追较蚴亲兓?，所以?duì)應(yīng)的觀測(cè)角是不同的.對(duì)于小柱體上下表面處的觀測(cè)角的計(jì)算，由觀測(cè)點(diǎn)的坐標(biāo)與小柱體上下表面處的坐標(biāo)，可以確定觀測(cè)視線LOS方向的矢量坐標(biāo)，然后由小柱體上下表面處的垂直方向的矢量坐標(biāo)，計(jì)算出小柱體上下表面處的觀測(cè)角μ的大小.

第四步，調(diào)用大氣模型MSISE－00計(jì)算每個(gè)小柱體的密度參數(shù).當(dāng)小柱體上下表面處的經(jīng)度、緯度坐標(biāo)確定之后，（Glat，Glon）連同表1中給定參數(shù)以及該大氣區(qū)域的上下邊界及分層數(shù)作為輸入?yún)?shù)，調(diào)用大氣模型MSISE－00計(jì)算出小柱體上下表面處的氮?dú)夥肿用芏圈袾2（z）與氧氣分子密度ρO2（z）.

第五步，計(jì)算每個(gè)小柱體的體輻射率.對(duì)于每個(gè)小柱體體輻射率的計(jì)算，我們調(diào)用AURIC體輻射率計(jì)算模塊來(lái)完成.與AURIC原算法不同的是，此處每個(gè)小柱體體輻射率的計(jì)算中均考慮了太陽(yáng)天頂角的變化.AURIC在計(jì)算沿觀測(cè)視線LOS方向的柱輻射率時(shí)將太陽(yáng)天頂角作為常數(shù)，這種處理方式只有當(dāng)觀測(cè)視線LOS方向?yàn)樘斓追较驎r(shí)是完全合理的.因此，對(duì)于每個(gè)小柱體體輻射率的計(jì)算，均調(diào)用AURIC體輻射率計(jì)算模塊來(lái)計(jì)算天底方向的對(duì)應(yīng)值，這樣不同天底方向上的小柱體所對(duì)應(yīng)的太陽(yáng)天頂角也是不同的.每個(gè)小柱體的體輻射率計(jì)算完成之后，輸入對(duì)應(yīng)的fλ數(shù)據(jù)，根據(jù)（2）式，計(jì)算出對(duì)應(yīng)的LBH帶的體輻射率jλ.

對(duì)于每個(gè)小柱體參數(shù)的確定均采用梯形積分法，即，每個(gè)小柱體的氮?dú)夥肿用芏圈袾2，氧氣分子密度ρO2，觀測(cè)角μ以及體輻射率jλ均為每個(gè)小柱體上下表面處對(duì)應(yīng)的平均值.

第六步，輻射傳輸計(jì)算.對(duì)于LBH帶的輻射傳輸過(guò)程，只考慮氧氣分子的吸收.已知氧氣分子密度及吸收截面數(shù)據(jù)，根據(jù)（6）式，可計(jì)算出每個(gè)小柱體傳輸至觀測(cè)點(diǎn)處的光學(xué)厚度tλ.

最后，將確定的各個(gè)參數(shù)代入（5）式進(jìn)行計(jì)算，就可以得到某一觀測(cè)視線LOS方向的LBH日輝柱輻射率.

綜上所述，我們將這種基于AURIC代碼改進(jìn)的LBH日輝柱輻射率計(jì)算方法稱(chēng)為RAURIC模式.RAURIC模式對(duì)于AURIC輻射算法的局限性主要有兩點(diǎn)改進(jìn)：一是增加了觀測(cè)方位角，從而在觀測(cè)角相同時(shí)可對(duì)不同的觀測(cè)方位進(jìn)行計(jì)算，解決了AURIC計(jì)算結(jié)果對(duì)稱(chēng)的局限性；二是考慮了沿觀測(cè)視線LOS方向上太陽(yáng)天頂角的變化，這對(duì)于LBH帶的輻射計(jì)算非常重要，因?yàn)長(zhǎng)BH輻射強(qiáng)度與太陽(yáng)天頂角密切相關(guān)（Huffman，1992），從而使得RAURIC模式適用于大視場(chǎng)觀測(cè)的情況.

3 LBH輻射算法的驗(yàn)證

為驗(yàn)證LBH輻射算法，本文應(yīng)用RAURIC模式計(jì)算了電離層140～180nm波段的LBH日輝輻射，計(jì)算結(jié)果包括：天底方向的柱輻射率；其他觀測(cè)方向的體輻射率與柱輻射率.同時(shí)，對(duì)RAURIC與AURIC的計(jì)算結(jié)果進(jìn)行了比較、分析與討論.

3.1 天底方向

天底方向柱輻射率的計(jì)算結(jié)果如圖5所示，左右兩欄分別對(duì)應(yīng)表1中所給定的兩組參數(shù)所描述的大氣模型，觀測(cè)點(diǎn)的海拔高度為830km，經(jīng)度坐標(biāo)自上而下依次為0°，20°，50°，緯度坐標(biāo)的范圍均為－80°～80°（采用地理坐標(biāo)系）.

由圖5可以看出，不論是在太陽(yáng)活動(dòng)與地磁活動(dòng)較強(qiáng)的情況下（左邊一欄），還是在太陽(yáng)活動(dòng)與地磁活動(dòng)較弱的情況下（右邊一欄），RAURIC與AURIC的LBH柱輻射率計(jì)算結(jié)果在天底方向上都具有很好的一致性.

3.2 其他觀測(cè)方向

其他觀測(cè)方向的柱輻射率計(jì)算結(jié)果如圖6所示，左右兩欄分別對(duì)應(yīng)Group 1與Group 2兩組參數(shù)所描述的大氣模型.觀測(cè)點(diǎn)的海拔高度為830km，緯度與經(jīng)度坐標(biāo)自上而下分別為（20°，20°），（35°，35°），（50°，50°），觀測(cè)角的變動(dòng)范圍均是0°～60°.本文采用的是球面幾何大氣模型，當(dāng)觀測(cè)角相同時(shí)（天底方向除外），還存在觀測(cè)方位的不同.因此，我們選取了東西南北四個(gè)方位的柱輻射率進(jìn)行計(jì)算，對(duì)應(yīng)的觀測(cè)方位角β分別為0°，180°，270°，90°.

從圖6中可以看出，在所有觀測(cè)位置，當(dāng)觀測(cè)角比較小時(shí)（大約20°），太陽(yáng)天頂角的變化比較小，RAURIC與AURIC的計(jì)算結(jié)果非常接近.但是，隨著觀測(cè)角的增大（大約40°），太陽(yáng)天頂角的變化也增大，RAURIC與AURIC的計(jì)算結(jié)果的區(qū)別越來(lái)越明顯.如圖6a，6b中，東、北方位的太陽(yáng)天頂角變化比西、南方位變化較大，所以與AURIC計(jì)算結(jié)果的區(qū)別就比較明顯；圖6c，6d中，RAURIC與AURIC計(jì)算結(jié)果的區(qū)別比圖6a，6b更明顯.隨著觀測(cè)角的進(jìn)一步增大，如圖6e，6f所示，四個(gè)方位的太陽(yáng)天頂角變化都非常顯著，RAURIC與AURIC的柱輻射率計(jì)算結(jié)果的差別就越來(lái)越大.圖5與圖6的結(jié)果反映出RAURIC模式與AURIC柱輻射率計(jì)算方法的區(qū)別，具體討論如下.

3.3 討論

電離層LBH日輝輻射是由光電子與氮?dú)夥肿优鲎布ぐl(fā)而產(chǎn)生的，LBH輻射與太陽(yáng)天頂角密切相關(guān)（Huffman，1992）.因此，在計(jì)算某一觀測(cè)視線LOS方向的LBH日輝柱輻射率時(shí)需要考慮太陽(yáng)天頂角的變化.為了說(shuō)明太陽(yáng)天頂角的變化對(duì)LBH日輝輻射的影響，我們選取觀測(cè)點(diǎn)位于海拔高度830km，緯度與經(jīng)度坐標(biāo)為（50°，50°）且觀測(cè)角為60°時(shí)，分別對(duì)東西南北四個(gè)方位的LBH體輻射率進(jìn)行了計(jì)算，如圖7所示.

由圖7可知，對(duì)LBH日輝體輻射率貢獻(xiàn)最大的海拔高度范圍大約是120～220km.同時(shí)，在此海拔高度范圍內(nèi)，AURIC與東西南北四個(gè)方位的LBH日輝體輻射率有明顯的區(qū)別，這是因?yàn)椋ㄔ诖烁叨确秶鷥?nèi)）四個(gè)方位的太陽(yáng)天頂角是不同的，詳見(jiàn)以下討論.

為了進(jìn)一步說(shuō)明太陽(yáng)天頂角對(duì)LBH日輝輻射的影響，我們選取固定海拔高度處來(lái)討論g因子與太陽(yáng)天頂角SZA的關(guān)系.如圖8所示，海拔高度分別為150km與180km處的g因子，太陽(yáng)天頂角的變化范圍為0°～90°，g因子通過(guò)（4）式進(jìn)行計(jì)算，圖8a，8b分別對(duì)應(yīng)Group 1與Group 2兩組參數(shù)所描述的大氣模型.

由圖8可知，在不同的海拔高度處，g因子隨太陽(yáng)天頂角的變化趨勢(shì)是一致的.同時(shí)，在太陽(yáng)天頂角小于40°時(shí)，g因子變化比較平緩；然而，當(dāng)太陽(yáng)天頂角大于40°時(shí)，g因子的變化就比較顯著了.接下來(lái)，我們說(shuō)明圖6e，6f中太陽(yáng)天頂角隨觀測(cè)角的變化情況，結(jié)果如圖9所示，其中縱坐標(biāo)為觀測(cè)視線LOS方向上下邊界處對(duì)應(yīng)的太陽(yáng)天頂角.

圖8 g因子與太陽(yáng)天頂角SZA的關(guān)系Fig.8 Relationship between g factor and SZA

圖9 觀測(cè)角與太陽(yáng)天頂角的關(guān)系Fig.9 Relationship between view angle and SZA

由圖9可以看出，AURIC將太陽(yáng)天頂角作為常數(shù)處理，因而隨觀測(cè)角的變化為直線；RAURIC模式考慮了太陽(yáng)天頂角的變化，當(dāng)觀測(cè)角比較小時(shí)，太陽(yáng)天頂角的變化亦不明顯，但是當(dāng)觀測(cè)角比較大時(shí)，太陽(yáng)天頂角在各個(gè)方位的變化就比較顯著了.當(dāng)觀測(cè)角為60°且觀測(cè)方位為東時(shí)，RAURIC與AURIC太陽(yáng)天頂角的差值大約是13°，而g因子相差近一個(gè)數(shù)量級(jí).因此，當(dāng)觀測(cè)視場(chǎng)比較小時(shí)，將太陽(yáng)天頂角作為常數(shù)處理是近似成立的；而當(dāng)觀測(cè)視場(chǎng)比較大時(shí)，不能將太陽(yáng)天頂角作為常數(shù)處理，就需要考慮太陽(yáng)天頂角的變化.

由以上討論可知，當(dāng)觀測(cè)方向?yàn)樘斓追较驎r(shí)，RAURIC與AURIC在天底方向上都將太陽(yáng)天頂角作為常數(shù)來(lái)處理，因此，二者對(duì)LBH柱輻射率計(jì)算結(jié)果是一致的（圖5）；在其他觀測(cè)方向上，RAURIC模式考慮了觀測(cè)視線LOS方向上太陽(yáng)天頂角的變化，當(dāng)觀測(cè)視場(chǎng)比較大時(shí)，RAURIC與AURIC計(jì)算結(jié)果的差別比較大（圖6），需要使用RAURIC模式進(jìn)行計(jì)算.

4 結(jié)論

本文針對(duì)大視場(chǎng)觀測(cè)模式，基于現(xiàn)有通用的大氣紫外輻射積分代碼AURIC模式，提出了一種改進(jìn)的電離層LBH日輝柱輻射率計(jì)算方法，即RAURIC模式.RAURIC模式對(duì)大視場(chǎng)且太陽(yáng)天頂角不大于90°時(shí)的體輻射率計(jì)算方法進(jìn)行了改進(jìn)，考慮了沿觀測(cè)視線方向上的太陽(yáng)天頂角變化對(duì)體輻射率的影響，并沒(méi)有修改AURIC模式中的光電子通量計(jì)算模塊.RAURIC模式在計(jì)算體輻射率時(shí)仍需要調(diào)用AURIC模式中原有的光電子通量計(jì)算模塊，而這一模塊僅僅適用于太陽(yáng)天頂角不大于90°的情況，因此RAURIC模式目前還只適用于太陽(yáng)天頂角不大于90°的電離層圓盤(pán)觀測(cè)時(shí)的LBH日輝輻射計(jì)算，而不適合太陽(yáng)天頂角大于90°的臨邊探測(cè)情況.

RAURIC模式為電離層LBH日輝圖像模擬技術(shù)與數(shù)據(jù)反演技術(shù)的研制奠定了基礎(chǔ)，從而為研究廣角極光成像儀的電離層LBH日輝圖像對(duì)于空間天氣事件（馬淑英等，2002；徐文耀，2009）的響應(yīng)情況以及極光圖像中的雜光去除技術(shù)（Strickland et al.，1994；Gladstone，1994）奠定了基礎(chǔ).

附錄A

如附圖1所示，地球半徑R＝6371.0km，觀測(cè)點(diǎn)A的海拔高度H＝830km.

附圖1 幾何示意圖Appendix Fig.1 Geometrical sketch map

（1）小視場(chǎng)對(duì)應(yīng)太陽(yáng)天頂角的變化

當(dāng)觀測(cè)點(diǎn)A的視場(chǎng)張角∠DAE為11.8°（DMSP／SSUSI和TIMED／GUVI兩臺(tái)儀器的瞬時(shí)掃描視場(chǎng)）時(shí)，∠OAE＝5.9°.觀測(cè)視線AD、AE與海拔高度為155km處的球面大氣相交時(shí)（典型的LBH帶投影高度），OD、OE對(duì)應(yīng)交點(diǎn)處的天底方向，因此D、E兩點(diǎn)之間的太陽(yáng)天頂角之差等于OD、OE之間的夾角.

在三角形ΔOAE中，根據(jù)正弦定理得

即D、E兩點(diǎn)之間的太陽(yáng)天頂角之差為1.2°.

（2）大視場(chǎng)對(duì)應(yīng)太陽(yáng)天頂角的變化

當(dāng)觀測(cè)點(diǎn)A的視場(chǎng)張角∠CAB為130°（FY－3氣象衛(wèi)星廣角極光成像儀的瞬時(shí)掃描視場(chǎng)）時(shí)，觀測(cè)視線AB、AC約與海拔高度為155km處的球面大氣相切；因?yàn)镃O、BO為C、B兩點(diǎn)處的天底方向，所以C、B兩點(diǎn)之間的太陽(yáng)天頂角之差等于CO、BO之間的夾角∠COB.計(jì)算過(guò)程如下：

因?yàn)橛^測(cè)視線AB、AC分別與BO、CO垂直，∠COB在四邊形ACOB中與∠CAB互為補(bǔ)角，所以∠COB≈50°，即C、B兩點(diǎn)之間的太陽(yáng)天頂角之差約為50°.

附錄B

垂直方位投影建立了地理經(jīng)緯度坐標(biāo)（Glat，Glon）與投影平面的水平方向及垂直方向坐標(biāo)（x，y）之間的轉(zhuǎn)換關(guān)系，以下是由投影坐標(biāo)（x，y）至經(jīng)緯度坐標(biāo)（Glat，Glon）的轉(zhuǎn)換方程（Snyder，1987），

H，φ，λ分別表示觀測(cè)點(diǎn)的海拔高度、緯度及經(jīng)度坐標(biāo)，R表示地球半徑.

致謝 作者對(duì)提供AURIC代碼的美國(guó)計(jì)算物理公司CPI表示衷心的感謝；對(duì)兩位審稿專(zhuān)家所提供的寶貴意見(jiàn)與建議表示衷心的感謝.

Ajello J M，Shemansky D E.1985.A reexamination of important N2cross sections by electron impact with application to the dayglow：the Lyman－Birge－Hopfield band system and NI（119.99nm）.J.Geophys.Res.，99（A10）：9845－9861，doi：10.1029／JA090iA10p09845.

Bush B C，Chakrabarti S.1995.A radiative transfer model using spherical geometry and partial frequency redistribution.J.Geophys.Res.，100（A10）：19627－19642，doi：10.1029／95JA01209.

Conway R R.1992.Self－Absorption of the N2Lyman－Birge－Hopfield bands in the far ultraviolet dayglow.J.Geophys.Res.，87（A2）：859－866，doi：10.1029／JA087iA02p00859.

Dashkevich Z V，Sergienko T I，Ivanov V E.1993.The Lyman－Birge－Hopfield bands in aurora.Planet.SpaceSci.，41（1）：81－87，doi：10.1016／0032－0633（93）90019－X.

Evans J S，Strickland D J，Paxton L J.1995.Satellite remote sensing of thermospheric O／N2and solar EUV，2.Data analysis.J.Geophys.Res.，100（A7）：12227－12233，doi：10.1029／95JA00573.

Eastes R W.2000.Emissions from the N2Lyman－Birge－Hopfield bands in the Earth′s atmosphere.Phys.Chem.Earth，25（5－6）：523－527，doi：10.1016／S1464－1917（00）00069－6.

Gladstone G R.1994.Simulations of DE 1UV airglow images.J.Geophys.Res.，99（A6）：11441－11448，doi：10.1029／93JA03525.

Huffman R E.1992.Atmospheric Ultraviolet Remote Sensing.London：Academic Press.

Link R，Strickland D J，Daniell R E.1992.AURIC airglow modules：phase 1development and application.SPIE，Ultraviolet Technology IV，1764：132－1441，doi：10.1117／12.140843.

Mcewen D J，Nicholls R W.1966.Intensity distribution of the Lyman－Birge－Hopfield band system of N2.Nature，209（5026）：902，doi：10.1038／209902a0.

Meier R R.1991.Ultraviolet spectroscopy and remote sensing of the upper atmosphere.SpaceSci.Rev.，58（1）：1－185，doi：10.1007／BF01206000.

Mende S B，Heetderks H，F(xiàn)rey H U，et al.2000.Far ultraviolet imaging from the IMAGE spacecraft.2.wideband FUV imaging.SpaceSci.Rev.，91（1）：271－285，doi：10.1023／A：1005227915363.

Ma S Y，Liu H X，Schlegel K.2002.A comparative study of magnetic storm effects on the ionosphere in the polar cap and auroral Oval－F－Region negative storm.ChineseJ.Geophys.（in Chinese），45（2）：160－169.

Oran E S，Strickland D J.1978.Photoelectron Flux in the Earth′s Ionospheric.Planet.SpaceSci.，26（1）：81－87，doi：10.1016／0032－0633（78）90056－9.

Paxton L J，Meng C I，F(xiàn)ountain G H，et al.1992.Special Sensor Ultraviolet Spectrographic Imager（SSUSI）：An instrument description.SPIE，1745：2－15，doi：10.1117／12.60595.

Paxton L J，Christensen A B，Humm D C，et al.1999.Global ultraviolet imager（GUVI）：measuring composition and energy inputs for the NASA Thermosphere Ionosphere Energetics and Dynamics（TIMED）mission.SPIE，3756：265－276，doi：10.1117／12.366380.

Snyder J P.1987.Map Projections—a Working Manual.Reston：U.S.Geological Survey.

Strickland D J，Evans J S，Bishop J，et al.1992.Atmospheric ultraviolet radiance integrated code（AURIC）：Current capabilities for rapidly modeling dayglow from the far UV to the near IR.SPIE，2831：184－199，doi：10.1117／12.257195.

Strickland D J，Cox R J，Barnes R P，et al.1994.Model for generating global images of emission from the thermosphere.Appl.Opt.，33（16）：3578－3594，doi：10.1364／AO.33.003578.

Strickland D J，Evans J S，Paxton L J.1995.Satellite remote sensing of thermospheric O／N2and solar EUV，1.Theory.J.Geophys.Res.，100（A7）：12217－12226，doi：10.1029／95JA00574.

Strickland D J，Bishop J，Evans J S，et al.1999.Atmospheric ultraviolet radiance integrated code（AURIC）：Theory，software architecture，inputs，and selected results.J.Quant.Spectrosc.Radiat.Transfer，62（6）：689－742，doi：10.1016／S0022－4073（98）00098－3.Xu W Y.2009.Variations of the auroral electrojet belt during

substorms.ChineseJ.Geophys.（in Chinese），52（3）：607－615.

附中文參考文獻(xiàn)

馬淑英，劉會(huì)欣，Schlegel K.2002.磁暴期間極光橢圓與極蓋區(qū)電離層效應(yīng)比較研究——F區(qū)負(fù)暴.地球物理學(xué)報(bào)，45（2）：160－169.

徐文耀.2009.亞暴期間極光電集流帶的變化.地球物理學(xué)報(bào)，52（3）：607－615.