非线性迭代物理法反演大气湿度廓线研究.pdfVIP

非线性迭代物理法反演大气湿度廓线 吴雪宝,张凤英 (国家卫星气象中心 北京100081) 摘 要 利用ATOVSAdvanced TIROS-N Operational Vertical Sounder资料可以反演晴空和有云情况下 的大气温度 湿度廓线等大气参数 本文简要介绍了非线性迭代物理反演法 并利用该方法对2002年7 月的NOAA-16极轨气象卫星ATOVS观测数据进行了大气湿度廓线的反演试验 试验结果表明 大气相对 湿度反演参数与无线电探空湿度观测资料的均方根误差在 15-23%之间 1 引言 利用卫星遥感资料反演大气温度 湿度的研究始于 70年代 随着卫星大气探测技术的不断发展以 及计算机运算速度的不断提高 反演方法也因此得到不断的改进 在过去的近 30年中 发展了许多卫 星大气参数反演方法 但归纳起来 这些方法可以分为两大类 即统计反演法和物理反演法 统计反演法是以大气参数廓线和对应的卫星观测辐射强度之间的回归分析为基础的 回归系数通常 由一组历史的卫星仪器观测辐射强度值和无线电探空资料确定 然后再将计算的回归系数应用于卫星的 [1] 实时测量值 基于统计回归分析的 特征向量统计反演法 是应用最广的一种统计反演方法 由于该 方法迫使大气参数反演解处于气候平均值附近 因此 可以保证反演解的真实性和避免不稳定解的出现 统计回归反演的另一个优点是无需计算权重函数 因此运算速度快 物理反演方法试图从卫星仪器辐射强度观测值直接推导出大气温度和湿度廓线 因此 需要给出温 度 湿度廓线初估场 并利用快速辐射传递正演模式计算出初估场对应的光谱通道辐射值 通过对初估 场的调整 使计算的光谱通道辐射值逼近实际观测值 并给出最终反演解 因此 物理反演方法的求解 过程实际上是一个迭代过程 该类方法需要计算权重函数 应用较广的物理反演方法是 同步物理反演 法 [2] 和在此基础上改进的 最小方差同步物理反演法 [3] 该方法的特点是可以同时反演出大气温 度 湿度廓线和地面温度 法国科学院动力实验室Chedin等人提出了 改进初估值反演 方法[4] 国 家卫星气象中心对区域 TOVSTIROS-N Operational Vertical Sounder资料业务处理系统进行过多 次更新换代 用改进的同步物理反演法取代统计回归法[5 6] 从 1998年开始 随着ATOVS的业务运行 人们开始发展新的大气参数反演方法[7,8,9] ,不同的反演方法采用的初估场 快速精确的正演辐射传递模 式 云检测方法等也各不相同 相对于温度反演而言 大气湿度反演更加困难 大气温度反演方面的文 章相对较多 因此本文不作讨论 而着眼在湿度反演这个难点上 我们采用非线性迭代物理法进行大气 湿度廓线反演试验 为了更好地获取云区的大气湿度垂直结构 在反演处理中充分利用不受云影响的微 波通道测值 2 ATOVS资料介绍 1998年5月 13日第五代美国极轨气象卫星NOAA-15发射成功 该系列的NOAA-16和NOAA-17也分 别于 1999年和2002年成功升空 星上装载有先进的大气垂直探测器ATOVS它由2个探测仪器组成 即20通道的高分辨率红外辐射探测器HIRS/3 20通道的先进微波探测器 AMSUThe Advanced Microwave Sounding Unit分别由微波温度探测器 AMSU-A和湿度探测器AMSU-B组成 新一代极轨气象卫星的成功发射标志着卫星大气探测水平进入了一 个新阶段 其最显著特点是增加了微波探测通道 可以极大地改善有云情况下的大气温湿探测能力 使 得原先以晴空大气探测为主 过渡到以全天候大气探测为主的阶段 AMSU-B是首次搭载在NOAA极轨卫 星上的微波湿度探测器,由英国气象局制造 共有5个通道 其中1个通道位于窗区 1个通道位于弱水 汽吸收线 3个通道位于水汽吸收带 主要用于大气湿度探测 尤其是云区 AMSU-B是一个扫描探测器 扫描宽度约为2250公里 扫描周期为8/3秒 扫描步数为90星下点分辨率约为 15公里 表 1给出了 AMSU-B探测通道光谱特征及其主要探