基于modis数据的广州市高分辨率气溶胶光学厚度时空变化

基于modis数据的广州市高分辨率气溶胶光学厚度时空变化

所谓总相桥，是指由总相桥和固体颗粒组成的多相体系。总相桥碎片的直径通常为0.00110m。根据研究，10m的颗粒可进入肺部，7m以下的颗粒可进入咽喉，但超过10m的颗粒可进入肺部，深入肺部，积累并添加到血液循环系统。这对人类的健康有直接影响，从而降低人体肺部的功能，以及呼吸系统疾病的加剧。此外，每年增加的死亡率与气和声的密度密切相关。总气音凝胶的吸收和散射影响太阳辐射，影响地气系统的辐射收支。总气中气凝胶在地球辐射收支中发挥着最重要的作用。另一方面，总气凝胶也可以用作云和雾的结合体，影响云的反射率，对地气系统的辐射平衡起到间接影响。然而，这种间接辐射强度的大小是不确定的。因此，总气位是气候系统研究最重要的不准确源。2004年的《东亚总气第六个实验》（亚欧性环境试验）试图研究中国地区的总气第六个实验。目前,已经有大量针对气溶胶光学特性的研究,美国NASA采用太阳分光光度计(CE-318)在全球建立了一个气溶胶自动地面观测网(AERONET);最近十几年,很多学者通过卫星遥感来研究气溶胶的特性.Levy等和Remer等在Kaufman的暗像元算法基础上开发了NASA的V5.2气溶胶反演算法,并利用MODIS传感器进行了气溶胶光学厚度的业务反演,为全球提供分辨率10km气溶胶产品.对于研究区域尺度和全球尺度的气溶胶污染物输送有重要应用价值.然而,对于城市地区的气溶胶来说,分辨率10km气溶胶产品的分辨率不能满足要求.高分辨率的卫星遥感反演,可以更好地认识城市气溶胶状况,李成才(2005)等利用MODIS资料遥感香港地区高分辨率气溶胶光学厚度,并与香港地区的地面污染物PM10质量浓度的变化进行了比较,结果表明气溶胶光学厚度产品可以用来描绘城市尺度的气溶胶污染分布.利用MODIS数据反演高分辨率的气溶胶光学厚度来研究城市气溶胶的时空动态具有可行性和应用价值.然而,这方面的研究较少.本文利用MODISL1B数据和NASA最新的V5.2气溶胶业务反演算法,对广州市进行高空间分辨率气溶胶光学厚度的反演,并分析广州市的气溶胶光学厚度时空变化特征.1数据收集1.1方面榆树MODISL1B数据从NASA的LADSWEB数据中心获得.该数据已经进行了辐射校正和几何校正.中分辨率成像光谱仪(MODIS)是EOS卫星TERRA和AQUA上搭载的一个重要传感器.上午星TERRA发射于1999年12月18日,大约当地时间上午10:30上行飞过赤道;下午星AQUA发射于2002年5月4日,大约当地时间下午1:30下行飞过赤道.提供了在可见光、近红外和红外共36个通道的全球观测,为反演有关陆地、云、气溶胶、水汽、臭氧、海色、浮游植物、生物地球化学等产品提供了丰富的信息.可见光通道1(660nm)和通道2(860nm)具有250m星下点的分辨率,可见和近红外的通道3~7具有500m的分辨率,其它的分辨率为1000km.MODIS扫描宽度2330km,覆盖全球只要1d.NASA利用MODIS数据进行气溶胶特性反演,发布了气溶胶光学厚度(AOT)和其他光学特性的10km的全球分布产品.1.2气溶胶光学特性参数地面数据利用法国CEMIL公司制造的自动跟踪扫描的太阳分光光度计(CE-318)进行观测.该仪器具有10个滤光片,中心波长为340,380,440,500,670,870,936,1020,1064nm,936nm是水汽的强吸收波段.它可以自动跟踪太阳进行太阳直接辐射测量、太阳等天顶角天空扫描、太阳主平面扫描和极化通道天空扫描.CE-318测得的太阳直接辐射数据可用来反演计算大气透过率、消光光学厚度、气溶胶光学厚度、大气水汽柱总量和臭氧总量.天空扫描数据可以反演大气气溶胶粒子尺度谱分布及气溶胶相函数.目前,利用多波段光度计遥感气溶胶光学厚度是目前气溶胶遥感手段中最准确的方法,通常被用来校验卫星遥感的结果.AERONET(AERosolRoboticNetwork)提供了大量的数据进行气溶胶光学特性的研究和卫星反演结果的检验.本次研究的观测仪器CE-318架设在中山大学生命科学学院实验楼顶,地理坐标为23.06667°N、113.39167°E,海拔高度大约为0.04km.位置如图1所示:2反演法是计算气溶胶的光学厚度2.1大气后向散射情况气溶胶遥感是通过观测到的大气顶辐亮度和表面双向反射率特性之间的关系来反演的,假设卫星观测的目标表面为均匀朗伯表面,不考虑气体吸收,卫星观测到的表观反射率(ρ*λ(θ0,θ,φ))可表示为:式中:θ为观测天顶角;θ0为太阳天顶角;φ为太阳光线的散射辐射的方位角;Fλ(θ0)为规一化地表反射率的下行辐射通量;Tλ(θ)为向上的总透过率;sλ为大气后向散射比;ρaλ(θ0,θ,φ)为程辐射反射率;ρsλ(θ0,θ,φ)为下垫面的反射率.在单次散射近似中,ρaλ(θ0,θ,φ)可用气溶胶单次散射相函数(Pa(θ0,θ,φ))和单次散射反照率(ω0)的关系来表示,则有:式中:ρmλ(θ0,θ,φ)为分子散射的程辐射;µ和µ0分别为观测(θ)和入射(θ0)的余弦cos(θ)和cos(θ0).卫星观测到的表观反射率ρ*λ(θ0,θ,φ)既是气溶胶光学厚度的函数,又是下垫面反射率的函数.如果知道了下垫面的反射率ρsλ(θ0,θ,φ),并根据不同地区的气溶胶特征,确定大气气溶胶模型,就可以得到气溶胶光学厚度;反之,如果已知气溶胶光学厚度和相应的大气参数,也可以反演下垫面反射率.Kaufman等研究发现,在清洁大气条件下,绿色植被表面的反射率在红、蓝和2.1µm通道相似,表明它们之间存在着一定的关系.由于可见光通道受气溶胶散射的影响,而2.1µm通道不受气溶胶的影响,反映了地面的特征,气溶胶下方的植被覆盖区地表的红、蓝通道的反射率可以从2.1µm通道的反射率估算出来.观测得到的红、兰通道的表观反射率和利用2.1µm通道估算出的地面反射率之间的差异就是气溶胶的影响.2.2日—CE-318气溶胶光学厚度的计算根据Beer-Bouguer-Lambert定律,在给定地点和一年中的给定时间测量的给定波长的太阳辐照度可表示为:求对数得:式中:I0(λ)为日地平均距离处大气上界太阳直接辐射辐照度;τtot(λ)大气总光学厚度;θ0为太阳天顶角.根据当地地方时的时角、地理经纬度、太阳赤纬可以计算出太阳高度角;(d0/d)2为日—地距离修正因子.采用近似公式:式中:J为一年中的第几天;m(θ0)为大气质量:大气总光学厚度[τtot(λ)]包括气溶胶光学厚度[τa(λ)]、瑞利分子散射光学厚度[τR(λ)],气体吸收光学厚度[τg(λ)].在不受水汽影响的波段大气总光学厚度表示为:瑞利分子散射光学厚度可用经验公式表示为式中:Z为海拔高度(km).气体吸收光学厚度可表示为:气体吸收光学厚度主要是臭氧、二氧化氮、二氧化硫和甲烷等各个分量的贡献的总和.利用CE-318观测结果对遥感反演结果进行验证时,由于波长不是一一对应的,所以运用埃斯屈朗方程对太阳分光光度计的结果进行插值.埃斯屈朗认为气溶胶光学厚度可表示为大气浊度参数β和波长指数α的函数,其中大气浊度参数β是波长等于1µm时的气溶胶光学厚度.埃斯屈朗公式可表示为:通过对CE-318观测的440,500,670,870nm处气溶胶光学厚度进行插值,生成与MODIS对应波段(550nm)的气溶胶光学厚度值进行比较.3结果3.1气溶胶光学厚度由于云对气溶胶光学厚度反演结果的精度影响很大,所以选取无云的MODISL1B数据进行反演.上午过境(10:30)TERRA卫星的MODIS数据32d,下午过境(13:30)AQUA卫星的MODIS数据38d.太阳分光光度计(CE-318)的气溶胶光学厚度用与卫星过境的时间对应的值,卫星遥感反演的气溶胶光学厚度取与地面观测点对应单个像元的值进行回归分析.由于(CE-318)没有和MODIS对应的550nm波长,通过埃斯屈朗方程进行插值,求得550nm波长处的气溶胶光学厚度.10:00~11:00的气溶胶光学厚度的平均值与TERRA卫星的MODIS反演的数据进行比较;13:00~14:00的气溶胶光学厚度的平均值与AQUA卫星的MODIS反演数据进行比较.图2对太阳分光光度计和卫星遥感反演的结果进行了回归分析,MODIS/TERRA和MODIS/AQUA与CE-318的回归方程分别为:τTERRA550=0.90*τCE-318550+0.14和τAQUA550=0.98*τCE-318550+0.14.相关系数分别为0.91和0.93,反演结果精度较高.与其他学者的比较研究结果相近.Levy等认为,在550nm处,MODIS/AERONET的回归方程近似于:TERRA/CE-318和AQUA/CE-318的气溶胶光学厚度的均方根误差(RMSE)分别为0.26和0.18.与李成才等利用MODIS资料遥感香港地区高分辨率气溶胶光学厚度的均方根误差(0.12)相比较大.然而,李成才等的研究认为:在植被密集的地表,均方根误差小于0.1,而在靠近海岸的地区均方根误差达0.3.所以,TERRACE-318和AQUA/CE-318均方根误差(RMSE)结果比较合理.用AQUA的遥感数据进行气溶胶光学厚度的反演精度略高于TERRA.TERRACE-318和AQUA/CE-318的回归分析表明,气溶胶较小时,MODIS反演结果高估了气溶胶光学厚度,反之亦然.这与其他研究结果相一致但是,MODIS的值普遍大于CE-318,说明本次研究MODIS反演的气溶胶光学厚度偏高,由于空间分辨率高,增大了气溶胶光学厚度反演的不确定性.3.2气溶胶光学厚度的日变化由于大气中的气溶胶变化较快,在不同的时刻气溶胶光学厚度的变化也比较大.图3和图4分别是2008~2009年冬季CE-318观测的气溶胶光学厚度和其相对应MODIS(TERRA、AQUA)反演的气溶胶光学厚度的变化,从2008年11月开始,气溶胶光学厚度逐渐增大;2008年12月的气溶胶光学厚度较低,平均大约为0.65,到2009年2月平均大约为1.35.总体上冬季的气溶胶光学厚度较大.TERRA的平均值为0.87,对应的CE-318为0.81,差值为0.06;AQUA的平均值为0.81;对应的CE-318为0.69,差值为0.12.遥感反演的气溶胶光学厚度略高于CE-318观测值.图5为广州市2008~2009年冬季平均气溶胶光学厚度的空间分布.从图5中可以看出气溶胶光学厚度空间差异显著,在0.1~1之间变化,呈现东北低西南高的特征.较小的气溶胶光学厚度分布在增城市、从化市.因为这些地区人为活动较小,且有大片森林覆盖,空气较为清洁,是气溶胶光学厚度较小的主要原因;广州市气溶胶光学厚度的空间分布下午观测的AQUA气溶胶光学厚度低于上午观测的TERRA.这和其他地区气溶胶光学厚度的日变化相似.原因可能是广州地区冬季早上大气水汽含量较大,出现了薄雾,但到中午时分薄雾逐渐散去,天空比较晴朗,大气状况稳定.所以,中午时分(13:00~14:00)的气溶胶光学厚度较小.通过计算图像的相似度来分析两个传感器反演结果在空间分布的差异,图像的相似度计算可转化为直方图的距离来计算,相似度的定量度量描述公式可表示为:式中:G,S分别为图像的直方图;N为像元样本数.gi,is为两个图像的像元.从图6和图7也可以看出气溶胶光学厚度的分布频率相似,呈现双峰分布,TERRA的峰值分别大约为0.37和0.72AQUA的峰值分别大约为0.29和0.61.两个反演结果的相似度为0.90,相似度较高,说明TERRA和AQUA的气溶胶光学厚度反演结果的空间分布很相似.反演结果的空间分布相似性也可以用相关系数来衡量,TERRA和AQUA平均气溶胶光学厚度的相关系数达到0.99.图8为气溶胶光学厚度的标准差分布TERRA卫星反演的气溶胶光学厚度的标准差集中在0.19~0.72之间;AQUA的标准差集中在0.05~0.35之间,AQUA的标准差整体低于TERRA也就是说气溶胶光学厚度较小的下午波动也较小;AQUA的标准差分布与平均值一致,呈东北低西南高;而TERRA却不然.这说明AQUA反演的气溶胶光学厚度在不同时间的空间分布相似,而TERRA在不同时间的空间分布变化较大.4气溶胶光学厚度反演4.1在广州市利用地面太阳光度计进行了气溶气光学厚度观测和利用MODISL1B数据进行了高空间分辨率(1km)气溶胶光学厚度的反演.通过利用地面多波段太阳光度计观测的气溶胶光学厚度对卫星反演结果进行验证,MODIS卫星遥感高空间分辨率气溶胶光学厚度达到了较高的精度,TERRA/CE-318的线性回归方程的斜率和截距分别为0.90和0.14,R为0.91;AQUA/CE-318的线性回归方程的斜率和截距分别为0.98和0.14,R为0.93;CE-318观测的10:00~11:00和13:00~14:00气溶胶光学厚度的平均值分别为0.81和0.69,对应的TERRA和AQUA反演的平均值分别为0.87和0.81.气溶胶较小时,MODIS反演结果高估了气溶胶光学厚度,反之亦然.遥感反演的气溶胶光学厚度略高于CE-318观测得到的气溶胶光学厚度;利用AQUA的遥感数据进行气溶胶光学厚度反演精度高于TERRA;AQUA的标准差整体低于TERRA,也就是说较小的气溶胶光学厚度的下午波动也较小.4.2利用MODIS卫星遥感反演高空间分辨率气溶胶光学厚度可以反映广州市气溶胶光学厚度的时空变化特征.2008年12月的气溶胶光学厚度较低,平均大约为0.65;随后气溶胶光学厚度逐渐增大,到2009年2月,气溶胶光学厚度平均大约为1.35.广州市气溶胶光学厚度空间差异显著,在0.1~1之间变化,呈东北低西南高的空间分布特征.即森林覆盖率比较高的地区气溶胶较低.本文利用NASA的V5.2气溶胶反演算法是Levy等和Remer等在Kaufman的暗像元算法的基础上发展起来的,是目前陆地上空气溶胶反演应用最为广泛的算法.该算法主要做了两个方面的改进:第一,可见光和中红外通道的反射率比值不再是常量,而是一个函数关系;第二,该方法考虑了植被指数(NDVI)的影响.进行了VIS/SWIR(可见光和中红外