大氣臨邊輻射傳輸快速計算研究

楊春平, 馬小莉， 曲兆俊， 徐振亞

(1.電子科技大學(xué) 光電信息學(xué)院 激光雷達(dá)實驗室, 成都　610054; 2. 中國空空導(dǎo)彈研究院， 河南 洛陽　471009)

?

大氣臨邊輻射傳輸快速計算研究

楊春平1, 馬小莉1， 曲兆俊2， 徐振亞2

(1.電子科技大學(xué) 光電信息學(xué)院 激光雷達(dá)實驗室, 成都610054; 2. 中國空空導(dǎo)彈研究院， 河南 洛陽471009)

摘要：對臨邊探測情形下大氣路徑透過率的計算方法進(jìn)行了研究。 利用球面大氣模型、 C-G近似方法、 光線追跡與迭代算法以及統(tǒng)計帶模式算法等建立了臨邊大氣透過率的快速計算模型。 計算分析了幾個常用遙感波段內(nèi)， 典型大氣條件、 典型切向高度下的大氣輻射傳輸特性。 最后， 將本文計算結(jié)果與CDI模型的計算結(jié)果進(jìn)行比較， 結(jié)果表明： 本文模型與CDI的計算結(jié)果相對差異在2%以內(nèi)。

關(guān)鍵詞：臨邊探測； 透過率； 光線追跡； 迭代法； 帶模式算法

0引言

臨邊遙感作為一種近年來興起的大氣探測技術(shù)， 由于具有空間覆蓋性好和垂直分辨率高的優(yōu)點， 受到大氣科學(xué)研究者的關(guān)注。 臨邊遙感是星載探測器采用以一定高度掠過地表的觀測方式， 視線路徑可長達(dá)上千公里。 目前德國不萊梅大學(xué)、 加拿大薩斯克切溫大學(xué)等對臨邊輻射傳輸建模取得了一些成果， 但模型也存在不足， 如計算光譜窄、 計算速度慢[1]。 目前國內(nèi)主要是利用SCIATRAN開展氣體含量反演算法研究[2-3]， 未見有對大氣臨邊輻射快速計算建模的相關(guān)報道。

針對上述問題， 基于球面大氣模型與帶模式算法建立了臨邊探測下大氣輻射傳輸快速計算模型。 首先， 基于HITRAN2008大氣分子光譜數(shù)據(jù)庫計算建立了0.2～20 μm波段大氣分子的吸收數(shù)據(jù)庫， 并建立了相應(yīng)的氣溶膠數(shù)據(jù)庫、 太陽光譜輻照度數(shù)據(jù)庫。 利用分層迭代算法和光線追跡算法提高了物質(zhì)吸收量的精度， 使用了帶模式算法和C-G近似算法計算臨邊路徑上的大氣透過率。 最后， 利用CDI模型的計算結(jié)果對本文模型進(jìn)行了校驗。

1大氣臨邊透過率的計算方法

在輻射計算中， 核心問題是大氣透過率的計算。 臨邊遙感下的大氣透過率計算涉及兩個關(guān)鍵問題， 一個是球面大氣中大氣分子的物質(zhì)吸收量計算， 另一個是臨邊非均勻路徑的等效均勻化計算。

1.1　臨邊球面大氣分層方法

由于臨邊探測時會受到地球曲率的影響， 從視線遠(yuǎn)端大氣頂處沿臨邊方向到達(dá)探測器， 整個視線路徑可長達(dá)數(shù)千公里， 視線路徑上的不同垂直高度處的大氣溫度、 壓力和各種大氣分子的含量隨高度分布明顯變化。 因此， 需要先對大氣進(jìn)行分層處理， 大氣按垂直高度在0～100 km范圍內(nèi)共被分為35層， 從路徑遠(yuǎn)端大氣頂處到近端大氣頂處共有70層。

臨邊切向高度決定了從切點高度一直到大氣頂處所需計算的大氣層數(shù)， 臨邊切向高度不同， 所需計算的大氣層數(shù)則不同[4]， 臨邊情形下大氣分層與切向高度之間的示意圖如圖1所示。

圖1　切向高度與大氣層數(shù)

1.2　大氣分子吸收數(shù)據(jù)的處理

基于HITRAN2008大氣分子譜線數(shù)據(jù)庫建立了0.2～20 μm波段大氣分子的吸收數(shù)據(jù)庫， 每種大氣分子的帶模式參數(shù)S/d，γc及1/d分別由下式計算得出：

(1)

(2)

(3)

式中：Δυ為光譜分辨率； Sj(T)為某種大氣分子第j條譜線在溫度T時， 光譜間隔內(nèi)的積分線強(qiáng)； 參考溫度Ts=296K； Qr， Qv分別為轉(zhuǎn)動配分函數(shù)和振動配分函數(shù)； h為普朗克常數(shù)； c為光速； k為玻爾茲曼常數(shù)； P為壓強(qiáng)； T0=273.15K； P0=101.325kPa; 計算CO2分子時x取3/4， 其他分子取1/2； γc為洛倫茲線寬。

1.3　層吸收量的計算方法

將每一層大氣的底高Hmin和頂高Hmax內(nèi)再進(jìn)一步細(xì)分為z1,z2, …，zN個高度， 新高度處的大氣廓線可利用插值算法得到。 則各層大氣內(nèi)的物質(zhì)吸收量為細(xì)分各層內(nèi)物質(zhì)吸收量的總和， 即

(4)

(5)

(6)

式中：ρ為某種分子的密度；i為大氣細(xì)分層序號； Δzi為高度zi+1和zi的高度差； ds， dz分別為兩點間的路徑微分距離和垂直高度微分距離。 由于臨邊路徑具有幾何對稱性， 左、 右半球中對稱大氣層內(nèi)的物質(zhì)吸收量是相等的。 因此， 可從切點層向左(或右)半球進(jìn)行分層迭代計算， 再結(jié)合光線追跡算法， 可同時得到細(xì)分各層大氣分子物質(zhì)吸收量與局地天頂角， 吸收量的計算流程[5]見圖2。

圖2　分層迭代算法計算流程

得到帶模式參數(shù)和各層大氣分子的物質(zhì)吸收量后， 利用C-G近似算法將臨邊非均勻路徑轉(zhuǎn)換為等效的均勻路徑， 即采用光學(xué)厚度加權(quán)平均方法來計算帶模式參數(shù)：

(7)

γc/d

(8)

γd/d

(9)

式中：γc，γd分別為洛倫茲展寬寬度和多普勒展寬寬度。

1.4　路徑透過率的計算方法

臨邊遙感下的大氣路徑透過率t為

t=tm·ta·tc

(10)

式中：tm，ta，tc分別為大氣分子、 氣溶膠和云的透過率。 將大氣分子透過率tm分為吸收透過率tabs和散射透過率tsca， 吸收透過率tabs又分為中心吸收透過率tcent、 尾翼吸收透過率ttail和典型分子的連續(xù)吸收透過率[6]， 則

(11)

式中：M為大氣分子種類， 計算中包括了H2O， CO2， O3， N2O， CO等12種氣體；tct1，tct2分別為H2O， N2的連續(xù)吸收透過率。

氣溶膠透過率ta和云透過率tc分別利用了氣溶膠和云參數(shù)數(shù)據(jù)庫， 數(shù)據(jù)庫取自本課題組前期的研究成果[7-8]。

2計算結(jié)果處理與模型校驗

計算中使用的六種大氣模式的廓線數(shù)據(jù)取自AFGL報告[7]。 計算了1976 U.S.大氣模式， 鄉(xiāng)村23 km氣溶膠， 0.2～20 μm波段不同切向高度(20 km， 40 km， 60 km)時的大氣透過率及各分子透過率如圖3所示。

由圖3(a)可以看出， 切向高度越低大氣透過率越小， 因為切向高度降低造成臨邊路徑距離增大，且密度隨高度滿足負(fù)指數(shù)關(guān)系， 切向高度越低大氣的消光能力越強(qiáng)。 由圖3(b)可以看出， 在0.2～20 μm波段內(nèi)存在眾多的氣體吸收帶。 其中， H2O的吸收帶主要分布于2.7 μm和6.3 μm波段； CO2的吸收帶主要分布于2.7 μm， 4.3 μm和15 μm波段； O3的吸收帶主要分布于9.6 μm和14 μm波段。

圖3　不同切向高度時的大氣透過率及各分子透過率

為了驗證本文所提出模型的可靠性， 利用CDI模型的計算結(jié)果進(jìn)行了對比與分析。 依據(jù)CDI模型在1976 U.S.大氣模式， 310～360 nm范圍內(nèi)5個切向高度(25 km， 30 km， 35 km， 40 km， 45 km)的大氣透過率計算結(jié)果及比較如圖4所示。

圖4　大氣透過率的計算結(jié)果與比較

由圖4可以看出，紫外波段的大氣透過率隨波長增加而逐漸增大，且切向高度越高，大氣透過率越大。在紫外光譜區(qū)內(nèi)，大氣對太

陽輻射的吸收主要是O3造成的。 在切向高度25 km以下時， 波長小于320 nm的大氣透過率幾乎都為零。 本文對大氣臨邊透過率的計算結(jié)果與CDI模型的計算結(jié)果符合很好， 兩者的相對差異在2%以內(nèi)。

3結(jié)論

根據(jù)文中臨邊遙感下大氣輻射傳輸?shù)膶捁庾V快速計算模型研究, 可以得到以下結(jié)論：

(1) 通過利用球面大氣模型與帶模式算法建立了臨邊遙感下快速大氣輻射傳輸計算模型， 解決了目前大氣臨邊輻射傳輸模型存在計算光譜窄、 計算速度慢等問題；

(2) 本文計算結(jié)果與CDI模型計算結(jié)果的相對差異在2%以內(nèi)， 驗證了本文模型的可靠性。

與逐線積分算法相比， 本文利用帶模式算法計算大氣透過率的優(yōu)點是計算速度快， 同時計算精度較高。 對大氣臨邊輻射計算建模的研究成果可用于大氣背景臨邊輻射亮度預(yù)估、 大氣環(huán)境監(jiān)測及氣體含量反演等領(lǐng)域。

參考文獻(xiàn)：

[1] Bourassa A E, Degenstein D A， Llewellyn E J. SASKTRAN: A Spherical Geometry Radiative Transfer Code for Efficient Estimation of Limb Scattered Sunlight[J]. Journal of Quantitative Spectroscopy & Radiative Transfer, 2008, 109(1)： 52-73.

[2] 齊瑾, 張鵬, 張文建, 等. 基于SCIATRAN模型的二氧化氮DOAS反演敏感性試驗[J]. 氣象學(xué)報, 2008， 66(3)： 396-404.

[3] 張恭正, 陳圣波, 王明常, 等. 基于SCIATRAN模型的Limb大氣輻射傳輸模擬[J]. 科學(xué)技術(shù)與工程, 2010， 10(6)： 1506-1509.

[4] Loughman R P, Griffioen E, Oikarinen L, et al. Comparison of Radiative Transfer Models for Limb-Viewing Scattered Sunlight Measurements[J].Journal of Geophysical Research:Atmospheres, 2004.

[5] 楊春平. 天空背景光譜特性建模及仿真[D]. 成都： 電子科技大學(xué), 2008.

[6] 孟雪琴, 吳健, 楊春平. 無云地球大氣背景輻射光譜亮度的計算模型[J]. 應(yīng)用光學(xué), 2009， 30(1)： 167-171.

[7] Anderson G P, Clough S A, Kneizys F X, et al. AFGL Atmospheric Constituent Profiles (0-120 km)[R]. Air Force Geophysics Laboratory, Hanscom, AFB, MA， 1986.

[8] 廖子君, 楊春平. 可見光波段卷云散射特性的研究[J]. 光學(xué)與光電技術(shù), 2011， 9(6)： 25-28.

Study on Rapid Calculating Method of Atmospheric

Limb Radiative Transmittance

Yang Chunping1, Ma Xiaoli1, Qu Zhaojun2, Xu Zhenya2

(1.Lidar Laboratory, School of Optoelectronic Information,University of Electronic Science and Technology of China,

Chengdu 610054, China; 2. China Airborne Missile Academy, Luoyang 471009, China)

Abstract:The calculation method of atmospheric radiative transmittance under limb detecting is studied. With the spherical atmosphere model, the C-G approximation way, the light line tracing way, the algorithms of iteration and band mode, a rapid model is established for calculating atmospheric limb radiative transmittance. The radiation properties in the conditions of several typical remote sensing bands, atmospheric modes, and tangent heights are calculated. Finally, the calculation results of this paper with CDI model are compared, and the results show that relative differences between the two models are within less than 2%.

Key words:limb detecting; transmittance; light line tracing; iteration; band mode algorithm

作者簡介：楊春平(1969-)， 男， 四川成都人， 博士， 副教授， 研究方向為目標(biāo)與環(huán)境特性研究。

基金項目：航空科學(xué)基金項目(20140180002)

收稿日期：2015-07-28

中圖分類號：P412

文獻(xiàn)標(biāo)識碼：A

文章編號：1673-5048(2015)06-0014-03