大气辐射与遥感-第七章2.

第七章大气辐射在遥感中的应用兰州大学大气科学学院专业必修课-《大气辐射与遥感》授课人：王天河2015·春季第七章大气辐射在遥感中的应用§7.1引言§7.2利用透射的太阳光进行遥感§7.2.1气溶胶光学厚度和尺度谱的确定§7.2.2确定臭氧总量§7.2.3临边消光技术§7.3应用反射的太阳光进行遥感§7.3.1卫星-太阳几何光学和理论基础§7.3.2臭氧的卫星遥感§7.3.3气溶胶的卫星遥感§7.3.4陆地表面的卫星遥感§7.3.5云的光学厚度和粒子尺度§7.4利用发射的红外辐射进行遥感§7.4.1地球表面温度的确定§7.4.2温度廓线的遥感§7.4.3水汽和痕量气体廓线的遥感§7.4.4云的红外遥感§7.4.5红外冷却率和地表通量的遥感§7.5利用发射的微波辐射进行遥感§7.5.1微波波谱和微波辐射传输§7.5.2由微波发射辐射确定降雨率和水汽§7.5.3微波探测器的温度反演§7.6利用激光和微波能量进行遥感§7.6.1后向散射方程§7.6.2激光雷达的差分吸收方法§7.6.3激光雷达的差分退偏方法§7.6.4用于研究云的毫米波雷达§7.4利用发射的红外辐射进行遥感mdIv(t,m)dt=Iv(t,m)-Bv(t)理论基础红外辐射传输基本方程向上辐亮度的解：其中红外谱区地表发射率接近于1。大气顶处卫星垂直对地（μ=1）测量的辐亮度：仪器仅能分辨有限的带宽，在某一波数间隔内的辐亮度按归一化形式表示：Iv(t,m)=Ivt*()exp-t*-tmæèçöø÷+Bv¢t()exp-¢t-tmæèçöø÷d¢tmtt*òIvt*()=evBvTs()Gvt()=exp-t()¶Gvt()¶t=-exp-t()Iv(0)=BvTs()Gvps()+BvTp()éëùû¶Gvp()¶pdpps0òIv=IvYv,v()dvv1v2òYv,v()dvv1v2ò响应函数的有效频谱间隔通常很小，以至于普朗克函数的变化微不足道，于是可以用表示不带来明显误差。气象卫星进行大气和地面遥感的基本原理：考虑响应函数的光谱透射比：天基遥感关键：从实测的红外辐亮度中提取大气状态（温度、温度廓线）和大气组成（水汽、痕量气体、云）的有关信息。Iv(0)=BvTs()Gvps()+BvTp()éëùû¶Gvp()¶pdpps0òGvp()=Yv,v()exp-1gkv¢p()q¢p()d¢p0pòéëêùûúdvv1v2òYv,v()dvv1v2òBvTs()（＊）热红外光谱中主要吸收光谱：CO2：15μm，4.3μmO3：9.6μmH2O：整个红外谱区，最重要6.3μm和波数小于500cm-1CH4：7.6μmN2O：7.9μmCFCs：大气窗区有吸收线大气窗区：800-1200cm-1，除9.6μm的臭氧吸收带外，大气气体的吸收表现为一极小值。§7.4.1地球表面温度的确定大气窗区向上的辐亮度与地表的发射辐射密切相关。为了确定地表温度，引入大气分裂窗技术，即利用两个通道上的实测值消去Ta项。在10.5～12.5μm大气窗区，Ta的变化幅度小于1k，同时大气窗区地表发射率变化很小。透射比主要是水汽的连续吸收形成的：两个通道的窗区方程：（7.4.11）Iv=BvTs()Gvps()+BvTa()1-Gvps()éëùûGvt()»exp-kvm()=1-kvmI1=B1Tb1()=B1Ts()G1+B1Ta()1-G1[]I2=B2Tb2()=B2Ts()G2+B2Ta()1-G2[]（7.4.13）（7.4.14）为了消去Ta项，利用Ta的泰勒级数展开T的普朗克函数：消去T-Ta：利用亮温Tb2和Ts代替T，并利用（7.4.14）分析得到：从(7.4.17)和(7.4.13)中消去B1(Ta)，得到分裂窗方程：在窗区的小光谱间隔内，两亮温间可建立一种局部线性关系：用10.9和12μm确定地表温度，夜间加上3.7μm增加反演精度。BiT()=BiTa()+¶Bi¶TT-Ta()B2T()=B2Ta()+¶B2/¶T¶B1/¶TB1T()-B1Ta()éëùû（7.4.15）（7.4.16）B1Tb2()=B1Ts()G2+B1Ta()1-G2[]（7.4.17）B1Ts()=B1Tb1()+hB1Tb1()-B1Tb2()éëùûh=1-G1G1-G2»k1k2-k1Ts=Tb1+hTb1-Tb2()（7.4.18）（7.4.20）§7.4.2温度廓线的遥感根据卫星热红外辐射观测推求大气温度廓线最早是由King(1956)提出的，证明了由辐亮度推求温度廓线的可行性。Kaplan(1959)证明了具有垂直分辨率的温度场可根据大气发射的谱分布导出，发展了上述遥感探测的概念。在光谱带翼区的观测可以感知大气深处，而在带中心的观测仅能感知大气的顶层。因此，通过适当选择一组不同的探测波数，就可以利用实测的辐亮度进行解译，最终求得大气得垂直温度分布。如果地表发射辐射对观测的辐亮度有重要贡献，则必须确定下垫面的温度。如果舍弃地表贡献项，则前向辐射传输方程可表示为：给出对应于不同波数和已知权函数的一系列实测辐亮度之后，温度廓线反演的基本问题是怎样求解函数由于普朗克函数随波数而变化，所以在反演中必须首先消除普朗克函数对波数的依赖关系。由于普朗克函数是光滑函数，而且在小光谱间隔内，可以用线性形式近似表示Iv(0)=BvTs()Gvps()+BvTp()éëùû¶Gvp()¶pdpps0òIv(0)=BvTp()éëùû¶Gvp()¶pdpps0òBvTp()éëùû（7.4.22）（7.4.23）BvT()=cvBvrT()+dv（＊）方程(7.4.22)可表示为即著名的第一类弗雷德霍姆方程。其中：由热红外发射观测值确定大气温度，发射源必须是一种含量相当丰富的已知气体，且其分布要比较均匀。否则，气体含量的不确定性将使我们不能根据观测结果来明确地确定温度。100km以下的地气系统中，有两种气体具有均匀丰富的含量，且在易于观测的光谱区具有发射带：(1)CO2，具有振转带谱线；(2)O2，具有微波自旋－转动带。gv=fp()Kvp()dpps0ògv=Iv-dvcvfp()=BvrTp()éëùûKvp()=¶Gvp()¶p（7.4.24a）（7.4.24b）(1)带中心时黑体温度逐渐下降，这与对流层温度随高度增加而降低有关。(2)波数690附近，温度出现极小值，这与较冷的对流层顶有关。(3)波数690之外，随波数下降，温度却又升高，这是由于平流层温度升高所致。(4)可据此选择一组探测波数，把对流层和平流层低层温度廓线大部分包含进去。云雨4号卫星上的红外干涉光谱仪观测到的15μm带附近按黑体温度给出的出射辐亮度谱对装载在NOAA2号卫星上的一组垂直温度廓线辐射仪(VTPR)计算的权函数廓线和透射比廓线的个例。权函数的每个峰值代表对向上辐亮度的最大贡献部分。668.5cm-1677.5cm-1695.0cm-1708.8cm-1725.0cm-1745.0cm-1云的去除在许多温度反演的情况中，主要的问题是探测辐射仪的视场内存在云，进行温度反演之前去除云的效应非常重要。考虑一对相邻的像素(扫描点)，像素1和2的辐亮度可表示为：如果相邻像素内探测通道的云光学性质相同，则相邻像素内的温度场也是相同的，于是有如果预先知道N＊值，则可根据以下方程确定晴空气柱辐亮度：（7.4.46）每个像素点的有效云量，即云量与发射率的乘积（7.4.47）（7.4.48）在辐射仪视场中消除云贡献的N＊方法须遵守的限定条件：(1)假定相邻都得云具有同样的高度、温度和光学性质；(2)有效云量η1和η2必须是不同的；(3)需要附加信息来根据N＊值确定晴空气柱辐亮度。为了确定N＊值，通常的办法是使用一个独立的微波通道：即注意：微波观测经常会受到低分辨率、降水云的污染和预估的地表发射率造成的误差的影响，在温度反演中检测和校正云的作用时，处理探测资料时需要许多经验方法和阈值检验。由6.3μm振转带遥感水汽对上式积分项进行分部积分如果已通过15μm和4.3μmCO2带反演了温度廓线，则未知数为透射比。如果选择水汽6.3μm振转带波数，则仅有的未知量为路径长度廓线u(p)。与温度的反演相比，推求气体的廓线的辐射传输方程更为复杂。§7.4.3水汽和痕量气体廓线的遥感Iv(0)=BvTs()Gvps()+BvTp()éëùû¶Gvp()¶pdpps0ò（＊）（7.4.51）（7.4.50）Gvp()=exp-kvp()dup()0pòéëêùûúdvDvDvòIv(0)=BvT0()éëùû-Gvp()¶Bvp()¶pdpps0ò理论上，温室气体(O3,N2O,CH4,CFCs)的浓度可以根据俯视观测的光谱仪来推求得到。然而：(1)臭氧在9.6μm带显现出振转谱线。由于浓度极大值位于平流层，目前还没有设计出反演使用的基于热红外辐射的俯视观测仪器。(2)N2O和CH4在7.9μm和7.6μm光谱区具有谱线结构。其他微量气体与此相同，而且由于混合比较小，这些气体浓度的遥感使用的是临边扫描技术。平行辐射仪或光谱仪，其视角是水平的，它接收来自较薄高度层的大气辐射，称为临边辐射。由临边观测获得痕量气体廓线的一般方法叫做临边探测。上图为临边探测的几何光路。大气临边扫描探测优点：(1)由于大气密度和压力随高度快速递减，发射辐射主要在切点之上几千米的范围。因此可获得固有的高垂直分辨率。(2)接受的所有辐射都来自大气，相比朝天底观测的仪器，不存在因下垫面变化而引起的变化。(3)沿水平路径有很大的不透明性。对确定中层大气的微量气体特别有用。(4)卫星的仪器视场方向可取在相对于卫星运动的任何方位角方向，并且能覆盖相当大的范围。缺点：沿光线路径上高云的干扰；较大水平范围上急剧变化大气状态的解译问题。考虑一个由云层构成的视场，云层位于温度为Tc的气压层pc上，云量为η，云上的发射率为εv。为了说明问题，仅考虑单色辐射传输，卫星观测的辐亮度可以写成：式中右端四项的意义：地表、云下大气、云层、云上大气发射辐射对大气顶处的贡献。§7.4.4云的红外遥感卫星地面视场Tc,pc,h,evˆIv=1-hev()BvTs()Gvps,pc()+Bvp()¶Gvp,pc()¶ppspcòdpéëêùûúGvpc,0()+hevBvTc()Gvpc,0()+Bvp()¶Gvp,0()¶ppc0òdp（7.4.55）单色透射比定义为：根据单色透射比的增加遵循指数运算的规则：对于晴空视场：从有云视场中减去晴空视场，即云的辐亮度信息：进行分部积分：Gvps,pc()´Gvpc,0()=Gvps,0()Gvp,pc()´Gvpc,0()=Gvp,0()Gvp1,p2()=exp-1gkvp()qp()dpp1p2òéëêùûúIvclr=BvTs()Gvps,0()+Bvp()¶Gvp,0()¶pps0òdp（7.4.57）ˆIv-Ivclr=-hevBvTs()Gvps,pc()-hevBvp()¶Gvp,0()¶ppspcòdp+hevBvTc()Gvpc,0()ˆIv-Ivclr=hevGvp,0()¶Bvp()¶ppspcòdp（7.4.58a）（7.4.58b）对云顶气压和发射率的CO2切片技术将方程(7.4.58b)应用到15μmCO2带的两个谱带上：即CO2切片技术，利用两个CO2通道在相同视场中观测的云信息的比值来确定云顶气压的想法是Smith&Platt(1978)提出的。一旦确定云顶气压，就可以由红外窗区通道计算有效发射率，在窗区里（令波数指数为3）ˆI1-I1clrˆI2-I2clr=G1p,0()¶B1p()¶ppspcòdpG2p,0()¶B2p()¶ppspcòdpˆI3=he3B3Tc()+1-he3()I3clr（7.4.59）he3=ˆI3-I3clrB3Tc()-I3clr（7.4.60a）（7.4.60b）CO2切片方法对于高空卷云是最适用的，这是由于权函数在15μm带的性质使然。Wylie等(1994)利用此方法，根据NOAA