20世纪90年代以来，随着城镇化的快速发展，诸如水土流失、草场退化、土壤盐渍化、土地荒漠化、森林资源危机、水资源短缺、生物多样性减少等各项生态问题也纷纷出现[1-3]，生态系统面临越来越严重的威胁，生态安全问题逐渐引起人们的重视，各界学者从不同的视角对生态安全展开研究[4-5]。其中，我国著名学者俞孔坚教授提出的“景观安全格局理论”为生态安全的研究做出了重要贡献。俞孔坚教授提出了通过“确定源—建立阻力面—根据阻力面判别安全格局”的步骤构建景观生态安全格局[6]，这一方法已被众多学者引用来研究景观生态安全格局[7-10]。源地是现存的乡土物种的栖息地，是物种扩散和维持的源点，具有内部同质性和向四周扩张或向“源”本身汇集的能力，是历史时期生态系统相对稳定的区域[6,11]，是构建景观生态安全格局的基础。因此，对于源地的识别，既应考虑斑块空间上的生态质量优良性，又应考虑斑块时间上的生态质量稳定性。目前的研究，多是从空间尺度进行评价，选取生态质量优良的斑块作为构建景观生态安全格局的源地，而从斑块动态变化的稳定性方面进行评价的较少。本研究从空间和时间两个尺度综合考虑，选取生态质量优良且变化稳定的斑块，作为景观生态安全格局构建的源地。 北戴河新区是我国首批国家智慧城市试点之一，是渤海湾著名的黄金旅游城市，其快速发展的旅游产业为当地经济发展做出了巨大贡献，然而，随之出现的一些生态环境问题阻滞了北戴河新区生态旅游的发展。因此，构建景观生态安全格局，完善生态旅游规划，意义重大。本研究以北戴河新区为例，在借鉴前人研究成果的基础上，尝试从空间和时间两个尺度综合考虑，对景观斑块质量进行评价，探讨景观生态安全格局源地的识别方法，为景观生态安全格局的合理构建奠定基础。 1 研究区概况与数据来源 1.1 研究区概况 北戴河新区位于秦皇岛市域滨海地区的西侧，在119°1′—119°24′E、39°25′—39°50′N，总面积514.79 km2。新区地势平坦，海拔较低，年平均气温11℃，属暖温带半湿润大陆性季风气候，降水丰沛，雨热同季。新区产业发展以旅游度假和农业生产为主，有连片的耕地和少量的园地等农业资源，有诸如沙滩海水浴场、海洋大漠、泻湖湿地和连绵林地等丰富的旅游资源，为发展生态农业和生态旅游产业提供了得天独厚的条件。然而，近些年来，人类活动程度的加剧，对北戴河新区的生态环境造成了一定破坏，耕地、水域等生态用地面积不断减少。为维持生态系统稳定运行，为北戴河新区提供一个良好的旅游业发展环境，必须采取措施保护生态用地。 1.2 数据来源 研究所需不同时期土地利用类型分布图由遥感影像解译而来。研究选取了1991年9月的landsat 5 TM影像、2001年9月的landsat 5 TM影像和2013年10月的landsat8 OLI_TIRS影像。此三期遥感影像云量少、成像质量好，有利于解译的精确性。此外，还获取了Google Earth 历史图像和北戴河新区2013年土地利用现状矢量数据，用于解译标志的建立。其他社会经济数据来自于北戴河新区统计年鉴和北戴河新区整治规划资料。 遥感影像的解译借助ENVI 4.7软件，通过波段叠加、几何校正、图像融合、图像镶嵌、图像裁剪和图像增强等手段对3期遥感影像进行预处理。然后，进行监督分类，根据研究所需数据类别，将研究区土地分为耕地、园地、林地、草地、水体、建设用地和未利用地。最后，对分类结果进行修正，并运用混淆矩阵和Kappa系数分别对1991年、2001年和2013年三期监督分类后的影像进行解译精度检验。经检验，三期研究区遥感影像分类结果的Kappa系数分别为0.862 6，0.863 8，0.870 3，解译精度均在80%以上，能够满足研究需要。将栅格图像转换为可以在ArcGIS 10.0软件中编辑的矢量数据，参考前人对生态用地概念和分类的研究[12-13]，提取耕地、园地、林地、草地、水体作为生态用地，进行源地的识别研究。 2 研究思路与方法 源地的选取通过景观斑块质量的评价来实现。景观斑块质量的优劣既表现在斑块自身的特征属性方面，又表现在斑块对外界环境的响应方面[14]。本研究从这两个方面构建景观斑块质量评价体系，采用层次分析法确定权重，建立评价模型，分别对北戴河新区1991年、2001年和2013年的生态用地景观斑块质量进行评价，并将结果划分为五个级别，定义级别数值越小，景观斑块质量越好。在此基础上，叠加三期景观斑块质量评价结果，选取三期都为1级和2级的斑块作为2013年景观生态安全格局的源地。 2.1 评价体系构建 研究基于维持景观过程完整性、维持生物多样性和维持生态系统稳定性的目的，从斑块自身特征属性和斑块对外界环境的响应两方面选取6项指标，运用层次分析法确定指标权重，构建景观斑块质量评价体系(表1)。 2.2 评价指标计算 本研究以斑块为评价单元，计算各指标分值，并采用自然间断点分级法将各指标分值划分为1—5五个级别数值，定义级别数值越小，斑块质量越高，以此实现指标的无量纲化。各指标的计算方法如下： (1) 斑块面积、分维数、边界破碎度的计算。斑块面积、分维数和边界破碎度的计算，借助Fragstats 3.4景观分析软件进行，并运用ArcGIS 10.0软件将计算结果转换为栅格格式。 表1 北戴河新区景观斑块质量评价指标体系与权重目标层准则层指标层指标层对目标层的综合权重指标生态学意义景观斑块质量评价斑块自身特征属性斑块面积分维数边界破碎度0..影响生物活动的扩散和生物多样性的保护,属正向指标反映斑块形状的复杂程度和人类活动对斑块的干扰程度,属负向指标景观破碎化程度的直接反映,属负向指标斑块对外界的响应斑块重要值生境质量生态服务价值0...2500斑块在维持景观连通性上的重要程度,属正向指标生态系统对人类活动的适宜程度,属正向指标生态系统服务的效用价值,是维持生物多样性能力的反映,属正向指标 (2) 斑块重要值的计算。斑块重要值通过景观连通性指数来实现，本研究选取可能连通性指数PC来进行斑块重要值dPC的衡量。可能连通性指数是一种可能性模型，反映斑块之间连通性的可能性与斑块之间的距离有关[15]。计算公式如下： 斑块的重要值指在此点处断裂(或者移除后)整个区域景观连通性的变化量，斑块重要值的计算公式如下： 式中：dPC为斑块重要值；I表示景观中所有斑块的整体指数值；Iremove是去除某单个斑块后剩余斑块的整体指数值。 斑块重要值的计算借助ArcGIS 10.0软件、Conefor Inputs for ArcGIS 10插件模块和Conefor Sensinode 2.2软件共同完成。先计算斑块节点数据和各斑块之间的距离数据，再计算斑块重要值，并将结果转换为栅格格式。 (3) 生境质量的计算。生境质量的评估借鉴由美国斯坦福大学、世界自然基金会和大自然保护协会联合开发的生态系统服务功能评估工具InVEST模型来实现。InVEST模型通过分析生境斑块受周围环境的威胁程度来实现生境质量的评估[16]，其执行需要土地覆盖图、各威胁因子的影响程度及最大影响距离、土地覆盖类型对各威胁因子的敏感度数据以及土地能够受到合法保护的程度数据。土地覆盖图来自遥感解译的30×30 m土地利用类型分布栅格图，耕地、园地、林地、草地和水体为生境数据，生境适宜度(Habitat)赋值为1，其他土地为非生境数据，生境适宜度(Habitat)赋值为0。威胁源包括高速公路、铁路、主要公路、其他建设用地和工矿用地五种，通过在土地利用类型分布图中提取得到。各威胁因子的影响程度和最大影响距离及土地覆盖类型对各威胁因子的敏感度数据参考相关文献[14,17-19]并结合研究区实际得到(表2，表3)。考虑到土地合法保护程度数据的难获取性，本研究不将此指标考虑在生境质量的评估指标体系中。在ArcGIS 10.0软件中添加InVEST 2.2.0工具，导入相应数据，评估斑块生境质量。 表2 北戴河新区生境威胁因子影响程度与最大影响距离名称影响程度最大影响距离/km高速公路0.50.75铁路0.40.50主要公路0.60.50其他建设用地0.72.00工矿用地1.03.00 表3 北戴河新区生境类型对各威胁因子的敏感度威胁因子生境类型高速公路铁路主要公路其他建设用地工矿用地耕地园地林地草地水体 (4) 生态服务价值的计算。生态服务价值的评估借鉴谢高地等人制定的中国生态系统单位面积生态服务价值当量表来实现[20]。依据此表，得到耕地、园地、林地、草地和水体的生态服务价值，其中，园地的生态服务价值取林地与草地的平均值(表4)。将结果转换为栅格格式。 表4 北戴河新区不同生态用地类型生态服务价值用地类型耕地园地林地草地水体生物多样性服务价值 2.3 评价模型构建与源地识别 综合评价在权重确定和指标无量纲化的基础上，按加权模型来运算(公式3)，结果为五个级别，级别数值越小，斑块质量越好。 式中：Fi为第i个评价单元的质量总分值；pij为第i个评价单元第j个指标的级别数值；wij为第i个评价单元第j个指标的综合权重。 叠加1991年、2001年和2013年三个时期的景观斑块质量综合评价结果，选取三期都为1级和2级的斑块，作为2013年景观生态安全格局的源地。 3 结果与分析 根据上述研究方法，分别计算1991年、2001年和2013年的各评价指标分值，并采用自然间断点分级法分为五个级别(图1—6)，然后计算3个时期的综合评价结果(图7)，最后叠加得到源地分布图(图8)。 图1 北戴河新区1991年、2001年和2013年斑块面积级别 图2 北戴河新区1991年、2001年和2013年分维数级别 图3 北戴河新区1991年、2001年和2013年边界破碎度级别 图4 北戴河新区1991年、2001年和2013年斑块重要值级别 根据景观斑块质量综合评价结果，结合各年各级别斑块的面积统计数据(表5)，得到其时间变化和空间分布上的特征。 图5 北戴河新区1991年、2001年和2013年生境质量级别 图6 北戴河新区1991年、2001年和2013年生态服务价值级别 图7 北戴河新区1991年、2001年和2013年景观斑块质量综合评价结果级别 图8 北戴河新区2013年景观生态安全格局源地动态识别 3.1 1991—2013年，北戴河新区景观斑块质量呈现先下降后上升的时间变化规律 从时间动态变化来看，1级斑块的面积1991年为53.24 km2，2001年为47.10 km2，较1991年下降了11.53%，2013年为62.26 km2，较2001年上升了35.19%，可见，1991—2013年，1级斑块的面积总体呈现先下降后上升的态势。2级斑块的面积1991年为67.04 km2，2001年为53.26 km2，较1991年下降了20.55%，2013年为47.02 km2，较2001年下降了11.72%，可见，1991—2013年，2级斑块的面积总体呈现下降的态势，但下降幅度逐渐降低。1级斑块和2级斑块的总面积1991年为120.28 km2，2001年为100.36 km2，较1991年下降了16.56%，2013年为109.28 km2，较2001年上升了8.89%，可见，1991—2013年，1级斑块和2级斑块的总面积总体呈现先下降后上升的态势。总体来说，1991—2001年，景观斑块质量呈现下降的态势，景观破碎化加剧，人为活动对景观的干扰加强，生态安全受到威胁；2001—2013年，景观斑块质量逐渐提高，主要是因为这12 a间，新区经过河道疏浚、沿海防护林培育、土地整治等有利于生态防护措施的实施，生态环境得到改善，总体斑块质量得到一定提升。 3.2 2013年，北戴河新区源地呈现向连片的水体及林地附近散布的空间分布规律 从空间分布来看，2013年源地的总面积为73.34 km2，占北戴河新区总面积的14.25%。其中，1级斑块的面积为37.44 km2，占新区总面积的7.27%；2级斑块的面积为35.90 km2，占新区总面积的6.97%。源地主要分布在水产局及周边、七里海及周边、刘台庄镇东南部以及国有林场南部地区，国有林场所占比例最大，约有73.07%。水产局以及七里海周边生态用地类型主要为水体，平均斑块面积较大，斑块连通性和生态系统服务价值较高，七里海属于黄金海岸国家级自然保护区的核心区范畴，水生动植物资源丰富，其中有国家二级保护动物文昌鱼。刘台庄镇东南部分布有面积较大的成片水田，景观连通性比较好；国有林场是国家级自然保护区，其南部地区有较大面积的水体资源和森林资源，动植物物种丰富，其中有国家二级保护植物狭叶瓶尔小草和世界濒危动物黑嘴鸥，生态系统服务价值非常高。应当加强景观生态安全格局源地的保护力度，降低人为活动对生态环境的不良干扰。 表5 北戴河新区1991年、2001年和2013年各级别景观斑块面积及其占新区总面积的比例级别年份1991年面积/km2比例/%2001年面积/km2比例/%2013年面积/km2比例/%1级级级.2826..8525.424级.2119..2819.875级 4 结论与讨论 本研究以北戴河新区为例，探讨了景观生态安全格局源地的动态识别方法，为更加合理地构建景观生态安全格局提供基础。研究从斑块自身特征属性和斑块对外界环境的响应两个方面选取指标，对北戴河新区景观斑块质量进行了动态评价和分级，在此基础上，选取质量良好且生态系统稳定的斑块作为景观生态安全格局构建的源地。源地分布在景观连通性好、生态服务价值高、生物多样性丰富、生态系统较稳定的区域，这些区域不仅对新区的生态安全格局有重要的影响，也对新区产业尤其旅游业的发展有重大影响，因此，应当加强对这些区域的保护，合理限制人为活动对其景观的干扰。 本研究从空间和时间两个层面识别源地，比单单考虑空间层面更能相对准确地定位源地的位置与规模。由于数据获取方面的限制性，本研究对某些一定程度上影响源地识别的指标如生境质量评估中的土地合法保护程度等未加考虑，今后的研究中，应当进一步探索调查源地识别的影响因素，完善指标体系，以更加准确地定位源地的位置与规模。