于天舒,郑拴宁,朱捷缘,唐明方,董仁才,王 岩

1 中国科学院城市环境研究所城市环境与健康重点实验室,厦门 361021 2 中国科学院大学,北京 100049 3 中国科学院生态环境研究中心,北京 100085 4 赤峰市自然资源局,赤峰 024000

随着社会经济的快速发展,人类活动强度大、自然资源开发利用不合理等多重压力的存在,对生态环境资源产生了巨大影响,导致土地退化、生态系统服务功能降低、城市无序扩张、景观破碎化、生物多样性降低、水污染等生态环境问题日益严峻,严重威胁区域生态安全,制约生态-社会-经济可持续发展[1-2]。所谓生态安全,有广义和狭义两种理解,广义上指的人的生活、健康、安乐、基本权利、生活保障来源、必要资源、社会秩序和人类适应环境变化的能力等方面不受威胁的状态,包括自然生态安全、经济生态安全和社会生态安全,组成一个复合人工生态安全系统；狭义上是指自然和半自然生态系统的安全,即生态系统的完整性和健康水平的整体反映[3]。目前,国内外学者对生态安全已有深入研究,包括生态安全概念内涵[4-8]、评价方法[9-11]、管理对策[12-13]等多个方面。跨流域调水工程水源区的生态安全问题作为调水工程成败的关键,也是诸多学者关注和研究的热点之一[14-17]。

南水北调工程起源于河南省南阳市丹江口水库,是我国规模宏大的跨流域调水工程,对于缓解我国华北地区水资源短缺问题和促进区域间协调发展具有重要作用[15]。同时,南水北调工程建设和跨流域调水也会改变周边地区土地利用和土地覆被状况[18],影响水源区水量和水质,威胁生态系统健康和生态安全。因此,研究调水工程水源区生态安全状况,对于保障水源区生态建设和供水安全,实现调水区域的生态安全和可持续发展至关重要。

目前,我国学者从不同视角对南水北调工程水源区生态安全进行了研究,如水源区水质评价[19-20]、景观格局变化特征分析[21]、生态安全评价[22-23]、生态补偿标准制定[24-26]等。尽管已有研究应用“压力-状态-响应”(Pressure-State-Response,PSR)模型对丹江口库区和南阳市生态安全进行评价[27-28],但仅选用统计数据对区域(行政区)尺度生态安全状况进行评估,并没有考虑土地利用类型、植被覆盖度等空间数据,评价结果难以准确反映水源区生态安全空间分布特征及邻近空间对其影响的关系。

因此,本研究以南水北调中线工程水源地南阳市为例,构建了基于PSR框架的生态安全评价指标体系,结合统计数据和遥感数据,对2000、2010、2015年南阳市、丹江口汇水区、干渠缓冲区生态安全时空变化特征进行分析,并进一步分析了南阳市生态安全变化驱动因素,最后提出相应的对策和建议。本研究能够有效识别南阳市生态安全水平较低的区域,厘清影响其生态安全的主要因素,为政府管理部门精准实施生态工程提供决策依据,对于保障南水北调中线工程高效安全运行和周边区域用水安全具有重要意义。

1 研究区概况

南阳市是南水北调中线工程水源地和渠首所在地,位于河南省西南部(东经110°58′—113°49′,北纬32°17′—33°48′,图1),市辖1市2区10县,总面积约为2.65万km2。南阳市东西北三面环山,南部为南阳盆地,地貌类型复杂多样,地处北亚热带向暖温带过渡地带,属典型的季风大陆半湿润气候,四季分明。市内水系发达,分属汉水流域、淮河流域、黄河流域。南水北调工程水源丹江口水库位于淅川县,工程干渠由市域西南部向东北部延伸,贯穿南阳市。2018年末,全市GDP达到3566.77亿元,年末总人口为1198.07万人,城镇化率达到46.23%。

图1 研究区域位置Fig.1 Location of the study area

2 数据来源和研究方法

2.1 数据来源

本研究主要涉及四种数据类型。生态系统类型数据来源于生态系统十年评估研究成果(http://www.ecosystem.csdb.cn/),包括林地、灌丛、草地、水域及湿地、耕地、城镇用地和未利用地7类用地类型,空间分辨率为30 m。数字高程模型(DEM)来源于地理空间数据云平台(http://www.gscloud.cn)提供的美国地质调查局(U.S.Geological Survey,USGS)SRTM数据,空间分辨率30 m,用于计算地形位指数。NDVI数据来源于美国地质调查局官网(https://lpdaac.usgs.gov)空间分辨率250 m,用于计算植被覆盖度指数[29]。统计数据来源于南阳市统计年鉴(2001年、2011年、2016年)以及河南省水资源公报(2000年、2010年、2015年),并使用ArcGIS软件实现栅格化。空间数据坐标系统投影为China_Lambert_Conformal_Conic,数据重采样后分辨率均为30 m×30 m。

2.2 生态安全评价指标体系构建

本研究结合水源区生态系统自身特点,遵循数据可靠性、易获取性、科学代表性、系统完整性、指标多样性与层次性等原则,构建了基于PSR框架的南阳市生态安全评价指标体系,包含1个目标层、3个准则层、16个指标层。生态压力指标主要用于反映人类生产生活活动对生态系统产生的干扰和破坏[30-31],生态状态指标主要用于反映土地利用强度和生态环境质量状况,生态响应指标主要用于反映针对生态问题的管理政策和治理措施,其中生态压力和生态响应指标为统计数据、生态状态指标为遥感数据。对要素指标进行标准化处理消除矢量单位的影响,然后采用层次分析法[32]结合文献调查及专家打分确定指标权重,具体指标、指标含义及权重如表1所示。

表1 南阳市生态安全评价指标体系Table 1 The evaluation index system used to assess the ecological security of Nanyang City

2.3 生态安全指数

生态安全指数(ESI)反映了指标因子对生态安全的影响程度和生态环境抵御胁迫的能力,指数越高表明生态安全水平越高。本研究通过ArcGIS 10.2的Raster Calculator工具计算获取2000—2015年南阳市生态安全指数。参考相关文献[27-28,35-36],并结合研究区实际情况,将ESI分为5个安全等级:低度安全(0.00—40.00)、较低安全(40.00—50.00)、一般安全(50.00—60.00)、较高安全(60.00—80.00)、高度安全(80.00—100.00)。生态安全指数计算公式如下所示:

式中,ESI为南阳市生态安全评价指数,ESI为指标Ci标准化后的指标值,Wi为指标的权重值。

3 结果

3.1 南阳市生态安全水平时空变化

2000、2010、2015年南阳市生态安全指数平均值分别为48.88、49.76、50.60,整体处于上升趋势；3个时期生态安全水平分别处于较低安全、较低安全、一般安全,生态安全状况有所改善。从空间分布上看(图2),2000、2010、2015年南阳市生态安全水平均呈现出西北高、东南低的趋势,整体处于较低水平；同时在西峡县、内乡县、桐柏县等县区内部也存在明显区域差异。2000—2015年南阳市生态安全水平变化具有明显的区域差异:西部地区ESI整体呈增加趋势,其中西峡县增幅较小,南阳市增幅明显,超过10%；东部地区ESI整体呈减小趋势,其中社旗县减小幅度较大。

图2 2000—2015年南阳市生态安全等级图Fig.2 The distribution of ecological security degrees in Nanyang City from 2000 to 2015

统计结果显示(图3),2000年南阳市低度安全、较低安全、一般安全、较高安全、高度安全的面积分别为380、13590、11590,10、0 km2；2010年南阳市低度安全、较低安全、一般安全、较高安全、高度安全的面积分别为1280、9380、14860、40、0 km2；2015年南阳市低度安全、较低安全、一般安全、较高安全、高度安全的面积分别为270、9920、15250,130、0 km2,3个时期生态安全水平整体较低,有待进一步改善。2000—2015年南阳市较低安全水平面积明显减少(27.01%),一般安全水平面积明显增加(31.58%),较高安全水平面积也呈现一定的增加趋势。

图3 南阳市生态安全等级面积分布Fig.3 Area distribution of ESI levels in Nanyang City

2000—2015年南阳市各行政单元生态安全水平具有明显的区域差异(图4),ESI均值从大到小依次为淅川县>方城县>南召县>市辖区>西峡县>内乡县>桐柏县>社旗县>邓州市>镇坪县>唐河县>新野县,其中生态安全水平处于一般安全的县区有6个,处于较低安全的县区有6个,研究时段内均未出现生态安全等级为高度安全、较高安全和低度安全的县区,生态安全水平为一般安全的县区数量有所增加。2000—2015年南阳市淅川县、内乡县、邓州市3个南水北调水源区以及南阳市辖区生态安全水平明显提高,社旗县生态安全水平呈降低趋势,而唐河县和新野县在2010年出现一个明显的转折。

图4 2000—2015年南阳市各行政单元平均生态安全指数Fig.4 Average ESI value of each administrative region unit in Nanyang City from 2000 to 2015

3.2 丹江口水库汇水区及干渠周边生态安全水平时空变化

2000、2010、2015年丹江口水库汇水区生态安全指数均值分别为51.52、52.80、53.70,整体处于上升趋势,生态安全水平有所提升。空间上(图5),丹江口水库汇水区生态安全水平整体呈现出南部高、北部低的趋势,具有明显的区域差异,特别是2000年北部地区大部分区域处于较低安全水平。2000—2015年丹江口水库汇水区生态安全水平整体均呈现明显的增加趋势,北部地区较低安全水平面积大幅度减少,南部地区逐渐出现较高安全水平面积分布。

图5 南阳市丹江口水库汇水区生态安全等级分布图Fig.5 The distribution of ecological security degrees of Danjiangkou Reservoir catchment area

为了进一步分析南水北调工程干渠生态安全状况,本研究以干渠两侧3 km为缓冲半径,提取得到南水北调干渠缓周边生态安全状况(图6)。2000、2010、2015年南水北调工程干渠缓冲区生态安全指数均值分别为47.98、49.57、51.02,整体处于上升趋势。南水北调工程干渠下游生态安全水平明显高于上游区域。

图6 南阳市南水北调中线工程水渠周边生态安全等级分布图Fig.6 The distribution of ecological security degrees within a 3 km buffer zone for the channel of the Middle Route of South-to-North Water Diversion Project in Nanyang City

3.3 南阳市生态安全影响因素

生态压力指数、生态状态指数、生态响应指数对2000—2015年南阳市生态安全指数贡献度增加量存在明显差异,生态相应指数的变化幅度最大,生态状态指数变化幅度最小。如图7所示,生态压力指数对生态安全指数的贡献度增加量均为负值,贡献度明显减少,其中市辖区、邓州市、新野县、社旗县的贡献度增加量低于-10%；生态状态指数对生态安全指数的贡献度增加量有正有负,以负值为主,但变化幅度整体相对较小,其中市辖区的贡献度增加量低于-10%；而生态响应指数对生态安全指数的贡献度增加量均为正值,且变化幅度均相对较大,其中市辖区、邓州市、西峡县、新野县、南召县、内乡县的贡献度增加量大于10%,仅有唐河县的贡献度增加量小于5%。

图7 南阳市各区县生态压力、状态、响应指数贡献度增加量 Fig.7 The contribution increases of ecological pressure, state and response index in Nanyang city

根据Pearson相关分析结果(表2),不同时期南阳市生态安全水平影响因素存在差异,且生态安全水平与生态状态指标的相关性整体较高,与生态响应指标的相关性较低。2000年,南阳市生态安全水平与人口密度、单位面积化肥施用量、生态弹性度、均匀度指数四个指标显著相关；2010年,南阳市生态安全水平与单位面积化肥施用量、生态弹性度、多样性指数、均匀度指数四个指标显著相关；2015年,南阳市生态安全水平与多样性指数、均匀度指数两个指标显著相关；2000—2015年生态安全水平与生态响应指标之间不存在显著相关关系。

表2 南阳市生态安全指数与影响因素相关分析结果Table 2 Correlation analysis results of ESI and influencing factors in Nanyang City

4 讨论

南阳市各地区生态环境条件、社会发展水平存在差异,导致区域内生态安全水平时空差异明显。空间上,南阳市西北部地区以林地、草地、灌丛等生态用地为主,同时还是丹江口水库汇水区所在地,生态质量高,生态压力小,因此生态安全水平相对较高；而在东南部地区以耕地、建设用地等人类活动强度较大的用地类型为主,生态干扰大,生态环境质量相对较低,因此生态安全水平处于较低水平。时间上,随着人们生态环境保护意识的增强,人为造成的生态环境破坏问题得到缓解,生态环境压力不断减小,生态环境状态持续改善,加上退耕还林、污染治理、生态恢复的重要措施的大力实施,使得区域内生态安全水平得到明显改善和提高。南阳市市辖区城镇用地的扩张给区域内生态环境造成巨大压力,但是市辖区经济发展水平高,科技投入、环保技术水平高,很大程度上弥补甚至超过了生态压力对生态安全水平的贡献度。

南阳市生态安全水平整体上得到了持续改善,但仍然处于相对较低的水平,2000、2010、2015年区域内生态安全等级主要为较低安全和一般安全,较高安全水平以上的面积均低于1%,因此还有很大的提升空间。尤其是在东南部地区的唐河、社旗、新野、镇坪等农业种植区,要严格控制农药化肥施用量,倡导使用有机肥,同时加强点源、面源污染治理力度,防治生态环境污染；同时加强农田道路两旁生态廊道体系建设和河渠两边生态环境保护,提高景观多样性,改善景观破碎化；合理控制城镇无序扩张,加强城市生态空间管理,对于城市绿地的数量、质量、空间布局进行合理规划,全面提升城市生态系统服务功能。而在西北部的淅川、西峡、内乡、南召等县区,基于良好的生态环境本底,继续强化森林植被保护和抚育,提升植被覆盖度；加强丹江口水库汇水区以及南水北调工程沿线生态保护与修复,保障南水北调工程供水安全,全面提升南阳市生态安全水平。

本研究从宏观角度对南水北调工程实施前和通水后,南阳市各区县、丹江口汇水区和干渠周边生态安全水平进行了分析。不仅需要社会、经济数据等作为生态响应数据,还需要能够反映自然生态系统状态的空间数据。与常规的仅用统计数据构建综合指标体系的方法相比,我们的研究结合了植被覆盖度、地形位指数等空间栅格数据,并将评价结果落实到空间栅格上,不仅不能够显示南阳市不同县区之间的差异,也能够显示出各个县区内部的区域差异,具有良好的可视化效果和实用性。

需要指出的是,为了便于在时间维度上进行生态安全水平的比较,本研究将2000年、2010年和2015年3期指标值进行统一标准化。数据标准化是相关研究中的基本问题,但针对不同时期开展评价时,目前尚无文献明确说明其标准化是每个时期单独进行还是把不同时期的数据合在一起进行。本研究认为,将2000年、2010年和2015年3期指标值进行统一标准化是合理的,这样得到的不同时期的结果可以进行比较。例如,假设某指标两个时期的指标值分别在[10,20]和[80,90]上均匀分布,若将每一时期的数据单独标准化,标准化后两个时期的指标值一样,而两者的实际情况却相距甚远,因此每个时期单独标准化不宜进行时间维度上的比较。

5 结论

2000—2015年南阳市生态安全水平持续改善,且呈现出东南低、西北高的趋势,这与区域内生态环境状况变化、生态工程建设力度、快速城镇化等多种自然和人为因素密切相关。随着人们生态环境保护意识的不断增强,生态保护措施实施力度不断加大,生态压力对南阳市生态安全水平的贡献程度被生态响应的贡献度所替代,促使生态安全状况不断向好。为了进一步提升区域生态安全水平,南阳市需要继续加强生态环境保护与治理力度,全面改善生态环境质量。

本研究基于栅格单元进行生态安全水平空间量化,能够更加准确地反映南阳市生态安全水平的区域差异,有助于精准实施生态环境保护治理工程措施。此外,在开展不同时间的评价时,为了便于评价结果在时间维度上可比,应该将不同时间的指标值统一进行标准化。