-
景观格局指空间格局,即大小和形状各异的景观要素在空间上的排列和组合,包括景观组成单元的类型、数目及空间分布与配置模式,景观格局是景观异质性的体现,又是各种生态过程在不同尺度上作用的结果。景观空间格局研究是进行景观功能和动态研究的基础,它影响着生物动态和物质交换等多方面的生态过程[1]。随着GIS技术和景观结构分析软件Fragstats的应用,以此为手段开展的景观空间结构及其动态变化研究覆盖到建设规划[2,3]、物种栖息地质量及生境利用[4-7]、区域景观质量[8,9]等众多研究领域,为研究者从大尺度掌握区域生态质量、制定宏观管理策略提供了有力的技术手段。
大熊猫是全球生物多样性保护的旗舰物种,目前只分布于我国的四川、陕西、甘肃部分地区[10]。景观格局对大熊猫种群的生存繁衍影响重大,随着社会化进程加快导致大熊猫栖息地景观格局改变,结合景观生态学原理的栖息地景观格局研究受到关注[10],如白文科等人将家域模型与景观格局分析结合,以大熊猫生物学特性评估栖息地的空间格局变化特征[11]。目前大熊猫栖息地景观格局研究方法侧重于景观与“3S”技术相结合,用景观格局指数进行定量化分析,研究集中于秦岭、卧龙等区域[10]。
岷山山系是四川省大熊猫种群及栖息地的主要分布区,目前关于岷山山系大熊猫栖息地的研究很少,陈家德等采用景观生态学的斑块理论分析了土地岭大熊猫栖息地走廊带的植被分布格局、斑块结构、廊道特征及主要影响因素,提出建立土地岭走廊带的关键措施[12];张超基于“3S”技术完成了平武县小河沟自然保护区潜在栖息地的判别与预测,圈定该保护区的扩建范围[13];余翔等对岷山山系南端大熊猫栖息地的研究认为龙溪虹口-白水河-九顶山一线的大熊猫栖息地存在破碎化现象,建议继续加强大熊猫栖息地保护,增进各相邻保护区的交流和数据共享[14]。目前尚无针对岷山山系大熊猫栖息地景观格局方面的报道,本文以四川境内岷山山系大熊猫栖息地的景观格局为研究对象,利用Arc GIS和Fragstats软件对景观整体以及北段、中段、南段3个分区的景观格局进行分析,揭示大熊猫栖息地的景观格局特征并提出保护对策,对指导大熊猫国家公园建设实践和岷山山系大熊猫栖息地保护与质量提升具有重要意义。
HTML
-
岷山山系大熊猫栖息地主要分布于四川境内,根据研究区的地形地貌划分为北段、中段、南段3个区域以开展景观格局整体特征和分区比较研究,分区结果见图1和表1。
研究分区 面积/hm2及比例/% 行政区域 分区简介 北段 438773.67(50.16) 若尔盖县,九寨沟县,青川县,平武县北部,松潘县北部 呈西北-东南走向,包括九寨沟、白河、王朗、唐家河、黄龙、小河沟、贡杠岭、勿角、龙滴水等已建自然保护区 中段 299528.83(34.24) 平武县中南部,松潘县南部,北川县中北部,茂县北部 呈由北向南走向,包括雪宝顶、小寨子沟、宝顶沟北部、白羊、片口等已建自然保护区 南段 136487.98(15.60) 北川县东南部,茂县东部,安县,绵竹市,什邡市,彭州市,都江堰市,汶川县东部 位于四川盆地向青藏高原的过渡地带,包括龙溪-虹口、白水河、九顶山、千佛山、宝顶沟南部、蓥华山等已建自然保护区。 合计 874790.48 13个县(市) Table 1. Regionalization of giant panda habitats in the Minshan Mountains, Sichuan province
-
利用地理信息系统软件Arc GIS 10.5平台,以研究区各县(市)近年的林地变更、森林二调等数据为基础,结合野外调查和研究区2019年卫星影像对林地数据进行归类和校准。参考《土地利用现状分类》(GB/T21010-2017)分类体系,根据研究区土地利用现状、大熊猫栖息地植被类型的生态功能,将研究区景观划分为针叶林、阔叶林、人工经济林、灌丛、草地、湿地、耕地、建设用地8种类型,制作研究区景观格局分布现状图(见图1)。其中:针叶林包括云杉、冷杉、落叶松、高山松、圆柏、铁杉、柏木等针叶树种为优势形成的植被景观;阔叶林包括白桦、红桦、糙皮桦、川滇高山栎、槭树、山杨、椴树、鹅耳枥、黑壳楠、青冈等多种常绿、落叶阔叶树种为优势形成的植被景观;人工经济林包括柳杉、杉木、刺槐、辐射松、日本落叶松、水杉等人工起源的造林地、经济林木和果木林地;灌丛包括白夹竹、箭竹、珍珠梅、柳、沙棘、杜鹃、蔷薇、悬钩子、小檗、马桑、胡颓子、高山栎等为优势形成的常绿、落叶灌丛景观;草地包括禾本科、莎草科等为优势形成的高山草甸和次生草丛群落;湿地包括湖泊、河流、溪沟、沼泽地等湿地景观;耕地包括区内的旱地、水田农业种植区;建设用地包括城镇、村庄、聚居点、公路、工矿设施等人工建筑区域。上述各景观类型根据生境自然度、大熊猫利用度、受干扰程度等因子合理归类划分,代表了研究区不同的生态价值和服务功能,彼此镶嵌分布形成研究区的景观基本格局。经抽样评估,本研究景观格局分布图抽样精度达到80%以上,由于研究区范围广、面积大,这一精度基本满足研究要求。
-
在Arc GIS 10.5中将shape格式景观图层转化成GeoTIFF Grid(.jpg)格式数据(转化的栅格大小为12 m×12 m)并导入景观结构定量分析软件Fragstats(v4.2.1),在景观类型(Class Level)和景观(Landscape Level)两个尺度水平计算景观结构指数(部分指数计算的搜索阈值为2 000 m)。本研究选择的景观格局指数[18]、计算公式、单位及意义参见相关文献[3,17,18]及Fragstats软件手册[19]。
2.1. 研究分区
2.2. 景观制图及建立分类体系
2.3. 景观结构指数选取及计算
-
研究区景观基本组成见表2,结果显示:(1)针叶林景观类型在研究区分布最广、面积最大,其AREA_MN(81.84 hm2·块−1)、ED(31.78 m·hm−2)和LPI(5.432%)3项指标与阔叶林相近,高于其他景观类型;PD(0.50个·km−2)处于中等偏低水平;阔叶林也是研究区的重要景观类型,除PLAND(32.41%)和NP(3453块)明显低于针叶林外,其他指标与针叶林基本相近,针叶林和阔叶林是研究区大熊猫栖息地的主体和主要活动区域。(2)人工经济林作为人工起源的景观类型在研究区PLAND仅4.22%,可见大熊猫栖息地内的人工林比例较低。(3)灌丛景观类型分布面积在研究区排第三位,PLAND为17.84%,但NP远高于其他景观类型,达15381块,导致PD和ED达研究区最高的1.76个·km−2和33.13 m·hm−2。(4)草地和湿地景观在研究区PLAND分别为2.71%和1.10%,是分布面积最小的自然景观类型,它们的NP相对较大,导致PD、ED、LPI都偏低。(5)耕地和建设用地属于人工起源性质,PLAND<1%,其他各项指标计算值也偏低,表明大熊猫栖息地自然度较高。研究区大熊猫栖息地内的针叶林、阔叶林、灌丛、草地、湿地5类自然景观类型的PLAND之和达94.94%,表明自然景观类型在研究区优势明显,其中大熊猫利用度最高的针叶林和阔叶林景观的PLAND之和为73.30%。
景观类型 斑块类型面积
CA /hm2斑块类型面积比
PLAND/%斑块数量
NP /个平均斑块面积
AREA_MN/(hm2·块−1)斑块密度
PD/(个·km−2)边缘密度
ED/(m·hm−2)最大斑块指数
LPI/%针叶林 357632.76 40.88 4370 81.84 0.50 31.78 5.432 阔叶林 283545.85 32.41 3453 82.12 0.39 23.76 7.035 人工经济林 36921.92 4.22 2655 13.91 0.30 6.71 0.081 灌丛 156074.67 17.84 15381 10.15 1.76 33.13 0.941 草地 23668.07 2.71 4300 5.50 0.49 6.22 0.125 湿地 9603.03 1.10 6656 1.44 0.76 5.42 0.177 耕地 5504.83 0.63 1638 3.36 0.19 2.16 0.020 建设用地 1839.36 0.21 10972 0.17 1.25 2.46 0.012 Table 2. Analysis of landscape components of the study area
研究区北段、中段、南段的景观基本组成见表3,结果显示:(1)针叶林的CA、PLAND、NP、AREA_MN在研究区北段达到最大值,在中段、南段依次递减,呈由北向南递减趋势;而PD、ED呈由北向南递增趋势,表明随着平均海拔的降低,针叶林景观的破碎化程度升高。(2)阔叶林在研究区北段、中段的CA、PLAND、AREA_MN和LPI相近;由于研究区南段的阔叶林分布面积大幅降低但斑块数量相对较多,导致南段分布的阔叶林AREA_MN明显降低,而PD和ED指标较研究区中北段明显上升。(3)北、中、南3区段的人工经济林AREA_MN相近,但研究区南段的人工经济林PLAND(7.85%)、PD(0.71个·km−2)和ED(14.39 m·hm−2)有所升高。(4)灌丛景观由北向南PLAND、AREA_MN、PD、ED值呈上升趋势。(5)草地、湿地、耕地、建设用地由于分布面积很小,在北、中、南3区分布没有明显变化。综上,针叶林、阔叶林、人工经济林和灌丛这4类面积较大的景观类型表现出由北向南破碎度增加的趋势,这与研究区南段平均海拔更低、人为活动干扰更强有关。
景观类型 景观分区 斑块类型面积
CA /hm2斑块类型面积比
PLAND /%斑块数量
NP/个平均斑块面积
AREA_MN/(hm2·块−1)斑块密度
PD/(个·km−2)边缘密度
ED/(m·hm−2)最大斑块指数
LPI/%针叶林 北段 206003.95 46.94 2060 100.00 0.47 29.24 8.272 中段 118337.41 39.50 1433 82.58 0.48 30.33 16.773 南段 40511.61 29.68 795 50.96 0.58 36.18 8.466 阔叶林 北段 139649.49 31.82 1379 101.27 0.31 22.48 14.035 中段 99295.04 33.15 859 115.59 0.29 18.16 15.945 南段 43476.83 31.85 1166 37.29 0.85 39.04 5.704 人工经济林 北段 15279.25 3.48 1083 14.11 0.25 5.85 0.144 中段 10931.27 3.65 575 19.01 0.19 4.48 0.222 南段 10712.82 7.85 973 11.01 0.71 14.39 0.667 灌丛 北段 55933.47 12.75 5901 9.48 1.34 24.99 0.739 中段 61042.31 20.38 5835 10.46 1.95 34.90 2.734 南段 33945.87 24.87 2882 11.78 2.11 46.86 2.097 草地 北段 14176.48 3.23 2863 4.95 0.65 7.91 0.121 中段 6604.98 2.20 819 8.06 0.27 3.37 0.366 南段 2117.20 1.55 352 6.01 0.26 5.21 0.119 湿地 北段 3293.04 0.75 2627 1.25 0.60 5.10 0.066 中段 797.34 0.27 3301 0.24 1.10 3.85 0.011 南段 5451.32 3.99 659 8.27 0.48 9.66 1.133 耕地 北段 3299.50 0.75 889 3.71 0.20 2.48 0.026 中段 2085.83 0.70 586 3.56 0.20 2.44 0.059 南段 118.58 0.09 127 0.93 0.09 0.52 0.007 建设用地 北段 1196.91 0.27 5019 0.24 1.14 2.93 0.011 中段 474.54 0.16 4384 0.11 1.46 2.08 0.035 南段 171.92 0.13 1544 0.11 1.13 1.80 0.006 Table 3. Analysis of landscape components in the north, middle and south section of the study area
-
研究区在景观类型水平的格局特征分析见表4。主要结果有:(1)从斑块形状和边缘褶皱程度综合分析,针叶林和阔叶林的斑块形状最复杂,其次为灌丛和人工经济林,建设用地SHAPE_MN和FRAC_MN值最低,表明其形状最简单规整,这与其人工起源性质有关。(2)针叶林(SPLIT=128.32)和阔叶林(SPLIT=110.74)的斑块分散程度最低;其次为灌丛、人工经济林、草地和湿地;建设用地分散程度最大(SPLIT=27363315.97)远高于其他景观类型,这与建设用地在研究区零散、孤立分布有关。由于研究区面积大,在设定阈值2 000 m时各景观类型的连接度指数都在1以内,连接度水平普遍偏低,与设定阈值较小有关。(3)针叶林和阔叶林的AREA_MN和PROX-MN均大于其他景观类型,表明它们的破碎化程度最低;其次是灌丛和人工经济林,破碎化程度处于中等水平;草地、湿地、耕地和建设用地的AREA_MN和PROX-MN都很小,表明它们的破碎化程度最高。(4)AI、IJI反映同类斑块聚集与分散的空间分布特征,针叶林、阔叶林、人工经济林、灌丛、草地的AI值均较高,表明这些景观类型的斑块分布相对集中;湿地、耕地、建设用地AI值较低、IJI值较高,表明这3类景观斑块在景观全域零散分布,且与其他景观的混合度较高。(5)由于湿地、耕地、建设用地斑块形状和边缘褶皱较简单,PAFRAC最接近1.5,表明这3类景观类型的稳定性较差,镶嵌结构简单;针叶林、阔叶林、人工经济林、灌丛、草地的PAFRAC值均远离1.5,表明这些景观类型的稳定性较高。
景观类型 平均形状指数
SHAPE_MN平均斑块分维数
FRAC_MN分离度指数
SPLIT连接度指数
CONNECT/%平均邻近度指数
PROX_MN/%聚集度指数
AI/%散布与并列指数
IJI/%周长面积分维数
PAFRAC针叶林 1.914 1.110 128.32 0.418 70383.86 97.56 61.51 1.277 阔叶林 1.942 1.110 110.74 0.405 40731.85 97.73 77.58 1.274 人工经济林 1.818 1.108 117431.41 0.699 262.16 95.09 65.70 1.253 灌丛 1.832 1.117 4883.88 0.332 1554.06 94.26 65.25 1.292 草地 1.795 1.117 210307.53 0.550 135.42 92.69 62.65 1.292 湿地 1.738 1.109 249911.04 0.834 60.41 85.08 79.47 1.479 耕地 1.738 1.108 4540926.87 0.824 50.80 89.48 89.33 1.315 建设用地 1.317 1.072 27363315.97 0.927 7.22 64.94 91.11 1.569 Table 4. Analysis of landscape pattern index at type level of the study area
研究区的针叶林、阔叶林作为大熊猫栖息地的主体和大熊猫痕迹点的主要分布区,它们的斑块形状、边缘褶皱较复杂,斑块分散和破碎程度低、聚集程度高,功能优于其他景观类型,对研究区的景观格局总体特征具有决定性作用。
-
研究区整体及各分段景观水平的格局指数分析见表5。主要结果有:(1)由北向南的分段梯度上,AREA_MN和CONTAG指数呈下降趋势,表明随着平均海拔降低、干扰强度升高各景观类型整体的聚集度降低,景观中连通度极高的优势斑块数量减少,景观破碎化程度加剧;(2)研究区各分段景观斑块边缘的平均褶皱程度无明显变化,FRAC-MN值在1.098—1.111之间波动;(3)研究区北段IJI=73.07值最高,表明北段各景观类型间的相邻程度和混合程度最高;(4)研究区南段的SPLIT=51.223,表明景观中同类型斑块更加分散,景观破碎度高于北段和中段,与CONTAG指数分析结果一致;(5)由北向南SHDI和SHEI指数依次升高,表明南段景观中的优势控制性组分减弱,景观组分更趋均衡,因为人类活动影响下人工起源的景观组分增多,所以由北向南景观多样性和均匀度都增加。
景观分区 平均斑块面积
AREA_MN/(hm2·块−1)平均斑块分维数
FRAC_MN散布与并列
指数IJI/%蔓延度指数
CONTAG连接度指数
CONNECT/%香浓多样性
指数SHDI香浓均匀地
指数SHEI分离度指数
SPLIT研究区-北段 20.11 1.106 73.07 64.67 1.01 1.299 0.625 27.106 研究区-中段 16.84 1.098 62.39 64.28 1.88 1.323 0.636 15.165 研究区-南段 16.06 1.111 66.31 58.60 2.68 1.479 0.711 51.223 研究区整体 17.70 1.104 69.71 62.74 0.56 1.364 0.656 58.666 Table 5. Analysis of landscape pattern index at landscape level of the study area