资源描述：
《modis卫星遥感北京地区气溶胶光学厚度及与地面光度计遥感的对比_毛节泰》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、.第13卷特刊应用气象学报Vol.13,Suppl.2002年1月JOURNALOFAPPLIEDMETEOROLOGICALSCIENCEJanuary2002MODIS卫星遥感北京地区气溶胶光学厚度及与地面光度计遥感的对比毛节泰1)李成才1)2)张军华1)3)刘晓阳1)刘启汉2)1)(北京大学物理学院大气科学系,北京100871)2)(香港科技大学海岸与大气研究中心)3)(DepartmentofPhysics,DalhousieUniversity,Canada)提要介绍MODIS卫星遥感气溶胶的方法,利用北京大学地面多波段太阳光度计的观测进行了对比,二者的相关性比较好。给出了描述北京地
2、区气溶胶光学厚度分布的几幅图片。卫星遥感对于更好地研究空气污染提供了一种新手段,卫星遥感的气溶胶光学厚度弥补了地面观测空间覆盖不足的缺陷。卫星遥感的气溶胶资料不仅对全球和区域气候研究而且对城市污染分析提供了丰富的资料。关键词:卫星遥感气溶胶光学厚度地面光度计引言对流层气溶胶粒子来源于工业活动、生物燃烧等人为源以及沙尘、海洋粒子等自然源。它们对太阳辐射的直接作用以及对云的微物理过程的影响、对大气化学过程的影响[1~5]分别表现为对气候变化的直接强迫和间接强迫。由于大多气溶胶只有几天的生命期[6,7]以及空间变化很大,利用卫星来观测全球气溶胶的分布和变化非常迫切。NASA发射的Terra卫星承载的
3、MODIS传感器具有36个通道,覆盖了紫外、可[8,9]见、近红外、红外等通道,为反演气溶胶和地表特征提供了丰富的信息。卫星接收到的辐射来源于地球大气的散射以及地表反射的复杂相互影响,卫星遥感陆地气溶胶长期..128应用气象学报13卷以来存在地表反射率和气溶胶光学特征两方面同时未知的难题。Kaufman等通过大量飞机试验发现对于植被密集的具有较低反照率的地表存在2.13m近红外通道反照率与0.47m、0.66m可见光通道反射率相关很好的结论,成功地运用于MODIS的气溶胶反演[10,11]。城市气溶胶光学厚度同时反映了大气污染的污浊程度,高分辨率的卫星遥感提供了监测城市大气污染的可能性。卫星遥
4、感弥补了一般地面观测难以反映污染物空间具体分布和变化趋向的不足。但是对于光学厚度反演过程中源于地表反照率和气溶胶模型带来的误差难以估计,卫星遥感需要同时有地面太阳光度计观测进行对比。NASA在全球几个洲的大陆部分和少数海岛海岸建立了一个包含大约70个地面太阳光度计的全球自动观[12,13]测网,为MODIS遥感气溶胶提供地面多通道遥感的对比资料。我们利用多通道太阳光度计在北京的连续观测,对照了MODIS气溶胶产品,二者的相关性比较好。对于中国大陆,各地区地表植被分布、气溶胶类型具有很大差异,有必要在不同代表性地区建立这种连续观测的地面站。1MODIS陆地气溶胶遥感[14,15]首先简单回顾一下
5、利用卫星遥感大气气溶胶的基本原理。卫星观测到的反射sat率R可以由下式表示:satSURFsatR=f(0,R)=L(0,,,0,0,)/0F0)(1)其中,F0为大气上界太阳辐射通量密度,0为整层大气光学厚度,(0,0)为入射光的方向,(,)是卫星观测方向,0,0,,分sat别为天顶角的余弦和方位角,L(,,,0,0)为卫星接收到的辐射,它满足辐射传输方程。SURFsat如果假设下垫面是反射率为R的均匀朗伯面,R可表示为:ATMATMSURFRsat(,,0,0)=RATM(,,0,0)+T(0)TSURF(ATM)R特刊毛节泰等:MODIS卫星遥感北京地区气溶胶光学厚度及与地面光度计遥感的
6、对比129ATMATMATMR为整层大气反射率,T(0)和T()分别表示从太阳到地面、从地面到SURFATM卫星大气层总的透过率(直射+漫射),系数1/(1-RR)代表地面和大气层多次sat散射的作用,因此卫星观测到的反射率R既是气溶胶光学厚度的函数,又是下垫面反SURF射率的函数,如果知道了下垫面的反射率R,并假定一定大气气溶胶模型,可以反演得到气溶胶光学厚度。辐射传输过程一般利用辐射传输模式来完成。业务的反演算法为节省计算时间采用一个含有多维变量的查算表(LUT)。NASA的MODIS气溶胶算法中,首先为了获得较高的信噪比,将250m和500m的晴空卫星资料合成为1km分辨率的资料,以下列
7、4个优先级进行蓝光0.47m和红光[10]0.66m通道地表反射率的确定:1)如果0.01R*2.10.05,则R0.47=R*2.1/4,R0.66=R*2.1/2;*2)如果R3.80.025,则R0.47=0.01,R0.66=0.02;3)如果0.01R*2.10.10,则R0.47=R*2.1/4,R0.66=R*2.1/2;4)如果0.01R*2.10.15,则R0.47=R*2.1/
