GPS可降水量与MODIS可降水量回归分析.doc

GPS可降水汽与MODIS可降水汽回归性分析 侯建国 杨成生 张勤 蒋光伟 王宏宇 （长安大学地质工程与测绘学院，陕西 西安，710054；） 摘要：大气水汽含量是气象以及大气科学研究中一个重要的参数，精确测定大气水汽含量的分布及变化对准确预报降水和灾害性天气具有重要意义。本文研究了两种空间遥感手段GPS和MODIS传感器获取大气水汽的原理，并利用2006年7月至2009年7月西安地区七期GPS观测数据和地面气象实测数据，研究了GPS PWV与MODIS PWV之间的回归关系，建立了MODIS反演水汽值的校正模型，为GPS遥感水汽、MODIS遥感水汽及其融合研究提供了参考价值。 关键词：GPS 中分辨率成像光谱20世纪90年代初发展起来的地基GPS技术利用GPS卫星信号延迟量来反演大气可降水量（Precipitable Water Vapour, PWV），且具有较高的精度（1-2mm）[3]，为探测大气水汽提供了一个有效的新途径，并成为各气象部门天气预报的主要数据源之一。然而，地基GPS系统在大范围内的建立，仍需要较长的一段时间，且GPS站点密度也很难详细反映大气的空间水汽变化信息。中分辨率成像光谱（Moderate Resolution Imaging SpectroradiometerMODIS)作为美国地球观测系统(EOS)系列卫星的最主要探测仪器之一，汇集了迄今为止最先进的空间星载传感器技术和遥感应用技术通过MODIS近红外反演含量 （1） 其中，，，分别表示大气折射常数；为特定的干空气气体常数；为干空气总质量密度；和分别为水汽局部气压和可压缩系数；和分别表示站点高度和对流层顶层高度。上式等式右边第一项表示天顶静力延迟（Zenith Hydrostatic Delay，ZHD）；等式右边第二项表示天顶湿延迟（Zenith Wet Delay，ZWD）。对于天顶静力延迟可利用站点信息和地面气压通过下式得到[6]： （2） 这里，单位为毫米；为地面气压，单位百帕；为站点纬度，单位为度；为站点高程，单位千米。天顶静力延迟一般为2.3米左右，其对地面气压测量的敏感度为[7]，也就是说的地面气压测量误差将会导致1的ZHD误差，实际上地面气压测量精度常常优于[6,7]。 因此，当GPS站点的精密坐标被确定以后，可利用GAMIT或BERNESE软件对获取的GPS数据进行解算，从而计算出GPS卫星信号的ZTD。同时，由GPS站点地面坐标信息和气压根据（2）式可高精度地获取ZHD。进而再由（1）式分析，从ZTD中减去ZHD，即可得到GPS信号的ZWD。需要注意的是，GPS数据记录的大气延迟是在截止高度角范围内卫星至GPS接收机的斜路径延迟，需经映射函数将斜路径延迟转换为天顶延迟。得到ZWD后，就可通过下式转化成PWV[6]： （3） 其中为液态水密度；为大气平均温度，单位为开尔文；为通用气体常数；为液态水摩尔质量。为一转换系数，通常取值范围在6.0到6.5之间，但有些地区可达到7.0[6-8]。2001年Scheuler等通过大量的研究进一步对（3）式进行简化，将表示为[9]： （4） 这样转换系数就变成大气平均温度的函数 3、MODIS水汽反演原理 MODIS传感器具有36个光学通道，分布在0.414μm的电磁波谱范围内。MODIS水汽反演则利用了位于0.8～1.3间的5个近红外通道，如表1所列。通道2和通道5是两个非吸收窗口，用于反演植被和云。位于0．936的18通道为强水汽吸收通道，在干燥环境下对水汽的反演很有效，而位于0．905的17通道为水汽弱吸收通道，在潮湿环境或太阳高度角很低的情况下对水汽的反演很有用[10,11]。 表1 MODIS上用于反演水汽的通道 MODIS通道号 位置/ 宽度/ 2 0.865 0.040 5 1.240 0.020 17 0.905 0.030 18 0.936 0.010 19 0.940 0.050 利用MODIS数据反演水汽是根据辐射传输方程来推导的，辐射传输方程简化如下： （5） 式中是波长；是传感器的测量辐射；是大气上界的太阳辐射；是整层大气透射比，它包含了太阳—下垫面—传感器路径上的大气中含水量的信息；是下垫面反射率； 是路径的散射辐射。将上式两边同除以，得： （6） 尽管反馈机制中包括通过大气到地面的后向散射过程，但因气溶胶光学厚度在近红外光谱区很小，反馈作用很小，所以(6)式右边第二项可以忽略。则上式简化为： （