论文部分内容阅读
痕量气体是大气重要的组成成分,它们受各种物理的、化学的、生物的、地球过程的作用参与生物地球化学的循环,对地球生命演化、大气化学发展及全球环境变化具有重要的影响作用。紫外辐射能破坏对生命至关重要的有机高分子,而臭氧可以吸收太阳紫外辐射,从而使生物器官免于其害;陆地生命的演化也被认为与臭氧密切相关。臭氧的光化学循环对维持大气温度垂直结构起着重要的作用,而且臭氧的减少还能指示人类对环境的破坏。NO2能对平流层臭氧光化学反应起到催化作用,还与光化学烟雾、酸雨、温室效应等密切相关。同时,由于近年工业的快速发展,大气中NO2浓度增长迅速,已经成为空气污染的一个重要指标。自上世纪80年代以来,南极平流层臭氧空洞的发现和大城市光化学烟雾污染的日益加剧,推进了大气痕量气体和区域环境污染问题研究,因此对大气痕量气体,尤其是臭氧和NO2的时空分布研究成为当前大气化学领域的一个前沿和热点问题。痕量气体的探测,经历了地基、探空气球、火箭及当前卫星等平台的发展,相对地基和探空仪等传统平台,利用卫星能对痕量气体进行一个全球范围、高空间分辨率地观测。传统的卫星大气探测几何方式,包括垂直观测地表的天底观测模式和指向太阳、月亮或星星的掩星观测模式。天底模式探测虽然能得到高空间分辨率的全球痕量气体总量分布,但是垂直分辨率非常有限;而掩星模式虽然能得到高垂直分辨率的痕量气体廓线,但是空间覆盖又比较差。自1981年美国国家航空航天局(National Aeronautics and Space Administration, NASA)第一次将临边扫描技术应用于太阳中间层探测卫星(Solar Mesosphere Explorer, SME)大气观测以来,这种新的卫星探测技术,即观测地表以上一系列切高上的大气临边辐射,开始在卫星大气遥感领域广泛应用。这种临边探测技术结合了天底和掩星两种技术的优点:提供如天底模式一样高空间分辨率全球覆盖的大气痕量气体分布,且能与掩星模式相媲美的高垂直分辨率的气体廓线。搭载在2002年3月1日发射升空的欧空局(EuropeanSpace Agency, ESA)环境卫星ENVISAT (ENVIronment SA Tellite)上的大气制图扫描成像吸收光谱仪(Scanning Imaging Absorption Spectrometer for Atmospheric CHartography, SCIAMACHY)仪器即具有临边探测能力,SCIAMACHY观测位于飞行方向前方3280km处地表至92km切高范围的临边辐射,垂直方向切高增加步长约为3.3km,工作波长为214-2386nm,光谱分辨率为0.24-1.48nm。其他具有临边探测能力的仪器有搭载在Odin卫星上的光学光谱仪与红外成像系统(Optical Spectrograph and InfraRed Imaging System, OSIRIS),于2001年2月20日发射升空,实际包括光学光谱仪(Optical Spectrograph, OS)和红外成像仪(InfraRed Imaging system, IRI)两个独立部件,分别探测7-70km范围的临边紫外-可见光散射辐射和7-100km范围的氧红外波段气辉辐射,视线切高步长约为2km。目前针对临边散射数据,特别是SCIAMACHY和OSIRIS数据的痕量气体反演研究,一般均是用波长配对或差分光学吸收光谱法(Differential Optical Absorption Spectroscopy, DOAS)应用于辐射处理,结合最优估计反演痕量气体数密度廓线,并且反演高度只限于平流层范围。本文将一种新型反演策略,即波长配对与加权乘代数重建技术(Weighted Multiplicative Algebraic Reconstruction Technique, WMART)联合SCIATRAN (radiative TRANsfer model for SCIAMACHY)正向模型,应用于SCIAMACHY临边散射辐射数据,以臭氧和NO2为例研究了大气痕量气体垂直廓线的反演。首先将临边辐射廓线正规化到参考切高,然后将吸收波长的正规化辐射廓线与参考波长的正规化辐射廓线组合生成反演向量,最后利用WMART结合SCIATRAN进行痕量气体反演,并开发了能批量处理SCIAMACHY多轨数据获取臭氧和NO2廓线的程序包SCIA_JLU,臭氧反演每轨需时约30分钟、NO2约15分钟。经过波长配对处理后,使得云、地表反照率和气溶胶对辐射值的影响降低一个数量级,而保持或提高了待反演成分的敏感性。根据临边辐射廓线、臭氧和NO2的辐射权函数及其反演向量权函数特征,选定臭氧Hartley和Chappuis吸收波段为臭氧反演波段,其中Hartley波段将用吸收较强的267.5 nm、286.5 mm、287.5 nm及305.1nm辐射分别与吸收较弱的307.5nm波长配对,Chappuis波段将吸收峰599nm辐射与两翼525nm和668nm辐射进行三波长组合,总体反演高度范围为10-68 km; NO2反演波段为420-450nm,并选定其中3组共6个波长进行组合,每一个组由吸收峰与两翼吸收较弱的3个波长组成,反演高度均为15-40km。将WMART与SCIATRAN辐射传输模型结合,确定了WMART不同切高上权因子和不同波段权因子,综合利用臭氧Hartley和Chappuis吸收波段SCIAMACHY临边散射辐射数据,反演了平流层到中间层低层大气臭氧数密度垂直廓线;并利用420-450nm波长范围多波段SCIAMACHY临边散射数据反演了大气平流层NO2数密度垂直廓线。为了分析加权乘代数重建反演误差来源及大小,对反演过程中一些不确定的参量进行了敏感性分析。通过正演模拟时引入参量的偏差,反演得到“不正确”的结果与没有引入偏差的“正确”反演结果进行对比,分析参量对反演结果的影响。研究中对切高定位、气溶胶、地表反照率和云这些主要误差源进行了误差分析,其中气溶胶参数包括影响较大的边界层能见度、平流层气溶胶负载和气溶胶消光系数,云参数包括敏感性较大的云高和光学厚度。结果表明,在参数可达到的最优精度条件下,整体误差最大的是切高定位偏差引起的、其次是气溶胶偏差引起误差、最小的误差是由地表反照率和云偏差导致的。然而,在平流层,不正确气溶胶引起的误差比切高定位偏差引起的误差更大,这主要由于平流层气溶胶负载和气溶胶消光系数两个参量导致,特别是平流层存在火山气溶胶时,会导致臭氧在平流层底层误差达42%左右,甚至会使NO2反演误差在15 km超过100%。臭氧和NO2反演结果受参量的影响基本随高度升高而减小,地表反照率和云的影响明显依赖太阳天顶角,而其他参量对反演结果的影响在不同天顶角下区别不大。相比之下,NO2反演结果受参量影响比臭氧更大一些。为了评价利用波长配对结合WMART算法反演的痕量气体廓线结果,将臭氧反演结果分别与Bremen大学提供的SCIAMACHY臭氧V2.3版本(称为BU臭氧)3天同步反演点数据、纬差±5°以内和经差±10°以内共3天Saskatchewan大学提供的OSIRIS臭氧V3.0版本和微波临边光谱仪(Microwave Limb Spectrometer, MLS)臭氧进行比较验证,同时将NO2反演结果分别与Bremen大学提供的SCIAMACHY NO2 V3.1版本(BU NO2)和Saskatchewan大学提供的OSIRIS NO2 V3.0版本进行了比较,NO2比较的时间及空间匹配条件与臭氧相同。发现臭氧反演结果与BU臭氧一致性最好,误差不超过15%,与OSIRIS对比最大偏差为20%,而与MLS比较仅在20-46km之间误差小于10%。对于N02,与BU NO2相对偏差小于10%,与OSIRIS NO2最大偏差在16%左右,且大部分高度范围内偏差均在10%以内。基于OSIRIS临边气辉辐射数据,尝试了中间层臭氧的反演。从层析成像反演技术在卫星临边探测中的应用出发,经过改进得到了用于卫星临边二维剖面反演的有效快速算法,利用IRI气辉辐射亮度数据反演了二维气辉体发射率垂直剖面。气辉辐射本质是激发态大气O2(a1Δg)回到基态发出光子的现象,衡量气辉辐射大小的体发射率与O2(a1Δg)浓度成正比。基于奇氧族光化学反应循环,建立了以大气连续方程为机理的O2(a1Δg)光化学模型,从而将臭氧和O2(a1Δg)定量地联系起来。通过由臭氧到体发射的正向模拟,及利用气辉基于层析成像反演的体发射相比,验证了光化学模型的正确性。基于光化学模型,最终分别结合剥洋葱和牛顿迭代法反演了50-90km高度范围的中间层臭氧二维垂直剖面。从反演结果与SCIAMACHY臭氧V2.3版本数据比较中发现,利用牛顿迭代法比剥洋葱法精度更高,平均偏差在12%以内,除个别高度外偏差基本在10%以内。臭氧剖面在高度上呈指数分布,但可惜地是,在80km附近臭氧第二峰值结构并没有明显地识别出来。最后受气辉辐射层析成像反演的启发,提出了一种新的临边散射探测技术与反演方案,即结合DOAS和层析成像反演技术,可同时得到沿轨道和垂直两个方向上的二维垂直剖面。根据新技术的实现要求,反演首先用DOAS技术对临边散射数据进行光谱拟合,得到痕量气体沿视线的柱量,然后以基于加权乘代数重建技术算法的层析成像技术反演大气网格中痕量气体数密度。并以N02为例进行了反演试验,结果表明,反演剖面结构和假设剖面在结构和数值大小两个方面一致性较高,剖面中的典型位置廓线差异在15%以内,且40km以下小于5%。通过该研究,应用波长配对与加权乘代数重建技术联合SCIATRAN这种新的反演策略从临边散射辐射中反演高垂直分辨率的痕量气体廓线,从精度、速度以及高度范围扩展等方面是可行和可靠的;利用临边气辉辐射数据通过层析成像技术能反演得到中间层臭氧剖面,且提出的DOAS-层析成像临边散射反演技术用于痕量气体二维垂直剖面反演是可行的。