气溶胶散射吸湿增长因子对霾强度变化的响应特征.pdf

生态与农村环境学报，（）：收稿日期：基金项目：四川省科技厅应用基础研发项目（）；国家重点研发计划（，）通信作者：佟景哲，米家媛，倪长健，等气溶胶散射吸湿增长因子对霾强度变化的响应特征生态与农村环境学报，（）：，（）：气溶胶散射吸湿增长因子对霾强度变化的响应特征佟景哲，米家媛，倪长健，蒋梦姣，杨寅山，陈云强（成都信息工程大学大气科学学院，四川成都；成都平原城市气象与环境四川省野外科学观测研究站，四川 成都；成都市环境保护科学研究院，四川 成都；四川省气象服务中心，四川 成都）摘要：基于成都市 年 月浊度计和黑碳仪的逐时观测数据，结合该时段同时次的能见度（）、相对湿度（）、二氧化氮（）以及细颗粒物（）浓度监测资料，利用光学综合法计算气溶胶散射吸湿增长因子，分析霾强度（轻微、轻度、中度和重度）变化对气溶胶散射吸湿增长因子的影响。结果表明：（）不同霾强度下的气溶胶散射吸湿增长因子均服从对数正态分布，该分布函数的数学期望与方差随霾强度的增加而增大。（）气溶胶吸湿性随霾强度的增加而增强，对轻微、轻度、中度和重度霾而言，强吸湿模态占比分别为 、和。（）分别构建了轻微、轻度、中度和重度霾强度下的二次多项式散射吸湿增长模型，对应全样本模拟值与实测值的决定系数（）分别为 、和 ，即气溶胶散射吸湿增长模型对样本的选择非常敏感，其泛化性能随霾强度的增加而增强。关键词：气溶胶散射吸湿增长因子；对数正态分布；强吸湿模态；二次多项式；霾强度中图分类号：文献标志码：文章编号：（），（，；，；，；，）：，（），（），（）（），（，）：（），（），（）（），：；在不同水汽条件下，单个气溶胶粒径、质量、密度、折射指数等微物理参数会发生变化，致使气溶胶粒子群宏观上的物理、化学及光学性质不断改变。气溶胶吸湿性不仅危害人体健康，也是大 生 态 与 农 村 环 境 学 报第 卷气环境、天气和气候过程的重要影响因子，在整个气溶胶科学研究中处于基础地位。研究发现，弱吸湿性气溶胶主要由含碳和尘埃颗粒构成，而强吸湿性气溶胶主要包含无机盐、被氧化的有机物和盐与有机物的内部混合物。气溶胶散射吸湿增长因子不仅表征气溶胶吸湿的光学效 应，也 是 计 算 气 溶 胶 辐 射 强 迫 的 重 要 参数。等利用 年实验观测数据，得到广州地区城市型、海洋型和城市海洋混合型气溶胶的散射吸湿增长因子分别为 、和 。等对韩国高山郡的气溶胶吸湿增长特性进行研究，得到沙尘型和生物质燃烧型气溶胶散射吸湿增长因子分别为 和 。程雅芳基于外场实验观测值，发现随着相对湿度从 增长到 ，散射系数可以增强 倍。另外，毛前军等观测不同粒径无机盐和碳质气溶胶粒子的吸湿增长发现，气溶胶粒径越小，其散射吸湿增长作用越明显。进一步研究指出，黑碳气溶胶通过吸附大量吸湿性物质（硫酸盐、硝酸盐和铵盐等），使黑碳气溶胶混合态、粒径以及化学组成显著变化，进而可能造成气溶胶吸湿光学效应的显著改变。张城语等基于相对湿度和黑碳质量浓度构建了成都地区气溶胶散射吸湿增长因子的双变量模型，显著提升了气溶胶散射吸湿增长因子的模拟精度。因此，气溶胶散射吸湿性增长因子的变化与气溶胶粒径分布、化学组成以及结构形态密切相关。霾是发生在大气近地层中以能见度降低为代表的灾害性天气。研究表明，我国大气污染具有高度复合污染特征，是一种新型的“霾化学”烟雾污染，二次气溶胶已成为大气气溶胶的主要部分，在重霾污染期间的占比可高达 左右。邱坚等研究发现，霾污染期间水溶性组分、的质量浓度显著增加，是非霾期间的 倍。吴丹等指出随霾污染的加重，粒径段颗粒物浓度迅速增加，细粒子占比增高。此外，相对湿度的增加可促进非均相化学反应发生，加速气态污染物向颗粒物转化，显著提升硫氧化率和氮 氧 化 率，从 而 影 响 二 次 气 溶 胶 浓 度 的 生成。由此可见，在霾的生消演化过程中气溶胶的理化性质也发生显著的改变，这必然会影响气溶胶吸湿的光学效应，其中的响应特征并不清晰。成都地区近年来社会经济发展迅猛，大气污染物排放量居高不下，受制于静稳和高湿的气象条件，该区域一直是四川盆地秋冬季污染最严重的区域之一。利用成都市 年 月浊度计、黑碳仪、环境颗粒物监测仪的逐时观测数据和该时段同时次大气能见度（）、相对湿度（）、二氧化氮（）监测资料，通过光学综合法计算了气溶胶散射吸湿增长因子，系统探究了该因子对霾强度变化的响应关系，据此深化对霾污染演化机理与气溶胶吸湿光学效应复杂性的认知。资料与方法 观测位置观测点位于成都市环境保护科学研究院综合大楼楼顶（，），距离地面高度，四周 内无高大建筑物，视野开阔。另外，观测点周围为集中居住区，范围内无明显工业大气污染源。数据获取和相关仪器浊度计和黑碳仪的逐时观测资料以及该时段同时次的环境气象监测数据（大气能见度、相对湿度、和 质量浓度）相关仪器介绍如下：（）浊度计（型，公司，澳大利亚）的观测波长为 ，采样频率为 次，切割头，检测范围 ，每 进行零点检查，零点漂移 ，每周用 气体进行跨度标定，通过内部温湿度传感器来控制浊度计内部加热系统，使得仪器内部腔室中气溶胶相对湿度控制在 以下，将其作为气溶胶的干燥状态。（）黑碳检测仪（型，公司，美国）获取等效黑碳质量浓度，有 个测量通道，波长分别为、和，数据采集频率为 次。黑碳检测仪采用 切割头，采样头与仪器连接中间增设硅胶管以减少水分对黑碳测量的影响。等效黑碳质量浓度统一用 表示。（）大气能见度由能见度仪（，公司，英国）进行监测，相对湿度由一体式气象站（，公司，德国）进行监测，质量浓度由化学发光、分析仪（，美国）监测，大气中细颗粒物 质量浓度由环境颗粒物监测仪（，公司，德国）实时测量，数据频率为 次。数据质量控制和筛选首先，将上述监测数据统一处理成小时均值数据；其次，剔除了出现降水、沙尘以及大风现象所在日的全部数据、仪器烘干后相对湿度仍大于 的 第 期 佟景哲等：气溶胶散射吸湿增长因子对霾强度变化的响应特征异常数据、超出界限值数据、连续无变化数据以及缺测数据；最后，根据 霾的观测和预报等级，在小时能见度数值低于 且环境相对湿度低于 时判断为霾；当环境相对湿度在 之间且 质量浓度高于 也判断为霾，能见度 时定义为重度霾，能见度在 时定义为中度霾，能见度在 时定义为轻度霾，能见度在 时定义为轻微霾，据此筛选获得匹配样本 个。气溶胶散射吸湿增长因子的计算 光学综合法大气消光系数代表光线在大气中传播单位距离时的相对衰减率，当对比感阈值（）为 时，在 波长处的环境大气消光系数与大气能见度的关系为。（）式（）中，为 波长处的环境大气消光系数，；为大气能见度，。环境条件下 波长处的大气消光系数的计算公式为 。（）式（）中，和 分别为环境条件下 波长处的气溶胶散射系数和吸收系数；和 分别为环境条件下 波长处大气的散射系数和吸收系数。采用 型浊度计观测 波长处干燥条件下的气溶胶散射系数。，。（）式（）中，为由，订正得到的 波长处干燥条件下气溶胶散射系数，；为成都市 波长指数，取值 。先用 型黑碳检测仪观测得到未经订正的 波长处等效黑碳质量浓度，然后利用黑碳质量 浓 度 反 演 和 波 长 处 的 吸 收系数。，（），。（）式（）（）中，为利用黑碳质量浓度反演 波长处的吸收系数，；为等效黑碳质量浓度，；为 波长处的吸收系数，。参照 的研究成果，对应 波长处的 一般取值为 。参照 等的计算方法，对应 波长处 的计算仅考虑 的吸收。（）式（）中，为 质量浓度，。利用光学综合法计算气溶胶散射吸湿增长因子（）：。（）函数参数化方案利用成都地区气溶胶粒径吸湿增长与散射吸湿增长的关系，构建以气溶胶粒径吸湿增长因子 为自变量的 函数参数化方案，计算气溶胶散射吸湿增长因子。()。（）气溶胶散射吸湿增长模型根据尹单丹等对幂函数、幂指数以及二次多项式 种气溶胶散射吸湿增长模型比对结果，二次多项式形式的散射吸湿增长模型可以最佳表征成都地区气溶胶散射吸湿增长因子随相对湿度的变化特征。（）。（）式（）中，为经验参数。对数正态分布函数对数正态分布在地质、水文、生物、材料、经济等研究领域已得到广泛的应用，其概率密度函数为（）（）（）。（）式（）中，为位置参数；为形状参数；为保证值参数。对应的数学期望（）、方差（）的计算如式（）（）所示，数学期望（）与方差（）决定了样本分布的均值与离散程度。，（）（）（）。（）结果与讨论 气溶胶散射吸湿增长因子分布对霾强度变化的响应 图 给出了成都地区不同霾强度（轻微、轻度、中度和重度）条件下气溶胶散射吸湿增长因子的频 生 态 与 农 村 环 境 学 报第 卷率直方图。由图 可见，不同霾强度条件下气溶胶散射吸湿增长因子均为右偏态分布，即相对正态分布曲线而言右侧偏长，左侧偏短；轻微、轻度、中度和重度霾下气溶胶散射吸湿增长因子的偏度分别为 、和 ，即气溶胶散射吸湿增长因子的右偏态程度随霾强度的增加而减小，说明气溶胶散射吸湿增长因子的散布范围随霾强度的增加而增大，样本的代表性也随之增强。图 种霾强度条件下气溶胶散射吸湿增长因子的频率 基于不同霾强度条件下气溶胶散射吸湿增长因子的右偏态分布特征，选取对数正态分布函数作为母体函数，据此对不同霾强度条件下气溶胶散射吸湿增长因子进行拟合分析。采用极大似然法进行参数估计，并借助 检验（检验）、相关系数（）以及均方根误差（）验证对数正态分布函数的适用性。检验、以及 的计算参见文献。不同霾强度条件下对数正态分布函数的拟合结果见表。表 种霾强度条件下气溶胶散射吸湿增长因子对数正态分布的拟合结果 霾强度参数估计值轻微 轻度 中度 重度 为位置参数，为形状参数，为保证值参数。为均方根误差，为数学期望，为方差。由表 可知，不同霾强度条件下的气溶胶散射吸湿增长因子均服从对数正态分布，通过了显著性 第 期 佟景哲等：气溶胶散射吸湿增长因子对霾强度变化的响应特征水平 的 检验（）；轻微、轻度、中度和重度霾强度条件下样本累积频率值与对应对数正态分布累积频率值的相关系数分别为 、和 ，对应均方根误差分别为、和 。由此可见，对数正态分布函数可以很好地表征不同霾强度条件下气溶胶散射吸湿增长因子的概率分布特征。进一步计算发现，轻微、轻度、中度和重度霾强度条件下对数正态分布函数的数学期望分别为 、和 ，数学期望随霾强度的增加而增大；在霾强度由中度上升至重度时增幅尤为明显，增幅为；另外，轻微、轻度、中度和重度霾强度条件下对数正态分布函数的方差分别为 、和 ，方差随霾强度的增加而增大；与数学期望类似，在霾强度由中度上升至重度时增幅尤为明显，增幅为。研究表明，霾强度的加重通常与高湿背景（值）有很好的对应关系，并且随霾强度增加，硫酸、硫酸铵以及硝酸铵这类二次水溶性无机粒子增多，小粒径粒子占比升高，由此导致气溶胶吸湿能力增强，表现为数学期望随霾强度的增加而增大，；另外，在静稳气象条件下，霾的演化与气溶胶多种理化过程有关，具有复杂性与不确定性等特点，这是方差随霾强度增加而增大的原因。由于环境气溶胶的化学组成与混合状态比较复杂，根据粒子吸湿性的不同，将其划分为不同吸湿模态，以反映气溶胶吸湿性能的强弱。为了进一步分析气溶胶散射吸湿增长因子对数正态分布参数变化的影响，基于文献的研究成果，按吸湿性的差异将气溶胶分为 种模态，即几乎不吸湿模态、弱吸湿模态和强吸湿模态，并根据可溶性体积比（）与气溶胶粒径吸湿增长因子（）给出了不同气溶胶吸湿模态的判定标准（表）。利用成都地区气溶胶粒径吸湿增长与散射吸湿增长的关系，得到气溶胶散射吸湿增长因子 在 种模态中的分布区间，在 时 为 几 乎 不 吸 湿 模 态，时为弱吸湿模态，时为强吸湿模态。基于气溶胶散射吸湿增长因子的对数正态分布函数，分别计算了不同霾强度条件下 种模态的出现概率，结果如图 所示。随霾强度的增加，几乎不吸湿模态和弱吸湿模态的占比总体出现降低趋势；在轻微和轻度霾强度下，强吸湿模态占比分别为 和 ，当霾强度上升至中度和重度时，强吸湿模态占比快速增长至 和 。表 种气溶胶吸湿模态的判定标准 指标几乎不吸湿模态弱吸湿模态强吸湿模态 为可溶性体积比，为气溶胶粒径吸湿增长因子，为气溶胶散射吸湿增长因子。图 种霾强度条件下各吸湿模态的占比 综上分析可见，成都地区气溶胶强吸湿模态占比随霾强度的增加而增大，这与气溶胶散射吸湿增长因子对数正态分布的数学期望随霾强度增加而增大是一致的，表明在霾的演化过程中气溶胶吸湿性与霾强度变化之间存在正反馈。不同霾强度条件下气溶胶散射吸湿性增长因子模型的适用性 由 节的分析可知，气溶胶散射吸湿增长因子的分布对霾强度的变化呈现出显著的响应特征，由此导致不同霾强度条件下气溶胶散射吸湿增长模型的差异。利用二次多项式散射吸湿增长模型对不同霾强度条件下的气溶胶散射吸湿增长因子进行拟合，相应参数和拟合结果如表 和图 所示。由图 可见，轻微、轻度、中度和重度霾强度下二次多项式散射吸湿增长模型的决定系数（）分别为、和 ，残差平方和（）分别为 、和 。这一分析结果表明，随霾强度的增加，对应气溶胶散射吸湿增长模型的拟合优度随之提升。由此可见，成都地区气溶胶散射吸湿增长模型对不同霾强度条件下建模样本的选择非常敏感，这无疑是模型不确定性的重要来源。生 态 与 农 村 环 境 学 报第 卷表 种霾强度条件下二次多项式散射吸湿增长模型的拟合参数 霾强度二次多项式散射吸湿增长模型拟合参数轻微 轻度 中度 重度 将对应的二次多项式散射吸湿增长模型对全样本进行模拟。将轻微、轻度、中度和重度霾强度下的二次多项式散射吸湿增长模型分别称之为模型、模型、模型 和模型，拟合效果如图 所示。实测值与模型、模型、模型 和模型 模拟值的决定系数（）分别为 、，残差平方和（）分别为 、，即模型的泛化性能也随霾强度的增加而增强。图 种霾强度条件下二次多项式散射吸湿增长模型的拟合结果 图 全样本的拟合结果 第 期 佟景哲等：气溶胶散射吸湿增长因子对霾强度变化的响应特征 综上可见，由于气溶胶散射吸湿增长因子对数正态分布的数学期望和方差随霾强度的增加而增大，故在中度和重度霾强度条件下气溶胶散射吸湿增长因子的代表性更强，能更好地表征气溶胶吸湿增长特性，这也是对应的二次多项式散射吸湿增长模型模拟和泛化性能更优的主要原因。从另一层面而言，中度和重度霾强度下气溶胶散射吸湿增长因子的占比不仅是衡量样本代表性的重要标准，也是决定研究区气溶胶散射吸湿增长模型适用性的关键因素。结论（）轻微、轻度、中度和重度霾强度下的气溶胶散射吸湿增长因子均服从对数正态分布，其数学期望与方差随霾强度的增加而增大。（）对轻微、轻度、中度和重度霾而言，强吸湿模态占比分别为 、和 ，气溶胶吸湿性随霾强度的增加而增强，表明在成都地区霾的演化过程中气溶胶吸湿性与霾强度变化之间存在正反馈。（）在轻微和轻度霾强度下构建的气溶胶散射吸湿增长二次多项式散射吸湿增长模型对应全样本模拟值和实测值的决定系数（）分别为 和 ，当霾强度上升至中度和重度时决定系数快速增长至 和 ，即成都地区气溶胶散射吸湿增长模型对建模样本的选择非常敏感，其泛化性能随霾强度的增加而增强。参考文献：，：，：，（）：，（）：，（）：，：，（）：，（）：，（）：，（）：，：，（）：，：，（）：程雅芳珠江三角洲新垦地区气溶胶的辐射特性：基于观测的模型 研 究 北 京：北 京 大 学，：，毛前军，王莉典型雾霾粒子及团簇的散射吸湿增长特性光学学报，（）：，（）：谭天怡，郭松，吴志军，等老化过程对大气黑碳颗粒物性质及其气候效应的影响科学通报，（）：，（）：，（）：张城语，倪长健，佟景哲，等成都地区气溶胶散射吸湿增长因子双变量模型中国环境科学，（）：，（）：张人禾，李强，张若楠 年 月中国东部持续性强雾霾天气产生的气象条件分析中国科学：地球科学，（）：，生 态 与 农 村 环 境 学 报第 卷 （），（）：，（）：楚碧武，马庆鑫，段凤魁，等大气“霾化学”：概念提出和研究展望化学进展，（）：，“”：，（）：张小曳，孙俊英，王亚强，等我国雾霾成因及其治理的思考科学通报，（）：，（）：，（）：邱坚，田苗苗，王晨波，等冬季南京城区气溶胶化学组分和吸湿性观测中国环境监测，（）：，（）：吴丹，曹双，汤莉莉，等南京北郊大气颗粒物的粒径分布及其影响因素分析环境科学，（）：，（）：，：刘凡，谭钦文，江霞，等成都市冬季相对湿度对颗粒物浓度和大气能见度的影响环境科学，（）：，（）：，：，霾的观测和预报等级 ，：杨寅山，倪长健，邓也，等成都市冬季大气消光系数及其组成的特征研究 环境科学学报，（）：，（）：李梅芳，叶芝祥基于太阳光度计的成都双流地区夏季气溶胶光学特性研究成都信息工程学院学报，（）：，（）：，：，（）：，（）：，：，（）：张智察，倪长健，张城语，等气溶胶粒径吸湿增长与散射吸湿增长的关系中国环境科学，（）：，（）：尹单丹，倪长健，邓也，等成都秋冬季气溶胶散射吸湿增长模型研究生态与农村环境学报，（）：，（）：臧红霞，郑华山，陈长茵三参数对数正态分布最大熵参数估计方法探讨水利科技与经济，（）：，（）：，（）：李琦，银燕，顾雪松，等南京夏季气溶胶吸湿增长因子和云凝结核的观测研究中国环境科学，（）：，（）：作者简介：佟景哲（），男，满族，吉林通化人，研究方向为大气物理学与大气环境。：（责任编辑：陈 昕）