第3章
CHAPTER 3


辐射大气透过率的工程


理论计算







辐射大气透过率对光电系统(跟踪系统、搜索系统、警戒系统、热成像系统等)的设计、性能、评价等具有重要影响。大气透过率对气象卫星光学遥感仪器的通道选择和气象卫星遥感的资料反演等民用工作也是至关重要的参数。而大气窗口在气象卫星成像遥感工作中更为重要,利用它们可得到表面温度、云图、监测火情和水情等。本章以红外辐射传输为重点、(延伸至可见光和近紫外),探讨大气中影响辐射传输的因素,在提供一些基本数据资料的前提下,研究这些因素对系统性能影响的计算方法——大气透过率的工程理论计算方法。
以往的计算方法,要么有一定的局限性,没有综合考虑高度修正、倾斜路程以及与大气衰减有关的因素,给系统设计和评价带来了一定的误差; 要么需要提供太多的参数(如大气辐射传输计算软件LOWTRAN、MODTRAN等),而且这些参数可能难以获得,给实际使用带来不便,尤其是对非大气辐射传输研究的专业人员的使用带来较大困难。为此建立有效、完整且满足精度要求、便于工程应用的大气光谱透过率、平均透过率、积分透过率的模型与计算方法是必要的。
本章首先在分析地球大气的组成与结构及其辐射吸收作用的基础上,阐述大气衰减与透过率、大气中辐射衰减的物理基础,并给出大气透过率数据表。其次,研究海平面上大气气体的分子吸收、不同高度时的分子吸收修正、大气分子与微粒的散射、与气象条件有关的衰减等问题。再次,提出平均透过率与积分透过率的计算方法,给出工程理论计算示例。最后,研讨几种大气辐射传输计算软件,并进行比较分析。
3.1地球大气的组成与结构及其辐射吸收作用
通过分析地球大气的组成以及大气层各层的特点,可以明晰辐射大气传输特性和规律,了解辐射吸收效应,尤其是对不同波段的具体大气窗口的情况,这将是光电系统的分析、设计和使用等诸多方面的基础。
3.1.1地球大气的组成
大气由多种气体——氮(N2)、氧(O2)、水蒸气(H2O)、二氧化碳(CO2)、甲烷(CH4)、一氧化碳(CO)、臭氧(O3)、一氧化二氮(N2O)(不稳定,遇光、湿或热变成二氧化氮NO2及一氧化氮NO,一氧化氮又变为二氧化氮)等,以及悬浮在大气之中的一些液体和固体粒子组成(这些悬浮粒子被称为“气溶胶”粒子)。大气的主要成分是氮(体积比约为78%)和氧(体积比约为20%)。水蒸气约占总量的1%左右,二氧化碳只占0.03%~0.05%。
在低层大气中,按气体组成百分比的变化程度,可分为不变气体和可变气体两部分。所谓不变气体,是指那些组成百分比不随时间、空间而变的气体。其中主要的有氮(N2)、氧(O2)、氩(Ar)等。所谓可变气体,是指组成百分比随时间、空间而变的气体。其中主要有水汽(H2O)、二氧化碳(CO2)和臭氧(O3)等。如果把大气中的水汽和气溶胶粒子除去,那么这样的大气常被称为干洁大气。在80km以下的干洁大气,其组成百分比基本保持不变。
表3.1列出了低层干洁大气的组成百分比的数字。可以看到,N2、O2及Ar的总含量约占整个大气的99.97%,其余的稀有气体如Ne、He、Kr、Xe、H2等,总和不超过0.03%。二氧化碳在表3.1中列为不变气体,但实际上自工业革命以来,因化石燃料的燃烧而一直在慢慢地增加。


表3.1海平面大气成分表



气体分子量容积百分比/%

氮N228.013478.084
氧O231.998820.9476
氩Ar39.9480.934
二氧化碳CO244.009950.0322①
氖Ne20.1830.001818
氦He4.00260.000524
氪Kr83.800.000114
氢H22.015940.00005
氙Xe131.300.0000087
甲烷CH416.0430.00016
一氧化二氮N2O440.000028
一氧化碳CO280.0000075

注: ①据报道,此值目前已达到0.034,甚至以上。

实测表明,在高度小于80km的大气层中,干洁大气组成百分比保持不变。若要计算80km以下某种成分的数密度n0时,可按下式计算
ni=FiN(z<80km)(3.1)
式中,Fi为第i种气体成分在海平面大气中的容积百分比; N为所研究高度上的大气总数密度。
3.1.2大气层结构
按照大气层内垂直方向的温度分布及运动特点,按国际通用的术语,可将整个大气层在铅垂方向分成五层。它们是对流层、平流层、中间层、暖层及散逸层(见图3.1)。现将各层的特点简述如下。


图3.1大气分层图



1. 对流层
对流层是大气层的最底层。它也是各层中最薄的一层。然而整个大气层质量的3/4和90%以上的水汽及主要的天气现象都集中在这一层内,因而此层也是和人类活动关系最密切的一层。由于不同纬度球表面受太阳辐射加热的程度不同,不同季节时地球表面的受辐射加热的情况也不同,与之相应,大气的垂直运动也会不同,因而对流层的厚度会随纬度和季节而变化。对流层的厚度在热带平均为17~18km,温带平均为10~12km; 在寒带则只有8~9km。夏季的厚度通常要大于冬季的厚度。
对流层的特征是: 它有强烈的对流运动,在高低层之间有质量和热交换,这就使近地面的水汽向上输送,形成云雨; 在层内的大气温度随高度升高而线性地降低,其温度递减率为0.65℃/100m。显然,对流层厚度大的低纬度区域的对流层顶部的温度要比高纬度区域的顶部温度低。
2. 平流层
平流层是从对流层顶至55km高度之间的大气层。平流层的特征是,它没有强烈的上下对流运动,因而气流平稳,远程的喷气式客机通常在此层内飞行; 在平流层的下部,温度随高度的变化很少,从30km左右的高度开始,温度随高度的增加而增加,到此层的顶层温度为-3℃~-17℃。这主要是臭氧对太阳辐射的强烈吸收而造成的。此层内水汽和尘埃的含量均很少,空气透明度很好。
3. 中间层
中间层是高度为55~85km的一层大气。此层和对流层有相似的特点,温度随高度而迅速下降,至此层顶时温度已降至-83℃~-113℃。此层内也有相当强烈的上下对流运动,故称为高空对流层。
在0~85km高度范围内,温度可以用7个连续的线性方程来描述,其形式为
T=Tb+Lb(H-Hb)(3.2)
式中,下脚标b的值为0~6,第一层为0,第二至第七层的b分别为1~6,Lb的值列在表3.2中。


表3.20~85km内温度廓线图各段的参考高度和梯度



b高度z/km温度梯度Lb/(K·km-1)

00-6.5
1110.0
2201.0
3322.8
4470.0
551-2.8
671-2.0
784.85—


在上述高度范围内的标准大气的温度高度廓线图如图3.2所示。


图3.2标准大气温度高度廓线图


4. 暖层
暖层是从中间层顶部向上直至800km高度之间的一层大气。此层的大气温度随高度的增加而迅速上升。到300km高度时,温度可接近1000℃。此层内的大气由于强烈的太阳紫外线和宇宙线的照射,存在着几个电离层,即D、E、F(包括F1、F2)层。这些电离层的高度及电子浓度如表3.3所示。


表3.3各电离层的特征



层次名称组成

高度/km最大电子浓度/(105个·cm-3)

白天夜间白天夜间


DO+2+e-60~90—0.15~2.00
EO+2+e-90~13090~1001.5~3.00.2
F1N+2+e-160~280—2.5~4.00
F2O++e-300~350280~30010~300.2


由于电离层的存在,使短波无线电的远距离传输成为可能。在1000km以下的整个大气层的温度廓线图可参见图3.3。


图3.31000km以下大气层的温度廓线图



5. 散逸层
散逸层自暖层向上,统称为散逸层。它是大气层和星际空间的过渡带。该层离地球较远,引力小,空气稀薄。某些高速运动的空气分子一旦被撞击出去后,就很难再回向地球,因而不断散逸到宇宙空间。近代卫星探测资料表明,大气层上界为2000~3000km。
3.1.3大气的辐射吸收作用
大气的吸收和散射与波长有关,即有明显的选择性。比如,对于波长范围为1~11μm、1.2~1.3μm、1.6~1.75μm、2.1~2.4μm、3.4~4.2μm、8~13μm等波段内的辐射,大气的吸收作用较小,有所谓“大气窗口”的称谓。表3.4列出了8个大气窗口的波段范围、透过率值及其物理机制,其中,透过率值为特定条件下的典型平均值。有效利用大气窗口,可以使系统的作用距离达到相对最大化。值得一提的是,在高空使用,对红外系统有着更好的条件,这是因为高空水蒸气及CO2很少的缘故。


表3.48个大气窗口的波段范围、透过率值及其物理机制



序号窗 口 名 称波 段 范 围透过率/%散射吸收机制γ为散射系数

1紫外、可见、近红外0.3~1.115μm70、95、80瑞利散射,γ∝λ-4,臭氧吸收
2近红外窗口1.4~1.9μm60~95米氏散射,γ∝λ-2,CO2吸收
3近红外窗口2.0~2.5μm80米氏散射,γ∝λ-2,H2O吸收
4中红外窗口3.5~5.0μm60~70米氏散射,γ∝λ-2,H2O吸收
5远(热)红外窗口8~14.0μm80米氏散射,γ∝λ-2,H2O吸收
6微波窗口1.0~1.8mm35~40瑞利散射,γ∝λ-4,可忽略
7微波窗口2.0~5.0mm50~70只在1.63mm和3.48mm处有水汽吸收带,其他波长无吸收损失
8微波窗口8.0~1000mm100无任何散射和吸收损失

CO2在2.7μm、4.3μm为中心的附近及11.4~20μm的区域有吸收带,在1.4μm、16μm、2.0μm、4.8μm、5.2μm、9.4μm、10.4μm处有弱的吸收带。水汽在1.87μm、27μm、6.27μm处有强吸收带,在0.94μm、1.1μm、1.38μm、5.2μm处有弱吸收带。N2O在45μm处有一个较强的吸收带,在2.9μm、4.05μm、7.7μm、8.6μm、17.1μm处还有弱吸收带。CO在4.6μm处有一个较强的吸收带,在2.3μm处有一个弱的吸收带。CH4在331μm、6.5μm、7.65μm处有吸收带。O3在9.6μm处有一个较强的吸收带,在4.7μm、89μm、14μm处有弱吸收带。
有资料把红外辐射划分为近红外、短波红外、中波红外、长波红外、甚长波红外等多种区间。由于大气会对红外辐射进行吸收,留下3个“窗口”波段,即1~3μm、3~5μm、8~14μm。在军事应用上将这3个波段称为近红外、中红外以及远红外波段,具体更小波段划分及其大气透过率、大气背景辐射强度的比较,如表3.5所示。


表3.5红外波段大气透过率与大气背景辐射强度定性分析



波长范围/μm波段大气透过率大气背景亮度

1.1~1.4J高夜晚很低
1.5~1.8H高非常低
2.0~2.4K高非常低
3.0~4.0L3.0~3.5μm: 中等

3.5~4.0μm: 高低
4.6~5.0M低高
7.5~14.5N8~9μm&10~12μm: 中等

其他波段: 低非常高
17~4017~25μm: Q
28~40μm: Z非常低非常高
330~370非常低低

大气透过率与观测地点、大气状况等有关。海拔较低的观测地点大气透过要差一些,在离地15km的高度范围内,只有部分窗口能透过天体或者空间目标的红外辐射。
有综合统计分析,吸收红外辐射的主要因素是水蒸气,它主要集中在2~3km大气层以下; 虽然CO2只占大气体积的0.03%~0.05%,但它是红外辐射衰减的另一重要原因,它在空气中的分布比水蒸气均匀。大气中O3含量很少,它主要位于10~40km空间范围中,特别是集中在20~30km高度上。而在20km以下的大气层中,O3对辐射的吸收影响是非常小的,只有在雷雨之后,大气中O3的含量才突然增大。因此对用在20km以下的红外系统,O3的衰减作用可以忽略不计。另外,在下层大气中,NO、NO2、CO等气体的吸收,通常可以忽略。
3.2大气衰减与透过率
大多数光电系统必须通过地球大气才能观察到目标,从设计者的角度来看,这是不利的。因为从目标来的辐射量在到达光电传感器之前,受到大气中某些气体的选择性吸收,大气中悬浮微粒的散射,同时还要经受大气某些特性剧烈变化的调制。对于干洁大气,服从瑞利的分子散射定律。但在通常情况下,实际大气中总含有较大的悬浮质点,特别是在近地层大气中,常常存在有雾、雨和尘粒等,其大小不仅与入射辐射的波长相仿,而且还可能远远超过波长。因此当辐射落在这些质点上时,将会发生散射、衍射及折射现象。
当大气分子或微粒遇到入射的光时,会受激而向四面八方发射出频率与入射光相同的光,这种现象就是光的散射。当有散射存在时,沿某一方向传输的光会因散射而使传输受到衰减。此种因散射产生的衰减将和吸收造成的衰减一样,服从朗伯比尔定律。通过大气而减弱的整个过程称为衰减。通过大气的透过率可以表达如下: 
τ=e-σx(3.3)
式中,σ称为衰减系数x是路程长度。大气对辐射能传播的衰减作用表现为散射损失、吸收损失和反射损失3种物理机制上,因此,大气光谱衰减系数σ可表示为
σ=α+γ+σ反射(3.4)
式中,α、γ、σ反射分别表示由大气分子吸收、散射、反射机制引起的大气衰减系数(km-1)。这里α是吸收系数,起因于大气分子的吸收; 而γ是散射系数,起因于气体分子、烟、雾等的散射。可想而知,α和γ二者均随波长而变化。其中,γ与辐射波长的λ4或λ2成正比(瑞利散射定律或米氏散射定律),即波长越小,散射和吸收损失越小。
实验测定,大气微粒包括水滴(云雾、降水)、冰粒和尘埃,它们的直径一般不超过100μm; 大气降水云层的粒子中包括有雨滴、冰粒、雪花和干湿冰雹,其直径均大于100μm,有的可达几毫米(如雨滴)、几厘米(如冰雹)。可见在绝大部分(除了冰雹)天气情况下,它们都远小于微波波长,尤其是对于波长8.0~1000mm的微波,不存在散射和吸收损失,大气透过率为100%。一般而言,对于红外、可见、紫外3个波段,它们的大气透过率大小的排序依次是τ红外>τ可见>τ紫外。
在大多数情况下,如果忽略σ反射,则衰减由吸收、散射因素造成。可得: σ=α+γ。
在红外波段,吸收比散射严重得多。雾和云是很强的散射物质,实际上对红外辐射也是不透过的。正因如此,红外系统受天气影响较严重,难以具有全天候性能,尤其是被动红外系统。
辐射的大气透过率取决于大气气体的性质和它的质点浓度、分布。因而也就取决于气象条件,特别是在大气中的水蒸气、其他气体和尘粒的含量不断变化,大气透过率也随之变化,同时还与海拔高度有关。理论研究只能对上述关系得出近似的结果,因此需要在各种气象条件下,对不同高度和不同波长进行大量的测量来弥补理论研究之不足,从实验数据中构成经验公式和半经验公式,从而可以近似求出各种大气条件下的大气透过率。因此,辐射的大气透过率取决于气象条件,而且随天气条件和高度而变。用理论方法只可求出近似的大气衰减和透过率,而在一定的天气条件和高度下,可用已推导出的经验公式求出精确的近似值。
3.3大气中辐射衰减的物理基础
光电(红外)系统用于观测、搜索、跟踪远距离目标时,目标辐射在到达系统的光学系统之前必须通过大气,并被衰减。在衰减的同时,因大气梯度和湍流引起空气折射率不均,使辐射发生畸变。此外,大气是固有辐射源。所有这些现象使远距离目标的图像质量变坏。大气对光电(红外)系统的影响基本表现在辐射衰减,红外辐射衰减与以下3种现象有关: 
(1) 大气气体分子的吸收; 
(2) 大气中分子、气溶胶、微粒的散射; 
(3) 因气象条件(云、雾、雨、雪)的衰减。
这样,在分析红外系统工作时必须全面考虑上述3种现象。但同时值得一提的是,当因气象衰减与天气条件有关时,气体分子和微粒的红外辐射吸收与散射经常发生。
大气传输的特征是,光谱透过率τa(λ)和衰减系数(消光系数)μ(λ)之间的关系可用布盖乐朗伯定律表示,即
τa(λ)=φe(λ,R)/φe(λ,0)=exp(-μ(λ)R)(3.5)
式中,R为目标与红外系统之间的距离; φe(λ,R)为距离R处(或红外系统处)目标或背景辐射通量光谱密度; φe(λ,0)为R=0时φe(λ,R)的特例; λ为波长。
平均透过率和平均衰减系数为
τ-a=1λ2-λ1∫λ2λ1τa(λ)dλ(3.6)

μ-=1λ2-λ1∫λ2λ1μ(λ)dλ(3.7)

式中,λ1~λ2为光谱范围。积分透过率为
τa=∫λ2λ1φe(λ,0)τa(λ)dλ∫λ2λ1φe(λ,0)dλ(3.8)
大气光谱透过率τa(λ)用下式确定: 
τa(λ)=τ1(λ)·τ2(λ)·τ3(λ)(3.9)
式中,τ1(λ)、τ2(λ)、τ3(λ)分别为被吸收、散射和因气象衰减制约的大气光谱透过率。
3.4大气透过率数据表
早期关于大气红外衰减方面的研究,主要集中在大气气体分子的带吸收研究,这一时期在实验室和野外都进行了大量的测量工作,限于当时的理论和实验手段,所得的大气红外辐射和传输特性数据没有足够高的光谱分辨率,限制了这些数据在实际工程中的应用。由波带模型法描述的大气吸收模型已被广泛地应用到工程计算中。但是,在用波带模型法计算中,需要知道谱线的参数,如线强、线宽、线距等,然而这些参数的获得并不容易。
很早以前,有些学者将实验确定的透过率数据,并考查了一些有用的测量,制成了以01μm为间隔、从0.3~7μm的一个很宽的吸收物质浓度范围内的光谱透过率表,表3.6是可降水分含量从0.1~1000mm的水蒸气光谱透过率表,表3.7是路程长度从0.1~1000km的二氧化碳光谱透过率表。后来,小哈得逊将这两个表格从7μm扩展到14μm,水蒸气的计算结果列在表3.8中,二氧化碳的计算结果列在表3.9中。这些表格是基于实验数据并用Elsasser波带模型插值而编制的,最后形成波长间隔为0.1μm的光谱透过率表。


表3.6海平面水平路程上的水蒸气的光谱透过率(0.3~6.9μm)



波长

/μm
水蒸气含量(可降水分毫米数)

0.10.20.51251020501002005001000

0.30.9800.9720.9550.9370.9110.8600.8020.7230.5740.4280.2630.0760.012
0.40.9800.9720.9550.9370.9110.8600.8020.7230.5740.4280.2630.0760.012
0.50.9860.9800.9680.9560.9370.9010.8610.8040.6950.5790.4330.2150.079
0.60.9900.9860.9770.9680.9550.9290.9000.8600.7790.6920.5750.3750.210
0.70.9910.9870.9800.9720.9600.9370.9100.8730.8000.7220.6150.4250.260
0.80.9890.9840.9750.9650.9500.9220.8910.8450.7580.6630.5390.3300.168
0.90.9650.9510.9220.8900.8440.7570.6610.5350.3260.1650.0500.0020
1.00.9900.9860.9770.9680.9550.9290.9000.8600.7790.6920.5750.3750.210
1.10.9700.9580.9320.9050.8660.7900.7070.5950.4060.2350.0930.0080
1.20.9800.9720.9550.9370.9110.8600.8020.7230.5740.4280.2630.0760.012
1.30.7260.6110.4320.2680.1160.0130000000
1.40.9300.9020.8440.7820.6950.5360.3810.2160.0640.005000
1.50.9970.9940.9910.9880.9820.9720.9600.9440.9110.8740.8230.7240.616
1.60.9980.9970.9960.9940.9910.9860.9800.9720.9560.9370.9110.8600.802
1.70.9980.9970.9960.9940.9910.9860.9800.9720.9560.9370.9110.8600.802
1.80.7920.7070.5550.4060.2390.0620.008000000
1.90.9600.9430.9110.8740.8220.7230.6170.4780.2620.1130.02400
2.00.9850.9790.9660.9530.9330.8940.8510.7900.6740.5520.4010.1840.006
2.10.9970.9940.9910.9880.9820.9720.9600.9440.9110.8740.8230.7240.616
2.20.9980.9970.9960.9940.9910.9860.9800.9720.9560.9370.9110.8600.802
2.30.9970.9940.9910.9880.9820.9720.9600.9440.9110.8740.8230.7240.616
2.40.9800.9720.9550.9370.9110.8600.8020.7230.5740.4280.2630.0760.012
2.50.9300.9020.8440.7820.6950.5360.3810.2160.0640.005000
2.60.6170.4790.2610.1100.00200000000
2.70.3610.1960.0400.004000000000
2.80.4530.2890.0920.0170.00100000000
2.90.6890.5710.3690.2050.0730.0050000000
3.00.8510.7900.6730.5520.4010.1840.0600.00800000
3.10.9000.8600.7790.6920.5740.3750.2100.0760.0050000
续表



波长

/μm
水蒸气含量(可降水分毫米数)

0.10.20.51251020501002005001000

3.20.9250.8940.8330.7660.6740.5060.3470.1840.0350.003000
3.30.9500.9300.8880.8430.7790.6580.5310.3770.1610.0480.00500
3.40.9730.9620.9390.9140.8800.8110.7350.6330.4480.2850.1300.0170.001
3.50.9880.9830.9730.9620.9460.9150.8810.8320.7360.6350.5020.2870.133
3.60.9940.9920.9870.9820.9730.9580.9470.9160.8660.8120.7380.5960.452
3.70.9970.9940.9910.9880.9820.9720.9600.9440.9110.8740.8230.7240.616
3.80.9980.9970.9950.9940.9910.9860.9800.9720.9560.9370.9110.8600.802
3.90.9980.9970.9950.9940.9910.9860.9800.9720.9560.9370.9110.8600.802
4.00.9970.9950.9930.9900.9870.9770.9700.9600.9300.9000.8700.7900.700
4.10.9770.9940.9910.9880.9820.9720.9600.9440.9110.8740.8230.7240.616
4.20.9940.9920.9870.9820.9730.9580.9470.9160.8660.8120.7380.5960.452
4.30.9910.9840.9750.9720.9500.9370.9100.8730.8000.7220.6150.4250.260
4.40.9800.9720.9550.9370.9110.8600.8020.7230.5740.4280.2630.0760.012
4.50.9700.9580.9320.9050.8660.7900.7070.5950.4000.2350.0930.0080
4.60.9800.9430.9110.8740.8220.7230.6170.4780.2620.1130.02400
4.70.9500.9300.8880.8430.7790.6580.5310.3770.1610.0480.00500
4.80.9400.9150.8660.8120.7360.5950.4520.2890.1170.0180.00100
4.90.9300.9020.8440.7820.6950.5360.3810.2160.0640.005000
5.00.9150.8800.8110.7360.6340.4510.2860.1320.0170000
5.10.8850.8390.7470.6490.5190.3080.1490.0410.0010000
5.20.8460.7840.6640.5390.3850.1690.0520.00600000
5.30.7920.7070.5550.4060.2390.0620.008000000
5.40.7260.6110.4320.2680.1160.0130000000
5.50.6170.4790.2610.1100.03500000000
5.60.4910.3310.1210.0290.00200000000
5.70.3610.1960.0400.004000000000
5.80.1410.0440.0010000000000
5.90.1410.0440.0010000000000
6.00.1800.0580.0030000000000
6.10.2600.1120.0120000000000
6.20.6520.5240.3130.1530.0430.0010000000
6.30.5520.4010.1820.0600.00800000000
6.40.3170.1570.0250.002000000000
6.50.1640.0490.0020000000000
6.60.1380.0420.0010000000000
6.70.3220.1620.0370.002000000000
6.80.3610.1960.0400.004000000000
6.90.4160.2500.0680.010000000000


表3.7海平面水平路程上的二氧化碳的光谱透过率(0.3~6.9μm)



波长

/μm

路程长度/km

0.10.20.51251020501002005001000

0.31111111111111
0.41111111111111
0.51111111111111
0.61111111111111
0.71111111111111
0.81111111111111
0.91111111111111
1.01111111111111
1.11111111111111
1.21111111111111
1.31110.9990.9990.9990.9980.9970.9960.9940.9920.9870.982
1.40.9960.9950.9920.9880.9840.9750.9640.9490.9190.8850.8380.7470.649
1.50.9990.9990.9980.9980.9970.9950.9930.9900.9840.9760.9670.9490.927
1.60.9960.9950.9920.9880.9840.9750.9640.9490.9190.8850.8380.7470.649
1.71110.9990.9990.9990.9980.9970.9960.9940.9920.9870.982
1.81111111111111
1.91110.9990.9990.9990.9980.9970.9960.9940.9920.9870.982
2.00.9780.9690.9510.9310.9030.8470.7850.6990.5410.3870.2210.0530.006
2.10.9980.9970.9960.9940.9920.9870.9820.9740.9590.9420.9190.8720.820
2.21111111111111
2.31111111111111
2.41111111111111
2.51111111111111
2.61111111111111
2.70.7990.7180.5690.4190.2530.0710.011000000
2.80.8710.8040.6950.5780.4320.2150.0790.01300000
2.90.9970.9950.9930.9900.9850.9770.9680.9540.9270.8980.8550.7720.683
3.01111111111111
3.11111111111111
3.21111111111111
3.31111111111111
3.41111111111111
3.51111111111111
3.61111111111111
3.71111111111111
3.81111111111111
3.91111111111111
4.00.9980.9970.9960.9940.9910.9860.9800.9710.9550.9370.9110.8590.802
4.10.9830.9750.9610.9440.9210.8760.8250.7550.6220.4850.3220.1180.027
4.20.6730.5510.4450.1820.0590.0030000000
4.30.0980.01600000000000
续表



波长

/μm

路程长度/km

0.10.20.51251020501002005001000

4.40.4810.3190.1150.0260.00200000000
4.50.9570.9490.9030.8630.8070.6990.5850.4390.2220.0840.01400
4.60.9950.9930.9890.9850.9780.9660.9510.9310.8910.8450.7830.6630.539
4.70.9950.9930.9890.9850.9780.9660.9510.9310.8910.8450.7830.6630.539
4.80.9760.9660.9450.9220.8910.8280.7590.6640.4920.3310.1690.0300.002
4.90.9750.9640.9430.9200.8860.8220.7500.6520.4680.3130.1530.0240.001
5.00.9990.9980.9970.9950.9940.9900.9860.9790.9680.9540.9350.8970.855
5.110.9990.9990.9980.9980.9960.9940.9920.9880.9840.9760.9610.946
5.20.9860.9800.9680.9550.9360.8990.8570.7990.6870.5690.4200.2030.072
5.30.9970.9950.9930.9890.9840.9760.9660.9510.9230.8910.8460.7600.666
5.41111111111111
5.51111111111111
5.61111111111111
5.71111111111111
5.81111111111111
5.91111111111111
6.01111111111111
6.11111111111111
6.21111111111111
6.31111111111111
6.41111111111111
6.51111111111111
6.61111111111111
6.71111111111111
6.81111111111111
6.91111111111111


表3.8海平面水平路程上的水蒸气的光谱透过率(7.0~13.9μm)



波长

/μm

水蒸气含量(可降水分毫米数)

0.20.5125102050100200

7.00.5690.2450.0600.004000000
7.10.7160.4330.1880.035000000
7.20.7820.5400.2920.0850.00200000
7.30.8490.6640.4410.1940.01700000
7.40.9220.8170.6660.4440.1320.0180000
7.50.9470.8740.7620.5820.2580.0660000
7.60.9220.8170.6660.4440.1320.0180000
7.70.9780.9440.8840.7960.5640.3280.1020.00300
7.80.9740.9370.8780.7710.5230.2730.0740.00200
7.90.9820.9590.9200.8420.6580.4330.1870.01500
续表



波长

/μm

水蒸气含量(可降水分毫米数)

0.20.5125102050100200

8.00.9900.9750.9510.9040.7770.6030.3650.0800.0060
8.10.9940.9860.9720.9450.8690.7540.5680.2440.0590.003
8.20.9930.9820.9640.9300.8340.6960.4840.1630.0270
8.30.9950.9880.9760.9530.8870.7860.6180.3000.0900.008
8.40.9950.9870.9750.9500.8800.7740.5990.2780.0770.006
8.50.9940.9860.9720.9440.8660.7500.5620.2370.0560.003
8.60.9960.9920.9820.9650.9150.8370.7020.4110.1690.029
8.70.9960.9920.9830.9660.9160.8390.7040.4160.1730.030
8.80.9970.9930.9830.9660.9170.8410.7070.4210.1770.031
8.90.9970.9920.9830.9660.9180.8430.7090.4250.1800.032
9.00.9970.9920.9840.9680.9210.8480.7190.4400.1930.037
9.10.9970.9920.9850.9700.9260.8580.7350.4640.2150.046
9.20.9970.9930.9850.9710.9290.8630.7440.4780.2280.052
9.30.9970.9930.9860.9720.9300.8670.7500.4890.2390.057
9.40.9970.9930.9860.9730.9330.8700.7560.4980.2480.061
9.50.9970.9930.9870.9730.9340.8730.7620.5070.2570.066
9.60.9970.9930.9870.9740.9360.8760.7660.5160.2650.070
9.70.9970.9930.9870.9740.9370.8780.7700.5210.2700.073
9.80.9970.9940.9870.9750.9380.8800.7730.5260.2770.077
9.90.9970.9940.9870.9750.9390.8820.7770.5320.2830.080
10.00.9980.9940.9880.9750.9400.8830.7800.5380.2890.083
10.10.9980.9940.9880.9750.9400.8830.7800.5380.2890.083
10.20.9980.9940.9880.9750.9400.8830.7800.5380.2890.083
10.30.9980.9940.9880.9760.9400.8840.7810.5400.2920.085
10.40.9980.9940.9880.9760.9410.8850.7820.5420.2940.086
10.50.9980.9940.9880.9760.9410.8860.7840.5440.2950.087
10.60.9980.9940.9880.9760.9420.8870.7860.5480.3000.089
10.70.9980.9940.9880.9760.9420.8870.7870.5500.3020.091
10.80.9980.9940.9880.9760.9410.8860.7840.5440.2950.087
10.90.9980.9940.9880.9760.9400.8840.7810.5400.2920.085
11.00.9980.9940.9880.9750.9400.8830.7790.5360.2870.082
11.10.9980.9940.9870.9750.9390.8820.7770.5320.2830.080
11.20.9970.9930.9860.9720.9310.8670.7500.4870.2370.056
11.30.9970.9920.9850.9700.9270.8590.7380.4670.2180.048
11.40.9970.9930.9860.9710.9300.8650.7480.4850.2350.055
11.50.9970.9930.9860.9720.9320.8680.7530.4930.2430.059
11.60.9970.9930.9870.9740.9350.8750.7650.5130.2620.069
11.70.9960.9900.9800.9610.9060.8200.6730.3720.1380.019
11.80.9970.9920.9820.9690.9250.8630.7330.4600.2120.045
11.90.9970.9930.9860.9720.9320.8690.7550.4950.2450.060
12.00.9970.9930.9870.9740.9370.8780.7700.5210.2700.073
续表



波长

/μm

水蒸气含量(可降水分毫米数)

0.20.5125102050100200

12.10.9970.9940.9870.9750.9380.8800.7730.5260.2770.077
12.20.9970.9940.9870.9750.9380.8800.7750.5280.2790.078
12.30.9970.9930.9870.9740.9370.8780.7700.5210.2700.073
12.40.9970.9930.9870.9740.9350.8740.7640.5110.2610.068
12.50.9970.9930.9860.9730.9330.8710.7590.5020.2520.063
12.60.9970.9930.9860.9720.9310.8680.7520.4910.2410.058
12.70.9970.9930.9850.9710.9290.8630.7440.4780.2280.052
12.80.9970.9920.9850.9700.9260.8580.7360.4660.2170.047
12.90.9970.9920.9840.9690.9240.8530.7280.4520.2040.041
13.00.9970.9920.9840.9670.9210.8460.7180.4370.1910.036
13.10.9960.9910.9830.9660.9180.8430.7090.4240.1800.032
13.20.9960.9910.9820.9650.9150.8370.7010.4110.1690.028
13.30.9960.9910.9820.9640.9120.8310.6900.3970.1530.025
13.40.9960.9900.9810.9620.9080.8250.6810.3820.1460.021
13.50.9960.9900.9800.9610.9050.8190.6700.3680.1360.019
13.60.9960.9900.9790.9590.9020.8130.6610.3550.1260.016
13.70.9960.9890.9790.9580.8980.8070.6510.3420.1170.014
13.80.9960.9890.9780.9560.8940.8000.6400.3280.1070.011
13.90.9950.9880.9770.9550.8910.7930.6290.3130.0980.010




表3.9海平面水平路程上的二氧化碳的光谱透过率(7.0~13.9μm)



波长

/μm

路程长度/km

0.20.5125102050100200

7.01111111111
7.11111111111
7.21111111111
7.31111111111
7.41111111111
7.51111111111
7.61111111111
7.71111111111
7.81111111111
7.91111111111
8.01111111111
8.11111111111
8.21111111111
8.31111111111
8.41111111111
8.51111111111
8.61111111111
续表



波长

/μm

路程长度/km

0.20.5125102050100200

8.71111111111
8.81111111111
8.91111111111
9.01111111111
9.1110.9990.9990.9980.9950.9910.9780.9550.914
9.2110.9990.9980.9950.9910.9820.9550.9130.834
9.30.9990.9970.9950.9900.9750.9510.9040.7760.6050.363
9.40.9930.9820.9650.9310.8370.7000.4910.1680.0280.001
9.50.9930.9830.9670.9350.8420.7150.5120.1870.0350.001
9.60.9960.9900.9800.9610.9060.8210.6750.3630.1400.029
9.70.9950.9860.9730.9470.8730.7610.5800.2560.0650.004
9.80.9970.9920.9840.9690.9240.8580.7300.4550.2060.043
9.90.9980.9950.9890.9790.9480.8970.8110.5850.3420.123
10.0110.9990.9970.9940.9890.9780.9450.8920.797
10.110.9990.9980.9960.9900.9800.9600.9020.8140.663
10.20.9970.9940.9880.9770.9430.8900.7920.5580.3120.097
10.30.9970.9940.9870.9750.9390.8810.7770.5320.2830.080
10.4110.9990.9980.9950.9910.9820.9550.9130.834
10.5110.9990.9980.9980.9950.9910.9780.9550.914
10.6110.9990.9990.9980.9950.9910.9780.9550.914
10.71110.9990.9990.9970.9950.9860.9730.947
10.8110.9990.9980.9980.9950.9910.9780.9550.914
10.910.9990.9990.9970.9930.9860.9730.9340.8720.761
11.010.9990.9990.9970.9930.9860.9730.9340.8720.761
11.110.9990.9980.9970.9920.9840.9690.9230.8550.726
11.210.9990.9980.9950.9890.9780.9550.8920.7960.633
11.30.9990.9990.9970.9940.9850.9710.9420.8620.7420.552
11.40.9990.9980.9970.9930.9830.9660.9340.8420.7090.503
11.50.9990.9980.9960.9920.9800.9600.9210.8140.6610.438
11.60.9990.9980.9950.9910.9770.9550.9120.7940.6320.399
11.70.9990.9980.9950.9910.9770.9550.9120.7940.6320.399
11.80.9990.9980.9970.9930.9830.9660.9340.8420.7090.503
11.910.9990.9980.9950.9890.9780.9550.8920.7960.633
12.0110.9990.9990.9970.9930.9860.9660.9340.872
12.1110.9990.9980.9980.9950.9910.9780.9550.914
12.2110.9990.9980.9980.9950.9910.9780.9550.914
12.30.9980.9950.9900.9810.9520.9070.8230.6140.3760.142
12.40.9940.9850.9700.9410.8590.7380.5450.2180.0480.002
12.50.9870.9680.9360.8770.7190.5170.2680.0370.0010
12.60.9800.9500.9030.8150.5990.3580.1290.00600
12.70.9960.9890.9790.9590.8990.8090.6540.3460.1200.015
续表



波长

/μm

路程长度/km

0.20.5125102050100200

12.80.9900.9740.9490.9010.7700.5920.3510.0720.0050
12.90.9850.9620.9250.8560.6770.4580.2100.02000
13.00.9910.9770.9550.9120.7940.6300.3970.0990.0100
13.10.9900.9740.9490.9000.7680.5920.3480.0710.0050
13.20.9780.9460.8950.8010.5750.3300.1090.00400
13.30.9520.8840.7820.6110.2920.0850.007000
13.40.9350.8460.7150.5120.1870.0350.001000
13.50.9010.7670.5930.3520.0700.0050000
13.60.9010.7920.6270.3510.0970.0090000
13.70.9160.8030.6440.4150.1100.0120000
13.80.8580.6810.4640.2150.02100000
13.90.7780.5340.2860.0820.00200000

上述表格的使用并不困难。例如,要想求得某一段水平路程上与水蒸气有关的透过率,那么可以根据已知的气象条件(温度、湿度),以及水平路程的长度来计算可凝结水量,再通过查表得到各波长上或0.1μm波长间隔的与水蒸气有关的透过率。同样,根据已知的水平路程,可以由表查得各个波长上或0.1μm波长间隔的与二氧化碳有关的透过率。
这些表只适用于海平面的水平路程上。在高空,由于大气压力和温度降低,光谱吸收线和光谱带的宽度都变窄了。可以预料,通过同样的路程时,吸收变小,所以大气的透过率就会增加。伴随温度的降低,也使透过率稍有增加,不过其影响很小,一般都忽略不计。如果作稍微简单的修正,这些表格就可用于高空。在高度为h的x水平路程长度上,其光谱透过率等于水平长度为x0的等效海平面路程上的透过率,x和x0有如下关系: 
x0=x(P/P0)k(3.10)
式中,P/P0是h高度的大气压力与海平面大气压力之比; k为常数,对水蒸气是0.5,对二氧化碳是1.5。值得一提的是,等效海平面路程是大气透过率计算时的一个重要概念,很显然,在具有相同透过率的情况下,高空的路程要比海平面的路程更长一些。如果要计算某一高度上的一段路程上的透过率,那么可以根据表3.10给出的高度修正因子,再由式(3.10)求出等效海平面路程,这样就可以计算不同高度的水平路程的透过率了。


表3.10用以折合成等效海平面路程的高度修正因子(P/P0)k的值



高度/km

高度修正因子
水蒸气二氧化碳高度/km

高度修正因子
水蒸气二氧化碳


0.3050.9810.9402.140.8690.660
0.6100.9610.8332.440.8520.620
0.9150.9420.8402.740.8350.580
1.220.9230.7743.050.8190.548
1.520.9040.7433.810.7900.494
1.830.8860.6994.570.7390.404
续表



高度/km

高度修正因子
水蒸气二氧化碳高度/km

高度修正因子
水蒸气二氧化碳

5.340.7140.36415.20.3480.042
6.100.6700.29918.30.2720.020
6.860.6430.26621.40.2140.010
7.620.6090.22624.40.1670.005
9.150.5520.16827.40.1340.002
10.70.4860.11530.50.1050.001
12.20.4410.085


3.5海平面上大气气体的分子吸收
水汽(H2O)、CO2分子产生最强的选择性红外辐射吸收,因此,综合透过率结果为水汽透过率τH2O(λ)和CO2分子透过率τCO2(λ)的乘积,即
τ1(λ)=τH2O(λ)·τCO2(λ)(3.11)
1. τH2O(λ)的确定
引用“可降水分”的概念,可降水分ω是由底面Sa和长度等于红外系统到目标距离R的圆柱体内大气含水汽凝结的水层厚度来量度的。也可表述为: 截面积为1cm2,长度等于全部辐射路程的水蒸气汽柱中所含水蒸气凝结成为液态水后的水后的水柱长度。重要的是,不要混淆给定厚度的可降水分的吸收和相同厚度的液体水的吸收。10mm厚的液体水层,在超过1.5μm的波段上,实际就不透过了,而对含有10mm可降水分的路程的透过率则超过60%。ω0为每km路程的可降水分(相对湿度Hr=100%时),可查表3.11,由绝对湿度Ha[大气中水蒸气的量,(g/m3)]而得,即
ω0=110ρHa(cm/km)(3.12)
式中,ρ为水的密度。


表3.11标准大气的绝对湿度

单位: g/m3


t/℃0123456789

-400.12000.10750.09620.08610.07690.06870.06120.05450.04850.0431
-300.3410.3080.2790.2520.2270.2050.18490.1660.1490.134
-200.8880.8100.7380.6720.6110.5560.5050.4580.4150.376
-102.1541.9711.8081.6581.5201.3931.2751.1661.0660.973
-04.844.474.133.823.523.252.992.762.542.33
04.845.185.555.946.356.797.257.748.268.81
109.3910.0010.6411.3312.0512.8113.6114.4515.3416.28
2017.318.319.420.521.723.024.325.727.228.7
3030.332.033.735.637.539.541.643.646.148.5
4051.053.656.359.262.165.268.471.875.378.9

ω=ω0R(3.13)

ω0也可查表3.12得到。


表3.12Hr=100% 时,不同温度下每千米大气中的可降水分厘米数

单位: cm/km


t/℃00.20.40.60.8

00.4860.4930.5000.5070.514
10.5210.5280.5350.5430.550
20.5570.5650.5730.5800.588
30.5960.6040.6120.6210.629
40.6370.6460.6550.6630.672
50.6810.6900.7000.7090.719
60.7280.7380.7480.7580.768
70.7780.7880.7980.8080.818
80.8280.8390.8510.8620.874
90.8850.8960.9070.9190.930
100.9410.9530.9650.9780.990
111.0021.0151.0281.0421.055
121.0681.0821.0951.1091.122
131.1361.1501.1651.1791.194
141.2081.2231.2381.2531.268
151.2831.2991.3161.3321.349
161.3651.3821.3991.4151.432
171.4491.4671.4851.5031.521
181.5391.5581.5721.5971.613
191.6321.6521.6721.6921.712
201.7321.7531.7731.7941.814
211.8351.8571.8971.9011.923
221.9451.9631.9912.0132.036
232.0592.0832.1082.1322.157
242.1812.2062.2312.2552.280
252.3052.3322.3592.3862.413
262.4402.4672.4952.5222.550
272.5772.6072.6362.6662.695
282.7252.7752.7852.8152.846
292.8762.9082.9412.9733.006
303.038
-300.046
-290.0500.0490.0480.0460.045
-280.0540.0530.0520.0520.051
-270.0590.0580.0570.0560.055
-260.0650.0640.0630.0610.060
-250.0700.0690.0680.0670.066
-240.0760.0750.0740.0720.071
-230.0840.0820.0810.0790.078
续表


t/℃00.20.40.60.8

-220.0910.0900.0880.0870.085
-210.0990.0970.0960.0940.093
-200.1080.1060.1040.1030.101
-190.1170.1150.1130.1120.110
-180.1270.1250.1230.1210.119
-170.1370.1350.1330.1310.129
-160.1490.1470.1440.1420.139
-150.1610.1590.1560.1540.151
-140.1740.1710.1690.1660.164
-130.1880.1850.1820.1800.177
-120.2030.2000.1970.1940.191
-110.2190.2160.2130.2090.206
-100.2370.2330.2300.2260.233
-90.2550.2510.2480.2410.241
-80.2740.2700.2660.2630.259
-70.2950.2910.2870.2820.278
-60.3180.3130.3090.3040.300
-50.3410.3360.3320.3270.323
-40.3670.3620.3570.3510.346
-30.3940.3890.3830.3780.372
-20.4230.4170.4110.4060.400
-10.4530.4470.4410.4350.429
-00.4860.4790.4730.4660.460

由已知空气温度tB,查表得到ω0,然后乘以实际空气相对湿度Hr,即得到此实际空气相对湿度下的可降水分ω,相应的可降水分为
ω=ω0HrR=ω0R·Hr(3.14)
由可降水分ω,查海平面水平路程上水汽光谱透过率表,即可求得不同波长时所对应的τH2O(λ)。如果从表中不能直接查到,可通过外推法或内插法求得。
2. τCO2(λ)的确定
研究实验证明,CO2的密度在大气近表层中实际保持不变,直到非常高的高空,CO2在大气中的浓度是常数,因而它在大气中的分布随时间的变化是很小的,系数τCO2(λ)只与辐射通过的距离有关。因此,由CO2的吸收的造成的辐射衰减,可以认为与气象条件无关,τCO2(λ)只与辐射通过的距离有关。表3.7和表3.9列出了海平面路程τCO2(λ)的值: 在03~6.9μm、海平面水平路程0.1~1000km,以及在7.0~13.9μm、海平面水平路程0.2~200km的光谱透过率。τCO2(λ)可以通过查表得到。
在3~5μm和8~14μm波段内,为求取与吸收损耗有关的平均大气透过率τ-a,J.M.Lloyd采用近似式: 
τ-a=0.8326-0.0277ω0(3.15)
为评价在水汽和CO2中与可降水分ω0有关的平均衰减系数,在10μm±0.1μm、10μm±2μm范围内,可采用经验公式: 
α-a(10±0.1)=0.0124ω0+0.0088(3.16)

α-a(10±2)=0.0338ω0+0.045(3.17)
式中,α-a的单位为km-1。
3.6不同高度时的分子吸收修正问题
对不同高度,由于气温和气压的不同,辐射的吸收是不同的; 同时,分子密度也是不同的。因此应从这两方面考虑分子吸收修正问题。
3.6.1吸收本领随高度而改变所引起的修正
由3.5节可知,水汽对辐射的吸收会随气温和气压而变,因此对于高空的情况需要进行修正。也就是说,修正时,只需用修正系数乘以该高度处的水平距离,就得到等效海平面距离,并以此等效海平面距离计算沉积水厚度(可降水分)。修正系数βH2O由下式确定: 
βH2O=PP01/2T0T1/4(3.18)
式中,P0、T0为海平面上的气压和气温; P、T为给定高度上的气压和气温。
由式(3.18)可知,温度的影响很小(≤4%),可以忽略不计,因此,一般取高度修正系数: 
βH2O=(P/P0)1/2(3.19)
假定用ωe表示辐射传输路程中按吸收本领折算成大气近地层水汽的等效可降水分的有效厚度,用ωH表示H高度下可降水分层的实际厚度,则
ωe=ωHβH2O(3.20)
βH2O可以通过查修正系数表得到,或者在实际应用中,得到具有足够精度的近似值,可由下式确定: 
βH2O=e-0.0654H(3.21)
式中,H的单位为km。
对于CO2,类似可得到下列关系式: 
βCO2=(P/P0)1.5≈e-0.19H(3.22)

Re=R′HβCO2(3.23)
式中,R′H为在高度H上辐射传输的距离(此处R′H未考虑CO2质量引起的修正); Re为按吸收本领折算成近地层的有效距离。
3.6.2分子密度随高度而改变所引起的修正
由于高度不同,引起水平、倾斜路程中水汽量的不同,以及空气压强和质量的不同,因而分子密度也是不同的。为此,因高度变化带来的对可降水分和路程长度的影响需要进行修正。
1. 水平、倾斜路程中水汽量变化引起的修正
由式(3.14)和式(3.20)得
ωH=ω0H′rR=ωR·H′r(3.24)
式中,R为给定高度处的水平路程; ω0H′r为相应高度处的可降水分。ω0H′r是与相应高度处的温度、湿度有关的,忽略气温变化的影响,ω0H′r与湿度有关,而湿度随高度的分布服从下面的定律(对于标准大气): 
Ha,H=Ha,oe-βH(3.25)
式中,Ha,H为高度H处的绝对湿度; Ha,o为近地处或海平面处的绝对湿度; β为045/km。
由相对湿度和热力学温度的定义可知: 
H′r=k·Ha,H(3.26)

Hr=k·Ha,0(3.27)
式中,H′r和Hr分别为高度H处、海平面处的相对湿度;  k为一定温度下饱和空气中的水蒸气质量,仅与温度有关,在忽略高度H处和海平面处的温度差异时,可知不同高度的k值相同。
由式(3.25)、式(3.26)和式(3.27)可得
H′r=Hre-βH=Hre-0.45H(3.28)
将式(3.28)代入式(3.24)得到
ωH=ω0R·Hre-0.45H(3.29)
综合考虑水汽的吸收本领和水汽量随高度的变化,结合式(3.20)和式(3.21),得到确定距海平面上高度为H的辐射沿水平传输路程中的可降水分的有效厚度的公式为
ωe=ω0R·Hre-0.45H·e-0.0654H=ω0R·Hre-0.5154H(3.30)
在倾斜路程中的可降水分的有效厚度可由下面的方法计算。根据式(3.30),位于高度h处的大气元层ds(见图3.4)中的可降水分有效厚度为
dωe=ω0Hre-0.5154Hds(3.31)
式中,h=scosγ,为高度的瞬时值; γ为地面(海平面)法线与辐射传输方向之间的夹角。



图3.4倾斜路程中可降水分的有效厚度计算图


在不考虑地面曲率的情况下,对式(3.31)积分可求得从高度H1到H2的倾斜程上大气中可降水分的有效总厚度为
ωe=ω0Hr∫s2s1e-0.5154scosγds=ω0Hre-0.5154H1-e-0.5154H20.5154cosγ(3.32)
式中,H1=s1cosγ,H2=s2cosγ。如果H1=0,H2=H,则有
ωe=ω0Hr1-e-0.5154H0.5154cosγ(3.33)
值得一提的是,令γ=0,就可以由式(3.32)求得垂直路程上大气中可降水分的有效总厚度。
2. 水平、倾斜路程中因空气压强和质量引起的修正
因为空气压强随高度的变化规律如下: 
PH=P0e-0.123H(3.34)
考虑到式(3.34),可按下式把距地面(海平面)高度为H的水平路程折算成近地水平路程: 
R′H=RHe-0.123H(3.35)
式中,RH为高度H上的辐射传输距离; R′H为等效近地水平路程。
在倾斜路程的情况下,与式(3.32)、式(3.33)类似,可得到折算成近地层辐射路程的公式,
在H1到H2范围内
R′H=e-0.123H1-e-0.123H20.123cosγ(3.36)
综合考虑CO2的吸收本领和质量随高度的变化,结合式(3.23)和式(3.35),得到折算成近地层的路程有效长度为
Re=RHe-0.123H·e-0.19H=RHe-0.313H(3.37)
在倾斜路程的情况下,按CO2的吸收本领折算成大气近地层的路程有效长度的计算公式,由式(3.37),类似得到折算成近地层的路程有效长度。
对于高度从H1到H2的范围
Re=e-0.313H1-e-0.313H20.313cosγ(3.38)
对于高度从0到H情况下,即H1=0,H2=H,则有
Re=1-e-0.313H0.313cosγ(3.39)
3.6.3纯吸收时的透过率计算程序
仅考虑纯吸收的透过率可按下述程序计算。
1. 对水蒸气
(1) 根据温度,查Hr=100%时,地面上每千米大气中的可降水分厘米数表,或根据表3.6、表3.7,求出地面上的可降水分,得到ω0。
(2) 根据高度进行辐射传输距离修正(吸收本领和大气本身密度随高度而减小所产生的影响)。
(3) 求出全路程的可降水分ωe,即ωe=ω0Hr·(修正以后得到的近地层有效距离)。亦即,根据不同的要求,由式(3.30)或式(3.32)计算ωe。
(4) 查海平面水平路程上水蒸气的光谱透过率表,得到仅考虑水蒸气时的大气透过率τH2O(λ)。
2. 对CO2
(1) 同上,根据高度进行路程距离修正,亦即,按不同的要求,由式(3.37)或式(3.38)计算成近地层的路程有效长度Re。
(2) 查海平面水平路程上的CO2光谱透过率表,得到仅考虑CO2时的大气透过率τCO2(λ)。
3. 连乘得到纯吸收时的透过率τ1(λ)
τ1(λ)=τH2O(λ)·τCO2(λ)(3.40)
3.7大气分子与微粒的散射
大气中传输的辐射通量,同样经受空气分子散射(大气分子散射和微粒散射),微粒散射即为仅存在于大气中的地球表面灰尘、烟雾、水滴、盐粒等不同粒子的散射。分子散射可以作较精确的计算,而微粒散射与大气状态有关。
由分子、气溶胶、雾、霾和云层引起的散射可用米氏散射理论来描述。米氏散射理论适用于比波长小很多的小微粒散射,即瑞利散射,并且也适用于微粒大小比波长大得多的微粒散射,即非选择性散射。若分子或微粒的尺寸小于波长,则是遵守λ-4律的瑞利选择散射; 对于尺寸比波长大的微粒可看作非选择散射。
由于分子热运动造成分子密度的局部涨落所引起的散射称为分子散射。此时有瑞利定律,即线性散射系数σλφ为
σλφ=8π3(n2-1)3λ4N(1+cos2φ)(3.41)
式中,σλφ为与波长为λ的辐射通量成φ角的方向上的线性散射系数; n为大气折射率; N为分子密度(cm-3)。
瑞利定律适用于散射质点的直径远小于波长的情况。对于直径不超过0.42μm的小质点,瑞利定律成立。对于红外波段来说,这种散射是很弱的,在大多数情况下,在每单位体积内仅散射原光束能量的10-6~10-7。因为瑞利散射与λ-4有关系,所以当波长超过大约2μm时,分子散射就没有意义了。结果是与吸收相比,因此微粒大小远小于波长的分子散射是没有意义的,而只有分子聚集才产生散射。
实际上大气并不是干净的,而含有各种液态或固态的杂质,因此它不服从瑞利定律,而服从米氏散射理论。在此情况下,它的线性散射系数为
σλ=2πNr2K(a)(3.42)

a=2πrλ(3.43)
式中,r为散射粒子的半径; K(a)为波长的函数,并由介质的折射率决定。对于水滴来说,K(a)的具体数值如表3.13所示。
在不同能见度的情况下,大气中所含微粒的情况如表3.14所示。计算σλ时,应根据当时当地的气象条件,选好N和a,则可按式(3.42)计算不同波长的线性散射系数,从而得到散射情况下的透过率。


表3.13K(a)函数表



aK(a)aK(a)aK(a)

10.0681.7151.35
20.0791.4201.35
30.82101.12251.15
41.45110.89301.09
51.8311.20.82351.06
61.95121.00401.03
71.90131.25


表3.14大气所含微粒情况



大 气 情 况大多数微粒直径/μm直径范围d/μm微粒数N/cm-3气象能见度/m

海烟161~451~52750
海雾16~321~451~1090~360
地烟0.20.2~1.41~200约9000
地雾(非工业区)3~181~401~300090~360
地雾(工业区)2~101~1001~8000约50
压层云2~101~70450~9000<50

由此可知,利用线性散射系数计算由于散射引起的大气衰减时,均需知道大气中悬浮粒子的材料、大小以及密度等详细资料,而这些资料又是很难确定和测量的。因此在气象学中,采用一种实验方法来处理散射问题,这就是如下所述的利用气象能见度求光谱线性散射系数的方法,即采用依标准气象能见度DV确定的实验数据来计算光谱透过率,这种方法在应用时既方便,又可靠。
气象能见度DV表征大气的模糊度,并且是白天能看见天空背景下水平方向上角尺度大于30′模糊物体的最大距离。它代表了大气的透射性能在可见光区的指定波长λ0处(通常取λ0=0.555μm或λ0=0.61μm),目标和背景之间对比减弱的程度。在这些波长处,大气的吸收为零,因而影响透射的原因将只是散射这一种因素。
表3.15列出了国际能见度等级和与之相应的能见度、透明系数及线性消光系数的值,以供参考。


表3.15国际能见度等级



等级能见度特征气象能见度DV/km观 察 条 件τ*μ



0
1
2很差

<0.05浓雾<10-34>78
0.05~0.2大雾、稠密的大雪10-8.519.5
0.2~0.5中雾、大雪10-3.47.8


3
4差
0.5~1薄雾、中雪0.023.9
1~2暴雨、中等薄雾或雪0.141.95


5
6中等
2~4大雨、小雾或小雪0.380.98
4~10中雨、很小的雾或雪0.680.39

7好10~20无沉积物或小雨0.820.195
8很好20~50无沉积物0.920.078
9非常好>50纯洁空气>0.92<0.078

线性消光系数μ=线性散射系数σ+线性吸收系数α

即透过率
τ=e-μR=e-(σ+α)R(3.44)
当距离R=1km时的透过率称为透过系数,用τ*表示,即
τ*=e-μ(3.45)
表3.16列出了国际能见度等级与散射系数的关系,可供参考。


表3.16国际能见度等级与散射系数关系表



等级气象能见度DV/km散射系数γ/km-1

0<0.05>78.2
10.05

0.278.2

19.6
20.2

0.519.6

7.82
30.5

17.82

3.91
41

23.91

1.96
52

41.96

0.954
64

100.954

0.391
710

200.391

0.196
820

500.196

0.078
9>500.078
纯净空气277km0.0141

注: 表中数值是由分子散射构成。

眼睛感知的最小对比度(阈值对比度)等于2%,因此,气象能见度DV就是目标对比度为1时,通过大气后感知的对比度为0.02的距离,即
KV(DV)=KV(0)exp(-αVDV)(3.46)
式中,KV(0)=1,KV(DV)=0.02。

由此得
αV=-1DVlnKV(DV)KV(0)=3.91DV(3.47)

式中,DV的单位为km; αV的单位为km-1。
由实测结果确定,在0.3~14μm区间衰减系数与散射αp~λ-q有关,因此,依赖大气分子微粒散射的光谱散射系数αp(λ)可用下式得到
αp(λ)=αV(λ0)λ0λq=3.91DVλ0λq(3.48)

式中,λ0取0.555μm(0.61μm); λ为红外辐射波长(μm)。对特别好的能见度(例如,DV>50km时),修正因子q=1.6; 对于中等能见度(DV=10km),q=1.3; 如果大气中的霾很厚,以致能见度很差(DV<6km)时,可取q=0.585D1/3V; 对于0.3~14μm区间,一般可取q=1.3。表3.17详细列出了不同气象能见度对应的修正因子。


表3.17不同能见度对应的修正因子



修正因子q气象能见度DV/km等级气 象 状 况

1.6>509极晴朗
1.36<DV<506~8晴朗
0.16DV+0.341<DV<64~6霾
DV-0.50.5<DV<13轻雾
0DV<0.5<3中、浓、厚雾

从而由式(3.48)求得纯粹由散射导致的透过率τ2(λ)为
τ2(λ)=exp(-αpR)=exp-3.91DVλ0λq·R(3.49)
式中,R为作用距离。
在有霾的情况下,人们根据多次观察数据,提出一个由于霾的散射所造成的大气透过率经验计算公式: 
τ′2= (0.998)ω(3.50)
式中,ω为可降水分(cm/km)。
3.8与气象条件有关的衰减
因为气象(雾、雨、雪)粒子尺寸通常比红外辐射波长大得多,所以根据米氏散射理论,这样的粒子产生非选择的辐射散射。
雾粒的尺寸各有不同。虽然辐射在雾中有衰减,但随波长的变化,比在大气分子和粒子散射时要弱些,而不比8~14μm透过窗口小,通常可见光区衰减比它大2~2.5倍。
对于小粒状雾,光谱衰减系数可按式(3.48)近似计算。
在可见光和红外光谱区,雨和雪的辐射衰减与雾的衰减有别,是非选择性的。因此,对于决定与其强度相关的雨、雪的衰减系数可采用在10.6μm波长得到的经验公式: 
α雨=0.66J0.66雨(km-1)(3.51)
α雪=6.5J0.7雪(km-1)(3.52)
式中,J雨、J雪分别为与气象条件有关的降雨、降雪强度(mm/h)。
在没有实测数据的情况下,可采用下面的降雨强度数值进行计算: 小雨,2.5mm/h; 中雨,12.5mm/h; 大雨,25mm/h。
由雨的衰减所导致的透过率为
τ′3=exp(-α雨·R)(3.53)
由雪的衰减所导致的透过率为
τ3(λ)=exp(-α雪·R)(3.54)
式中,R为作用距离。
3.9平均透过率与积分透过率的计算方法
对于非连续的离散型数据,可采用将积分形式变为求和形式的方法来计算平均透过率与积分透过率。
3.9.1平均透过率的计算方法
如上所述,首先求出τ1(λ)、τ2(λ)、τ3(λ),由式(3.9)即可求出大气光谱透过率τa(λ),然后将平均透过率的积分形式变为求和形式,即由式(3.6)变为下式: 
τ-a=1Δλτa(λ1)×12dλ+τa(λ1+dλ)dλ+…+τa(λ1+(n-1)·dλ)dλ+12τa(λ2)dλ

=dλΔλ12(τa(λ1)+τa(λ2))+∑n-1i=1τa(λ1+i·dλ)
(3.55)
式中,λ1~λ2为光谱范围; Δλ=λ1~λ2; dλ为光谱间隔,亦即求和间隔; n=Δλ/dλ,n应为正整数,n越大,τ-a的精确程度越高。
将有关数据代入式(3.55),可以求出平均透过率τ-a。
3.9.2积分透过率的计算方法
积分透过率的计算方法如下: 
(1) 用式(3.9)求出波段内的大气光谱透过率τa(λ)。
(2) 确定大气积分透过率。
把式(3.8)中的积分形式变为求和形式,即用下式确定大气积分透过率: 
τa=∑N-1i=0φe(λ1+i·dλ)τa(λ1+i·dλ)∑N-1i=0φe(λ1+i·dλ)(3.56)

式中,dλ为求和间隔; N=(λ2-λ1)/dλ。
目标辐射通量由下式确定(对朗伯灰体辐射源): 
e(λ)=ε0Me(λ,T)S0(3.57)
式中,ε0为灰体辐射系数; Me(λ,T)为绝对黑体的辐射出射度; S0为目标辐射的面积; T为目标热力学温度。
将式(3.57)代入式(3.56),得
τa=∑N-1i=0Me(λ1+i·dλ)τa(λ1+i·dλ)∑N-1i=0Me(λ1+i·dλ)(3.58)
(3) Me(λ,T)的计算。
在红外系统中,Me(λ,T)的常用形式为
Me(λ)=c1/[λ5(exp(c2/λT)-1)](3.59)
式中,c1为第一辐射常量; c2为第二辐射常量; λ为辐射波长(μm)。
(4) 求τa。
由Me(λ1+i·dλ)和τa(λ1+i·dλ),代入式(3.58)而求得τa。
3.10计算示例
例3.1气象条件如下: 水平路程距离D=2km,气象能见度DV=20km(在λ0=055μm处),空气温度tB=20℃,空气相对湿度Hr=80%,无雾和雨。计算λ1=4μm和λ2=10μm时的大气光谱透过率。
解: 
(1) 求τH2O(λ)。
根据tB=20℃,查表3.12得到空气相对湿度Hr=100%时的可降水分,ω0(100%)=17.3mm/km,由式(3.12)可知,相对湿度Hr=80%、路程距离D=2km的可降水分为
ω=ω0DHr=27.7mm
从表3.6和表3.8用插值法或外推法可分别求出: 
τH2O(λ1)=0.952,τH2O(λ2)=0.718
(2) 求τCO2(λ)。
路程距离D=2km时,从表3.7和表3.9中查得
τCO2(λ1)=0.991,τCO2(λ2)=0.997
(3) 求τ2(λ1)、τ2(λ2)。
根据式(3.49)求出纯粹由散射导致的透过率为
τ2(λ1)=0.972,τ2(λ2)=0.991
(4) 求τa(λ1)、τa(λ2)。
由于没有气象衰减,所以τ3(λ)=1,因此由式(3.9)得到
τa(λ1)=0.952×0.991×0.972=0.917
τa(λ2)=0.718×0.997×0.991=0.709
例3.2气象条件如下: 水平路程距离D=2km,气象能见度DV=10km(在λ0=055μm处),空气温度tB=5℃,空气相对湿度Hr=85%,小雨强度为1mm/h。计算8~14μm范围内的大气光谱透过率。
解: (1) 求τH2O(λ)。
根据tB=5℃,查表3.12得到空气相对湿度Hr=100%时的可降水分,ω0(100%)=68mm/km,由式(3.12)可知,相对湿度Hr=85%、路程距离D=2km的可降水分为
ω=ω0DHr=11.6mm
由表3.8用插值法或外推法,可分别求出8~14μm范围内的τH2O(λ),结果列于表3.18。
(2) 求τCO2(λ)。
由表3.9,对路程距离D=2km时,可确定8~14μm范围内的τCO2(λ),结果列于表3.18。


表3.18例3.2的大气光谱透过率



λ/μmτH2O(λ)τCO2(λ)τ2(λ)τa(λ)

8.00.56510.9760.147
8.40.74610.9780.195
8.80.82010.9790.214
9.60.8580.9610.9810.216
10.00.8670.9970.9820.226
10.40.8690.9980.9830.228
10.60.8710.9990.9830.228
11.00.8660.9970.9840.227
11.80.8420.9930.9860.220
12.60.8490.8150.9870.182
13.00.8260.9120.9870.199
13.60.7890.3510.9880.073
13.80.7740.2150.9880.044

(3) 求τ2(λ)。
根据式(3.49)求出纯粹由散射导致的透过率τ2(λ),结果列于表3.18。
(4) 求τ′3(λ)。
由式(3.51)得a雨=0.66×10.66=0.66km-1,由式(3.53)计算出小雨的衰减所导致的非选择气象透过率τ′3(λ)=e-2×0.66=0.267。
(5) 求τa(λ)。
由式(3.9)可计算出8~14μm范围内的大气光谱透过率τa(λ),结果列于表3.18,按照此表,可做出相应的大气光谱透过率曲线。为了说明计算方法,这里仅列出了有限的几组数据。

例3.3气象条件如下: 水平路程距离D=3km,气象能见度DV=15km(在λ0=055μm处),空气温度tB=27℃,空气相对湿度Hr=70%,无雾和雨。目标是温度to=27℃的朗伯灰体辐射源。计算在8~14μm范围内的大气积分透过率。
解: 与例3.2类似,首先用式(3.9)确定大气光谱透过率,然后由式(3.58)计算大气积分透过率。
(1) 求τH2O(λ)。
根据tB=27℃,查表3.12得到空气相对湿度Hr=100%时的可降水分,ω0(100%)=25.8mm/km,由式(3.14)可知,相对湿度Hr=70%、路程距离D=3km的可降水分为
ω=ω0DHr=54mm
由表3.8用插值法或外推法,可分别求出8~14μm范围内的τH2O(λ),结果列于表3.19。
(2) 求τCO2(λ)。
由表3.9,对路程距离D=3km时,用外推法确定8~14μm范围内的τCO2(λ),结果列于表3.19。
(3) 求τ2(λ)。
根据式(3.49)求出纯粹由散射导致的透过率τ2(λ),结果列于表3.19。
用式(3.9)计算出大气光谱透过率τa(λ),结果列于表3.19。


表3.19例3.3的大气光谱透过率



λ/μmτH2O(λ)τCO2(λ)τ2(λ)τa(λ)

8.00.0710.980.07

8.40.2610.980.25
8.80.4010.980.39
9.60.490.940.980.45
10.00.510.990.980.49
10.40.520.990.980.50
10.60.5210.980.51
11.00.510.990.980.50
11.80.440.990.990.43
12.60.470.740.990.34
13.00.410.870.990.35
13.60.340.270.990.09
13.80.310.150.990.05

(4) 计算大气积分透过率。
取求和间隔Δλ=0.5μm,由式(3.58)和式(3.59)计算大气积分透过率。计算得到
τa=0.39
例3.4气象条件如下: 水平路程距离D=1.8km,高度H=2km,气象能见度DV=13.8km(在λ0=0.61μm处),空气温度tB=2℃,空气相对湿度Hr=50%。计算2.0~25μm波段内的平均大气透过率。
解: (1) 求τH2O(λ)。
根据tB=2℃,查表3.12得到空气相对湿度Hr=100%时的可降水分,ω0(100%)=0557cm/km,由式(3.30)可知,相对湿度Hr=50%、D=1.8km、H=2km上的可降水分为
ωe=ω0R·Hre-0.515H=0.178cm≈0.2cm
由表3.6可查出2.0~2.5μm范围内的τH2O(λ),结果列于表3.20。
(2) 求τCO2(λ)。
由式(3.37)可得到有效海平面距离为
Re=RHe-0.313H=0.962km
查表3.7,τCO2(λ)结果列于表3.20。
(3) 求τ2(λ)。
根据式(3.49)求出纯粹由散射导致的透过率τ2(λ)。由于DV=13.8km,故取q=1.3。因为这个光谱带很窄,由散射导致的透过率随波长变化较慢,故取光谱带的中心波长λ=225μm处的τ2(λ)作为平均的τ2(λ),得
τ2(2.25)=91%


表3.20例3.4的大气光谱透过率和平均大气透过率



λ/μmτH2O(λ)τCO2(λ)τ2(λ)τ-a


2.00.9330.931
2.10.9820.994
2.20.9911.000
2.30.9821.000
2.40.9111.000
2.50.6951.000
0.91
0.845

(4) 求τ-a。
由式(3.55)可计算出2.0~2.5μm范围内的平均大气光谱透过率τ-a。取Δλ=05μm,而光谱间隔dλ=0.1μm。得
τ-a=0.845
这里所述的计算方法,仅需要查阅有关的基本数据表,就可以顺利进行在多种情况下对大气透过率的较准确计算,尤其是对大气窗口的透过率计算,并且便于工程应用,还可以在计算机上编制软件进行计算。
3.11几种大气辐射传输计算软件应用比较分析
当前国内外已有多种大气辐射传输的模式与算法,常见的有LOWTRAN(Low Resolution Transmission)、FASCODE(Fast high resolution Code)、MODTRAN(Moderate resolution Transmission)、PCMODWIN、DISORT(Discrete Ordinate Method)大气辐射传输计算软件。此外还有6S(Second Simulation of Satellite Signal in the Solar Spectrum)软件、SBDART(Santa Barbara Disort Atmospheric Radiative Transfer)软件、中国科学院安徽光机所的通用大气辐射传输软件CART(Combined Atmospheric Radiative Transfer)等。
大气辐射传输软件主要包括大气分子吸收谱线数据库[如HITRAN(HIghresolution TRANsmission)]、逐线积分程序[如LBLRTM(Linebyline Radiative Transfer Mode)]、谱带吸收透过率程序(如MODTRAN、6S)、辐射传输方程(多次散射)求解程序(如CART、DISORT、SBDART),以及矢量(偏振)辐射传输方程求解程序[如6SV1(6S的矢量版)]等。从实用角度看,常用算法软件主要分为两大类: 一类是以MODTRAN为代表的“应用型”软件; 另一类是以DISORT为代表的“数学型”软件。
就低分辨率(分辨率≥20cm-1)大气辐射传输模式而言,美国空军地球物理实验室(AFGL)开发的低分辨率透过率计算模式(LOWTRAN)程序是公认的有效而方便的大气效应计算软件,是一种FORTRAN的计算机程序。LOWTRAN从1970年提出至今已7个版本。LOWTRAN从一开始就是为应用而建,在40多年的发展过程中不断扩充和修订基础资料,改进算法,增加可计算的辐射传输结果,从原意义上的“低分辨率大气透过率计算模式”扩展到目前能导出复杂天气条件下多种辐射传输量的“低分辨率大气辐射传输计算模式”,提供了许多新的应用可能性,已被国际上许多应用专家广泛应用于各自的实际问题。它的主要用途是为了军事和遥感的工程应用,以20cm-1的光谱分辨率的单参数带模式计算0~50000cm-1的大气透过率、大气背景辐射、单次散射的阳光和月光辐射亮度,亦即用于计算低频谱分辨率(20cm-1)系统给定大气路径的平均透过率和路程辐射亮度。LOWTRAN可计算从紫外到微波(0.2μm~∞波段)的大气传输问题。早期的LOWTRAN主要用于计算光在大气中传输的透过率。后来的版本可用于计算给定倾斜路程的大气透过率、大气辐射及背景由太阳光和月光散射的辐射。应用于: 大气透过率计算、背景辐射计算、探测几何路径、大气折射、吸收气体含量、非对称因子的光谱分布。特点及优势: 计算迅速、结构灵活多变,可作为下层大气的辅助工具、地表面战术系统的辅助工具等。
LOWTRAN 7版本是一款功能非常强大的大气透过率和背景辐射计算软件,一个单参数带模式的低分辨率大气透过率和背景辐射计算软件包,它能够有效帮助用户快速地计算出大气透过率和背景辐射,软件规模小巧,操作简单,支持大气温度、气压、密度的垂直廓线、水汽、臭氧、甲烷、一氧化碳和一氧化二氮等大气模式输入,是工程计算人员必备的一个辅助工具。包含了大气分子的吸收和散射、水汽吸收、气溶胶的散射和吸收、大气背景辐射,日光或月光的单次散射和地表反射、直接大气辐射以及日光、大气热辐射的多次散射等。大气模式设立了热带、中纬度夏季/冬季、近北极夏季/冬季、1976年美国标准大气及自定义模式选择。气溶胶消光扩充为城市型、乡村型、海洋型、对流层和平流层等多种模式,并考虑了对风速的依赖关系,建立了雾、雨和卷云的模型。由于该版本软件引入了多次散射模型,从实际计算效果看,多次散射的影响可能会比较明显,会损失图像信息的细节,这种现象在一些气象条件较差的情况下会表现出来。
LOWTRAN 7大致可分为三部分。
(1) 大气模式输入。包括大气温度、气压、密度的垂直廓线、水汽、臭氧、甲烷、一氧化碳和一氧化二氮的混合比垂直廓线及其他13种微量气体的垂直廓线,城乡大气气溶胶、雾、沙尘、火山灰、云、雨的廓线和辐射参数,如消光系数、吸收系数、非对称因子的光谱分布及地外太阳光谱。
(2) 探测几何路径、大气折射及吸收气体含量。
(3) 光谱透过率计算及大气太阳背景辐射计算(包括或不包括多次散射)。LOWTRAN 7共有5个主输入卡。LOWTRAN考虑因素较全面,只要给定温度、气压、水汽含量、气溶胶模型、能见度距离、辐射波长范围、路径长度和类型(水平或斜程)等,就能得到光谱透过率等结果。在一些仿真软件系统中,LOWTRAN模块已成为一个重要的组成部分。
LOWTRAN从1989年的版本7以后就发展成中光谱分辨率大气辐射传输模式(MODTRAN)。其目的在于改进LOWTRAN光谱分辨率,它将光谱的半宽度由LOWTRAN的20cm-1减小到2cm-1; 它的主要改进包括发展了一种2cm-1光谱分辨率的分子吸收的算法和更新了对分子吸收的气压温度关系的处理,同时维持LOWTRAN 7的基本程序和使用结构; 程序以卡片的形式来进行参数设置,操作起来清晰简洁。
MODTRAN是目前流行的红外辐射传输计算模型。应用FORTRAN语言编写源代码,设计了MODTRAN软件图形界面。针对FORTRAN语言计算效率高而图形功能弱、Visual Basic(VB)计算效率低而图形功能强的特点,用VB和FORTRAN两种语言混合编程,实现FORTRAN计算程序资源的再利用。应用于: 计算吸收物质的路径透过率、计算吸收物质的大气辐射率、计算吸收物质单次(多次)散射的太阳/月亮辐射率、计算中等光谱分辨率、计算大气透过率。特点及优势为: 模式选择性强; 辐射过程几乎考虑了大气中所有大气分子的吸收、散射和气溶胶、云的吸收和散射效应; 数据输入、结果输出方便等。
与LOWTRAN相比,MODTRAN不但提高了光谱分辨率,而且包括了多次散射辐射传输精确算法——离散纵标法,对有散射大气的辐射传输如太阳短波辐射,比LOWTRAN中的近似算法有更高的精度和更大的灵活性。该软件在以下各方面有广泛的应用: 
(1) 计算不同地面反射率、太阳天顶角、云类型、气溶胶条件和臭氧分布情况下的紫外、可见光波段的辐射; 
(2) 计算平流层上层的辐射场,并与气球辐射探测做比较; 
(3) 计算地表辐射场,并与地基辐射观测作对比; 
(4) 对紫外和可见波段,作不同近似条件下的辐射传输方程的精度评估。
MODTRAN 3版本使用说明与早期的LOWTRAN 7版本类似,但是加入了一些新的参数,增加了一张新输入卡片,对另外两张卡片也作了一些小改动。新参数主要是为了控制: 
(1) 第二个多次散射选项(根据DISORT多流算法); 
(2) 新的带可选的三角光滑函数的高分辨率太阳光谱辐照度(根据全部计算的辐照度); 
(3) 对CO2混合比的“步”更新。
目前版本为MODTRAN 5。
与LOWTRAN、MODTRAN同步发展的是FASCODE。FASCODE提供了标准大气模式、大量的气溶胶模式(包括雾、沙漠尘埃和海洋气溶胶等)、水和冰云模式。它假定大气为球面分层结构,其最佳分层已经达到辐射的出射度或透过率计算中的特定精度。它具有比LOWTRAN更高的光谱分辨率,相应的计算量也大一些。FASCODE适用于研究精确的单色波长和激光大气传输问题。对一般光电系统的分析和设计来说,LOWTRAN已具有足够的精度。
PCMODWIN程序是以LOWTRAN模式系列为基础而发展的Windows环境下运行的程序包,PCMODWIN的内核是MODTRAN 3,最初发行年份为1997年,包含LOWTRAN 7低等分辨率和MODTRAN中等分辨率辐射模式,可以进行大气透过率、辐射传输计算。该程序具有以下特点: 
(1) 模式选择性强,可任意选择LOWTRAN 7或MODTRAN辐射传输模式,在选择的模式下,可以计算吸收物质的路径透过率、大气辐射率、单次(多次) 散射的太阳/月亮辐射率和直接透过的太阳辐射等。
(2) 辐射过程几乎考虑了大气中所有大气分子的吸收、散射和气溶胶、云的吸收和散射效应。
(3) 数据输入、结果输出方便,既可以应用现有标准模式大气和模式气溶胶、云等,又可以由用户直接输入观测或指定资料进行模式计算,当某些观测数据(如CH4、NxO等)缺省时,这种输入方式对辐射传输计算特别有效。
(4) 该程序在Windows环境下操作,操作界面简便直观,输出结果丰富,既可以在Windows窗口下直观地看到模式计算结果的辐射图解,还可以按设计要求的透过函数进行模拟计算。当计算所得大气透过率或辐射随波数起伏较大时,可以进行过滤平滑,得到以ASCII码形式给出的结果,便于以后进行有关分析。
DISORT是一种能处理多次散射的辐射传输计算方法,用于散射吸收大气的单色辐射传输计算,是用FORTRAN语言编写的按照结构化程序设计的辐射传输计算专用软件。DISORT软件包最早是由Stamnes等人于1988年开发的,采用离散坐标法求解辐射传输方程,给出了完全稳定的解析解,已经成为公认的辐射传输精确算法的实用软件包。该软件可求解从紫外到微波段、垂直非均匀、各向异性并含热源的平面平行介质中的辐射传输问题,计算包括热辐射、散射、吸收、下边界双向反射和发射等物理过程。后来DISORT软件包又加入了大气模式参数,以及消光系数、像函数、光学厚度计算等多个子程序,使其成为使用方便的实用软件(UVSPEC)。UVSPEC是针对紫外和可见光波段(176~850nm)设计的,光谱分辨率为1nm,计算0~100km高度的散射、直射辐射通量和辐射强度。与LOWTRAN 7比较,离散坐标法虽是一种精确的辐射传输计算方法,但它的计算精度与所采用的流数有关,而且在视角很大的方向上,它考虑的是平面平行大气,没有考虑地球曲率和大气折射的影响,可能有较大的误差。