生态安全这一思想最早出现于20世纪40年代的土地功能与土地健康评价当中, 其概念构筑在环境安全的理论基础之上^[1]。1989年国际应用系统分析研究所(International Institute for Applied Systems Analysis, IASA)在阐释全球生态环境问题时首次提出“生态安全”这一概念^[2]。随着生态安全理论框架的日趋成熟, 这一概念愈加明确：生态安全指一个区域内维持人类生存与发展所需的生态系统服务处于少受甚至不受威胁的状态, 使生态环境既能满足人类与生物群落生存与发展的需要, 又能保证生态环境本身的生态过程不受损害, 使其与社会经济处在可持续发展的状态^[3]。由此可见, 生态系统服务与生态安全有着紧密联系, 将生态系统服务与生态安全相结合, 有助于为生态环境管理提供依据^[4]。

传统的生态安全评价研究一般从生态系统的结构和功能角度去构建指标体系, 而将人类视为外部要素不加以考虑^[5]。然而在现实情况中, 人类早已成为生态安全格局中不可或缺的部分^[6]。基于以上情形, 生态系统服务成为耦合人类与其周围生物与非生物环境的有效纽带, 生态系统服务的可持续性成为人类直观感知生态是否安全的直接途径^[7-9]。因此, 生态系统服务理论的融入为生态安全评估提供了正向的工具, 将人类对生态系统服务的需求与生态系统所能提供相应服务的能力进行耦合, 成为目前生态安全评价中的新视角^[10-11]。

基于生态系统服务理论构建生态安全评价指标体系应同时考虑生态系统正常运行与人类社会的可持续发展^[3]。彭保发等^[10]通过辨识生态系统服务的供给水平与需求热点, 综合考虑生态系统服务的传输网络, 基于压力-状态-响应(Pressure-State-Response, PSR)模型^[12-13]构建了耦合生态系统服务的生态安全评价机制。在该机制中, 人类活动作为驱动因子, 对生态系统服务的获取强度对生态系统造成压力, 改变了生态系统的结构与功能, 生态系统服务的传输能力也受到影响, 最终生态系统服务流量对此作出响应。该框架能够评价生态系统服务传输网络的流量、稳定性与可持续性, 有助于揭示人类活动和生态安全之间的相互作用, 可以作为生态安全评价的基本框架。

黑河流域是中国第二大内陆河流域, 随着经济社会的快速发展, 黑河中游绿洲的耕地和建设用地不断扩张, 相应的林地、草地不断萎缩^[14-15]。黑河中游的绿洲扩张影响了生态系统服务的供给与传输, 进一步影响了区域的生态安全格局。针对黑河中游的绿洲扩张情景, 基于耦合生态系统服务理论的PSR模型框架^[10], 构建了黑河中游生态安全评价指标体系, 评估区域生态安全变化, 为制定生态管理策略、支撑生态管理权衡与决策等提供了支撑^[16]。

1 研究区概况

本文拟选择典型荒漠-绿洲黑河中游为研究区, 黑河中游位于河西走廊中段, 南部为祁连山脉, 北部紧邻合黎、龙首两山脉(97°20′—102°12′ E, 38°28′—39°50′ N), 主要包括甘肃省张掖市的甘州区、高台县及临泽县(图 1)。本区域海拔介于1252—3609 m之间, 大部分属于黑河的山前冲洪积扇, 水分条件相对较好, 形成了独特的荒漠-绿洲景观。黑河中游受大陆性气候和青藏高原的共同影响, 属中温带甘-蒙气候区, 年平均气温介于-1.1—9.5 ℃之间。由于远离水汽输送通道, 该处大部分地区年均降水量小于200 mm, 年均蒸发量超过1600 mm。黑河中游紧邻荒漠, 土壤类型以荒漠土与灰漠土为主(主要分布在绿洲和荒漠过渡地带)。黑河中游生长的天然植被主要为耐旱、耐盐碱的乔木、小灌木与半灌木等。中游绿洲农业区则具有显著的人工绿洲特征, 植被以人工灌溉栽培农作物和防风沙的行道林木等为主。该区域长期受到人类生产生活的影响和改造, 区域生态安全、用水调配和环境保护等关乎国计民生的问题始终较为突出。作为我国生态环境比较脆弱的地区之一, 黑河中游成为诸多学者关注的焦点区域^[17]。

2 研究方法与数据来源 2.1 研究方法 2.1.1 生态安全评价方法

构建PSR指标体系。PSR模型包括目标层、准则层与指标层3个部分：目标层即为区域生态安全, 准则层包括压力、状态及响应3个部分, 指标层则为指标体系中的各个具体指标。在参考相关研究^{[10, 18-20]}的基础上, 基于生态系统服务理论构建了黑河中游生态安全PSR指标体系：压力指标表示生态环境出现问题的原因, 侧重描述外界对生态系统的干扰及人类对生态系统服务的获取强度^[10], 例如伴随着经济发展对地下水等资源的过度消耗所导致的生态退化问题, 本文从社会发展、环境污染、资源消耗等方面展开指标遴选：其中社会发展侧重于描述人口增长及经济发展对生态系统所产生的压力, 环境污染侧重于描述农业生产活动对生态系统所造成的污染, 资源消耗侧重于描述区域生产生活对水资源的消耗程度；状态指标主要表示在各种压力状态下, 生态系统当前的结构与功能, 如景观指数、景观类型等所发生的变化；响应指标则主要表示经济系统、社会系统及生态系统所采取的应对生态环境变化的具体措施：其中经济系统侧重描述通过增加第三产业比重从而减少资源需求等措施, 社会系统描述政府及农户采取节水灌溉、转移劳动力等降低资源消耗的措施, 生态系统则根据生态系统服务价值变化反映生态系统的响应状况。同时依据科学性、可操作性、代表性、全面性及可比性等原则, 构建了PSR指标体系(图 2)。

指标无量纲化。由于各个指标量级不同, 因此本文采用极差标准化方法^[19]对指标进行无量纲化处理, 对于正向指标：

(1)

对于负向指标：

(2)

式中, X_ij为指标标准化值；x_ij为第i年第j个指标的值；max x_j与min x_j为第j个指标的最大值和最小值。

确定指标权重。指标赋权是生态安全评价中的重要内容, 指标赋权法包括主观赋权法与客观赋权法, 其中熵值法是一种常用的客观赋权法。熵是一种对不确定性的定量, 信息量越大, 不确定性越小, 熵值也越小。一般而言, 指标离散程度高, 熵值就大, 反之亦然。本文采用熵权法^[21]对指标权重进行赋值, 通过计算各个指标的特征比重、熵值, 进而确定各指标权重：

(3)

(4)

(5)

式中, m为研究期；n为指标个数；S_ij为指标的特征比重；e_j为指标熵值；w_j为指标权重。

计算生态安全指数。生态安全指数表征研究区在压力、状态、响应方面的生态安全水平, 指数越高, 则生态安全水平越高, 本文采用综合指数法计算生态安全指数, 计算公式如下^[18]：

(6)

(7)

式中, ESI_i为P(压力)、S(状态)、R(响应)3个层次的生态安全指数；ESI为生态安全综合指数；X_ij为指标标准化值；w_j为第j项指标权重。

障碍因子识别。障碍度模型可以识别出对生态安全产生障碍的主要影响因子, 障碍度计算公式如下^[6]：

(8)

式中, F_i为指标的障碍度, w_i为指标的权重, X_i为指标的标准化值。

2.1.2 其他评价方法

景观格局指数。景观格局指数是计算生态安全指数的重要指标, 选取景观破碎度和景观多样性指标衡量生态系统状态, 计算公式如下^[22-24]：

(9)

式中, C为景观破碎度, N为斑块数；A为总面积。

(10)

式中, H为景观多样性, P_i为景观类型i所占面积的比例, m为景观类型总数。

生态系统服务价值核算方法。采用效益转移法^[25]评估生态系统服务价值, 该方法通过定义每种土地利用/覆盖类型所对应的单位面积生态系统服务价值, 并将其分别乘以每种土地利用/覆盖的面积来估算生态系统服务价值。由于该方法不是所用到的主要方法, 因此不对其加以展开, 具体方法见谢高地等^[25]的文章(已根据研究区气候、经济等实际情况, 对谢高地等^[25]提出的中国单位面积生态系统服务价值当量进行了区域校正)。

2.2 数据来源 2.2.1 农户调查数据

基于参与式农村评估法, 对黑河中游干流流经区域(张掖市甘州区、临泽县及高台县)进行了半结构式农户问卷调查。首先于2018年7月对3个区县的节水政策进行了初步了解, 并对农户的生计方式及对节水政策的感知进行了初步调查, 在此基础上明确了问卷调查的主体内容。2018年10月开展了张掖市甘州区、临泽县及高台县农户对节水政策的适应与感知调查, 问卷主要通过对比节水政策前后生计策略的转变分析农户对节水政策的适应策略, 并通过农户对节水政策的感知调查评估现行政策的合理性, 为下一步政策的改进提供借鉴。正式调研覆盖3个区、县共计77个行政村(图 3), 取得问卷276份, 其中包含有效问卷259份。通过将有效问卷输入Excel中实现问卷数据电子化, 运用SPSS等软件对数据进行了统计与分析, 遴选出生态安全评价所需指标, 并为生态安全对策与建议的提出提供了基础。

2.2.2 其他数据

气象数据、土壤保持数据来源于中国科学院资源环境科学数据中心；经济数据来源于2001—2015年张掖市统计年鉴、全国农产品成本收益资料汇编等；水利数据主要来源于甘肃省水资源公报；土地利用数据来自地理国情监测云平台；中国单位面积生态系统服务价值当量来自谢高地等^[25]的文章。

3 结果与分析 3.1 黑河中游生态安全评价指标体系构建

水资源是维系黑河中游绿洲生态系统稳定与社会经济发展的命脉, 水分充足则绿洲扩张, 水分短缺则绿洲萎缩, 总之水资源的数量与质量决定了绿洲规模的大小^[26-27]。其次, 土壤是农业活动的基石, 虽然黑河流域拥有广袤的土地面积, 但是土地沙化、土壤质量退化等问题突出, 直接影响了土地生产力水平, 限制了干旱区人类生存空间的拓展^[28]。因此, 黑河中游水土资源的演变对区域生态安全维持有至关重要的作用, 本文即以水土资源为侧重点, 构建了黑河中游生态安全评价指标体系：由于中游人类活动频繁, 因此压力系统主要考虑社会发展、环境污染及资源消耗所产生的压力；状态指标仍侧重描述生态系统当前的结构与功能；响应指标则遴选了能够表征经济-社会-生态复合系统的应对策略的指标, 并在参考相关文献^[18-20]的基础上构建黑河中游生态安全评价指标体系见表 1：

3.2 PSR评价指标权重的确定

根据公式(1)—(5)计算了黑河中游生态安全评价指标的权重。由表 2, 在黑河中游生态安全指标体系中, 高效节水灌溉面积的权重最高, 近年来为保障国家西部生态屏障的安全, 黑河中游实施了大规模高效节水项目, 2001—2015年高效节水灌溉面积增长了55.17%, 因此该指标对生态安全的影响程度最高；其次, 生态用水量的权重为0.0617, 在生态脆弱的黑河中游地区, 生态用水量的多少与生态系统的稳定性高度相关, 因此生态用水量的权重较高；此外, 地下水开采率与第三产业比重的权重均高于0.60, 表明其对区域生态安全的影响程度也较高。

3.3 生态安全评价

根据公式(6)及表 2计算了黑河中游生态安全压力、状态、响应指数(图 4), 结合相关文献^{[18, 29]}与研究区实际情况, 建立生态安全等级划分标准(表 3)。结果表明：2001—2015年黑河中游压力指数呈现波动下降趋势：2001年压力指数为0.2615, 2008年降为0.1644, 2015年进一步降为0.0569。该指数越来越小, 表明生态系统所面临压力越来越大, 压力指数对生态安全的贡献率越来越低；2001—2015年状态指数呈现先增长后减少的趋势：2001年状态指数为0.1082, 2008年增至0.1313, 2015年降至0.1236；2001—2015年响应指数越来越高：2001年响应指数为0.0051, 2008年增至0.1590, 2015年高达0.4083, 响应指数增长表明经济-社会-生态复合系统采取了相应策略应对区域生态环境问题。

根据公式(7)计算了黑河中游生态安全综合指数(图 4), 结果表明2001年生态安全综合指数为0.3748, 压力指数对综合指数贡献率最高, 表明此时生态系统面临压力最小；2008年生态安全综合指数增至0.4548, 压力、状态与响应指数的贡献率均为30%左右, 表明此时生态系统面临压力有所增长, 同时开始采取积极措施提高生态安全水平；2015年生态安全综合指数进一步增长至0.5888, 响应指数贡献率最高, 表明经济-社会-生态复合系统采取了一系列措施显著提高生态安全水平。根据表 3, 2001年生态安全等级为“较不安全”, 2008年生态安全等级升为“临界安全”, 虽然2015年生态安全综合指数增加, 但生态安全等级仍处于“临界安全”状态。

3.4 障碍因子诊断

根据公式(8)计算了各个指标对生态安全的障碍度, 列出了2001、2008、2015年障碍度最高的5个因子(表 4), 结果表明2001年影响生态安全的障碍因子主要是响应层指标, 高效节水灌溉面积是最主要的障碍因子, 其次是生态用水量；2008年影响生态安全的障碍因子主要是响应层和压力层指标, 高效节水灌溉面积仍然是障碍度最高的因子；2015年影响生态安全的障碍因子主要是压力层指标, 随着人类生产生活对水资源的需求增长, 地下水开采率的障碍度最高, 化肥施用量也对生态安全产生了显著影响。以上表明随着响应水平的提高, 2001—2015年障碍因子逐渐由响应层指标转为压力层指标, 降低水资源消耗量和环境污染水平成为提高生态安全水平的重要途径。

4 讨论与建议

基于生态系统服务与生态安全格局间的内在联系, 根据人类社会与生态系统间以生态系统服务流动为纽带形成的反馈机制, 构建了基于生态系统服务理论的PSR模型。一方面为生态安全评价提供综合的量化指标体系提供了逻辑基础, 拓展了人们在生态安全评估方法领域的研究视野；另一方面该模型耦合了人类活动与区域生态安全, 利于深入理解黑河中游生态安全维持机制, 拓展人们对生态安全评价、生态安全维持与提高的认识。

黑河中游是流域生产生活的主要承载区, 具有较高的农业活动强度。近年来生态系统面临的环境污染和资源消耗问题越来越大：农业生产过程中化肥、农药、薄膜施用量显著增长, 2001—2015年化肥施用量增长22.84%；此外, 生产生活对水资源的需求量越来越高, 2001—2015年地下水开采量增长一倍有余, 随着生态系统面临的压力越来越大, 压力指数持续降低。在人类活动干扰下, 中游绿洲生态系统结构和功能略有降低, 植被退化、景观破碎度增加, 虽然在湿地保护工程实施背景下, 水域面积有所增长, 但状态指数仍呈现微弱的降低趋势。与此同时, 经济-社会-生态复合系统采取了相应策略应对生态环境问题：随着绿洲经济发展及技术水平的提高, 政府对高效节水灌溉工程投入越来越大, 2001—2015年高效节水灌溉面积增长了55.17%, 在一定程度上改善了缺水状况。2001—2015年GDP增长四倍有余, 经济发展促进了产业和劳动力转移, 湿地保护工程的实施也促进生态系统服务价值增长, 响应指数呈现增长趋势, 是生态安全水平提高的主要原因。

总体来看, 水资源是制约黑河中游绿洲发展的命脉, 水资源丰富与否决定了绿洲规模的大小。与水资源利用相关的节水灌溉面积、生态需水量、地下水开采率等均是影响区域生态安全程度较高的因子。维持区域水资源供应平衡, 实现水资源可持续利用是提高该区域生态安全水平的核心。在较短时间尺度上, 水资源供应格局呈现稳定趋势, 因此, 应从降低水资源消耗量、提高水资源利用效率角度入手提高区域生态安全水平。基于生态安全评价结果, 结合农户调查问卷, 提出了以下提高生态安全水平的建议：

(1) 坚持生态保护优先。黑河中游绿洲是全国重要的商品粮生产基地, 2001—2015年黑河中游耕地面积增长20.40%。荒地开垦、乱砍滥伐导致植被退化, 2001—2015年林地、草地面积分别减少19.97%、3.68%。此外, 耕地扩张及大水漫灌加剧了水资源消耗, 2001—2015年地下水开采量增加, 超越了生态系统的承载力, 影响了区域生态安全。因此, 为实现黑河中游生态经济可持续发展, 必须坚持生态保护优先原则, 恢复与重建黑河中游退化生态系统, 保障天然绿洲与人工绿洲协调发展。

(2) 调整种植和产业结构。2015年黑河中游绿洲农业用水量占比高达90%以上, 但农业产值占比不足30%, 水资源投入产出比例较低。因此, 首先需降低高耗水作物比例, 在发展制种玉米的基础上, 进一步发展大田反季蔬菜、花卉和酿酒葡萄种植等, 构建“节水型与高效益”的新型农业体系。问卷调查结果显示虽然2001—2015年经济作物种植比例有所提升, 但农户认为温室成本较高, 提高相关补贴标准、加大扶持力度有助于提高农户调整种植结构的积极性。此外, 应进一步推进产业结构优化, 发展生态产业, 提高水资源利用效率。

(3) 发展高效节水灌溉技术。2001—2015年张掖市微灌和低压管灌面积分别增长34767 hm²和43100 hm², 是生态安全水平增长的重要原因, 但仍存在覆盖率低、节水工程建设标准低等问题。问卷调查结果显示, 高达40.96%的农户认为滴灌等节水设施灌溉效果差, 导致粮食减产, 收入降低。而且农户认为滴灌水量过小, 耕地表层土壤干旱, 加重了风沙危害, 以上因素影响了农户参与节水灌溉的积极性。因此, 发展适于研究区的高效节水技术, 实现节水、增收双赢才能真正提高节水灌溉效益, 解决水资源短缺问题。

(4) 提高公众生态保护意识。问卷调查结果显示69.65%的农户认为研究区水资源不短缺, 84.05%的农户认为水资源污染不严重, 76.34%的农户认为经济发展比生态保护更重要, 表明农户对水资源短缺、农业面源污染了解较少, 生态保护意识淡薄。虽然86.82%的农户注意节约用水, 但主要原因是农户认为水资源费用定价高, 农户倾向通过节约用水减少开支, 而非自发节水行为。因此, 亟需加强生态保护宣传, 提高公众水资源节约和保护意识, 实现生态经济可持续发展。

本文运用熵值法对PSR模型的指标权重进行计算, 该方法根据各指标所提供信息量大小来确定权重, 克服了主观赋值法存在的评价结果不稳定的现象, 在一定程度上提高了评价质量^{[20, 30]}。但是PSR模型在指标遴选方面存在一定程度的重叠问题：以评价大城市生态安全为例, 压力层常选指标为表征生态系统服务获取强度的“总人口”, 而状态层常选指标为表征生态系统结构的“林地面积”。实际上人口数量增长越快, 城市扩张对林地的侵占程度也越高, 因此“总人口”与“林地面积”存在重复性, 导致生态安全评价结果出现偏差, 降低了将评价结果应用于生态环境管理中的可行性。因此, 如何科学遴选独立的PSR指标体系将是下一步生态安全评价研究的重点。

5 结论

本文基于PSR模型, 耦合生态系统服务理论, 构建了黑河中游生态安全评价指标体系：其中压力指标主要表征生态环境出现问题的原因, 即外界对生态系统所造成压力及人类对生态系统服务的获取强度；状态指标主要表征在各种压力状态下, 生态系统当前的结构与功能所发生的变化；响应指标则主要表征经济系统、社会系统及生态系统所采取的应对生态环境变化的具体措施。

黑河中游生态安全评价结果表明：2001—2015年化肥施用量与灌溉用水增加导致生态系统所面临压力越来越大, 但由于经济-社会-生态复合系统采取了相应策略应对区域生态环境问题, 黑河中游绿洲生态安全水平略有提高, 生态安全综合指数由0.3748增至0.5888, 生态安全等级由“较不安全”达到“临界安全”状态。

虽然黑河中游生态安全呈现波动增长趋势, 但2015年生态安全等级仍处于“临界安全”状态, 且稳定性较差, 因此黑河中游生态安全水平有待进一步提高。在以上研究的基础上, 结合问卷调查结果, 围绕干旱区绿洲发展命脉—水资源的可持续利用提出了坚持生态保护优先、调整种植和产业结构、发展高效节水灌溉技术、提高公众生态保护意识等提高研究区生态安全水平的相关对策与建议。