一种定量化研究城市森林对PM2.5净化作用的方法技术

本发明专利技术属于森林管理技术领域，涉及一种基于遥感影像数据和气象站点及森林调查资料数据相融合的森林生态功能定量化研究方法。从城市区域尺度出发，依据MODIS遥感影像反演得到PM2.5平均浓度空间分布图，使用空间统计分析方法，定量分析各影响因子间的直接和交互作用大小。该方法主要由生成PM2.5平均浓度空间分布图、PM2.5平均浓度数据和各生态因子数据的空间匹配融合、选择合适分析阈值识别冷热点区域、锁定主要影响因子并改造4大阶段组成。本方法数据获取途径较简单且成本较低，反演操作简单具有推广性，而且对各影响因子影响力的定量化结果为城市森林管理者改善城市森林布局提供了指导思想和科学基础。

A quantitative study of the purifying effect of urban forest on PM2.5

The invention belongs to the forest management technology field, and relates to a quantitative research method of forest ecological function based on remote sensing image data and meteorological station and forest survey data. Based on the urban scale, the spatial distribution map of PM2.5 mean concentration is obtained based on MODIS remote sensing image inversion, and the direct and interaction size of each influencing factor is quantitatively analyzed by spatial statistical analysis. The method is mainly composed of spatial matching of PM2.5 average concentration, PM2.5 average concentration data and spatial matching of different ecological factors, selecting suitable threshold to identify cold hot spots, locking main influence factors, and reforming 4 stages. The method of data acquisition is simple and low cost. The inversion operation is simple and has generalization. Moreover, the quantitative results of influence factors provide guidance and scientific basis for urban forest managers to improve urban forest layout.

【技术实现步骤摘要】
一种定量化研究城市森林对PM2.5净化作用的方法
本专利技术属于森林管理
，尤其涉及一种基于遥感影像对地观测数据和气象站点及森林调查资料的地面观测数据相融合的森林生态功能定量化研究方法。
技术介绍
PM2.5颗粒物在地球大气中含量虽然很少，但富含大量有毒、有害物质或成为病毒和细菌的载体，在大气中的停留时间长、输送距离远，人体本身无法过滤和防护，对慢性疾病的发病率和死亡率影响很大，会引发包括哮喘、支气管炎、肺癌和心血管病等疾病。近年来，国内外对近地面颗粒物浓度的时空分布进行了广泛的监测、预测和反演等研究，主要可以分为2种，一种是依据机制过程模型对区域大气质量进行预测预报，另一种是基于遥感技术反演与PM2.5有物理机制联系的气溶胶光学厚度，然后建立其统计关系的方法。第一类方法可用于预测预报，但由于输入参数要求高，且较为复杂，如若简化则精准度受影响较大。第二种方法由于其成本较低，且具有时空连续性，可有助于研究PM2.5的治理工作。该方法主要有仅考虑AOD和PM2.5浓度之间关系的简单的线性模型、考虑多种环境气象因子的多元线性回归模型、人工智能模型和广义加法模型。简单线性模型对数据要求量小，易于操作，而多源线性回归模型则能大幅提高拟合精度，人工智能模型虽然对提高精度有一定作用，但其结构复杂构造困难，广义加法模型是线性模型非参数化的扩展，可用于预测模拟。对PM2.5的进行防治，除了控制源头污染，还要需要利用现有资源对其进行控制。而现有研究已表明城市森林对PM2.5及其他大气颗粒物有明显的调控作用，表现在不仅能覆盖地表减少大气颗粒物的来源，而且能通过叶面吸附和移除大气颗粒物，加速其沉降过程，通常有停留、附着和黏附3种方式。但是，目前开展的许多针对PM2.5的研究，主要是从样地或植株的小尺度，而且，对城市森林生态系统减少大气颗粒物的研究主要集中在总悬浮颗粒物(TSP)和PM10等大粒径粉尘上，对净化PM2.5的定量化研究成果还很少。
技术实现思路
针对现有的技术问题和研究领域的空白，本专利技术的目的在于提供一种可以从区域尺度、快速稳定地定量化分析一定时间段内城市森林对PM2.5的滞尘调节作用方法，并探明各影响因子贡献力大小，为城市森林优化布局提供理论支持和指导方针。为了实现上述目的，本专利技术采用如下技术方案：1.根据MODIS数据生成当天500m高分辨率气溶胶光学厚度空间分布图。2.通过高度订正和湿度订正后，得到干燥后气溶胶的消光系数（）空间分布图。3.建立与卫星过境时间地面PM2.5小时观测值的线性函数关系，从而得该城市PM2.5浓度平均值空间分布图。4.计算森林资源调查数据图形库中每块森林小班斑块上的平均PM2.5浓度。5.使用局部Gi*统计方法，在合适的分析阈值下，识别整个城市和森林斑块的冷点区和热点区。6.使用地理探测器模型分别对森林斑块冷点和热点区进行影响因子贡献度分析。7.通过各影响因子间的直接作用和交互作用，确定主要影响因素，并结合城市规划蓝图，在城市热点区对相应影响因子进行改造。进一步，所述步骤1）包括以下步骤：a)收集包含研究区的MODIS遥感影像数据、地面监测站点小时整点的PM2.5浓度、相对湿度、能见度、风速和风向数据。b)进行AOD浓度空间反演计算，空间辨别率达到500米。进一步，所述步骤2）包括以下步骤：a)使用地面气象站点为中心周围5公里范围内的非零AOD均值和气象站点卫星过境时间观测的能见度（VIS）数据，通过公式（1）和（2）估计当天的标高（H），并且在研究区域内取同样的标高值，进而将MODIS反演的AOD空间分布订正为地面消光系数（kex,0）的空间分布。高度订正依据公式如下：AOD=×H（1）（2）b)假设研究区域内取一样的相对湿度值，利用气象站点卫星过境时间观测的相对湿度对地面消光系数空间分布进行订正。当观测相对湿度小于40%时不进行订正，大于40%时进行订正，依据的公式如下：（3）kex,dry是干气溶胶消光系数，RH为观测的相对湿度。进一步，所述步骤3）包括以下步骤：a)选择气象站点上空属于无云遮盖状态时的地面观测PM2.5数据与其所在像元的遥感反演所得的AOD数据建立线性函数关系方程式。b)建立函数关系使用的数据根据卫星过境时的定时风向和风速进行分类后生成多个函数关系式，风向的分类包括北风NW、NNW、N、NNE；东风NE、ENE、E、ESE；南风SE、SSE、S、SSW和西风SW、WSW、W、WNW，风速的分类包括>1.5米/秒、1.5-7.0米/秒、大于7.0米/秒。c)利用所得的不同函数关系式，根据分类结果，对每日干气溶胶消光系数空间分布图进行修订，得到若干张PM2.5浓度空间分布图。d)将步骤c）所得分布图中每个相同空间位置的栅格值进行平均，得到研究区研究时间段PM2.5浓度平均值空间分布图。进一步，所述步骤4）具体为：将步骤3）中所获取的数据与森林资源调查数据图像库进行空间匹配，匹配的原则是以相同投影坐标系为准，然后依据森林资源调查数据库中已经划定的城市森林小班斑块，计算每个斑块PM2.5平均浓度。进一步，所述步骤6）包括以下步骤：a)收集各影响因子空间分布数据与森林资源调查数据图像库进行空间匹配，匹配的原则仍是以相同投影坐标系为准。b)依据森林资源调查数据库中已经划定的城市森林小班斑块，计算每个斑块各因子的平均值。c)将各斑块的连续变量数值按照行业标准进行划分，对于无行业标准连续变量，可采用自然间断点分类法（NaturalBreaks(Jenks)method）或等量分类法（Quantilemethod）将数据分成5-7类，对于已有的类型变量和PM2.5浓度数据则保持不变。d)将各个斑块分类后的各影响因子和PM2.5浓度数据导入地理探测器模型进行分析。附图说明图1示出了示范例1中干气溶胶消光系数和PM2.5浓度观测值回归分析结果。图2示出了范例1中不同阈值距离下冷热点区面积和显著斑块周围最小斑块数量变化趋势。具体实施方式：（一）以福建省泉州地区的晋江市为例。晋江市属于南亚热带季风气候，兼具海洋性气候特征，四季较为不明显，5-10月为夏季，年均气温在18.3–21.3°C。城市森林呈现减少和破碎化趋势，2014年斑块数约为6000个，森林总面积为9566公顷。其具体实施方式如下：1.收集2014年空间分辨率为250米包含晋江市的AQUA-MODIS数据，泉州市（清源山、万安、涂山街和津头埔）和厦门市（鼓浪屿、湖里、洪文、溪东）各4个地面监测站点逐小时的PM2.5浓度和晋江市南北2个气象站点数据，包含相对湿度、能见度、风速和风向。2.利用MODIS数据进行AOD反演计算，得到分辨率为500米的AOD空间分布图。3.高度订正和湿度订正后，通过AOD空间分布图得到干燥后气溶胶的消光系数（）空间分布图。4.由于晋江属于沿海城市有云天气较多，且气象站点较少，达到可建立关系原则的数据量较少，所以没有对风速和风向进行分类，而是将所有达到无云遮盖气象站点要求的当天数据与进行函数关系式的建立，得到简单线性方程y=46.75x+19.34，式中x为，单位为km-1，x为PM2.5浓度，单位为μg/m3。样本量为95，拟合R平方为0.50。依据该方程，得到晋江市2014年本文档来自技高网...

【技术保护点】
一种利用地面观测和对地观测数据，针对改善PM2.5浓度的城市森林生态功能定量化研究方法，其特征在于该方法包括以下步骤：1）使用MODIS遥感数据对城市区域某研究时间段内每天550nm波长下的气溶胶光学厚度（AOD）进行反演后，再分类建立

【技术特征摘要】
1.一种利用地面观测和对地观测数据，针对改善PM2.5浓度的城市森林生态功能定量化研究方法，其特征在于该方法包括以下步骤：1）使用MODIS遥感数据对城市区域某研究时间段内每天550nm波长下的气溶胶光学厚度（AOD）进行反演后，再分类建立与卫星过境时间观测的地面PM2.5小时观测值的线性函数关系，从而得到该城市PM2.5浓度平均值空间分布图，空间分辨率为500m；2）将遥感反演所得空间分布数据和地面调查数据融合，计算森林资源调查数据图形库中每块森林小班斑块上的平均PM2.5浓度；3）依据步骤2）的结果，使用局部Gi*统计方法，在合适的分析阈值下，识别整个城市和森林斑块的冷点区和热点区；4）使用地理探测器模型分别对森林斑块冷点和热点区进行影响因子贡献度分析；5）通过各影响因子间的直接作用和交互作用，确定主要影响因素，并结合城市规划蓝图，在城市热点区对相应影响因子进行改造。2.根据权利要求1所述的方法，其特征在于步骤1）使用的函数关系建立的假设前提应为气象站点上空属于无云遮盖状态。3.根据权利要求1所述的方法，其特征在于步骤1）建立函数关系使用的数据为卫星过境时气象站点观测的定时风向和风速，将其分类后可生成多个函数关系式，并应用于以气象站...

【专利技术属性】
技术研发人员：左舒翟，任引，石光明，
申请(专利权)人：中国科学院城市环境研究所，
类型：发明
国别省市：福建,35