面向生态系统服务供需的开都-孔雀河流域生态安全格局研究

2023-06-08闫豫疆李建贵李均力

干旱区研究 2023年5期

闫豫疆，李建贵，李均力，蒋腾

（1.新疆农业大学经济管理学院，新疆乌鲁木齐 830052；2.新疆农业大学林学与风景园林学院，新疆乌鲁木齐 830052；3.中国科学院新疆生态与地理研究所，荒漠与绿洲生态国家重点实验室，新疆乌鲁木齐 830011；4.新疆林业科学院，新疆乌鲁木齐 830018）

生态系统是人类生存的基础条件，在社会经济快速发展和城市化进程不断推进的条件下，由生态系统服务供需失衡造成的生态安全问题日益突出。生态系统空间格局通过影响生态过程进而影响生态系统服务。为此，通过构建区域生态安全格局，优化生态网络，提高景观的连通性，实现对区域生态的有效调控和生态系统服务的供需均衡［1］。西北干旱内陆流域水资源缺乏且分布不均，生态安全脆弱，水土资源的过度开发加剧了生态安全风险，一定程度上限制了生态环境与经济社会的可持续发展。近年来，国家相继出台一系列的山水林田湖草沙冰生态保护修复工程，如孔雀河下游生态输水，就是为了恢复下游生态系统安全屏障。因此，从生态系统服务供需的角度构建流域生态系统安全格局，对提升和优化生态系统服务具有重要作用。

国外学者对生态安全格局的优化研究可追溯到20 世纪初［2］，主要围绕生态网络、绿色基础设施等方面展开，其技术体系相对成熟［1］。国内学者俞孔坚［3］最早在1990 年提出了景观生态安全格局概念，近年来相关学者逐步形成了“源地识别-阻力面构建-廊道提取-构建安全格局”的基本范式［4］。区域生态安全格局的构建方法有最小累积阻力模型、电路模型［5］和斑块重力模型，最小累积阻力模型应用最为广泛，其能较好地反映景观格局变化和生态过程演变的相互作用关系，源地识别主要从生境重要性［6］、生态敏感性［7］、景观连通性［8］等判定，而从生态系统服务供给和需求空间之间的联系考量较少［9］。从研究区域来看，国内大多学者多集中在经济发达、生态环境质量良好的区域，而对经济欠发达及生态环境脆弱的西北生态功能区的相关研究关注较少［9-10］，尤其是以干旱区独立流域为研究区域的较少。

开都-孔雀河流域位于新疆塔里木盆地北缘，流域内有我国最大的内陆淡水湖博斯腾湖，湖泊生态系统安全与否对中下游生态及社会经济生活有重要影响［11-13］。流域以博斯腾湖为节点，上游的开都河流域水源充足、人类活动较少，生态系统处于较好的状态，而下游的孔雀河流域由于水土资源开发活动强烈，造成湿地萎缩、草地退化等生态安全问题。从全流域角度评估生态系统服务供需，构建和优化生态安全格局，可在有限水资源条件下实现对流域的整体性保护。因此，本文选择食物供给、生境质量、碳固持和防风固沙4 项生态系统服务指标，采用热点分析法识别重要的生态供给源地，以人口密度、地均GDP、夜间灯光亮度和土地利用强度综合指标，为流域需求源地，利用最小累积阻力模型构建研究区生态安全格局，结合生态环境特征，提出开都-孔雀河流域生态安全格局的优化方案，为保障区域生态安全提供参考。

1 数据与方法

1.1 研究区概况

开都-孔雀河流域位于新疆巴音郭楞蒙古自治州境内，地处塔里木盆地东北部、塔克拉玛干沙漠东北缘，流域包括和静、和硕、焉耆、博湖、库尔勒市、尉犁和若羌等地，总面积约9.35×104km2（图1）。流域水源补给主要来自于开都河上游高山地带的降水和冰雪融水，流经焉耆盆地、库尔勒和尉犁县，是绿洲农业生产、城市发展及居民生活的重要水源供给，“山水林田湖草沙冰”等生态要素构成一个典型的内陆干旱生态系统。流域属于温带大陆性干旱气候，日照时间长、降水量小、潜在蒸散量大，且流域上游西北方向地势高，下游东南方向地势低，气象差异明显，年降水量从上游的300～500 mm 逐渐下降为下游的50～100 mm［14-16］。流域下游特别是孔雀河中下游地区的植被覆盖率较低，生态环境极为脆弱。

图1 开都-孔雀河流域示意图Fig.1 Sketch map of the Kaidu-Kongque River Basin

1.2 数据来源

开都-孔雀河流域2020年土地利用数据来源于地理空间数据云平台（http://www.gscloud.cn/）的GF-1 WFV 多光谱数据，空间分辨率为16 m，对其进行辐射校正、几何校正、波段融合、图像镶嵌、裁剪、分类提取等处理。DEM 数据来源于地理空间数据云（http://www.gscloud.cn/），空间分辨率为30 m；夜间灯光数据来源于资源环境科学与数据中心（http://www.resdc.cn）的中国长时间序列逐年人造夜间灯光数据集（1984—2020）；植被归一化植被指数来源于数据共享服务系统（https://data.casearth.cn/）；植被净初级生产力来源于地理监测云平台（http://www.dsac.cn/）；碳排放数据来源于中国碳核算数据库（https://www.ceads.net.cn/）；国内生产总值（GDP）数据来源于资源环境科学与数据中心（http://www.res⁃dc.cn/DOI）的中国GDP 空间分布公里网格数据集；人口空间分布公里网格数据集和全国第七次人口普查数据从WordPop（https://hub.worldpop.org/）和资源环境科学与数据中心获取，以及土壤数据来源于国家地球系统科学数据中心（http://www.geodata.cn）的中国土壤数据集（v1.1）（2009）；其他社会经济数据来源于新疆统计年鉴、巴音郭楞蒙古自治州和相关县市的统计年鉴。为保证数据的连续性和精准性，将数据坐标系统一为WGS-1984-UTM-zone-45，并将所有空间数据重采样为1 km×1 km分辨率的栅格图层。

1.3 研究方法

本文选取食物供给、碳固持、生境质量和防风固沙4项生态系统服务，从供给、需求两个层面对各项生态系统服务进行定量分析，借助ArcGIS中Spa⁃tial Statistics工具对研究区进行热点分析，从而确定生态供给源地和需求源地。运用最小累积阻力模型，建立生态阻力面，根据成本距离和路径分析，生成最小累积成本路径，并选择生态廊道之间交汇处的生态功能区作为生态节点。在此基础上，按照“识别生态源-构建阻力面-提取生态廊道”模式［17］，综合生态源地、廊道和节点等要素构建开都-孔雀河流域生态安全格局。

1.3.1 生态源地的识别生态源地是指区域中生态系统服务供给/需求量较大或供给种类较多的斑块，具有一定空间扩展性和连续性景观组分［2］。本文针对研究区域生态本底状况，选择食物供给、碳固持、生境质量和防风固沙4项生态系统服务指标进行定量评估，选取单位面积生态系统服务供给/需求量大于流域平均值的斑块，作为生态系统服务的重要或重点区域，再通过热点分析，选择空间上连续、斑块面积大于50 km2的区域，作为流域内生态系统服务供给/需求源地。

（1）生态系统服务供给量评估

食物供给服务的计算按区域归一化植被指数（Normalized Difference Vegetation Index，NDVI）比值法，得到各栅格单元食物供给能力［1］；碳固持和生境质量服务采用当前应用较为广泛的InVEST 3.10.2模型，其中的碳固持模块（Carbon Storage and Se⁃questration）和生境质量模块（Habitat Quality）分别用于评估生态系统供给服务［10］；利用修正的土壤风蚀方程（Revised Wind Erosion Equation，RWEQ）估算研究区防风固沙服务供给量［11］。具体计算方法如表1所示。（2）生态系统服务需求量评估

表1 生态系统服务供给量计算方法Tab.1 Calculation of ecosystem service provisioning

生态系统服务需求是指在一定时空范围内被人类实际利用的自然资源和服务，需求量的大小反映人类社会对自然生态环境的作用强度，以及对生态系统服务的需求程度［4］。本文选取人口密度、地均GDP、夜间灯光指数和土地利用强度4 项指标反映研究区生态系统服务的需求程度。由于人口、GDP和夜间灯光指数的单位不统一，在空间分布上存在显著差异，本文首先采用离差标准化方法计算各生态系统服务需求量，再利用自然对数法计算总体生态系统服务需求量，进而消除由单位和数据量的差异造成的局部剧烈波动［12］。计算公式如下：

式中：x为标准化后的生态系统服务需求量；xi为第i个单元的生态系统服务需求量；xmax为流域最大值，xmin为流域最小值；Yi表示研究单元i对生态系统服务的总需求；x1、x2、x3和x4分别表示土地利用强度、人口密度（人·km-2）、地均GDP（万元·km-2）和平均夜间灯光指数（Average Nighttime Light Index，ANLI）。

（3）热点分析

利用ArcGIS软件的热点分析工具，识别研究区生态系统服务在空间中的高值区（热点区域）与低值区（冷点区域）的集聚程度，得到具有显著统计学意义的正的Gi*（或Z）数值，Z值得分越高，高值（热点）的空间聚类就越明显［18］。将热点分析运用到生态系统服务研究中，其输出结果呈现表面连续性，表征景观的连通性，便于区域生态系统更好的管理［19］。其计算公式如下：

式中：wij为格网i与格网j之间的空间权重矩阵；xj为格网j的属性值；n为总格网数。其中：

在我国西北干旱地区，水资源作为非常重要的约束因子，很大程度上决定了生态环境和土地覆被类型，因此，本文将该地区的湖泊和湿地作为供给源地［20］。

1.3.2 生态阻力面构建累积阻力模型主要描述“源地”生态系统服务流运行的空间态势［8］，开都-孔雀河流域地貌复杂，从上游到下游分别为山间盆地、高山峡谷、焉耆盆地和塔里木沉降盆地，地形差异较大以及土地覆被类型的变化都是生态源地向外扩散的阻力来源［21-22］。干旱区绿洲生态系统受自然和人为因素的影响都比较大，选择景观类型、高程、坡度和地均GDP作为阻力因子［23-25］，其阻力分量及其系数设置见表2所示。生态源地与目标单元之间的最小累积阻力计算公式如下：

表2 景观累计阻力分量及其阻力系数设置Tab.2 Cumulative resistance components of the landscape and their resistance factor settings

式中：MCR 为最小累积阻力；Dij为物质流从源地j到源地i的欧氏距离；Ri为空间某一景观单元i的阻力值；Σ 表示源地i与源地j之间穿越所有单元阻力的累积；fmin为景观单元对不同的源取累积阻力最小值；m和n分别为供给源地和需求源地总数。

1.3.3 生态廊道和生态节点识别

（1）生态廊道识别

生态廊道是发挥景观生态系统服务的重要途径，其连通性是衡量廊道结构的关键指标［26-27］。根据上述计算的生态阻力面，利用GIS中的Distance分析工具，以生态供给源地为中心，基于路径成本分析工具生成成本回溯链和距离，进而得到生态供给源地到需求源地的最小阻力路径，即生态廊道。本文考虑到干旱区水资源的重要性，将河流、水系作为流域的关键廊道［28-29］。

（2）重力模型

重力模型用于定量计算源与目标间相互作用力大小，利用力越大，则代表廊道越重要，计算公式如下：

式中：Gij是供给源地i和需求源地j的相互作用，其大小反映供给和需求源地之间潜在生态廊道的重要程度；Ni和Nj是i和j两个斑块的权重值；Dij是i和j两个斑块之间潜在廊道阻力的标准化值；Pi为斑块i阻力值；Pj为斑块j阻力值；Si是斑块i的面积；Sj是斑块j的面积；Lij是i和j两个斑块之间的廊道阻力值；Lmax是所有廊道累计阻力的最大值。

（3）生态节点识别

生态节点是区域生态体系或景观格局中重要的区域，对区域生态安全起着控制点的作用。本文采用生态廊道的交点作为生态节点，将关键廊道与潜在廊道间的交点作为一级生态节点，潜在廊道间的交点作为二级生态节点［29］。

2 结果与分析

2.1 生态系统服务供给源地空间格局

通过对2020 年开都-孔雀河流域食物供给、碳固持、生境质量和防风固沙生态服务供给量进行评估（图2），食物供给高的区域主要分布在中游耕地较为密集的区域，其次是开都河上游的草原；碳固持服务与食物供给服务类似，且植被覆盖度越高，碳固持服务能力越强；生境质量服务供给指数较高的区域主要为大面积草地或水域，主要分布在开都河上游草原区与博斯腾湖湖区，而建设用地、耕地外围的生境质量较差；防风固沙服务供给高值区主要分布于地形起伏较大的丘陵山区，结合气象因子、土壤质地和植被覆盖度综合分析可知，开都河上游山区降雨丰富，植被覆盖度高，防风固沙量功能较高，但也存在部分草场退化、沙化等功能退化的现象。

图2 开都-孔雀河流域生态系统供给服务空间分布Fig.2 Spatial distribution of ecosystem supply services in Kaidu-Kongque River Basin

从研究区生态系统供给服务的冷热点分布可知（图2e），热点区同样分布于开都河上游草原区及中部耕作区。去除面积小于50 km2的斑块，最终确定开都-孔雀河流域生态供给源地14 个，总面积为20068.4 km2，占研究区面积的21.5%（图2f）。

2.2 流域生态系统服务需求源地空间特征

开都-孔雀河流域人口分布、夜间灯光指数、地均GDP 和土地利用强度空间分异显著（图3），前三者的高值区域主要分布在各县市城区以及各县政府驻地所在乡镇，这些地区整体生态压力较高，对生态系统服务需求较大，综合各项需求指标，去除面积较小的斑块，得到全流域生态需求的空间格局（图3e）。识别生态需求源地斑块共9 个，合计面积为4327.5 km2，占流域总面积4.6%。总体而言，流域生态系统服务需求集中分布于中游区域，高值聚集区为焉耆盆地和库尔勒市-尉犁县。生态系统服务低需求区域主要分布在流域的上游山区和下游荒漠区，该区域地理位置较为偏僻，人类活动干扰较小，故而对生态系统服务需求相对较小。

图3 开都-孔雀河流域生态系统需求服务空间分布Fig.3 Spatial distribution of ecosystem demand services in Kaidu-Kongque River Basin

2.3 生态廊道和生态节点构建

2.3.1 阻力面构建生态系统能量流动和物质循环遇到的阻力值越小，物质和能量在景观单元之间的流动就越容易，所需提供的生态系统服务就越多［1］。流域生态阻力低值区主要分布于上游巴音布鲁克国家级自然保护区、中游博斯腾湖大小湖区、焉耆盆地绿洲农业区和库尔勒-尉犁农业区；高值区主要分布在植被覆盖少的山区以及受人类影响大的城镇地区（图4）。

图4 开都-孔雀河流域生态扩张累计阻力面Fig.4 Cumulative resistance surface of ecological expansion in the Kaidu-Kongque River Basin

2.3.2 生态廊道识别利用最小阻力模型计算研究区生态系统服务供给源地与需求源地之间的最小成本路径，共得到生态廊道126条，然后利用重力模型计算出各生态廊道的重力矩阵（表3），相互作用力强度大于1 作为重要廊道，共计17 条，总长度654.7 km。生态廊道在开都河流域有两条轴线，一是从西北向东南延伸的轴线，二是从北部山区向南延伸的轴线，在焉耆盆地较为错综复杂；孔雀河流域生态廊道相对集中，分布在库尔勒-尉犁县区域。

表3 生态廊道重力值统计Tab.3 Statistical of gravity values for ecological corridors

2.3.3 生态节点识别生态节点实际上是处于生态廊道沿线的特定位置，相对生态源地面积较小的生态斑块［23］，通过生态节点提取，开都-孔雀河流域有生态节点65 个，其中，重点生态节点24 个，如图5所示。

图5 开都-孔雀河流域生态廊道和生态节点分布Fig.5 Distribution of ecological corridors and ecological nodes in the Kaidu-Kongque River Basin

2.4 生态安全格局构建

开都-孔雀河流域具有完整的“山水林田湖草沙冰”生态系统结构，且空间分布差异较明显。在原有生态廊道结构的基础上，结合水资源分布以自然和人工廊道构建起的研究区生态安全格局，以生态功能区划为基础背景，综合考虑生态基底、自然保护区位置及各生态安全要素分布，构建了“两核心、两片区、三横四纵多节点”的流域生态安全格局（图6）。

图6 开都-孔雀河流域“源地-廊道-节点”生态网络格局Fig.6 Source-corridor-node ecological network pattern in the Kaidu-Kongque River Basin

（1）两核心、两片区是核心生态保育区和重点农业发展区，其中，两个核心生态保育区是巴音布鲁克国家级自然保护区、巩乃斯国家森林公园和博斯腾湖国家湿地公园，这些区域保障着整个流域的水资源，丰富景观类型，提高生态系统完整性，其生态保护修复策略主要是增加生态廊道，扩大自然或半自然生态绿地比例，增强生境间的连通度。重点农业发展区包括焉耆盆地农业发展区和库尔勒市-尉犁县农业发展区，该区域一方面应注重农业产业结构调整以提高生态效益和提升生态系统服务；二是应充分利用闲置土地，深度挖掘农村建设用地潜力，提升闲置土地的综合效益。

（2）三横四纵多节点，三横为开都河流域生态保护与修复发展廊道和北侧重水源涵养与生物多样性保护修复廊道，以及铁门关至库尔勒市农业发展廊道；四纵为尤尔都斯盆地生物多样性保护区重点生态廊道、黄水沟重点生态廊道、孔雀河重点生态廊道以及北部源地间流通缓冲廊道；“多节点”位于生态系统服务需求和供给交汇处，以及生态廊道分布密集的焉耆盆地、库尔勒市。在构建“绿色斑块、生态廊道、湿地、农田”镶嵌的绿地生态网络系统，突出绿地网络的社会-生态复合功能，利用人工和自然生态系统，促进实施山水路林村综合治理，进而实现生态系统的良性循环，改善区域内生物栖息地环境。

3 讨论

开都-孔雀河流域是“山水林田湖草沙冰”一体化保护与系统治理的典型区域，研究该区域的生态系统服务的供需及生态状况对维护流域生态安全、促进流域经济社会可持续发展都具有重要意义。本文从生态系统服务供给和需求入手，识别流域生态系统供给和需求空间，依托生态源地、廊道和节点构建了流域生态安全格局，选择了供给/需求量大于50%的区域作为生态服务源地的重要区域，从各项服务的空间分布来看，其结果满足精度要求，但未考虑不同生态系统服务对不同区域的重要程度，忽略了区域生态系统服务的主体功能。

与以往的生态安全格局构建研究相比［18,29］，本文更注重以整个流域为研究对象，生态系统服务的供给和需求在相对独立的空间内形成闭环，更有利于供需关系的研究。开都-孔雀河流域中下游区域经济发展和人口数量较为聚集，存在高供给高需求和低供给高需求的状态，生活和生产空间重合，人工干预对生态系统服务的流通影响较大。基于此，本文构建了“两核、两片区、三横四纵多节点”多要素的生态安全网络结构，并提出相应的保护和优化生态系统的策略，其难点在于打通网络的连通性。开都-孔雀河流域道路、人工建筑、围栏等人工要素往往是生态廊道的主要阻断物，而由于缺少相应的数据未将其纳入生态网络结构体系建设，虽然研究结果具有一定的可操作性，但缺乏对小尺度的详细探讨，今后有待更深一步的研究。