黄河下游城市群生态系统健康时空演变特征

谷瑞丽多玲花邹自力黎德华

(1.东华理工大学 江西省软科学研究培育基地资源与环境战略研究中心,江西 南昌330013;2.东华理工大学测绘工程学院,江西 南昌330013;3.自然资源部 环鄱阳湖区域矿山环境监测与治理重点实验室,江西 南昌330013)

生态系统健康(ecosystem health)指一个生态系统所具有的稳定性和可持续性,即在时间上具有维持其组织结构、自我调节和对胁迫的恢复能力[1]。伴随着城市化进程的加快,人类活动影响加剧了生态系统结构的变化,从而对生态系统构成严重威胁,暴露出诸如大气污染、水土流失、城市热岛、景观结构遭到破坏、生境质量受到威胁等生态环境问题[2]。《黄河流域生态保护和高质量发展规划纲要》提出,沿黄河省区要毫不动摇地走生态优先、绿色发展为导向的现代化道路,加快构建人与自然和谐共生的现代绿水青山体系。黄河流域生态保护和高质量发展关乎中华民族伟大复兴,推动黄河流域生态保护建设,对促进黄河流域重要经济带的发展具有极其重要的战略意义,因此对黄河下游城市群的生态系统健康水平进行评估,有助于黄河流域生态保护和经济社会可持续发展目标的实现。鉴于黄河流域高质量发展中存在的困境及其重要地位,黄河流域的生态系统健康问题已经成为专家学者们关注的热点。牛明香等[3]以黄河河口区为研究区域,河口及与之相毗邻的陆域、海域作为整体生态系统,根据压力—状态—响应(pressurestate-response,PSR)模型建立评价体系,结合综合指数法对河口区生态系统健康水平进行测度。徐辉等[4]围绕结构、功能、过程、发展4个层面构建生态系统健康评价指标体系,以模糊物元模型为基础,从时空维度评价2010—2018年生态系统健康水平。刘玒玒等[5]选择关中5个城市为研究区,基于自然、社会和经济因素选择20项指标建立生态系统健康评价框架,运用模糊综合评价法评价关中5市生态系统健康,为城市生态环境保护和经济发展提供理论和现实依据。目前,黄河流域生态系统健康研究主要有3点不足:①以PSR模型指标体系为例,现有的生态系统健康评价指标体系适用空间尺度较窄,社会响应指标具有滞后性,量化结果容易受到主观影响[6]。②评价指标体系纷繁复杂,既有自然指标,也有社会和经济等指标,一定程度上脱离了生态系统健康评价的“生态”主旨[7]。③在研究尺度上,黄河流域已有的生态系统健康评估研究是以行政区为研究单元的,不能精细反映研究区内部生态系统健康水平的差异[8]。随着人类生活质量的改善,人们逐渐重视生态系统为人类福祉带来的直接或间接的利益,学者们将生态系统服务作为一个重要指标放进生态系统健康评价体系之中,提出更能反映自然生态环境对人居环境影响的VORS(vigor-organization-resilience-services)生态系统健康评价模型[9-11]。

黄河流域下游地区由于开发强度高,人与自然关系日益紧张,在当前快速城镇化背景下,人类活动导致下游地区生态系统脆弱性升高、稳定性较差,生态环境受到威胁。因此,本研究以黄河下游城市群为研究对象,针对黄河下游城市群特点,从生态系统健康内涵出发,以VORS模型建立评价指标,从格网尺度探讨黄河下游城市群生态系统健康水平及其时空分布特征,以期为城市群高质量发展提供有效参考。

1 研究区域与研究方法

1.1 研究区概况

迄今为止,黄河流域发展形成了7大城市群,构成“3+4”空间格局,即3个区域级城市群和4个地区性城市群[12]。山东半岛城市群和中原城市群,这两大区域级城市群位于华北平原,同属于黄河下游城市群[13]。总人口2亿,占黄河流域人口总数的47%,地区生产总值为13.46万亿元,占流域生产总值的56%[14]。两大城市群在空间位置上紧密相邻,经济相互依存,产业分工上相辅相成。参照已有研究成果对黄河下游地区的划分,选取中原城市群的核心区域,即河南省所辖地级市和山东半岛城市群作为研究对象[15]。山东半岛城市群和中原城市群覆盖了黄河下游大部分区域,在一定程度上能够体现出整个黄河下游的概况。

1.2 数据来源

2000,2010和2020年土地利用数据来源于中国科学院资源环境科学与数据中心(http:∥www.resdc.cn),2000年和2010年数据生产制作主要使用Landsat-TM/ETM遥感影像数据,2020年主要使用Landsat 8遥感影像数据,通过人机交互目视解译生成,空间分辨率为1 km×1 km,数据精度在88%以上[16],满足研究需要。NDVI(normalized difference vegetation index)采用NASA官网MODIS(moderate resolution imaging spectroradiometer)产品(MOD13Q1),通过MRT转换工具将hdf格式转换为tiff格式,空间分辨率为250 m,时间分辨率为16 d。社会经济数据来自2000,2010和2020年《全国农产品成本收益资料汇编》《河南省统计年鉴》《山东省统计年鉴》。为了避免按照行政区划指标均匀分布问题,研究尺度不再采用行政区划尺度,而是采用网格尺度,本研究所设定的网格尺度是10 km×10 km。

1.3 研究方法

1.3.1 生态系统健康评估框架构建 评估黄河下游城市群生态系统健康水平,对制定合理有效的生态系统健康保护政策具有十分重要的意义。本研究采用VORS模型(图1),选择不同指标分别测度生态系统活力、生态系统组织力、生态系统恢复力和生态系统服务价值,以此对黄河下游城市群生态系统健康进行动态评估。

图1 生态系统健康评价体系框架

本研究构建的VORS模型是VOR模型的扩展,将生态系统服务作为一个重要指标纳进生态系统健康评价体系之中,更好地将生态系统与人类社会联系起来,以借助不同指标衡量生态系统在功能、结构、过程和服务等方面的健康状况[11]。本研究主要从“目标层—因子层—指标层”构建VORS模型,目标层是指生态系统健康评估,可以揭示研究区生态系统健康水平,活力、组织力、恢复力和服务构成了因子层,因子层下的特定指标,即植被覆盖指数、景观格局指数(景观连通性、景观异质性)、生境质量和生态系统服务价值构成了指标层,指标的变化趋势可以揭示生态系统健康的过程,为政府精确有效地制定政策提供依据。本文通过构建VORS模型,定量化分析生态系统健康水平的时序演变、空间特征以及空间关联性,从而揭示黄河下游城市群生态系统健康水平的发展态势与空间规律。

(1)生态系统活力(ecosystem vigor)通常通过植被覆盖度或植被净初级生产力(net primary productivity,NPP)来反映生态系统的代谢和初级生产力[17]。本研究参考生态环境状况评价技术规范(HJ192-2015),以植被覆盖指数来表征生态系统活力。植被覆盖指数(vegetation coverage index)是评价区域植被覆盖的程度,利用评价区域单位面积归一化植被指数(NDVI)表示,具体公式为[18]:

式中:Pi表示5—9月像元NDVI月最大值的均值,NDVI采用NASA官网MOD13数据,空间分辨率为250 m;n代表区域像元数;Aveg代表植被覆盖指数的归一化系数,参考值为0.012 1。

(2)生态系统组织力(ecosystem organization)被用来反映生态系统各要素之间的相互作用,通常用以表达生态系统稳定性,本研究选择景观异质性和景观连通性来定量,并赋予不同权重计算组织力[19-20]。景观连通性主要由整体景观连通和生境连通性共同决定,分别使用蔓延度(CONTAG)和斑块凝聚度指数(COHESION)对其进行量化。景观异质性具体使用香农多样性指数(SHDI)来表征。参考以往研究[11],整体景观的异质性和连通性的权重赋值为0.35,生态系统组织力(O)计算公式为:

(3)生态系统恢复力(ecosystem resilience)也被称为生态系统弹力,主要是指系统在遭受到干扰或威胁后恢复到原状的能力。本研究采用国际上广泛使用的In VEST模型生境质量模块来表征[21]。生境质量所涉及到的参数主要参考以往的研究[22-24],并在此基础上结合黄河下游城市群区域特征进行修正。生境退化度和生境质量的计算公式为[25]:

式中:Dxj表示生境退化程度;Wr代表生态威胁因子的权重;ry表示威胁因子强度;irxy代表威胁因子r的最大影响距离;βx代表示生境抗干扰水平;Sjr表示生境对威胁因子的敏感度。

式中:Qxj为土地利用类型j上栅格x的生境质量;Hxj为土地利用类型j中x的生境适宜性;k为半饱和常数。生境质量值越大其恢复力越好。

(4)生态系统服务(ecosystem services)是指人类从生态系统获得支持服务、供给服务、文化服务、调节服务等对人类生活有贡献的相关服务或产品[26]。本研究采用当量因子法,依据Costanza建立的生态系统服务评价模型,参考谢高地[27]改进后的生态系统服务系数来修正研究区生态系统服务系数,1个生态服务价值当量因子的经济价值为单位面积粮食产量的1/7[28]。本研究基于2000,2010和2020年主要粮食作物的产量、播种面积和平均价格计算得到黄河下游城市群单位生态系统服务的经济价值为1 449.07元/hm2,研究区的生态系统服务价值计算公式[29-30]:

式中:ESV代表生态系统服务价值;j代表生态系统服务类型;Ai代表第i类土地利用类型的面积(hm2);VCi代表第i类土地利用类型单位面积生态系统服务价值〔元/(hm2·a)〕;ECj代表某类土地利用类型第j项生态系统服务价值当量;k代表生态系统服务类型数量;Ea代表1单位生态系统服务的经济价值〔元/(hm2·a)〕。

1.3.2 生态系统健康指数测度 由于各项评价指标的量纲各不相同,其实际数量也存在很大不同,所以本研究必须先将各指标标准化才能对生态系统健康进行评价。本研究所使用的方法是极值归一化法。计算公式为:

式中:Xi是第i年指标的实际值;min(Xi)是最小值;max(Xi)是最大值;Pi是标准值。

参考袁毛宁[11]对广州市生态系统健康的研成果,基于VORS框架建立评价模型:

式中:EHI是生态系统健康指数;V,O,R,S分别为活力、组织力、恢复力与生态系统服务。

1.3.3 空间自相关 空间自相关是用来度量附近观测值之间是否存在空间关联性。本研究采用全局莫兰指数(Moran’sI)分析整个研究区空间关联性。计算公式为[31]:

式中:n是空间上研究单元总数;xi,xj分别表示区域i和j的生态系统健康值;Wij表示区域i,j的空间权重矩阵;S2为所有要素权重的集合;ˉx表示平均生态系统健康值。

为了进一步分析研究区生态系统健康空间分布格局特征,运用Getis-OrdG*i指数测度局域关联特征,识别具有统计显著性的冷点区和热点区。计算公式为[32]:

式中:G*i(d)为局部空间自相关值;xi为研究单元i的观测值;Wij(d)为空间权重矩阵。

2 结果与分析

2.1 生态系统健康评价指标分析

根据黄河下游城市群生态系统健康评价模型,对2000,2010和2020年活力、组织力、恢复力和生态系统服务4个因子分析(如图2—5所示),展示了生态系统健康4个评价指标标准化后的空间格局。2000,2010和2020年黄河下游城市群生态系统活力的平均值分别为0.598,0.687和0.664,生态系统活力先上升后下降,总体上呈明显上升的趋势。生态系统活力高值区主要分布在中原城市群西部,低值则分布在山东半岛城市群北部沿海区域。由于植被覆盖指数受气温、年降水量、太阳辐射等多种因素综合影响,纬度更低的中原城市群的生态系统活力普遍高于山东半岛城市群。

图2 黄河下游城市群生态系统活力时空分布

2000,2010和2020年黄河下游城市群生态系统组织力的平均值分别为0.601,0.512和0.514,生态系统组织力在2010年前呈快速下降趋势,2010年后相对稳定。研究期内黄河下游城市群生态系统组织力的变化趋势是先下降后保持稳定,总体上呈明显下降的趋势,其中以研究区西部地区变化较为明显。生态系统组织力的高值主要分布在山东半岛城市群地区,而低值则集中在中原城市群,总体上呈现出东部高,西部低的分布格局。

图3 黄河下游城市群生态系统组织力时空分布

2000,2010和2020年黄河下游城市群生态系统恢复力的平均值分别为0.429,0.425和0.416,生态系统恢复力呈现下降的趋势。中原城市群西部地区生态系统恢复力的下降趋势尤为明显。中原城市群生态系统恢复力普遍高于山东半岛城市群,总体呈现出西部区域生态系统恢复力高,中部区域生态系统恢复力低,东部区域生态系统恢复力处于中等水平的分布特征。

图4 黄河下游城市群生态系统恢复力时空分布

2000,2010和2020年黄河下游城市群生态系统服务的平均值分别为0.207,0.213和0.125,生态系统服务呈现明显下降的趋势。近10 a间呈现明显下降的趋势,主要是由耕地、草地以及林地大面积流失所导致的。生态系统服务的高值区集中在中原城市群西部和山东半岛城市群南部,而山东半岛城市群西部和中原城市群东部区域则属于生态系统服务价值低值区。

图5 黄河下游城市群生态系统服务时空分布

2.2 生态系统健康时空格局及动态变化

综合考虑活力、组织力、恢复力、服务4个因子,得到黄河下游城市群生态系统健康时空分布图(图6)。根据以往的研究[33],将生态系统健康指数划分为5个等级,其等级分为差(0～0.2)、较差(0.2～0.4)、一般(0.4～0.6)、较好(0.6～0.8)、良好(0.8～1)。研究区内生态系统健康等级为较差的占主导地位,部分地区的生态系统健康等级为一般,少数地区的生态系统健康等级为差,本研究不存在生态系统健康等级为较好和良好的地区,表明该地区的生态系统健康是脆弱的。平均值代表了研究区生态系统健康的一般分布水平,标准差揭示了数据的离散程度,标准差越大表示数据越趋于离散。2000,2010和2020年黄河下游城市群生态系统健康发生了不同程度的变化,生态系统健康指数平均值从2000年的0.325下降到2020的0.320,总体上呈现略微降低的趋势,标准差从2000年的0.106上升到2020年的0.112,表明生态系统健康水平分布更加离散,即黄河下游城市群生态系统健康等级为一般的地区得到进一步保护,等级为差的地区进一步恶化,整体生态系统健康状况退化。从空间分布来看,2000,2010和2020年黄河下游城市群生态系统健康水平具有显著的空间异质性,各等级空间分布集中,整体生态系统健康水平脆弱。从研究区中部延伸到东北部的大部分区域,生态系统健康等级为较差且占绝大多数,占比高达约74%。由于该区域是耕地集中分布区,容易受到人类活动的严重干扰,故随着城镇化进程加快,建设用地增加导致生境斑块破碎化、连通性差并进一步导致该区域生态系统健康状况恶化。生态系统健康等级为一般则主要分布在中原城市群的西南部以及山东半岛城市群中部北部小部分地区。该区域林草地面积广布,景观多样性较为丰富,人类活动干扰比较少,同时结构相对完整,生态系统良好。总而言之,随着城市群的快速扩张以及人类活动对生态系统的干扰,生态问题越发突出,生态系统健康水平也随之下降。

图6 黄河下游城市群生态系统健康时空分布

2.3 生态系统健康空间集聚特征

为了进一步探究生态系统健康空间集聚分布情况,绘制了生态系统健康热点和冷点分布图(图7)。2000—2020年黄河下游城市群生态系统健康全局莫兰指数分别为0.512,0.522,0.543。由此可见,20 a间生态系统健康水平空间依赖性显著,生态系统健康分布在整体上具有正向的空间集聚特征,即生态系统健康指数高的区域空间集聚,而生态系统健康指数低的区域相互临近。生态系统健康热点区主要分布在中原城市群西南部以及山东半岛城市群中部及东北部,该区具有较高的植被覆盖度、生境质量及抗干扰能力,因此出现高值聚集。研究期间生态系统健康冷点区面积所占比例上升,热点区面积则略微下降,生态系统健康状况不容乐观。

图7 黄河下游城市群生态系统健康热点和冷点分布

3 讨论

2000—2020年,生态系统健康状况总体上呈现略微降低的趋势。一方面,黄河下游城市群耕地和林地面积持续减少,同时受到气候变化和人类干扰的影响,生态状况极其脆弱,生态系统健康面临巨大压力,另一方面,随着城市化进程的加快,人口规模的增长,建设用地规模快速扩张,对生态系统的干扰逐渐增强,从而导致生态系统健康水平下降[33]。在空间分布上,研究区中部延伸到东北部的大部分区域,生态系统健康水平较差,该区域是耕地集中分布区,容易受到人类活动的严重干扰,随着城镇化进程加快,建设用地增加导致生境斑块破碎化、连通性差并进一步导致该区域生态系统健康状况恶化[34]。因此,黄河下游地区开发建设要加强生态环境保护,禁止盲目扩大建设用地,实现人与自然的和谐相处,有效保护和修复生态环境,推动黄河下游城市群健康可持续发展。

本研究也存在一些局限之处,如评价指标体系的量化过于依赖于土地利用数据,造成结果相对单一,不能很好地反映社会经济等因素对生态系统健康的影响,在一定程度上影响了生态系统健康的评价结果,下一步可以在基于土地利用相关指标的基础上,增加社会经济、生态文化健康等指标。另外由于研究尺度较大,数据分辨率可能会对研究结果产生一定程度的偏差,为了更加精准的测算生态系统健康水平状况,要使用更加精确地土地利用数据。综上所述,未来可以通过提高土地利用分辨率和丰富评价指标体系两方面增强评价结果的科学性。

4 结论

(1)2000—2020年,研究区内生态系统活力、组织力、恢复力和生态系统服务均发生了不同程度的变化。生态系统活力从2000年的0.598上升到2020年的0.664,总体上呈明显上升的趋势;2000,2010和2020年生态系统组织力分别为0.601,0.512和0.514,总体上呈明显下降的趋势,但2010年后相对稳定;生态系统恢复力从2000年的0.429下降到2020年的0.416,生态系统服务从2000年的0.207下降到2020年的0.125,基本维持在较低水平,在一定程度上降低了生态系统健康水平。

(2)2000—2020年,生态系统健康指数平均值从2000年的0.325下降到2020年的0.320,标准差从2000年的0.106上升到2020年的0.112,黄河下游城市群生态系统健康水平总体上呈现略微降低的趋势,黄河下游城市群生态系统健康等级为一般的地区得到进一步保护,等级为差的地区进一步恶化。黄河下游城市群生态系统健康水平具有空间异质性,各等级空间分布集中,生态系统健康程度表现为较差的占主导地位,所占比例高达约74%,分布在研究区中部延伸到东北部的大部分区域,整体生态系统健康水平脆弱。

(3)2000—2020年,黄河下游城市群生态系统健康空间依赖性显著,2000,2010和2020年全局莫兰指数分别为0.512,0.522和0.543,生态系统健康分布在整体上具有正向的空间集聚效应。生态系统健康热点区主要分布在中原城市群西部及南部以及山东半岛城市群中部及东北部。