环渤海地区也称作环渤海经济圈，被辽东半岛、山东半岛和华北平原形成“C”字形所环抱，它主要包括辽宁的大连、丹东、锦州、营口、盘锦、葫芦岛，河北的秦皇岛、唐山、沧州，山东的青岛、东营、烟台、潍坊、威海、日照、滨州，江苏的连云港和天津共18个城市，分布着我国重要的海上货运进出港口群，与全世界众多国家保持良好的贸易往来. 针对环勃海经济圈城市，有学者^[19]提出“环湾保护、拥湾发展”的发展布局，以及环渤海地区和京津冀地区协同发展、创新驱动经济增长、进行生态环境修复和改善等的思路^[20].2013年是“一带一路”倡议的重要时期，也是环渤海地区的经济发展和生态环境均衡发展的重要转折点. 因此对环渤海地区2014-2018年的大气环境脆弱性进行研究，探索该区域大气环境污染、人口和自然压力、绿色发展、经济发展系统的状况以及政策对地区发展的影响，可以为该区域后续发展提供政策建议.

本研究基础数据来源于2014-2018年各省(市)环境质量公报、国家统计年鉴、省(市)级统计年鉴、国民经济和社会发展统计公报等；通过对基础数据进行计算得到二氧化硫排放强度、烟粉尘排放强度、环保经费投入强度、单位GDP能耗、人均绿地面积、第二产业增加值占GDP比重和科学技术经费投入强度相关数据.

本研究采用熵值法计算评价指标的权重，基于暴露性-敏感性-适应能力模型构建评价指标体系，得到环渤海地区大气环境脆弱性的最终评价结果；利用ArcGIS 10.5和地理探测器进行大气环境脆弱性空间特征探索.

(1) 权重计算. 熵值法的本质是一种客观赋权方法，是根据各项指标的观测数据提供的信息量的大小来确定指标权重的方法. 熵是热力学名词，在信息论中又称为平均信息量. 观测数据提供的信息量越大，不确定性就越小，熵值也越小，对综合评价影响越大；反之信息量越小，不确定性越大，熵值越大，对综合评价影响越小^[21]. 熵值法确定权重的计算公式为：

式中：P_ij为第i个评价对象在第j个评价指标下的特征比重；x_ij为第i个评价对象第j个指标的观测值. 由于P_ij在计算熵值时不能为负，对观测数据进行调整，需要在同一时期不同评价对象的某指标加上一个较小的值使评价对象这一指标的最小值大于零. 本研究利用以下公式对人口自然增长率进行了调整，将指标的熵值计算影响降到较小，

式中，k>0，e_j>0. e_j为第j项指标的熵值，定义k=1/lnn.

式中：g_j为指标的差异系数，g_j越大，指标作用越大；w_j为指标j在某一年份的权重，最终指标权重为5年该指标权重的均值.

(2) 大气环境脆弱性综合指数. 正向指标的指标观测值越大，大气环境脆弱性综合指数越高，脆弱性越强；负向指标的指标观测值越大，大气环境脆弱性综合指数越低，脆弱性越弱. 因此首先需要对指标进行预处理，使指标无量纲化后对大气环境脆弱性进行综合评价，计算公式如下：

正向指标标准化：

负向指标标准化：

大气环境脆弱性综合指数(AEVI)：

式中：x_ij为第i个评价对象的第j个评价指标观测值；x_ijmin为最小观测值；x_ijmax为最大观测值.

空间分异性是地理现象的基本特点之一，地理探测器是探索研究对象空间分异性的工具. 本研究主要使用地理探测器中的因子探测和交互作用探测^[22]进行计算，计算之前对数据进行预处理，对指标进行分级并将数值量转化为类型量.

(1) 因子探测. 该探测用以探测大气环境脆弱性Y的空间分异性以及影响因子X对Y空间分异性的影响力，用q值进行度量，具体公式如下：

式中：h为变量的分层，即X的分级；N_h和N分别为层的单元数；σ_h²和σ²是h和全区Y值的方差；SSW为层内方差和；SST为全区总方差. q的区间为[0, 1].

(2) 交互作用探测. 该探测用以识别不同影响因子交互叠加时对研究对象空间分异性的影响力.

一个合理的评价指标体系能够对评价对象进行全面分析，得到较好的评价结果，最大程度地反映研究区域存在的问题. 暴露性-敏感性-适应能力(ESA)模型在脆弱性评价的研究上具有广泛应用性，得到了国内外学者的普遍认可. 该模型是在人地耦合系统的基础上包含生态环境和社会经济，从而对评价对象脆弱性的内在机制进行深度解析^[23]. 暴露性指研究对象面对压力时的响应程度；敏感性指研究对象受到压力时的易改变程度；适应能力指评价对象处理不利因素的能力以及从中恢复的能力.

基于ESA模型，本研究选择大气环境污染作为一级指标研究大气环境的暴露程度. 根据大气环境质量标准和工业污染排放，选择SO₂，NO₂，PM_2.5，PM₁₀ 4种主要的大气污染物和AQI指数(空气质量指数，对空气质量进行定量描述)作为二级指标，反映大气环境污染状况. 工业二氧化硫排放强度和工业烟粉尘排放强度是指单位GDP的污染物排放量，反映大气污染物排放和经济的协调发展关系，对脆弱性的增强有正向影响，因此选择其作为大气环境的二级指标.

选择人口和自然压力为一级指标研究大气环境的敏感性. 人口会对大气环境造成压力，这种压力超过对人口的承载力，系统便难以进行自我调节，不利于可持续性的发展. 除此之外，由于气候变暖，大气温度的升高对人类社会和地球环境影响巨大，海平面上升，雾霾积累无法扩散，严重影响人类身心健康，因此选择人口死亡率、人口自然增长率、人口密度、平均温度作为人口和自然压力的二级指标.

选择绿色发展和经济发展为一级指标研究大气环境的适应能力. 绿色发展是将环境保护作为内在要素协调经济、社会、环境协调发展的一种新理念，依据2016年国家发改委印发的《绿色发展指标体系》考核内容，结合本文研究对象，从对大气环境脆弱性有直接影响或间接影响的角度，选择环保经费投入强度、能源消费总量、单位GDP能耗和人均绿地面积作为绿色发展的二级指标. 经济发展和生态环境相互作用，选择经济发展作为一级指标也是探索经济发展和大气环境发展内在联系的关键，结合环渤海地区海上贸易的特点，选择第二产业占GDP比重、科学技术经费投入强度、民用车拥有量和水运货运量作为二级指标研究地区经济发展. 评价指标体系如表 1所示.

大气环境脆弱性分级标准如表 2所示，综合研究区特点及已有研究成果^[24-25]，本研究按照0~1将其均等划分为5个等级.

根据熵值法权重计算结果，工业二氧化硫排放强度、工业烟粉尘排放强度、能源消费总量、民用车拥有量和水运货运量是指标体系中权重较大的前五项，其权重值分别为0.070 5，0.121 0，0.097 9，0.069 6，0.193 0，环渤海地区2014年和2018年该5项指标空间分布如图 2所示，指标分级采用等间距分级法. 由图 2可知，工业二氧化硫排放强度和工业烟粉尘排放强度的空间分布相似，西部城市明显高于东部城市，2018年数值大幅低于2014年，说明环渤海地区对于工业二氧化硫和工业烟粉尘的排放进行了有效控制，表明经济增长对大气环境的消耗呈减弱的趋势. 能源消费总量和民用车拥有量在2014年和2018年的空间分布上没有明显变化，环渤海西南地区的这2项指标明显高于东北地区，河北作为工业大省，环渤海经济圈范围内的河北城市能源消费总量略高于山东相关城市. 从权重计算结果来看，水运货运量的权重明显高于其他指标，环渤海经济圈的沿海特性，使水运货运量成为导致大气环境脆弱性的主要因素，水运会带来船舶硫氧化物排放、碳排放等问题，从图 2-e可以发现，环渤海地区天津、大连、烟台为主要货运港口，指标数据显示，2018年山东的青岛、威海以及江苏的连云港货运量明显提高，这是在国家“一带一路”倡议和“十三五”规划的推动下港口大力发展的成果.

对辽宁省、河北省、山东省和北京市的大气环境脆弱性进行综合评价，结果如图 3-a所示. 整体来看环渤海经济区的三省一市大气环境脆弱性呈波动下降的趋势，环渤海地区大气环境脆弱性在2014-2018年间由强到弱依次为2014年辽宁省、河北省、北京市、山东省，2018年辽宁省、河北省、山东省、北京市. 辽宁省大气环境为环渤海的三省一市中最脆弱地区，2016年脆弱性达到最弱后阶段性增强再下降，河北省在2017年大气环境脆弱性阶段性增强后下降，山东省大气环境脆弱性在2014-2016年增强后逐渐下降，北京市大气环境脆弱性波动性较小，处于长期平稳的减弱状态.

从单系统对大气环境脆弱性的影响(图 3-b)发现，辽宁省、河北省、山东省经济发展系统的脆弱性指数波动较大，对大气环境脆弱性影响最为显著，也是导致2018年3省份大气环境脆弱性较强的原因. 其次，国家对大气环境的重视和一系列的防治措施，使大气环境污染状况有了非常明显的改善，从图中可以看出，大气环境污染系统的脆弱性随时间变化下降，这也是令三省一市大气环境脆弱性下降的重要因素，人口和自然压力系统脆弱性较为平稳，辽宁省和北京市有下降的趋势，山东省和河北省指数则变化较小，主要是由于人口和自然变化是一个长期的过程，很难在短时间内有大范围变化. 图 3-b所示的绿色发展系统显示，山东省的绿色发展推行较好，对于大气环境脆弱变化具有积极作用，而辽宁省、河北省、北京市的绿色发展系统脆弱性增强，因此今后发展仍需重点关注地区绿色发展，以达到一个更好的大气生态环境水平.

(1) 环渤海地区省市评价结果. 对2014-2018年环渤海地区城市进行综合评价发现(表 3)，环渤海城市群大气环境脆弱性整体呈下降趋势，大气环境在面对外界冲击时的承受力增强，仅连云港和青岛2018年的大气环境脆弱性高于2014年. 环渤海地区北部的辽宁省部分，2014年大连、营口、盘锦、葫芦岛的大气环境脆弱性低于辽宁省整体水平，在辽宁省区域范围内呈现出相对较好的大气环境，2018年仅有盘锦脆弱性高于辽宁省，表明辽宁省沿海城市大气环境改善效果显著. 环渤海地区中部的河北省部分，唐山2014-2018年大气环境脆弱性变化小，对河北省脆弱性下降贡献小，秦皇岛与沧州大气环境脆弱性下降明显，2018年均低于河北省整体水平. 环渤海地区南部山东省域内的沿海城市，2014年大气环境脆弱性均高于山东省平均水平，到2018年，仅有青岛大气环境脆弱性增强，山东省整体以及其他沿海城市大气环境脆弱性下降，但青岛、烟台、威海、滨州4市高于省域平均水平，因此需要对这4个城市的大气环境状况进行把控.

环渤海城市群不同等级大气环境脆弱性面积比例如图 4所示. 2014年中度脆弱占比最大为50.21%，重度脆弱地区占比为7.67%；2015年重度脆弱性地区消失，中度脆弱性面积增加，较重脆弱性面积减小，同时出现了较轻脆弱性地区(占比5.18%)；2016年和2017年中度和较重脆弱性面积减小，较轻脆弱性面积增加；2018年较重脆弱性地区完全消失，出现轻度脆弱性地区(占比7.67%)，较轻脆弱地区占比最大达到50.78%.

(2) 重要港口城市单系统评价结果. 环渤海经济圈以大连为始辐射整个东北地区，以天津为始辐射京津冀地区，以烟台、青岛为始辐射中原地区，以连云港为始辐射珠三角地区，整体经济圈向南又辐射到朝鲜和韩国. “一带一路”倡议加强了环渤海经济区与其他地区的合作，而作为连接和起始的重要节点城市，经济发展下的环境发展显得尤为重要. 如图 5所示，天津、大连、青岛、烟台和连云港大气环境脆弱性单系统评价结果发现，5个城市的大气环境污染系统脆弱性逐年稳步下降(图 5-a)，只有连云港在2018年上涨明显. 对比大气污染系统各指标数据发现，连云港2018年二氧化硫排放强度比2017年有所增加，导致环境空气质量有所下降，水运货运量所带来的船舶硫氧化物排放是导致二氧化硫排放强度增加的重要因素之一，因此，在发展港口经济的同时要提升对大气环境污染的控制能力，包括污染源排放和污染治理技术的提高. 5个城市人口和自然压力系统脆弱性在2014-2018年处于一个波动性的状态(图 5-b)，没有明显的上升或下降趋势，因为短时间内人口和年平均温度的变化都是较为稳定的，青岛的人口和自然压力略有增长趋势，其主要是因为青岛人口密度的增加，提高了大气环境系统的压力. 除了天津外的其他4个城市的经济发展系统脆弱性在2014-2018年涨幅较大(图 5-c)，大连、青岛、烟台、连云港4市2018年民用车拥有量和水运货运量较2014年均有增长，是导致经济发展系统脆弱性发生变化的主要原因. 经济系统和大气环境污染系统联系紧密，相关部门应增加5座港口城市的科学技术投入强度，发展新兴技术产业，实现大气环境和区域经济的协调发展，提升与其他地区进行合作的环境基础. 绿色发展系统脆弱性(图 5-d)显示，除了大连外，其他4市都有上涨趋势，绿色发展很早就已经被提出，但真正作为国家重点发展理念去实行的时间较短，在国家自然资源日益紧缺的趋势下，能源消耗和单位GDP能耗都是需要重点关注的，对于前者应大力开发重点港口城市的新能源技术，例如风能、潮汐能等，对于后者应当实现能源的最大利用，不仅要保护和节约能源，更要依据环渤海地区的经济发展模式进行技术创新，优化产业结构.

(1) 大气环境脆弱性时空分布. 利用ArcGIS 10.5对环渤海地区城市大气环境脆弱性综合评价结果进行分级，对环渤海城市群大气环境脆弱性进行空间分类，得到相似的城市群. 分级结果显示(图 6)，2014-2018年，环渤海地区整体大气环境脆弱性呈逐渐减小的趋势，西部城市脆弱性减小幅度较大，2014-2018年山东半岛大气环境脆弱性从与辽宁地区相对平衡到高于辽宁地区，尽管地方政府已推进了环渤海地区一体化政策，但环渤海城市群大气环境发展仍存在不平衡的问题.

(2) 大气环境脆弱性空间分异分析. 利用地理探测器(Geodetector)对环渤海地区18个城市2014年和2018年大气环境脆弱性空间分异性进行因子探测和交互作用探测. 首先将指标的数值量转化为类型量，利用自然间断点分级法将2014年和2018年指标离散化为5个等级，将大气环境脆弱性和处理后的指标导入地理探测器生成结果.

因子探测结果显示(表 4)，2014年和2018年因子驱动力变化较大，但人均绿地面积和能源消费总量保持较高驱动力，q值分别为2014年的0.569 61(排名1)和0.411 37(排名1)，2018年的0.600 59(排名3)和0.398 33(排名2)，表明人均绿地面积和能源消费总量对大气环境脆弱性的空间分异性影响较大. 二氧化硫年均浓度、环保经费投入强度和单位GDP能耗在2014年对大气环境脆弱性空间分异影响力较大，但到2018年影响力下降明显. 第二产业增加值比重和民用车拥有量2014-2018年间对大气环境脆弱性的空间分异影响力增强. 分析结果显示，环渤海地区城市大气环境脆弱性空间分异性由绿色发展系统主导转化为经济发展系统主导.

因子交互探测结果显示，任意2种因子的交互作用均会增强对大气环境脆弱性的空间分异性，部分交互作用探测结果如表 5所示. 选取2014年和2018年因子交互作用较强的前19项进行探测，可以看出，人均绿地面积出现次数最多，与其他因子的交互作用强，尤其是与能源消费总量的交互作用最强，从因子分析角度出发，发现这2个因子在2014年和2018年均为大气环境脆弱性空间分异的主要影响指标，进行交互作用后的空间分异影响力更强. 2014年工业烟粉尘排放强度和人口自然增长率的交互作用最强，q值达0.997 84，2018年人口自然增长率和人均绿地面积的交互作用最强，q值为0.971 19. 因子探测器结果表明，人口自然增长率因子本身对大气环境脆弱性空间分异影响较小，2014年和2018年q值排名分别为15和17，但与人均绿地面积和工业烟粉尘排放强度的交互作用后对大气环境脆弱性的空间分异的发挥作用更强.

本研究采用暴露性-敏感性-适应能力模型构建了包含大气环境污染、人口和自然压力、绿色发展、经济发展4个子系统的环渤海地区大气环境脆弱性评价指标体系，采用客观赋权的熵值法进行指标权重的计算，将指标观测值进行标准化后进行各系统和整体的综合评价. 结果发现：①以省域和直辖市进行划分，2014年环渤海地区大气环境脆弱性由强到弱依次为辽宁省、河北省、北京市、山东省，2018年为辽宁省、河北省、山东省、北京市，三省一市大气环境脆弱性均呈下降趋势. ②从大气环境污染、人口和自然压力、绿色发展、经济发展4个系统对环渤海地区各省份沿海城市与所属省份进行比较分析发现，城市大气环境污染系统的脆弱性指数均下降，但辽宁省、河北省、山东省部分城市在2018年仍高于省平均水平；人口和自然压力系统的脆弱性指数没有明显变化；绿色发展系统的脆弱性指数波动性较强，但都处于一个较高的状态，城市间差异较大，除了秦皇岛，辽宁省、河北省沿海城市均低于省域水平，而山东省沿海城市均高于省域水平；经济发展系统的脆弱性指数波动性增强，辽宁省的大连、盘锦、葫芦岛，河北省的唐山，山东省的青岛、烟台、威海、滨州高于省域水平. ③从空间特征进行分析，2014-2018年环渤海地区整体大气环境脆弱性从四周脆弱性较强转变为中心脆弱性较强，地理探测器结果表明人均绿地面积和能源消费总量对大气环境脆弱性的空间分异性影响力最强，任意两因子的交互作用均增强地区空间分异性.

针对环渤海地区大气环境脆弱性评价结果，本研究提出以下建议：①在保证环渤海地区发展的同时，各级政府要明确环境保护职责，多角度考虑大气污染问题，除了工业生产排放污染、车辆排放污染外，面对日益增长的水上货运量，更多的船舶运输，制定针对性的控制政策，并加大对海上大气污染排放监测，颁布法规政策(如排放责任具体到各企业)，划定污染排放红线等. ②各城市对自己的优势进行定位，深入贯彻绿色发展理念，结合环渤海地区的地理位置，打造具有标签性的城市经济. ③减少一次能源消耗量，根据沿海城市独特的地理环境，加大力度开发风能、潮汐能以及水力资源等；减少一次能源对环境造成的污染，加大力度改进工业技术以提高能源利用效率，增大人均绿地面积，加大环保投入强度，最终增强大气环境的适应能力，降低大气环境脆弱性. 环渤海地区是“一带一路”沿线的必经之路，也是对外贸易的重要区域，未来应投入更大的力度加强港口建设和经济发展，更多关注经济发展所带来的环境问题.