本文以生态系统健康理论为基础,构建了以景观格局指数为主的活力–组织力–弹力–生态系统服务的生态健康评价体系,来研究三江源湿地的生态健康性。采用层次分析法(AHP)确定评价指标及各斑块类型的权重,结合计算出的景观格局指数值,利用综合指数法评价了三江源地区湿地环境的生态健康状况。通过对比2018年与1988年的景观指数值表明三江源湿地景观破碎度增大,景观异质性增强,景观的形状复杂化和不规则化增大。 Based on the theory of ecosystem health, this paper constructs an ecological health evaluation system based on landscape pattern index-based vitality-organization-elasticity-ecosystem services, and studies the ecological health of Sanjiangyuan wetlands. Using the hierarchical analysis method (AHP) to determine the evaluation index and the weight of each plaque type, combined with the calculated landscape pattern index value, the ecological health status of wetland environment in Sanjiangyuan area was evaluated by the comprehensive index method. By comparing the landscape index values in 2018 and 1988, it is shown that the landscape fragmentation of Sanjiangyuan wetland severity increases, the landscape heterogeneity increases, the shape of the landscape is complicated and irregularly increases.
常 敏1,2,3,4
1陕西省土地工程建设集团有限责任公司,陕西 西安
2自然资源部退化及未利用土地整治工程重点实验室,陕西 西安
3陕西地建土地工程技术研究院有限责任公司,陕西 西安
4陕西省土地整治工程技术研究中心,陕西 西安
收稿日期:2019年9月3日;录用日期:2019年9月18日;发布日期:2019年9月25日
本文以生态系统健康理论为基础,构建了以景观格局指数为主的活力–组织力–弹力–生态系统服务的生态健康评价体系,来研究三江源湿地的生态健康性。采用层次分析法(AHP)确定评价指标及各斑块类型的权重,结合计算出的景观格局指数值,利用综合指数法评价了三江源地区湿地环境的生态健康状况。通过对比2018年与1988年的景观指数值表明三江源湿地景观破碎度增大,景观异质性增强,景观的形状复杂化和不规则化增大。
关键词 :湿地,三江源,景观格局指数,生态系统健康,AHP
Copyright © 2019 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
湿地是地球上非常重要的资源,也是很重要的生态系统。它具有蓄洪防旱、降解污染、调节气候等诸多生态功能,被誉为“人类摇篮”、“人类基因库”、和“地球之肾” [
80年代,Schaeffer首次提出了有关度量生态健康的问题 [
本研究以黄河的源头玛多县为例,通过研究三江源湿地的活力、组织力、弹力和生态系统服务,更清楚地认识到三江源湿地生态建设本质特征,从而为湿地恢复与保护、湿地生态环境规划和管理提供科学依据,对促进区域的可持续发展有重要的现实意义。
三江源保护区气候属青藏高原气候系统,在青藏高原强大的高原下垫面和周围大致均匀环境场,巨大的地理空间孕育出了一个独特的气候单元。全年平均气温一般在−6℃~4℃之间,极端最低气温−48℃,极端最高气温28℃。年平均降水量一般在250~705 mm之间,年蒸发量相对较大,一般在730~1700 mm之间。总的气候特征是年温差小,日温差大,日照时间长,辐射强烈,风沙大,植物生长期短,绝大部分地区无绝对无霜期,自然灾害频繁。
源区水资源十分丰富,水质良好。区内河流纵横,湖泊湿地星罗棋布,是世界上河流最多、高海拔地区湿地面积最大的保护区,被誉为“中华水塔”。同时也是我国重要的湿地分布区,是世界上分布海拔最高的湿地 [
1) 层次分析法(AHP法)
由于每个指标对生态系统健康的贡献值是不同的,因此需要定量的确定指标因子的权重。指标因子的权重值发生变化会直接影响到评价结果,科学定量地给湿地生态系统健康评价指标进行权重赋值显得尤为重要。层次分析法(Analytical Hieraehy Process,简称AHP)是现在常用的一种用来定量分析指标权重的一种方法。这种方法的特点是:所需的定量化数据较少,对问题的本质以及所涉及的因素和内在关系分析的比较透彻;思路简单明了,将使用者的思维过程数量化、条理化以便于计算。同时也存在一些缺点:不能为决策提供更加新颖的方法,而且其所需要的定量数据较少,定性分析偏多,缺乏一定的客观性 [
2) 综合指数法
综合指数法将各项指标转化为同度量的个体指数,便于将各项指标综合起来,以综合指数为综合评比排序的依据。各项指标的权数是根据其重要程度决定的,体现了各项指标在综合值中作用的大小。综合指数法的基本思路则是利用层次分析法计算的权重和通过fragstats计算取得的数值进行累乘,然后相加,最后计算出指标的综合评价指数。
以玛多县为例,空间数据包括2018年的湿地遥感影像数据,利用Landsat卫星影像,条带号为134-35、134-36、133-35、133-36、133-37、135-35、135-36。遥感影像均为ETM数据。本文遥感影像均为5月至10月份,云量小于10,比较适合本文研究。
遥感影像的解译是指从遥感图像上获取目标地物的信息的过程,有目视解译和计算机解译两种常见的方法 [
图1. 分类影像图
湿地的生态系统功能是由系统功能的转运动态以及景观结构来支持的,故景观格局在一定程度上对湿地生态系统健康非常重要。景观格局指数是在空间邻接关系的基础上提出来的,因此可以作为系统健康评价的比较适宜的指标。本研究从三江源区湿地生态系统自身出发,以景观格局指数为主体建立了三江源湿地生态系统健康评价指标体系,如图2所示。
图2. 三江源湿地生态系统评价指标体系(APH结构)
1) 层次分析法确定指标权重
将各层之间的每项因子的联系标出,分层建立比较判断矩阵。矩阵由每一对因子的相对重要性比值构成。求出该层的相对权重,然后逐层推算求出全部因子的权重系数。矩阵中 b i j 为两因子的相对重要标度, b i j 表示对于Ak而言,取因子重要性程度对比 b i / b j 的比值。一般采用1、3、5、7来作为标度,其含义依次为:两个因素一样重要、一个因素比另一个较重要、重要、重要的多;2、4、6则为两相似判断的中间值,表示折中。而各自的倒数则具有类似相反的意义。矩阵计算后要进行一致性检验,一致性指标C.I. (Consitency Index)计算公式如下(1),它的值越大,表明判断矩阵偏离完全一致性的程度越大,C.I.的值越小,表明判断矩阵越接近于完全一致性。对于多阶判断矩阵,引入平均随机一致性指标R.I. (Random Index),表1给出1~8阶正互反RI为随机一致性比例。判断矩阵一致性指标C.I.与同阶平均随机一致性指标R.I.之比称为随机一致性比率C.R. (Consistency Ratio)。当C.R. < 0.10时,便认为判断矩阵具有可以接受的一致性。当C.R. ≥ 0.10时,就需要调整和修正判断矩阵,使其满足C.R. < 0.10,从而具有满意的一致性。
n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
---|---|---|---|---|---|---|---|---|
RI | 0 | 0 | 0.58 | 0.90 | 1.12 | 1.24 | 1.32 | 1.41 |
表1. 1~8阶正互反RI一致性指标RI
CI = λ max − n n − 1 (1)
首先对准则层和指标层各个要素构建判断矩阵。各层次中要素的排序(相对重要程度),亦即权重计算,实际上可归结为利用和积法计算判断矩阵的特征向量,再利用公式(2)计算最大特征值。通过计算整理并检验,得到如下权重结果。以表2~6为各目标层参评因子的判断矩阵及其计算结果。
B1 | P1 (NDVI) | P6 (多样性指数) | P2 (景观总面积) | Wi | |
---|---|---|---|---|---|
P1 | 1 | 3 | 7 | 0.6434 | |
P6 | 1/3 | 1 | 5 | 0.2828 | |
P2 | 1/7 | 1/5 | 1 | 0.0738 |
表2. 准则层B1 (活力层)与其相关因子P的判断矩阵
备注:经检验λmax = 3.0657,其一致性指标CI = 0.0329, RI = 0.58, CR = 0.0567 < 0.1。
B2 | P5 (蔓延度指数) | P4 (聚合度指数) | P3 (分形维数) | P6 (多样性指数) | P7 (均匀度指数) | P9 (丰度) | P2 (景观总面积) | Wi |
---|---|---|---|---|---|---|---|---|
P5 | 1 | 2 | 4 | 5 | 6 | 7 | 9 | 0.3738 |
P4 | 1/2 | 1 | 3 | 4 | 5 | 6 | 8 | 0.2641 |
P3 | 1/4 | 1/3 | 1 | 2 | 3 | 4 | 6 | 0.1366 |
P6 | 1/5 | 1/4 | 1/2 | 1 | 2 | 3 | 5 | 0.0938 |
P7 | 1/6 | 1/5 | 1/3 | 1/2 | 1 | 2 | 4 | 0.0638 |
P9 | 1/7 | 1/6 | 1/4 | 1/3 | 1/2 | 1 | 3 | 0.0440 |
P2 | 1/9 | 1/8 | 1/6 | 1/5 | 1/4 | 1/3 | 1 | 0.0239 |
表3. 准则层B2 (组织力)与其相关因子P的判断矩阵
备注:经检验λmax = 7.3144,其一致性指标CI = 0.0524, RI = 1.32, CR = 0.0397 < 0.1。
B3 | P6 (多样性指数) | P1 (NDVI) | P4 (聚合度指数) | P7 (均匀度指数) | P8 (斑块类型面积) | Wi |
---|---|---|---|---|---|---|
P6 | 1 | 3 | 5 | 7 | 9 | 0.5028 |
P1 | 1/3 | 1 | 3 | 5 | 7 | 0.2602 |
P4 | 1/5 | 1/3 | 1 | 3 | 5 | 0.1344 |
P7 | 1/7 | 1/5 | 1/3 | 1 | 3 | 0.0678 |
P8 | 1/9 | 1/7 | 1/5 | 1/3 | 1 | 0.0348 |
表4. 准则层B3与其相关因子P的判断矩阵
备注:经检验λmax = 5.2426,其一致性指标CI = 0.0607, RI = 1.12, CR = 0.0542 < 0.1。
B4 | P6 (多样性指数) | P9 (丰度) | P3 (分形维数) | Wi |
---|---|---|---|---|
P6 | 1 | 3 | 5 | 0.6334 |
P9 | 1/3 | 1 | 3 | 0.2605 |
P3 | 1/5 | 1/3 | 1 | 0.1061 |
表5. 准则层B4与其相关因子P的判断矩阵
备注:经检验λmax = 3.0386,其一致性指标CI = 0.0193, RI = 0.58, CR = 0.0333 < 0.1。
A | B1 (活力) | P1 (组织力) | P4 (弹力) | P7 (生态系统服务) | Wi |
---|---|---|---|---|---|
B1 | 1 | 3 | 5 | 7 | 0.5578 |
B2 | 1/3 | 1 | 3 | 5 | 0.2634 |
B3 | 1/5 | 1/3 | 1 | 3 | 0.1219 |
B4 | 1/7 | 1/5 | 1/3 | 1 | 0.0569 |
表6. 目标层A与准则层B的判断矩阵
备注:经检验λmax = 4.1185,其一致性指标CI = 0.0395, RI = 0.90, CR = 0.0439 < 0.1。
λ max = ∑ i = 1 i = n B W i n W i (2)
最终各参评因子的总权重由公式(3)算得出(表7)。
A | B1 | B2 | B3 | B4 | 总权重Wi |
---|---|---|---|---|---|
0.5578 | 0.2634 | 0.1219 | 0.0569 | ||
NDVI P1 | 0.6434 | - | 0.2602 | - | 0.3906 |
景观总面积P2 | 0.2828 | 0.0239 | - | - | 0.1640 |
分形维数P3 | - | 0.1366 | - | 0.1061 | 0.0420 |
聚合度指数P4 | - | 0.2641 | 0.1344 | - | 0.0859 |
蔓延度指数P5 | - | 0.3738 | - | - | 0.0985 |
多样性指数P6 | 0.0738 | 0.0938 | 0.5028 | 0.6334 | 0.1633 |
均匀度指数P7 | - | 0.0638 | 0.0678 | - | 0.0251 |
斑块类型面积P8 | - | - | 0.0348 | - | 0.0042 |
丰度P9 | - | 0.0440 | - | 0.2605 | 0.0264 |
表7. 最终各评价指标的总权重
W i = ∑ j = 1 3 B j ⋅ P i (3)
B j ——表示B层对于A层各相关因子权重;
P i ——表示P层对于B层各相关因子权重。
2) 层次分析法确定斑块类型贡献度
根据全国第二次湿地普查,本研究将湿地分为四种斑块类型:分别是湖泊、河流、沼泽化草甸、植被四类。沼泽化草甸是区域中最具有生物多样性的草地生态系统,所以它的贡献度最大。湖泊和河流作为水体一方面直接为人类提供动植物资源,另一方面在维持生物多样性、生态平衡、涵养水源、净化水质、调节小气候、保持水土、补充地下水均起到重要作用。湖泊是源头,河流是分支,因此湖泊的贡献度是大的,接下来是河流。植被在土壤形成和保护还有原材料价值上比较占优势,同时也对气体调节、水源涵养和生物多样性有所贡献,但不及前三个类型。本研究依旧运用层次分析法计算各斑块类型的权重值,如表8。
A | 沼泽化草甸 | 湖泊 | 河流 | 植被 | Wi |
---|---|---|---|---|---|
沼泽化草甸 | 1 | 3 | 5 | 8 | 0.5578 |
湖泊 | 1/3 | 1 | 3 | 6 | 0.2658 |
河流 | 1/5 | 1/3 | 1 | 4 | 0.1284 |
植被 | 1/8 | 1/6 | 1/4 | 1 | 0.0480 |
表8. 目标层A与准则层的判断矩阵
备注:经检验λmax = 4.0602,其一致性指标CI = 0.0201, RI = 0.90, CR = 0.0223 < 0.1。
生态健康性随时间、空间的变化而改变,与自然因子和人为因子相关,自然因素和人为因素中的某个因子在不同时间和空间范围内作用强度不同,在一定的时空背景下,每一因素的改变,或二者影响的叠加,都会导致生态系统的变化,对湿地生态系统健康性产生不同的影响,导致生态要素和生态系统的结构和功能变化,进而影响湿地生态系统健康性。所以,健康性的评价应是在一时间段内的变化情况,是一个动态的评价,可随时间的推移,呈现一定的波动。
利用综合指数法计算三江源区湿地生态系统健康状态综合指数E。其他学者是利用专家打分确定其评分值并做归一化处理后分别乘其权重值之和得到总分值用来计算综合健康值 [
E = ∑ i = 1 n V i ⋅ W i (4)
对于景观水平的景观指数,表6中值可直接运用;而对于斑块类型以及斑块水平的景观指数,则需先用表4和表5的值与各斑块类型的权重(见表7)相乘得到加权值再得以运用。其结果可见表9。
景观指数 | 景观类型面积(CA) | 聚合度指数(AI) | 分形维数(PARA) | 归一化植被指数(NDVI) |
---|---|---|---|---|
加权值 | 53267 | 50.3537 | 113.4977 | 0.1547 |
表9. 部分景观格局指数加权值
用公式3-5计算三江源区玛多县生态健康值结果,可得CEL = 31,921.39023。
对于三江源湿地生态健康评价,本研究选用1988年的值作为对照,因为1988年湿地生态状况较为良好,可供研究参考。
同法计算出1988年湿地景观指数及生态系统健康状态综合指数E = 31,729.2171,1988年的景观指数数据与2018年数据对比见表10。
景观类型面积 | 聚合度指数 | 分形维数 | 景观总面积 | 蔓延度指数 | 丰度 | 香农多样性 | 香农均匀度 | 归一化植被指数 | |
---|---|---|---|---|---|---|---|---|---|
1988年 | 50,263 | 50.9718 | 112.446 | 192,078 | 75.5302 | 4 | 0.5363 | 0.409 | 0.0997 |
2018年 | 53,267 | 50.3537 | 113.4977 | 193,176 | 75.5576 | 4 | 0.5738 | 0.4139 | 0.1547 |
表10. 1988年与2018年数据对比表
景观的异质性是一个描述景观格局的非常重要的属性,主要变现在空间分布的异质性和景观类型的多样性这两个方面。空间分布异质性常可以用景观的聚合度指数来描述,而景观多样性则可以用景观的多样性和均匀度指数来说明。从表10可以看出,研究区的景观多样性增大,由1988年的0.5363增大到2018年的0.5738这种变化趋势表明湿地景观中的优势景观优势越来越弱,景观的多样性呈上升趋势;景观均匀度指数由0.409升高到0.4139这说明湿地景观的均匀度指数也在不断升高。蔓延度指数则由1988年的75.5302下降到2018年的75.5576。湿地景观的多样性指数和均匀度指数的增大,以及蔓延度指数的减小,说明了湿地景观的异质性程度正在增加,整体景观中各景观类型所占的比例趋于平衡和近似,优势景观类型受到削弱。景观多样性指数和均匀度指数的增加能够在一定程度上维持和提高湿地的生态稳定性,使其生态系统的抗干扰能力有所提高。不过聚合度指数的减少,使得湿地景观的斑块连接度有所降低,分离性加大,景观在空间上的分布呈松散趋势。这样将会导致湿地整体景观的能量流动以及物质交换受到障碍,湿地生物种群的繁衍和迁移也会受到影响。湿地景观的生态系统过程和生态功能在一定程度上受斑块形状的影响。总体来讲,各个景观指数的值变化不大,而且也基本都在健康的范围内波动,说明湿地的健康状况良好。通过综合健康值的增加可知,湿地整体的健康度还有所提高。
本论文以玛多县为例,针对三江源区湿地生态环境特征及整体功能特征进行分析,据此建立评价指标体系,采用层次分析法与综合指数法结合的评价方法,对三江源区湿地的健康性进行了初步研究,结果表明:
1) 三江源湿地生态环境特性,从活力、组织力、弹力、生态系统服务四个方面确定了NDVI、景观总面积、分形维数、聚合度指数、蔓延度指数、多样性指数、均匀度指数、斑块类型面积、丰度这9个影响因子,作为三江源湿地生态健康评价指标。
2) 本论文通过运用层次分析法确定各个评价指标的权重,得到各指标的总权重值:NDVI (0.3906);景观总面积(0.1640);多样性指数(0.1633);蔓延度指数(0.0985);聚合度指数(0.0859);分形维数(0.0420);丰度(0.0264);均匀度指数(0.0251);斑块类型面积(0.0042);这表明影响三江源湿地生态环境健康性的主要指标有NDVI、景观总面积和多样性指数。
3) 对三江源湿地景观格局进行了分析和讨论:1988年和2018年这两年中,湖区湿地景观破碎度增大,景观异质性增强,景观的形状复杂化和不规则化。
常 敏. 三江源湿地生态健康状况的评价Analyzing and Assessing on the Ecological Environment Health of Wetlands in Sanjiangyuan[J]. 可持续发展, 2019, 09(04): 411-419. https://doi.org/10.12677/SD.2019.94072