偏差原则在基于卫星红外资料反演大气廓线中的研究论文.pdfVIP

偏差原则在基于卫星红外资料 反演大气廓线中的研究 赵强 杨世植 陈新兵 (中国科学院安徽光学精密机械研究所 合肥 230031) 摘要遥感方程或大气辐射传输方程是属于第一类非线性Fredholm积分方程，这类方程的显著特征是不适定 的，尤其此类问题解的不稳定性，增加了反演的难度。因此需要研究控制反演不稳定性和提高反演精度的有 效方法。为了使解稳定，先滑参数作为约束因子是越要的。结合牛顿非线性迭代法反演大气廓线，利用偏 差原则来最优选择光滑参数．y，即在同步反演大气廓线的同时对1采用分步选代的原则。最后利用MODIS红 外资料进行反演试验，反演结果表明采用偏差原则选取1明显优于经验法。 关键词遥感偏差原则 牛顿非线性迭代 MODIS 引言 自60年代美国发射第一颗气象卫星以来，世界各国都十分重视气象卫星的发展，它在监视和预告天气 变化以及大气物理科学研究中都很有重要的意义。大气遥感很重要的一步就是将定标、定位后的卫星观测 数据转化为各种气象要素，例如大气温度、湿度及大气各种气体成分分布等，这一步的工作通常称为反演。 卫星遥感大气中，大气参数反演是一个较为复杂而又颇受重视的数学物理问题，在用牛顿非线性迭代法反 演大气廓线中，需要一个矩阵H。来对反演进行控制，H。可以是一个先验误差协方差矩阵的逆矩阵，如果卫 星观测量及先验误差协方差矩阵皆为高斯分布，则最小方差解就变成典型的最大似然解，如果取H～=1 I(I为单位矩阵)则有此得到的解为最小信息解，1称为光滑参数。由于^y的选取依赖于卫星观测量、观 测误差、通道的光谱特征及反演中大气初始状态的选取等诸多因素的影响，因此参数_y的选取十分困难。 题是参数1的选择，不合理的选择1导致近似解不收敛。本文采用偏差原则，比较客观、合理、有效地选 择参数^y及考虑到问题的非线性性。在同步反演大气廓线的I一时对_y采用分步迭代的原则，并且应用到 MODIS红外资料反演中，计算结果表明，采用偏差原则选择参数．y的方法明显优于经验法。 l反演方法及偏差原则 如果忽略大气散射的影响，则地气系统到达卫星探测器的红外辐射量为 3日幽’+R‘ (1) R=￡母，，一Ij4肋(0，p)+(1--8)Ij 其中R是光谱辐射量，P是气压，e是表面发射率，B是普朗克函数，T是大气透过率函数，下标S表 示地表，T。=T，2，f。方程(1)右端第一项表示地表发射辐射贡献，第二项表示大气发射辐射贡献，第三 项表示反射大气辐射贡献。第四项表示红外区的反射太阳辐射贡献。由于在红外区常有8。一l故可略去反射 大气辐射贡献，当入4．Oum，可忽略R’ (1)式可表示为大气温度廓线、湿度廓线、地表温度等大气参数的非线性算子方程： Y=F(X) (2) 对(2)式线性处理，我们采用一阶泰勒展开村(x)=Y(X0)+F，·(x—x。) 其中F’是Y对x的一阶导数矩阵。 假定有N个通道同时观测大气，定义Y=(Y。，Y：…，Y。)’，y为第i通道的辐射量或亮度温度。x=(T。 表温度相当于求解方程： 292 (3) I，(彳)：Ip。一，(z)02+llx—X。If2：极小 即[1t“-f∞订一I“一f似仆(^一A∥川A—I】l。)=rain (4) 其中Y…是卫星观测值向量，E是观测误差协方差矩阵，它包括仪器的噪声和模式误差，H是x．1的误差协 方差矩阵，H～=1I，上标T定义为矩阵转置。(4)式解表示为： 片“)=爿。+(F。7E‘F+y，)～．F一．E“}…一F(X。)J (5) 解x(1)与1的选择有关，对于不同的．y导致不同的解。现在的问题是：如何去确定一个最优的光 滑参数^，，在选择1的过程中又要体现方程(1)的非线性性，以达到最优的解。本文采用偏差原则来选择 ．y： 2 ㈣Y“)】_y…卜o