基于MODIS数据近红外波段反射大气水汽含量技术流程.docVIP

基于MODIS数据近红外波段反演大气水汽含量的技术流程 前言 水汽是大气中重要的气象参数，其含量随时间和空间变化很大。水汽体积混合比的变化范围一般在0.1%~3%之间，是大气温度变化范围内唯一可以发生相变的成分。水汽的分布、传输和季节变化对于研究水循环、全球气候变化、天气预报、遥感大气校正等具有重要意义。缺乏精确、稳定、长期的水汽数据记录，是阻碍深入研究气候系统中水汽影响的主要原因。全球气球探空观测数据是离散的水汽数据,需进行空间插值扩展后，才能应用到气候模式中。而“以点代面”会使模式模拟结果产生较大误差。卫星遥感反演的面状分布特点,使得卫星资料估计大气总水汽含量的分布,受到了广泛的重视。 基本原理 卫星传感器所接收的入瞳辐射亮度的简化形式可以表示为（Hansen and Travis, 1974; Franer and Kaufman, 1985）： （1） 其中为波长，为传感器所获得的入瞳辐射亮度，为大气顶层太阳辐射亮度，为总的大气透过率，是指辐射从太阳到达地表再从地表到达传感器所经过的大气路径的总的透过率，是指地表的二向反射率，是指程辐射。式（1）忽略了光子在地表的上的多次反射，即假设光子在地表上只反射一次。式（1）右边第一项可以看作是地表直接反射的太阳辐射项，用表示。而将定义为表观反射率，用表示。因此式（1）又可以改写为： （2） 总透过率包含了在太阳-地表-传感器路径上的总的水汽含量信息。为已知量。在1附近，可以忽略瑞利散射项，对程辐射的主要贡献是气溶胶的散射。同时在1附近，相对来说是非常小的一部分。由于大多数气溶胶和水汽均位于大气底层的约2km处，因此水汽也会对气溶胶的单次或多次散射造成影响。因此，也包含有水汽信息（Gao and Goeta, 1990）。这里我们假定在气溶胶光学厚度非常小的时候可以忽略。这样的假设可以让我们在反演气柱水汽含量时不必去考虑单次或多次散射效应。因此，在忽略的情况下，可以将式（2）改写为 （3） 图1（Gao and Goeta, 1990）显示了典型土壤的光谱反射率曲线，关于在1微米附近的地标反射率，我们从图1中可以看出在0.85-1.25微米之间，这几种地物基本上随着波长的增加而线性的增加。只有在富含铁的土壤和矿物中这种趋势不明显。由于各种地物的地表反射率不相同，想要从单一水汽吸收波段来确定水汽的透过率困难比较大。然而，如果地物的反射率在水汽吸收波段和大气窗口波段之间变化不明显，即假定地表反射率在这两个波段是恒定不变的，那么，我们可以利用这两个波段的比值来确定相应水汽吸收波段的透过率。如式（4）所示（Kaufman and Gao, 1992）： （4） 如果地表反射率随波长线性的变化，那么可以再加上一个大气窗口波段，利用3波段比值来确定水汽吸收波段的透过率，如式（5）所示（Kaufman and Gao, 1992）： （5） 其中，C1的值为0.8，C2的值为0.2。实际上，式（5）中的分母是在没有水汽的情况下利用0.85微米和1.24微米两个波段估计的0.94微米波段的反射率。 图 1 五种典型土壤反射率光谱曲线（Condit, 1970; Stoner and baumgartner, 1980）：（1）有机质含量占主导，纹理较细腻；（2）受有机质含量影响，纹理较粗糙；（3）铁含量占主导的红色土壤；（4，5）富含铁和有机质的土壤。 算法描述 2.1 建立近红外水汽反演查找表 利用MODTRAN模拟800-1300nm的大气透过率曲线，然后利用MODIS的0.865微米、0.905微米、0.936微米、0.940微米和1.24微米的光谱响应函数处理MODTRAN输出的大气透过率曲线以获得这几个MODIS波段的模拟透过率值，利用式（4）或式（5）所述的公式计算波段比值。在MODTRAN的输入参数中可以指定不同的水汽含量来获得不同水汽含量下的透过率曲线，也就可以计算响应水汽含量下的波段比值，这样就可以建立一个波段比值随着水汽含量变化而变化的查找表（查找表包括水汽含量、两波段比值、三波段比值）。水汽含量变化范围为：0~6.77cm, step = 0.01cm。下面是一个计算透过率曲线的MODTRAN输入参数的例子： TS 2 2 0 0 1 1 1 1 1 1 1 0 0 0.000 0.00 ff 2f 0 365.00000 g0.01 0 f t t f 0.000 DATA/BM