近年来大气环境质量问题日益严重,乌鲁木齐市各地区空气质量差、能见度低等问题严重影响到了人们的生活。AOD受地表人类活动以及下垫面类型的影响而产生了复杂多样性以及不均匀等性质。通过大量数据研究气溶胶光学厚度的时空分布,有助于掌握AOD的分布规律,对大气颗粒物污染的监测和预测具有十分重要的作用。本文利用MODIS1 B级数据MOD021 KM日数据反演出2016年初乌鲁木齐市一次严重灰霾天气过程中的AOD时空动态情况,初步分析了这次天气过程中大气颗粒物的来源与去向。 In recent years, the problem of atmospheric environmental quality has become increasingly serious. The problems of poor air quality and low visibility in various regions of Urumqi have seriously affected people’s lives. AOD is subject to complex diversity and inhomogeneity due to the influence of surface human activities and the type of underlying surface. Studying the spatiotemporal distribution of aerosol optical thickness through a large amount of data helps to grasp the distri-bution law of AOD and plays an important role in monitoring and predicting atmospheric particu-late pollution. In this study, using MODIS1 B-level data MOD021 KM daily data to reflect the AOD spatiotemporal dynamics during a severe haze weather in Urumqi in early 2016, combined with meteorological data and using the backward airflow trajectory analysis model, the pollutants spread over Urumqi during this period. The process was analyzed, and the source and destination of atmospheric particulate matter during the weather process were analyzed.
张明远
江西理工大学,建筑与测绘工程学院,江西 赣州
收稿日期:2019年6月5日;录用日期:2019年6月20日;发布日期:2019年6月27日
近年来大气环境质量问题日益严重,乌鲁木齐市各地区空气质量差、能见度低等问题严重影响到了人们的生活。AOD受地表人类活动以及下垫面类型的影响而产生了复杂多样性以及不均匀等性质。通过大量数据研究气溶胶光学厚度的时空分布,有助于掌握AOD的分布规律,对大气颗粒物污染的监测和预测具有十分重要的作用。本文利用MODIS1 B级数据MOD021 KM日数据反演出2016年初乌鲁木齐市一次严重灰霾天气过程中的AOD时空动态情况,初步分析了这次天气过程中大气颗粒物的来源与去向。
关键词 :MODIS,遥感,气溶胶光学厚度
Copyright © 2019 by author and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
大气中悬浮的直径在0.001~100 μm的固体小颗粒、液滴以及有机物等众多物质共同构成的稳定混合物都被称作大气气溶胶。大气气溶胶由众多复杂物质构成,同时这些物质又来源于不同的地方,也是构成陆地-大气-海洋系统的一个非常重要的部分 [
中国区域的气溶胶时空分布的研究区域主要都是经济发达或者经济发展迅速的城市,对于远离海洋、经济发展相对来说较为缓慢的中国西北地区的气溶胶时空分布研究不多。国家对于西北地区经济发展一直都很重视,大力扶持西北地区经济发展的同时也很关注经济发展带来的各种环境问题。“十三五”规划提出深入实施西部大开发,并且在经济迅速发展的同时也要不断强化生态环境保护,通过培育城市群的发展带动西部经济发展。工业发展带动经济的发展,同时也带来了例如能见度下降、可吸入悬浮颗粒等许多环境问题,严重影响到了人们的生活及健康。因此,在发展经济的同时也必须注重环境的保护。在本次前期的资料及文献的相关阅读中,发现西北地区大气悬浮颗粒物监测站点数量较少,地面监测数据时间上不连续。由于气溶胶研究较少涉及到西北地区,因此有必要利用遥感技术研究西北地区AOD的时空分布特征,这有利于西北地区经济和环境的可持续发展。
国内外学者首次注意到气溶胶是在1999年,在亚洲南部区域的天空上出现了一片厚度大致为3 KM的棕色云 [
MODIS 1B级1KM地球数据产品:包含MODIS波段8至36的校准地球观测值观测结果,空间分辨率为1千米。它还包含来自MODIS频段1到7的校准数据,每个频段聚合到1 km分辨率。从MODIS双面扫描镜的一侧扫描的地球视图部分扫描出沿轨道10个像素的区域,通过1354个像素(或帧)交叉轨道。在最低点,这些像素是1 × 1公里。每帧包括来自36个频带中的每一个的10个像素(沿着轨道)。两个250米波段的每个像素具有16个子帧(数字样本),并且五个500米波段的每个像素具有4个子帧。此外,频段13和14输出高增益和低增益值。因此,每个帧由830个数字(DN)组成。每个值由带,检测器,子帧和镜像侧(BDSM)索引,并且每个值都是独立校准的。
遥感检测领域之中越来越多的运用到了AOD反演技术通过反演获得AOD产品对于研究十分的方便,气溶胶对辐射的散射及吸收作用会减弱辐射强度甚至是改变辐射的基本性质。对气溶胶进行反演其实就是测量辐射的变化情况。将地球表面假设为伯朗面,卫星观测的反射率为:
ρ ∗ = ρ α + ρ 1 − ρ ⋅ s ⋅ ( θ s ) ( θ v ) (2.1)
在此公式之中 ρ α 为大气中悬浮着的气溶胶在大气层中的大气分子的反射率,后面一项则用来表示气溶胶与大气和地表共同的反射率, ρ 为地表反射率,S为大气球面反照率,T( θ s ),T( θ v )上下层大气透过率, θ s 为太阳天顶角、 θ v 为卫星天顶角进行气溶胶反演之前一定要首先确定气溶胶的类型以及地表反射率数据在反演过程之中必不可少,本研究当中提到的反演方法主要就是暗像元法又被称为浓密植被法,其原理是植被密集区域的影响之下,传感器接收地表反射率会受到一定的影响特别是在是在可见光蓝(0.459~0.479 μm)、红(0.620~0.670 μm)两个波段,Kaufinan等 [
ρ 470 nm = 0.32 ρ 2100 nm (2.2)
蓝波与短波红外波长的系数为0.320,其关系式如下
ρ 660 nm = 0.55 ρ 2100 nm (2.3)
地面红波段和蓝波段的反射率与短波红外的反射率的关系分别为
ρ B = ρ 2.1 4
ρ R = ρ 2.1 2
通过红蓝波段与短波红外反射率的关系可以推测出气溶胶类型。基于大陆气溶胶查表(LUTs)反演,以红波与蓝波通道系数的比值为基础,确定气溶胶类型是否为粉尘;non-sand尘埃的类型、地理分布的方法被用来确定硫酸气溶胶类型,生物燃烧或其他混合类型,然后是AOD的红色和蓝色波段气溶胶模型修正达到陆地气溶胶光学厚度的反演。
通过红蓝波段与短波红外反射率的关系可以推测出气溶胶类型。基于大陆气溶胶查表(LUTs)反演,以红波与蓝波通道系数的比值为基础,确定气溶胶类型是否为粉尘;non-sand尘埃的类型、地理分布的方法被用来确定硫酸气溶胶类型,生物燃烧或其他混合类型,然后是AOD的红色和蓝色波段气溶胶模型修正达到陆地气溶胶光学厚度的反演。
本章反演使用到的方法为暗像元法也叫暗目标法。气溶胶的影响是气候模拟中最大的不确定性来源之一。气溶胶在空间上的时间不同,可能导致云微物理的变化,从而影响云的辐射特性和气候。将暗目标(DT)气溶胶检索算法应用于多光谱卫星传感器数据,并推导出气溶胶特性,包括陆地和海洋上的气溶胶光学厚度(AOD),以及海洋上的光谱AOD和气溶胶尺寸参数。DT检索产品用于开发全球和区域气溶胶气候学,研究气溶胶与云的相互作用,以及空气质量评估和预报。
DT算法已经应用于MODIS传感器的数据,从1999年的Terra卫星开始,然后是2001年的Aqua。实际上有两个独立且不同的“暗目标” (DT)算法。第一个用于检索海洋上的气溶胶信息(可见和较长波长的黑暗),第二个用于植被/黑暗污染的土地。自Terra推出以来,这些算法不断完善和改进。更新的算法定期应用于原始数据以生成新的数据产品集合。
尽管最初开发用于MODIS传感器,但这些算法可应用于任何测量反射率的传感器,该传感器测量覆盖可见光(VIS),近红外(NIR)和短波红外(SWIR)的适当波段。目前,DT检索算法正被移植到其他传感器上,包括SUOMI-NPP上的VIIRS,AHI,EPIC以及通常在ER2上飞行的机载传感器eMAS,并将在不久的将来应用于ABI (GOES-R-16)。
利用MODIS L1B数据反演气溶胶的流程如图1。
1) MODIS数据使用MODISL1B 1KM (MOD021KM)数据。操作应严格按照说明书进行,需要注意的地方主要是数据分辨率为1 KM,数据合成时反射率在上,发射率在下部。
2) 在对于乌鲁木齐市的气溶胶反演过程当中,主要就是使用ENVI5.1对MOD02数据进行辐射校正、几何校正、角度校正后将经过辐射校正和反射校正后的影像数据进行合成,将角度数据校正完合成之后还要经过BAND MATH工具处理,最后将经过云处理插件处理的影像数据、查找表和角度数据载入气溶胶反演工具当中,就可得到反演结果。
3) mod02气溶胶数据反演过程当中使用到的反演算法是经典的暗像元法也被称为浓密植被法。
4) 气溶胶反演工具中需要的查询表是从网上找到的适用于处理冬季影像数据其主要原理就是通过IDL调用6S大气辐射传输模型。另外还用到了两个ENVI插件,modis Aersol inversion和modis cloud。
2016年乌鲁木齐市空气质量问题还很严重,虽然不断的施行环境保护方面的措施,空气质量情况良好的时间占全年仅有百分之六十左右。全年大气污染主要由颗粒物引起。PM2.5是污染的主要原因。冬季空气质量污染最为严重,且优良率仅为4.88%。这主要是因为与冬季其他季节相比,乌鲁木齐的气候及地理特征都不利于污染物的扩散。整个城市呈现出由北向南递增的格局。乌鲁木齐市山谷风盛行白天风速明显大于夜间,因此,水塔山可以阻断污染物向南方的传播和扩散。市中心很高,周围很低,除自然来源外主要是由于燃烧燃料和大气颗粒物污染物来自城市汽车尾气。
反演结果如图2所示。乌鲁木齐AOD在城市高于农村整体分布,高值区主要集中在乌鲁木齐市天山区,沙区,水磨沟区,头屯河,米东区的市中心位置,达坂城区,乌鲁木齐市郊区位置低AOD区域,从郊区到市中心AOD的增长加快。AOD的高值区域是局部聚集现象,主要集中在天山区和沙依巴克区的市区。大气悬浮的各种污染物不易扩散,在此地区AOD的不断积累促进了城市AOD的增长。
图1. 反演流程图
图2. 反演结果
乌鲁木齐地处中部天山北麓,准噶尔盆地南缘。乌鲁木齐是新疆维吾尔族自治区首府,是自治区的经济、政治、科技、金融、文化中心和交通枢纽。乌鲁木齐三面环山。乌鲁木齐北部与准噶尔盆地相连,准噶尔盆地大致是向北开放的敞开状。乌鲁木齐市地形地貌复杂,山谷气候复杂,经常受到东南风的袭击。乌鲁木齐特殊的地形和气象条件不利于大气污染物的扩散。此外,乌鲁木齐作为工业城市,已基本形成以煤炭、钢铁、陶瓷、化工、建材为主导的产业经济体系,能源结构以燃煤为主,环境压力较大。
据历史气象资料统计可知,2016年1月23日至2016年2月9日乌鲁木齐市出现了一次严重污染过程如图3所示,重污染过程期间AQI指数均超过200,1月的28、29及2月的8、9号,出现五次AQI指数超过300为严重污染。
图3. AQI变化
由表1可以知道在这次严重的天气污染过程AQI值与6种大气污染颗粒物的相关系数,分别为0.960 (PM10),0.999 (PM2.5)、−0.194 (SO2)、0.096 (NO2)、0.246 (CO)和0.057 (O3),其值从大到小排列为PM2.5 > PM10 > NO2 > CO > O3 > SO2。在此次污染过程中,PM2.5与PM10和AQI的相关性很高,尤其是AQI与PM2.5相关系数达到0.9,由此表可看出在乌鲁木齐的这次天气过程当中空气中的可吸入颗粒物PM2.5和PM10是主要的污染物。
PM2.5 | PM10 | SO2 | NO2 | CO | O3_8 h | |
---|---|---|---|---|---|---|
相关系数 | 0.9 | 0.96 | −0.19 | 0.69 | 0.25 | 0.06 |
表1. 与AQI的相关系数
由表2可以看出,在本次污染过程中,AQI,PM10,PM2.5,O3和气温显著相关,相关系数从大到小排列为r(O3) > r(AQI) > r( PM2.5) > r(PM10)。NO2和风速负相关,风速较大有利于NO2的扩算。
相关系数 | 平均气温 | 平均风速 | 平均相对温度 | 降水量 |
---|---|---|---|---|
AQI | 0.609 | 0.110 | −0.276 | 0.151 |
PM10 | 0.521 | 0.139 | −0.357 | 0.182 |
PM2.5 | 0.601 | 0.099 | −0.278 | 0.151 |
SO2 | 0.156 | 0.207 | 0.060 | 0.100 |
NO2 | −0.339 | −0.611 | −0.124 | −0.435 |
CO | 0.075 | −0.384 | −0.453 | −0.563 |
O3 | 0.778 | 0.465 | −0.276 | −0.021 |
表2. 相关性分析
Hysplit-4模型是NOAA空气资源实验室和澳大利亚在过去20当中年联合开发的专业模型,主要是通过大气运动情况来预测出大气污染物的扩散情况。该模型具有较为完整的运动、扩散和沉积模型,能够处理多个气象输入场、多个物理过程和不同类型的污染物排放源。Hysplit-4气流轨迹分析模型分为前向和后向两种分析方式可以用来前推或者后推出气流轨迹的变化情况。在线地址为 http://ready.arl.noaa.gov/HYSPLIT_traj.php。使用到的数据集是NCEP GDAS把全球1˚ × 1˚数据插值到正形投影的地图上需要注意的是使用的时间为UTC时间,该数据每月都会更新,2005年起每七天数据为一组。
图4. 后向气流轨迹分析
选择2016年1月28日上午10点对应UTC为2点,地点设置为新疆维吾尔自治区省会乌鲁木齐,向后回推120 h的模拟轨迹图。在图中分别有红绿蓝三条轨迹线路,其中红线为500米高度轨迹,蓝线为1500米高度轨迹,绿线为2000米高度轨迹。结果如图4所示,乌鲁木齐地区近地面风场受地形影响严重,在乌鲁木齐附近可能发生山谷风,地面高度1500 m以上,地形对气象影响很微弱;由后向气流轨迹模型可知,乌鲁木齐主要受三种气流控制,即吐鲁番盆地的东南气流、来自准噶尔盆地的西北气流以及从天山西南附近的气流影响,大气污染物会在乌鲁木齐滞留。这一现象的主要原因是白天,天山斜率接收太阳辐射、空气温度升高,热空气在山坡上上升,沿着山坡寒冷的空气上升到山顶,从顶部流向准格尔盆地和吐鲁番盆地。该盆地抑制了盆地内强上升气流的形成,而乌鲁木齐南部吐鲁番盆地和西北部准噶尔盆地则以弱上升或弱下沉气流为主,天山、乌鲁木齐山谷和盆地均是如此。它们之间形成热循环,直接导致上层难以进入,部利于大气污染物的扩散。
据气象资料显示,1月24日至28日夜间,天空云层较多,乌鲁木齐附近最高风速仅为1.3 m/s,风力很弱影响大气污染物的扩散。
2016年乌鲁木齐市一次严重灰霾天气情况过程表明,乌鲁木齐市AOD整体分布趋势为市区高于郊区,乌鲁木齐市经济的发展以及大量重工业产生的污染对该市AOD增长起促进作用,另外其地形特征会阻碍大气污染物的扩散容易使得污染物产生堆积的情况出现。与其他季节相比,在冬季乌鲁木齐大气中会出现逆温层,逆温层较厚且维持时间较长也会阻碍大气中污染物的扩散。整个城市呈现南北递增的态势。主要原因是乌鲁木齐市城区河谷风盛行,白天河谷风速明显高于夜间的山风,因此水塔山在一定程度上可以阻断污染物向南方的传播和扩散。
张明远. 乌鲁木齐市一次灰霾天气过程中的AOD时空动态分析Spatial and Temporal Distribution of Aerosol Optical Thickness in Urumqi and Its Influencing Factors[J]. 环境保护前沿, 2019, 09(03): 486-495. https://doi.org/10.12677/AEP.2019.93068
https://doi.org/10.1029/95GL01660
https://doi.org/10.1126/science.173.3992.138
https://doi.org/10.1175/1520-0442(1995)008<0773:EOAPOT>2.0.CO;2
https://doi.org/10.1126/science.287.5459.1763
https://doi.org/10.1109/36.701075
https://doi.org/10.1029/2005GL024724
https://doi.org/10.1029/2009GL041828
https://doi.org/10.1080/01431161.2012.700135
https://doi.org/10.1016/j.atmosenv.2013.08.024
https://doi.org/10.5194/acp-15-11653-2015
https://doi.org/10.1007/s12145-016-0277-z