基于Kolmogorov-Zurbenko滤波法分析2015—2021年京津冀地区大气颗粒物变化趋势.pdf

生态与农村环境学报，（）：收稿日期：基金项目：江苏省自然科学基金（）；国家环境保护城市大气复合污染成因与防治重点实验室开放基金（）通信作者：占青，张运江，陈红，等基于 滤波法分析 年京津冀地区大气颗粒物变化趋势生态与农村环境学报，（）：，（）：基于 滤波法分析 年京津冀地区大气颗粒物变化趋势占 青，张运江，陈 红，张可馨，盖鑫磊（南京信息工程大学环境科学与工程学院 江苏省大气环境监测与污染控制高技术研究重点实验室 大气环境与装备技术协同创新中心，江苏 南京）摘要：京津冀地区是我国近年针对大气颗粒物污染治理的重点区域。采用（）滤波分析 年京津冀地区 个城市细颗粒物（）、可吸入颗粒物（）和粗颗粒物（）的时间序列。结果表明，和 年均浓度总体呈下降趋势，其下降速率分别为 和 。滤波分析结果表明，、和 的短期分量、季节分量和长期分量之间相互独立。短期分量贡献占比（大于）最高，其次为季节分量，说明原始时间序列的波动主要受控于污染源排放以及气象条件的短期变化和季节变化。针对 滤波分解出的浓度序列，建立逐步多元回归线性模型，定量评估气象和减排对变化趋势的贡献。逐城分析结果表明，秦皇岛市气象条件对 浓度的影响（）最大，衡水市减排措施对 浓度的影响（）最大；秦皇岛市气象条件对 浓度的影响（）最大，廊坊市减排措施对 浓度的影响（）最大；沧州市气象条件对 浓度的影响（）最大，廊坊市减排措施对 浓度的影响（）最大。整体分析结果表明，减排对 年京津冀大气颗粒物污染改善起着至关重要的作用。关键词：大气颗粒物；滤波法；逐步多元线性回归；变化趋势；京津冀地区中图分类号：；文献标志码：文章编号：（），（，）：，（）（），（），（），（），（），生 态 与 农 村 环 境 学 报第 卷 （）（），（）（），（），：；大气颗粒物是影响人体健康和环境质量的主要污染物，按照颗粒物粒径进行划分，可分为细颗粒物（）、可吸入颗粒物（）和粗颗粒物（）。我国城市规模扩大、工业化迅速扩展以及长期积累的能源结构矛盾，使得区域大气状况日益恶化。为了统筹推进区域空气环境治理工作，年国家开始实施大气污染防治行动计划打赢蓝天保卫战三年行动计划。北京、天津、河北三地在积极响应国家号召的同时也在积极探讨区域环境污染的协调机制，以降低环境中有机物排放量和改善三地空气环境质量；和 年连续发布京津冀及周边地区秋冬季大气污染综合治理攻坚行动方案；年出台天津市大气污染防治条例。京津冀作为国家经济发展中心，同时也是连接中原和华北平原的主要通道，通过分析京津冀大气环境质量和污染物空间分布特征，对实现区域空气质量联合控制起到至关重要的作用。自 年以来，我国各大城市空气环境监测中心陆续公布各大城市大气污染指数，并发表了一批关于京津冀、长三角、珠三角、成渝和汾渭平原等区域的研究成果。基于大气化学传输模型和污染源排放清单数据，近年我国学者研究发现减排行动是主导 年京津冀地区大气颗粒物污染改善的重要因素，贡献达。采用大气化学传输模式评估排放对大气颗粒物污染贡献的优点之一是可以定量掌握格点化的大气污染物排放总量与空气质量变化的直接响应关系。然而，基于大气化学传输模式的方法需要最新的排放清单数据作为模型输入项。因此，近年来很多学者选用基于观测数据的统计学方法来评估气象条件和排放对大气污染的贡献。基于观测数据的统计学方法可以直接采用观测数据进行评估，而不依赖于排放数据，有助于快速掌握排放、气象条件对污染物变化趋势的贡献，目前该方法已被国内外学者广泛使用。例如，等基于逐步多元线性回归模型（）量化气象条件对中国各地 的贡献趋势。在国外，该方法还在美国亚利桑那州图森地区和韩国得到应用，研究这些地区 和颗粒物的污染趋势、滤波分量与气象参数之间的相关性，以及耦合其他方法开展 污染预测研究。但这些研究基本集中在 年及之前，而近年来减排对大气颗粒物的影响尚不清楚，对 尤其是 变化趋势的了解也较少。虽然目前京津冀地区大气环境质量得到一定改善，但区域颗粒物污染状况依然不容乐观。因此，以 年京津冀地区典型城市大气颗粒物地面观测数据和气象要素数据为基础，基于观测数据的统计学方法，即（）滤波法，从时间、空间角度剖析区域大气颗粒物（包括 、和）污染特征及变化趋势，同时定量评估气象条件和排放对变化趋势的贡献，为制定京津冀地区大气颗粒物减排对策提供科学参考。材料与方法 研究区域及数据来源选取京津冀地区 个城市作为研究对象，即保定、北京、沧州、承德、邯郸、衡水、廊坊、秦皇岛、石家庄、唐山、天津、邢台和张家口。采用 滤波方法，选取 年 月 日至 年 月 日 个城市 、和 数据进行分析。空气质量数据来自于全国大气环境质量监控网格，包括、和 数据。相应气象数据来自于欧洲中期天气预报中心（）全球再分析数据，使用的气象要素包括：经纬风（和）、边 界 层 高 度（）、相 对 湿 度（）、总降水量（）、太阳辐射（）、距地面 温度（）、气压（）、云覆盖（）和 百帕经纬风和垂直输送变量（、和）。研究方法 滤波法 滤波方法计算过程简便，且无需对缺失值进行特殊处理，已广泛用于北京、天津、河北等地的大气污染指数、和颗粒物的相关分析。滤波器通常被应用于消除高频噪声、将低频成分与原 第 期 占 青等：基于 滤波法分析 年京津冀地区大气颗粒物变化趋势始信号分开以及消除季节效应 个典型场合。在此基础上，空气质量数据的时间序列可以表示为（）（）（）（）。（）式（）中，（）为原始序列；（）为长期分量；（）为季节分量；（）为短期分量。长期趋势是由总排放、污染物传输、气候、政策或经济方面的改变引起，季节分量是由太阳角度的变化引起，短期分量可以解释为天气和污染物排放的短期波动。滤波器是一种通过移动平均的重复迭代产生的低通滤波器。将（，）滤波器（表示窗口长度为 和迭代次数为 的滤波器）的移动平均定义为 。（）式（）中，为经 滤波后的时间序列；为 在滤波时其两端的滑动窗口长度；为滑动窗口长度，；为序列的时间间隔；为滑动窗口变量；为输入时间序列。迭代时，上一次传递的输出将成为下一次传递的输入。通过调整窗口长度和迭代次数，可以控制不同运动尺度的滤波。对于小于 天的周期，采用以下准则来确定滤波器有效宽度：。（）因此，（，）滤波器（表示窗口长度为 和迭代次数为 的滤波器）将去除（）的周期。方法可用来确定滤波器有效宽度，即滤波器权值约为 时的窗口大小，而不是截止周期，这是滤波器传递函数的一个性质。在以往 滤波 器 研 究 中，有 效 宽 度 和 截 止 时 间 都 被 使用，。采用（，）过滤器过滤后的气象和空气质量时间序列（分别为 和）被称为基线分量。将基线空气质量定义为长期和季节分量的总和，计算公式为（，）（）（）。（）每年的周期和更小的时间尺度可以通过选择更大的窗口来移除。（，）滤波器可以通过创建包含周期为 （）或 的周期的时间序列来从数据中提取长期分量，其结果是原始空气质量时间序列的长期成分：（）（，）。（）根据长期分量，可以计算出季节和短期分量。数据的季节分量是由空气质量和气象资料的基线序列减去长期成分而得到：（）（，）（，）。（）短期时间序列是由原始时间序列减去基线而得到：（）（）（，）。（）逐步多元回归模型通过对浓度时间序列的短期成分（）和基线成分（）的多元线性回归模型进行研究，结果表明：（）（）（），（）（）（）（）。（）式（）（）中，（）和（）分别为短期分量回归残差和基线分量回归残差；（）为短期分量回归的第个气象因子；（）为基线分量回归的第个气象因子；、和 为回归系数。污染物长期分量在均值附近的波动主要有 个方面原因：一是污染物排放的变动，二是天气状况的改变。在进行长期成分研究时，应将其与大气和污染物的环境因子区分开。所以，采用逐步多元线性回归模型是一个较好的选择。采用 建立各个时期大气污染物与气象条件之间的对应关系模型。浓度时间序列回归的总残差（）表达式为（）（）（）（）（）（）（）。（）式（）表示去除气象条件影响后，仅考虑排放对污染物浓度变化的影响。通过（，）对（）进行分析，可以得到受排放变化影响的污染物浓度的长期变化分量（）。将该变化分量与长期分量均值（）相加，可得到去除气象影响，仅受排放影响的污染物浓度长期分量时间序列（），通常也被称为经气象调整的污染物浓度长期分量时间序列：（）（）（）。（）评价指标（）为协方差总和与总方差之比，值越小，分解效果就越好，其计算公式为 （）（）（，）（，）（，）（）。（）式（）中，（）为 （）、（）和（）时间序列总方差；（）为污染物短期分量、季节分量和长期分量的总方差；（，）、（，）和（，）分别为长期分量与季节分量、长期分量与短期分量、季节分量与短期分量的协方差。生 态 与 农 村 环 境 学 报第 卷气象条件贡献率计算公式为 。（）式（）中，为气象条件贡献率；和 分别为原始数据趋势斜率和滤波后趋势斜率。结果与分析 总体污染概况京津冀地区西侧为太行山脉，北侧为燕山山脉，东临渤海湾，地形呈现西北高、东南低的特征，因复杂地形产生的山谷风是造成该地区颗粒物污染的重要原因。已有调查结果显示，京津冀地区主要受到北京西侧太行山东侧的西南向传输路径的影响，而西南部和东南部的传输路径为北京地区典型的污染物传输路径，在西南部的传播通量最大，有时甚至超过。同时，京津冀地区极易出现严重的沙尘暴。京津冀及其周围区域位于华北大草原，受到西北风沙的影响较大，而风沙不但增加 浓度，还增加 浓度。图 为京津冀地区 年 、和 年际变化趋势及平均浓度分布，京津冀 个城市 浓度均呈下降趋势，其中 浓度高值区以河北省南部地区（邯郸、邢台、石家庄和保定）为主，该区域 浓度下降趋势也尤为明显，下降幅度介于 之间；河北省南部地区（张家口、承德和秦皇岛）浓度最低，该区域 浓度下降幅度约为。个城市 浓度下降幅度介于 之间，其分布特征与 浓度分布特征一致，浓度高值区也以河北省南部地区（邯郸、邢台、石家庄和保定）为主，该区域 浓度下降幅度介于 之间。个城市 浓度下降幅度介于 之间，其分布特征与 和 浓度分布特征一致，高值区以河北省南部地区（邯郸、邢台和石家庄）为主，其下降趋势介于 之间。图 年京津冀地区 、和 年际变化趋势及平均浓度分布特征 ，（）图 为京津冀地区 、年均浓度和不同程度污染时段（优、良、轻度污染、中度污染和重度污染）占比。如图 所示，京津冀地区 年均浓度从 年的 迅速下降至 年的 ，年均下降约 ；年均浓度从 年的 迅速下降至 年 的 ，年 均 下 降 约 。第 期 占 青等：基于 滤波法分析 年京津冀地区大气颗粒物变化趋势根据 环境空气质量标准 和 环境空气质量指数（）技术规定（试行）评价标准，进一步分析颗粒物浓度发现，年 和 污染（即 平均浓度、平均浓度）天数分别为（占全年有效数据天数的）和 （占全年的），其中，和 重度污染（即 平均浓度、平均浓度）天数分别为（占全年的）和（占全年的不到）。年，和 污染天数分别下降为（占全年的）和（占全年的），其中，和 重度污染天数分别为 和 。有关研究结果表明，京津冀城市群颗粒物浓度分布不是完全随机的，而是表现出空间聚集性，城市群相似的高能耗工业结构布局是影响区域颗粒物污染的重要因素，因此，调整产业结构、统筹工业布局是降低京津冀区域颗粒物污染的主要措施之一。图 年京津冀地区 和 年平均浓度和各污染时段占比 典型城市 、和 浓度时间序列分解图 为采用 滤波器得到的 、和 时间序列各分量。在 滤波过程中，由于移动平均的重复迭代，时间序列开始和结束部分数据被截断。更长的窗口长度会导致额外数据丢失。因此，长期趋势分量比短期时间序列有更多数据被截断。需要注意的是，如果不应用时间序列分离方法，长期分量的相对较小的变化范围将很难从原始时间序列中分离出来。如图 所示，、和 长期分量均略有波动，但总体呈降低趋势，其中 、和 浓度分别降低 、和 。短期分量变化幅度明显，尤其冬天变化幅度最大。在短时间内，排放源变化不会太大，所以短期分量的剧烈变动可看作是气候的短期变动。季节分量表现出明显周期性变化，主要表现为 月和 月次年 月的波动。冬季大气颗粒物浓度高，可能是受到冬季取暖的影响。冬季供暖能源大多为一次能源，燃烧会产生大量颗粒物污染。除 生 态 与 农 村 环 境 学 报第 卷受供暖影响外，还受到冬季气象条件影响，例如，降雪对颗粒物具有高达 的去除率。另外，年 月 日至 年 月 日为停工期（暂停项目包括水利工程、土石工程、各种道路施工和房屋拆除），同时 年初新冠疫情爆发，综合各种情况，京津冀地区颗粒物浓度在冬季虽有所增加，但在一定程度上得到有效控制，未出现长期严重污染现象。在污染源一定的条件下，由于各季节气象因子不同，导致各地区颗粒物浓度存在明显季节性变化，如，秋季较弱的冷空气活动和低地面风速导致大气边界层降低，从而使得大气低层水汽积累，形成重污染天气。在春季，当风速大于 级时，强风使空气中颗粒物浓度急剧上升，同时由于重污染天气，颗粒物浓度也随之升高。夏季降水对颗粒物的去除效果显著，但夏季高温高湿天气也会使空气中污染物扩散条件出现不同程度恶化，从而使空气中颗粒物浓度增加，因此总体上夏季空气中颗粒物浓度变化不大。图 京津冀 个城市 、和 浓度时间序列分解结果 ，通过分析短期分量、季节分量和长期分量对原始序列的影响，可以得到各分量对原始序列的贡献率。在理想条件下，个成分是彼此独立的，所以原始序列变化应等于 个成分的变化之和。图 为京津冀地区 、和 各分量方差对原始序列方差的贡献，其中未知项为原始序列方差与 个分量方差和的差值。如图 所示，京津冀地区各城市 、和 的短期、长期和季节 个分量的总和均超过，且 值在 之间（表），说明 个分量基本上都符合独立性，证明 滤波法适用于京津冀城市群 个城市 、和 浓度的序列分解。个分量中，短期分量贡献最大，均在 以上，季节分量次之，长期分量最少。由此可以看出，原始时间序列的波动主要与短期和季节影响有关，即与污染源排放和气象条件的短期和季节变化有关。因此，为了更好地研究长期变化趋势，必须将长期分量与原始序列分离。减排措施与气象条件对颗粒物长期分量变化趋势的影响 图 为 年京津冀 个城市经气象调整后 、和 浓度长期分量序列以及气象条件对 、和 浓度长期分量 第 期 占 青等：基于 滤波法分析 年京津冀地区大气颗粒物变化趋势的影响。如图 所示，京津冀 个城市 、和 长期分量变化趋势在调整后与调整前基本一致，只在细节上有一些差异，长期分量与经气象调整后长期分量都呈现波动下降趋势，说明长期分量变化趋势主要与污染物排放变化有关。为了得到气象条件对长期分量的影响，且能更直观表现气象条件对 浓度的影响，将长期分量与经气象调整后长期分量做差值。当差值为正时，表明污染源导致的长期分量小于污染源和气象因素影响的长期分量，气象条件对污染状况改善是有利的；当差值为负时，表明由污染源导致的长期分量大于污染源和气象因素影响的长期分量，表明气象条件对环境改善不利。如图 所示，京津冀 个城市气象条件对 、和 长期分量的影响均在“”上下波动，这表明气象条件在某些时段有利，在有些时段则为不利。年京津冀 个城市气象条件对 、和 长期分量的影响呈上升趋势，说明该时段气象条件有利于污染状况改善，但是至 年底，气象条件的影响呈现下降趋势，且逐渐趋于“”以下，说明这时气象条件从有利于转为不利于大气污染状况改善。图 京津冀 个城市 、和 分量贡献 ，表 年京津冀 个城市污染物评价指标（），污染物 值保定北京沧州承德邯郸衡水廊坊秦皇岛石家庄唐山天津邢台张家口 生 态 与 农 村 环 境 学 报第 卷图 年京津冀 个城市经气象调整后 、和 浓度长期分量序列以及气象条件对 、和 浓度长期分量的影响 ，表 为京津冀城市群减排和气象条件分别对、和 浓度变化趋势的具体贡献。由于京津冀各城市气象条件和大气环境污染防治措施实施效果的不同，京津冀各城市气象条件和减排措施对 、和 长期分量改变的贡献也不同。就 而言，气象条件的贡献介于 之 间，减 排 措 施 的 贡 献 介 于 之间，其中，秦皇岛市气象条件对 的影响最大，占比为 ，衡水市减排措施影响最大，占比为 。就 而言，气象条件的贡献介于 之 间，减 排 措 施 的 贡 献 介 于 之间，其中，秦皇岛市气象条件影响最大，占比为 ，廊坊市减排措施影响最大，占比为。就 而言，气象条件的贡献介于 之 间，减 排 措 施 的 贡 献 介 于 之间，其中，沧州市气象条件影响最大，占比为 ，廊坊市减排措施影响最大，占比为 。为了验证 滤波模型的分离成果，等将基于 模型分离出的结果与基于 滤波模型分离出的结果进行对比，发现 种方法的结果相似：基于 滤波模型的结果显示，气象条件和减排措施对北京 浓度减少的贡献分别为 和；基于 模型结果显示，气象条件和减排措施对 浓度变化的贡献分别为 和。这 个模型都表明减排措施对大气颗粒物污染的改善至关重要。此外，张小曳等和 等研究结果也与之类似。这表明笔者研究结果具有一定可信度。结论基于 滤波法评估京津冀 个城市大气、和 变化趋势特征，以及气象条件和减排对变化趋势的贡献，主要结论如下：（）京津冀区域 年均浓度从 年的 迅速下降至 年的 ，年均下降 ；年均浓度从 年的 迅速下降至 年的 ，年均下降 。和 污染天数比例从 年的 和 分别下降至 年的 和。、和 浓度高值区主要分布在河北省南部邯郸、邢台、石家庄和保定等地，其变化趋势也最为显著。（）采用 滤波法对京津冀城市群 个城市大气 、和 浓度时间序列进行分解，发现短期分量贡献最大，表明原始序列的波动主要是受污染物排放量、天气和季节变化的影响。京津冀城市群 个城市大气 、和 浓度长期分量呈下降趋势。基于 滤波分析结果建立多元气象因子回归模型，得到气象条件和减排措施对颗粒物长期分量变化的影响规律。结果表明，京津冀地区长期分量变化趋势主要与污染物排放变化有关，说明 至 年京津冀地区实施的大气污染治理措施取得良好成效，同时，气象条件对 第 期 占 青等：基于 滤波法分析 年京津冀地区大气颗粒物变化趋势大气环境改善也存在一定影响。表 年京津冀 个城市减排措施和气象条件对颗粒物浓度的贡献 污染物减排措施气象条件排序贡献率 城市排序贡献率 城市 衡水 秦皇岛 沧州 保定 承德 张家口 廊坊 石家庄 天津 邢台 邯郸 北京 唐山 唐山 北京 邯郸 邢台 天津 石家庄 廊坊 张家口 承德 保定 沧州 秦皇岛 衡水 廊坊 秦皇岛 衡水 保定 天津 石家庄 邯郸 沧州 张家口 唐山 邢台 北京 承德 承德 北京 邢台 唐山 张家口 沧州 邯郸 石家庄 天津 保定 衡水 秦皇岛 廊坊 廊坊 沧州 张家口 北京 天津 石家庄 衡水 唐山 承德 保定 邯郸 秦皇岛 邢台 邢台 秦皇岛 邯郸 保定 承德 石家庄 衡水 唐山 天津 北京 张家口 沧州 廊坊（）年减排措施对京津冀地区大气颗粒物污染状况改善起着至关重要作用。就 而言，气 象 条 件 和 减 排 措 施 的 贡 献 分 别 介 于 和 之间，其中，秦皇岛市气象条件对 的影响最大，为 ，衡水市减排措施影响最大，为 。就 而言，气象条件和减排措施的贡献分别介于 和 之间，其中，秦皇岛市气象条件影响最大，占比 ，廊坊市减排措施影响最大，占比。就 而言，气象条件和减排 措 施 的 贡 献 分 别 介 于 和 之间，其中，沧州市气象条件影响最大，为 ，廊 坊 市 减 排 措 施 影 响 最 大，为 。参考文献：宋晓辉，杜亮亮，李剑东，等邯郸采暖期环境污染特征及空气质量预报方法研究生态与农村环境学报，（）：，（）：马国霞，於方，张衍燊，等大气污染防治行动计划实施效果评估及其对我国人均预期寿命的影响环境科学研究，（）：，（）：潘晓滨，曹媛“双碳”背景下京津冀大气污染协同治理研究资源节约与环保，（）：，“”，（）：张小曳，徐祥德，丁一汇，等 年气象条件变化对中国重点地区 质量浓度下降的影响中国科学：地球科学，（）：，（），（）：王跃思，李文杰，高文康，等 年中国重点区域颗粒物质量浓度和化学成分变化趋势中国科学：地球科学，（）：，（），（）：，（）：，（）：，生 态 与 农 村 环 境 学 报第 卷 ：，（）：，（）：，：（），（）：，：，：白鹤鸣，师华定，高庆先，等基于气象调整的京津冀典型城市空气污染指数序列重建生态与农村环境学报，（）：，（）：，（），：，（）：，：，（）：，（）：，：国家质量监督检验检疫总局，中国国家标准化管理委员会环境空气质量标准：北京：中国环境科学出版社，：张洁琼，王雅倩，高爽，等不同时间尺度气象要素与空气污染关系的 滤波研究中国环境科学，（）：，（）：赵域圻，杨婷，王自发，等基于 滤波的京津冀 年大气污染治理效果分析气候与环境研究，（）：，（）：田秀霞，田婷婷，周志博，等基于 滤波的臭氧评估研究环境科学与管理，（）：，（）：余益军，孟晓艳，王振，等京津冀地区城市臭氧污染趋势及原因探讨 环境科学，（）：，（）：秦人洁，张洁琼，王雅倩，等基于 滤波法的河北省 和浓度不同时间尺度分析研究环境科学学报，（）：，（）：，（）：，（）：，（）：，（）：董群北京山谷风环流及降雨对 浓度的影响北京：中国气象科学研究院，：，肖致美，徐虹，李鹏，等京津冀区域重污染期间 垂直分布及输送环境科学，（）：，（）：李鹏，肖致美，唐邈，等天津典型沙尘天气过程颗粒物污染特征生态与农村环境学报，（）：，（）：魏巍，皮冬勤，晏平仲，等 年春季华北地区一次典型沙尘重污染天气过程研究环境科学学报，（）：，（）：第 期 占 青等：基于 滤波法分析 年京津冀地区大气颗粒物变化趋势 李雪梅，许东明京津冀城市群 的空间分布及相关性分析生态与农村环境学报，（）：，（）：李星，赵文吉，熊秋林，等北京采暖季 水溶性无机离子污染特征及其影响因素生态环境学报，（）：，（）：黄亚林，刘超，曾克峰，等 年武汉市 的时空分布特征及其与气象条件的关系生态环境学报，（）：，（）：潘本锋，赵熠琳，李健军，等气象因素对大气中 的去除效应分析环境科技，（）：，（）：杨复沫，贺克斌，马永亮，等北京 浓度的变化特征及其与、的关系中国环境科学，（）：，（）：许万智，赵凤生，张寅超，等北京市区秋季气溶胶粒子浓度与特性参数的观测研究安全与环境学报，（）：，（）：樊守彬，田刚，秦建平，等道路降尘与扬尘 排放的关系研究环境科学与技术，（）：，（）：马启翔，刘赵梅，杜利劳西安南郊地区夏季降雨和风对可吸入颗粒物数浓度粒径分布的影响气象与环境学报，（）：，（）：古金霞，陈泽鑫，李琳希，等京津冀区域颗粒物污染特征研究南开大学学报（自然科学版），（）：，（）：作者简介：占青（），女，江苏盐城人，主要研究方向为大气环境。：（责任编辑：王昌群）