20世纪初, 美国学者Clements将Ecotone这一术语引入到生态学中^[1], 之后衍生出生态脆弱性的概念, 但直到现在, 生态脆弱性仍未有被大众接受的概念或定义。综合现有的研究可总结出：生态脆弱性是当区域生态环境受到外界干扰时, 所表现出的抗干扰能力弱, 受干扰后恢复能力低, 且发生转变后难以恢复的性质特点。全球变暖、冰川融化、干旱、洪水、污染以及生物多样性丧失等环境问题使得生态环境异常敏感, 生态脆弱区逐步扩张, 开展生态脆弱性研究是相当必要的。

国内外对生态脆弱性的研究已从初期的定性描述逐渐演变为普遍的定量评价^[2-5], 研究尺度精细化显著增强^[6-7], 研究领域也在不断拓宽, 如林草交错带^[8]、流域^[9]、海岸带^[10]等, 但对于河口三角洲这一特殊区域的生态脆弱性研究相对较少。目前生态脆弱性评价的方法众多, 常用的包括模糊隶属度评价法、层次分析法、灰色关联评价法和主成分分析法等^[11-14], 各个方法都有自身的缺点, 如层次分析法受人为主观思想的影响较重, 灰色关联度评价法的复杂计算过程被大多数专家诟病。评价指标是整个研究的核心之一, 学者们建立了多种指标选取框架, 如“压力-状态-响应”模型^[15]、“暴露度-敏感度-适应性”模型^[16]以及“驱动力-压力-状态-影响-响应-管理”模型^[17]等, 随着自然环境和社会发展中各种科学问题的提出, 生态脆弱性指标的全面化也在不断完善。

黄河三角洲是由黄河长年冲积形成的, 该地区拥有丰富的自然景观类型, 具有重要的生物多样性保护和生态屏障作用, 但其充裕的土地和石油资源也带来了严重的人类活动干扰, 再加上频繁的自然灾害影响, 出现了土壤退化、土地浪费、水体污染和海岸侵蚀等问题, 使该地区生态环境受到严重威胁, 所以对该区域进行生态脆弱性评估对区域生态保护和管理具有重要的科学意义和现实意义。目前也有部分学者对黄河三角洲的生态脆弱性做了评价研究, 如德国宇航局的Michel通过实地调查和问卷咨询的方式开展研究, 对农村环境在受到内外干扰情况下的脆弱性状况进行定性描述, 但缺少相关的定量数据支持^[18]；其他相关研究多是完全利用层次分析法进行评价^[19-21], 在指标的权重设定中受人为主观因素的影响过重。

本文拟采用模糊评价方法与层次分析法相结合对黄河三角洲生态脆弱性进行评价, 也称为模糊层次分析法, 将模糊三角函数作为层次分析法中指标重要性比较的赋值方式, 可以降低人为主观作用太强的缺点^[22-23], 而层次分析法为模糊评价法选取指标时提供了更清晰的思路和逻辑^[11], 此外, 本研究利用模糊隶属度作为指标和评价结果的定量化分类依据, 削弱了以往布尔运算等级界线太过严格的缺点。本文研究目的一方面是验证方法的可行性, 另一方面根据研究区的生态脆弱性状况, 从空间和定量统计方面做原因分析, 为黄河三角洲的可持续发展提供必要支撑。

1 研究区与数据 1.1 研究区概况

黄河三角洲位于中国山东省东北部, 属东营市辖区, 坐标区间为37°22′—38°04′N和118°14′—119°05′E, 总面积约为5062.59 km², 东部和北部与渤海相邻(图 1)。研究区地形平缓, 呈西高东低, 南高北低趋势, 高程在0.00 m至12.50 m之间, 受黄河改道和冲积作用, 微地貌类型多, 包括岗阶地、河成高地、低洼地、河滩地、平地以及滩涂地等。黄河三角洲属于温带大陆性季风气候, 降雨集中于6—9月, 年均蒸发量远大于降水量, 分别为1885.00 mm和537.40 mm；研究区土壤主要为潮土, 盐渍化重, 砂粒含量较高, 土壤质量普遍较低。研究区自然植被以草本植物为多, 主要是耐盐性的芦苇、柽柳、白茅和翅碱蓬等, 人工植被以农作物为主。研究区东部和北部设有两个自然保护区—黄河口自然保护区和一千二自然保护区；研究区人为活动日趋严重, 农业开垦范围不断扩张, 工业建设以及石油资源的开发增加了原有自然景观的破碎度。

1.2 评价指标系统建立

生态脆弱性评价指标和因子选取的合理性和科学性, 决定着评价过程的可行性以及评价结果的可靠性, 所选指标要能够充分反映研究区面临的主要生态环境问题。目前黄河三角洲面临的主要生态环境问题可大致总结为：海洋潮汐和风暴潮导致的海水入侵造成陆地原有生态环境的破坏；地下水埋深较浅, 地下水矿化度普遍较高, 土壤盐碱化严重；人类活动极大改变了自然景观, 也导致了土壤的次生盐渍化；农业化肥和农药的使用, 工业的排污造成了严重的面源和点源污染, 加剧了研究区的生态脆弱性。

根据以上的生态环境问题, 本研究参照“压力-状态-响应”模型框架从地下水状况、土壤、土地利用状态、地形地貌、植被覆盖、社会经济、气候和海洋影响8个角度选取了21个评价指标, 由于本文是将层次分析法与模糊理论相结合, 该理论未将各指标划分至具体类别(压力、状态、响应)中, 详细体系如图 2和表 1所示。

2 数据来源研究方法 2.1 数据获取与处理

地下水数据：课题组近些年在黄河三角洲埋设了16口地下水井用于研究, 本研究获取了2014年每口井的平均地下水位和平均地下水矿化度, 利用普通克里格插值得到两者在整个研究区的分布状况；

海洋影响：本研究收集了1984年至2014年间13期海岸线分布特征, 统计整个海岸每个位置的摆动频率, 作为海岸侵蚀系数；距海岸距离以2014年海岸线为基准, 利用欧氏距离法计算研究区各位置与海岸线的距离。

气候要素：从研究区内部和周边地区选取了13个气象站点, 并从中国气象数据网(http://data.cma.cn/site/index.html)获取了各站点的累年日均气温和累年日均降雨量, 并得到相应的年均≥10℃积温和年均降雨量, 进而获得干燥度(K)^[24], 利用普通克里格插值得到各参数在研究区的分布状况。

社会经济统计数据：查询2015年中国统计年鉴, 获取研究区内行政区的人口数量和GDP, 参考廖顺宝等人^[25-26]提出的人口空间化方法, 以“人口-居住地-影响要素”模式获取研究区的人口密度分布状况, 继而得到GDP密度分布。

植被覆盖度：本研究以2014年9月份Landsat8 OLS数据为基础, 首先获取研究区的NDVI指数, 并利用以下公式计算得到研究区植被覆盖度：

本研究中将NDVI_veg和NDVI_soil分别设定为研究区内的NDVI最大值和正的最小值, NDVI小于等于0的区域直接设定其覆盖度为0。

其他数据：土壤含盐量、土壤质量和2014年土地利用状况引用自吴春生对黄河三角洲土壤质量的研究成果^[27-29]；研究区DEM、微地貌数据、土壤类型和土壤质地取自中国科学院地理科学与资源研究所资源环境数据中心(http://www.resdc.cn/)。土地垦殖率、人类干扰指数、道路密度和水渠网密度数据源均为研究区2014年土地利用现状。所有数据均投影至统一坐标系下, 并重采样为30 m空间分辨率。

2.2 模糊层次分析法

首先, 参照层次分析法建立指标对比矩阵, 根据有关黄河三角洲研究的专家给出的建议, 多次修改对比矩阵, 并通过一致性检验, 参考Kahraman^[30]设定的模糊重要性对比赋值规则对原对比数值进行替代(表 2)。

其次, 计算指标权重。假设任意一个指标i与另一个指标j的重要性对比结果可以表示为(l_ij, m_ij, u_ij), 则指标i的累积扩展值可以表示为：

整个矩阵的累积扩展值表示为：

指标i的综合扩展值为：

l, m, u为模糊三角函数3个顶点的横坐标, n为指标个数。对于任意两个综合扩展值S_i和S_j, S_i ≥ S_j的概率大小可表示为：

设d_i′为S_i与其他综合扩展值比较后的最小值, 即：

对w′标准化即可获得各指标的权重：

在计算过程中会产生V(S_i≥S_k) =0, 导致d_i′=0, 致使最终权重不合理, 在此情况下对指标对比矩阵中的元素进行标准化处理, 即可避免最终权重出现0的情况^[31]。

2.3 模糊逻辑理论模型

模糊逻辑模型即是用模糊隶属度的方式来表达某个指标隶属于某个级别的概率, 而不是将指标指定为某一固定的等级^[32], 常用的隶属度函数为：

其中0＜MF_xi ≤1, MF_xi 为指标i的隶属度值, x_i为i的值, d为i的过渡区间的宽度, 一般选取隶属度值等于0.5和1时的指标值的差作为d的值, b为隶属度等于1时的指标值^[33]。

对于每个指标需要根据研究目的设定其适宜范围, 当超出该范围时其隶属度值全为0或1, 范围之内的则利用隶属度函数来计算具体的隶属度值。而在确定适宜范围是往往需要借助于已有的研究成果、文献资料以及相关部门设定的标准规范等。

2.4 生态脆弱性评价方法

研究区最终的生态脆弱性评价利用加权求和的方法, 即：

EVI即为生态脆弱性值, w_i为各指标的权重值, A_i为各指标的隶属度。

3 研究结果

依据模糊层次分析法建立过程, 各一级指标的重要性比较矩阵见表 3。

由于二级指标较多, 不再逐个展示, 参照权重计算的方法和过程, 各指标的权重如下表 4所示。

从表 4中看出, 一级指标中的土壤条件权重为0.31, 比其他各指标权重都大, 说明土壤属性特征对于黄河三角洲生态安全稳定具有重要的作用, 而由于研究区地形平缓, 无较大起伏, 所以地形地貌的影响相对较小, 其所占权重也相应较小。

依据模糊逻辑模型构建过程, 经查阅与各指标有关的文献资料、书籍和已有的研究结果^{[22, 24, 34-36]}, 并结合研究区实际情况, 确定了各指标的脆弱性适宜范围、各指标的趋向性以及各指标在模糊函数中的b值和d值。定量指标详情如下表 5所示：由于定性指标无数值型界线, 无法利用模糊函数计算隶属度, 故根据相关资料指定各定性指标的隶属度, 详情如表 6和表 7所示。

根据上述的适宜性范围并参照模糊隶属度函数获取各定量指标在研究区内的生态脆弱性隶属度空间分布状况, 结合各指标的权重, 利用加权求和模型获取研究区整体的生态脆弱性隶属度, 并利用自然断点法将脆弱性结果分成6个级别, 各级别空间分布及面积统计如下图 3和下表 8所示。

从图 3中可以看出, 在空间位置上越靠近海岸, 脆弱性越高, 以黄河入海口周边和研究区西北沿海最高, 向内陆逐渐降低, 研究区中部黄河和刁口河交叉处周边脆弱性最低, 整体规律与预计结果相符合, 与其他研究成果相比, 整体结果也更为合理。

从统计结果看(表 8), 重度脆弱区分布面积最大, 范围最广, 面积达到1113.25 km², 占总面积的22.03%, 相对来说研究区北部沿海更多, 尤其刁口河东侧面积较广；其次是极度脆弱区, 面积为952.85 km², 该区域分布更靠近海岸线, 除刁口河入海口两侧和一千二自然保护区沿海外, 其他沿海地区基本上为极度脆弱生态环境；中度、轻度、微度和不脆弱区被极度脆弱区和重度脆弱区包围, 且四种类型整体上也是由外向内脆弱度逐渐降低, 中度脆弱区和轻度脆弱区处于外围, 两者呈镶嵌分布, 不脆弱区与微度脆弱区在最内侧, 也呈零星交叉分布, 不脆弱区面积最小为464.17 km²。

为了更好的对研究区生态脆弱性做原因分析, 根据各指标设定结果, 土地状况所占权重仅次于土壤条件, 故将评价结果与土地利用现状进行叠加, 重点分析几个主要土地利用类型的生态脆弱性特点, 结果如下表 9所示。表 9显示, 农田面积最大, 其主要分布在前3个脆弱性级别, 以微度脆弱区面积最大为405.06 km², 其次是不脆弱区和轻度脆弱区, 农田在重度脆弱区和极度脆弱区也有少量分布。盐碱地同样在6个脆弱性级别区均有分布, 相对于农田来说, 除了不脆弱区, 在其他级别面积分布较为均匀, 轻度脆弱区和中度脆弱区面积基本相等, 重度脆弱区面积最大为290.10 km²。根据整个研究区脆弱性分布可以看出, 沿海滩涂基本上都属于极度脆弱类型, 面积达到359.50 km², 而在一千二保护区北侧的滩涂则属于重度脆弱区。居民地以分布在轻度脆弱区和中度脆弱区为主, 草地和林地在各级别内分布面积较均匀, 内陆滩涂和园地总面积较少。

4 分析和讨论

从评价结果中生态脆弱性各级别的空间分布状况来看, 空间格局较为合理, 与预期设想和实际的自然条件影响相符, 也证明了该方法的合理性和适用性, 可以为其他类似研究提供借鉴和依据。与前期黄河三角洲生态环境脆弱性研究成果相比^[19-21], 模糊理论的加入使得评价过程中人为主观要素影响得到减弱, 隶属度函数的使用不但实现了定量化评价, 而且在评价等级划分上更加灵活, 评价结果更为详细且在空间上的连续性强。

从评价指标上看, 表 4显示一级指标中土壤条件所占权重最大, 其中又以土壤含盐量和土壤质量最为重要, 所以生态脆弱性的空间分布受土壤含盐量和土壤质量的空间分布影响较重, 从吴春生对土壤质量和土壤含盐量的研究成果看(图 4), 两者与生态脆弱性评价结果在空间上具有相似性, 土壤质量与生态脆弱性的空间相关性达到-0.55；在土壤质量较差的区域, 生态脆弱性高, 如黄河入海口周边, 研究区的西北部和东南部等, 而在刁口河与黄河交叉处周边土壤质量较高, 土壤含盐量低, 故其为不脆弱区。所以可以从如何提高土壤质量和降低土壤含盐量方面采取措施, 来降低区域生态环境的脆弱度, 如减少对耐盐植被的破坏, 提高区域的植被覆盖, 从而减少因蒸发而产生的表土聚盐现象, 植被的枯枝落叶也能够提高土壤的肥力。

另外与生态脆弱性有较大关系的是土地状况和地下水条件。土地利用状况中人类干扰指数影响占比重较大, 尤其是农田的开发和人工湿地(养殖和盐田)的建设对生态环境的改变较大, 上述土地利用状况与生态脆弱性结果叠加统计结果以及下图 5中主要土地利用类型的空间分布状况显示, 一部分农田处于生态脆弱性较高的区域, 这不利于生态环境的维持和保护, 尤其当土壤肥力不足, 作物产量降低时, 该区域往往被弃耕, 进而产生次生盐渍化, 处于高生态脆弱性等级的盐碱地中有部分即为此类土地；另外, 养殖和盐田开发破坏了原有生态环境的演化方向, 加重了生态脆弱性程度。所以避免在脆弱性高的区域进行人为开发, 对于保护区域生态环境有重要的作用。

图 5中地下水埋深的空间分布状况显示其与生态脆弱性的空间分布格局也存在一定的相似性, 两者的空间相关性达到-0.74, 地下水矿化度较高, 其埋深越浅, 地表返盐越严重, 长期的低植被覆盖率也使得净初级生产力较低, 当受到外界干扰时, 易改变原有状态, 且难以恢复。

此外, 对比研究区沿海区域的脆弱性空间差异可以看出, 在设有防护措施或修复工程的位置, 生态脆弱性相对较低, 如从黄河口自然保护区北部至刁口河河口西岸的沿海岸段, 分别建有不同类型的防护堤, 有效阻挡了海水和潮汐向陆地的侵入, 保障了植被的生长, 从而提高了生态环境的稳定性；而黄河入海口及其周边区域由于要为动植物提供海陆交互的栖息地环境, 故保留原有自然状态, 减少人为影响则是正确的处理方式；所以根据不同的岸段类型, 采取相应的保护措施对于黄河三角洲生态稳定性具有积极作用。

与其他利用模糊层次分析法对生态脆弱性进行评价的研究相比, 方法运用上存在很大差异, 相互之间缺少对比, 无法得知所选的指标对比函数以及权重计算方法的优劣性, 这在以后的研究中还需要进一步探索；另外本研究只是针对2014年的黄河三角洲生态脆弱性做了评估, 所以在做原因分析时只能从现状上讨论, 总体上论据较为牵强, 说服力不足；为更好的分析黄河三角洲生态脆弱性空间分布特征, 还需要利用多期数据进行对比, 动态的展示空间上的变化, 但由于前期的数据不足, 如土壤质量等, 只能期待后期数据的完整性, 而且还要注意本研究中选取的指标在后期研究中的适用性等。

5 结论

本研究利用模糊层次分析法和模糊隶属度理论对黄河三角洲生态脆弱性进行了评价, 评价结果显示生态脆弱性级别分布具有一定的规律, 从沿海到内陆脆弱性逐渐减弱, 其中以黄河入海口周边、东南沿海以及西北沿海等地区脆弱性最高, 而处于内陆的黄河和刁口河沿岸为不脆弱区域, 空间分布格局整体上合理可信, 也证实了方法的可用性；土壤含盐量作为黄河三角洲土壤的特殊指标, 其对生态脆弱性影响较大, 以土壤养分为主要内涵的土壤质量对其影响更大, 空间相关性达到-0.55；地下水埋深直接控制着土壤含盐量的大小, 其与生态脆弱性的空间相关性大小为-0.74；另外, 人为活动如农田耕作、盐田和养殖开发以及石油开采对研究区生态环境产生了不利影响, 而沿海堤坝的建设则为生态环境免受海水侵扰提供了屏障, 有效保护了生态环境。为更好的分析黄河三角洲生态脆弱性原因以及提供保护生态环境的方法, 未来要着眼于生态脆弱性的动态研究, 利用多期结果进行对比。