神农架林区景观格局时空演变及其驱动力分析

曹嘉铄，邓政宇，胡远东，吴 妍

（1. 东北林业大学 园林学院，黑龙江 哈尔滨 150040；2. 东北林业大学 黑龙江省寒地园林植物种质资源开发与景观生态修复重点实验室，黑龙江 哈尔滨 150040）

景观格局反映景观要素在空间结构上的特征，可以从不同尺度讨论要素间的相互作用与生态效应[1]，是景观异质性的外在表征。目前，国内外学者对于景观格局的研究主要分为2个方面，即景观格局的空间质异性与时间异质性研究[2]。景观异质性在空间上的表达往往是复杂而连续的，因此以遥感(RS)和地理信息系统(GIS)等方法结合景观格局指数进行讨论，有助于实现格局演变的可视化。在景观生态学中，认为土地即景观[1]。据以往研究来看，土地利用变化可以反映景观格局演变过程，将土地利用同生态环境与自然资源系统相结合进行探讨，成为国内外研究的热点问题[3−4]。伴随3S技术的快速发展，土地利用的空间分异研究取得了长足的进展，土地利用结构分析已应用到城市及山区用地空间扩展[5]、景观组分在地形梯度上的分异[6]、湿地格局分异[7]等诸多方面。神农架林区位于中国地势二、三阶梯过渡地带，是中国重要的气候敏感区和生态环境脆弱地带[8]。20世纪六七十年代间，由于不合理开发导致区域生态环境遭严重破坏，林区政府开始探索绿色可持续发展模式[9]。本研究将基于多源数据，对神农架林区土地利用结构的组成及其时空演变特征与驱动力进行探讨，以探究20世纪90年代以来森林保护工程的成效，同时为林区的生态系统健康评价研究提供参考。

1 研究区域概况与数据来源

1.1 研究区域概况

神农架林区地处鄂西北山区，是中国唯一以“林区”命名的行政区划，总面积约3 215.80 km2。林区地处 31°15′～31°57′N，109°56′～110°58′E，西南部以近东西向的山脉为主，其中最高峰神农顶高 3 105.4 m，是华中地区最高点[10]。神农架林区是中国东西植物区系的交汇地带和南北植物区系的过渡地带，其亚热带森林生态系统是全球中纬度地区唯一保存完好的地区，拥有世界上最为丰富的生物多样性[11]，是具有全球意义的基因库，被称为“绿色奇迹”，林区内包括森林、灌丛、草甸和湿地等各类生态系统，具有调节气候、涵养水源、保持水土等重要生态环境服务价值[12]。近年来，以旅游业为主体的第三产业发展迅速，2018年全区旅游收入达57.29亿元，林区内建设用地规模大幅扩张，生态风险有所升高。

1.2 数据来源与处理

以1996、2007、2018年Landsat遥感影像为数据基础，影像空间分辨率为30 m，选取每个年份植被状态良好的2景影像，影像的云量均低于2%；相关规划资料详见表1。

表1 数据资料来源统计Table 1 Data source statistics

在ENVI 5.3.1软件中对遥感影像进行预处理，通过对遥感影像的目视判读和对研究区域的实地调查，采用最大似然法对遥感影像进行监督分类。根据研究目的，结合国家标准GB/T 21 010−2017《土地利用现状分类》，将神农架林区土地划分为林地、草地、建设用地、耕地和水域等五大类。利用Erdas Imagine 2015软件计算解译精度，结果显示3个年份Kappa系数都达85%以上，符合研究需要。

2 研究方法

2.1 土地利用转移矩阵

土地利用转移矩阵用于描述不同土地利用类型之间的转移情况，能反映研究时间段内土地利用类型转移的数量与方向[13−15]。其通用公式为

式(1)中：Aij为土地面积；n为土地利用类型数；i、j为转移前后的土地利用类型。

2.2 景观格局分布重心模型

景观格局分布的重心模型用于研究景观格局在空间上的分布情况，而多年份数据的叠加则可以反映景观格局重心在不同时刻的迁移变化，是过程研究中重要的表征指标[16−17]。

重心坐标计算公式为：

重心距离迁移公式为：

式(2)～(4)中：Xt、Yt为研究区域内某一景观类型在时间t的重心坐标；Dti为t时刻研究区域内第i个斑块的面积；Xi和Yi表示研究区域的坐标；D为某一景观类型的重心在t1和t2时间段内移动的距离。

2.3 景观格局指数

景观格局指数是描述景观变化的一种定量研究方法，是景观格局结构组成与空间配置信息的集中反映[18−20]。选取8个对于景观结构变化较为敏感的指数：斑块数(NP)、最大斑块指数(LPI)、散布与并列指数(IJI)、聚合度指数(AI)、分离度指数(SPLIT)、景观分割指数(DIVISION)、内聚力指数(COHESION)与香农多样性指数(SHDI)，采用移动窗口法[21−22]，在Fragstats软件中每次移动1个栅格，以900 m大小的移动窗口，计算得到研究区域各指数的空间分布，分析神农架林区景观格局破碎化及多样性特征。

3 结果与分析

3.1 土地利用结构变化特征

从图1和表2可以看出：林地是研究区域最主要的用地类型，占林区总面积的90%以上，22 a间，林地面积累计增加39.34 km2，森林覆盖度呈现增长趋势；草地与耕地面积总体上呈现减少趋势；建设用地面积持续增加，近10 a间增速比较明显，2018年研究区域内建设用地面积与1996年相比增加43.55 km2，这与该地区城镇化建设和旅游的大力发展有关；水域面积累计增加6.75 km2。

图1 1996、2007、2018年神农架林区土地覆被类型分类示意图Figure 1 Land cover types in Shennongjia Forestry District basin in 1996, 2007, 2018

表2 1996、2007、2018 年神农架林区土地利用结构表Table 2 Land use structure in Shennongjia Forestry District basin in 1996, 2007, 2018

3.2 土地利用转移过程

由表3可知：1996−2007年，林地转移为其他类型景观的面积变化量最大，主要转移为耕地和草地，转出面积分别为49.29和23.47 km2，而转移为林地最多的景观类型也是耕地和草地，由此可知林地与耕地、草地间相互转移关系密切，90%以上的新增建设用地由林地和耕地转化而来。

表3 1996−2007 年、2007−2018 年、1996−2018 年神农架林区土地利用转移矩阵Table 3 Land use transfer matrix in Shennongjia Forestry District of 1996−2007, 2007−2018 and 1996−2018

2007−2018年，林地仍然是转移面积最大的景观类型，主要转移为耕地和建设用地，转移面积分别为17.38和15.31 km2，表明神农架林区在城镇化建设及旅游业发展过程中，仍存在以占用林地作为代价进行开发建设的现象；林地与耕地、草地间仍保持着密切的相互转移关系，而林地的转出面积减少37.05 km2，转入面积增加27.66 km2，林地的总面积增加趋势显著。

3.3 景观格局分布重心转移变化

通过ArcGIS计算，得到1996、2007和2018年神农架林区各个用地类型景观格局分布重心转移的距离、速度和方向(表4)。1996−2018年，研究区域内水域和建设用地重心迁移明显，水域重心在研究时段内持续向西南方向迁移，累计迁移11 745.64 m，其中，1996−2007年迁移速度较快，为829.11 m·a−1，主要是由于大九湖湿地恢复与水库修建导致的；建设用地重心向西南方向迁移10 752.54 m，说明神农架林区的建设向西南发展，尤其是2016年进入国家公园发展的新时代，其范围划定也与建设用地重心的迁移方向相印证；林地具有向东扩张趋势，其重心在研究时间内向东南方向移动5 594.52 m。

表4 景观类型转移方向与距离Table 4 Landscape type transfer direction and distance

3.4 景观格局指数变化分析

从同年份的同类型景观比较结果(表5)来看：林地斑块个数最少，最大斑块指数与内聚力指数均较高，景观分割指数小于0.2，说明林地斑块的完整性较好，景观破碎化程度低，林地为研究区域的景观基质，而建设用地和水域等类型则以斑块或廊道的形式镶嵌其中。

表5 神农架林区 1996−2018 年景观格局指数变化Table 5 Change of landscape pattern index from 1996 to 2018

从不同年份的同类型景观比较结果(表5)来看：林地的景观破碎化程度先增后减，草地、耕地与水域的景观破碎化程度与分散程度在研究时段内呈增加趋势；人类活动使建设用地不断增加并向外扩张，使其破碎化程度在总体上呈上升趋势，但随着城镇化的推进，建设用地景观类型内部也出现了一定程度的整合。

由图2可见：从整体上看，神农架林区景观的破碎度有所降低，趋于整合；从局部来看，东北部区域与中南部区域的破碎化程度加剧，建设用地等人工斑块的扩张导致部分区域景观破碎度升高，而景观多样性也随之呈现扩散趋势；而林区多年来封山育林、人工造林、退耕还林等保护工作的开展，使得部分重点保护区域内林地面积持续增加，在地表覆被类型中占据主导地位，因此这些区域内景观多样性会呈现一定的降低趋势。

图2 1996−2018年神农架林区景观破碎化与多样性空间分布示意图Figure 2 Landscape fragmentation and diversity index of Shennongjia Forestry District from 1996 to 2018

3.5 景观格局演变的人为驱动因素分析

3.5.1 经济因素 通过对神农架林区1996−2018年经济发展数据(图3)与景观格局演变特征的分析对比，可知在开发历史较久，受经济增长与人类活动干扰较大的地区，景观格局演变较为显著，如东北部的松柏镇是神农架林区政府的驻地所在，而中南部的木鱼镇是神农架自然保护区管理局的驻地，也是神农架旅游接待服务中心，这些区域的景观破碎化程度呈现出较为明显的升高趋势。近年来旅游业更是逐步成为神农架林区的支柱性产业，2018年林区的旅游总收入接近57.3亿元，约为地区生产总值的2倍。与此同时林区大力开展的旅游配套设施建设及房地产项目开发，使得人工斑块数量显著增加。22 a间林区建设用地面积骤增43.56 km2，所占比例由0.11%提升至1.45%。因此，经济因素对神农架林区景观格局的演变具有重要影响。

图3 1996−2018年神农架林区地区生产总值与年旅游总收入趋势图Figure 3 Trend of GDP and total tourism revenue change in Shennongjia Forestry District from 1996 to 2018

3.5.2 人口因素 由图4可见：神农架林区人口数量出现先增后减的趋势；神农坛、天生桥、官门山、大九湖在内的几大主要景区游客接待量逐年上涨，于2018年突破1 587.5万人次。为提高游客接待能力而进行的景区开发，使得区域景观类型发生显著变化，1996−2018年，林地累计向建设用地转移25.45 km2，局部空间与生态环境所承受的压力随之增大。林区于2014年开展“美丽乡村”建设，以村庄布局“减量、聚居”为总体要求，注重农村土地资源综合整治，因此耕地及建设用地中的住宅用地分布有所整合，但旅游服务用地的增加仍然造成了景观破碎化程度的升高。

图4 1996−2018年神农架林区本地人口数量与年游客接待量趋势图Figure 4 Local population and annual tourist reception change in Shennongjia Forestry District from 1996 to 2018

3.5.3 政策因素 土地利用结构及变化深受区域政策影响[23]。20世纪60−80年代，采取毁林开荒的开发方式追求经济效益，神农架林区的生态环境遭到严重破坏，森林覆盖率骤降至63.5%。神农架林区于2000年全面启动天然林资源保护工程和退耕还林还草工程[24]，研究时段内耕地向林地、草地转移面积超过65 km2；2003年确立“保护第一，科学规划，合理开发，永续利用”的发展方针，促进生态环境逐步改善，森林覆盖率达90%以上；2011年发布《神农架旅游总体规划》，加快推进神农顶、官门山和天生桥景区配套设施建设，区域内建设用地面积显著增长，同时也与建设用地重心向西南迁移相印证；2016年7月，神农架被正式列入世界自然遗产名录，同年11月，神农架国家公园管理局正式成立，标志着神农架的保护与管理工作步入国家公园新时代[9]。1996−2018年，林区内水域面积累计增加6.75 km2，其中湿地的生态恢复成果得益于大九湖国家湿地公园的建立，湿地科学监测与保护工作在政策的引导下逐渐步入正轨[25−26]。

3.6 景观格局演变的自然驱动因素分析

相比于社会经济因素，自然因子的作用在短时间内并非十分显著，但从长远时间看，气候等因子的变化将对区域景观格局产生重大影响。1996−2018年，神农架林区的林地面积显著增加，林地主要由耕地和草地转化而来，除去政策的引导作用，主要依靠自然更新因素完成转化。近年来，神农架林区年降水量及年平均气温均呈较为微弱的增加趋势[27]，小气候条件优越，加之林区内土壤有机质含量丰富，为耕地向林地的演替提供了有利条件。

4 讨论

1996−2018年，神农架林区森林覆盖度总体上呈增长趋势，这与前人研究的结果一致[28]。刘家琰等[29]利用SPOT-VEGETATION数据对林区1998−2013年植被覆盖度格局变化进行分析，发现植被覆盖度呈较好增长趋势，整体变化率为1.45%。水域景观面积显著增加，破碎度升高，湿地生态恢复稳中向好，这主要得益于湿地生态保护政策的有效实施，人为因素在过程中起到了决定性作用[24]，在保护优先的政策指导下，林区逐步走上良性循环的绿色发展道路。同时，也有学者对自然因素进行探究，张子一等[30]讨论了个别年份受极端气候与虫害等自然因素，导致常绿植物受灾使其覆盖度有所下降。

城镇化的持续推进与旅游产业的迅速发展，造成了人工斑块对林地、草地等自然景观类型的侵占，景区内旅游基础设施的兴建，对生态环境造成严重的干扰，增大生态压力。建设用地的扩张促使区域内景观破碎度与多样性的增加，生态空间与社会经济空间融合程度不断加深。

由于旅游业的发展需求，连通林区内外的交通建设逐年推进。截至2018年林区内公路里程数达2 048.5 km，其中各大景区范围内公路网建设较为密集，造成林区景观局部破碎化[31]，导致原本连续的生态廊道网络出现生态间隙[32]，阻碍生态流的正常流动，降低了生物迁移的安全性。神农架机场位于林区中西部的红坪镇，于2014年正式通航，机场的落成使得区域内建设用地与耕地面积骤增，分割并部分取代了原有大面积草地斑块，强烈的人工干扰可能会降低区域生态系统的稳定性。

5 结论

1996−2018年，神农架林区土地利用结构变化显著，森林覆盖度总体上呈现增长趋势，湿地生态恢复稳中向好。随着神农架林区对于生态发展政策方针的不断调整与优化，森林保护工程成效日益显现，同时也在逐渐改变生态系统的结构和功能。

景观破碎度在整个研究区域内呈现降低趋势，表明神农架林区景观总体上趋于整合，但局部破碎化程度加剧，林区东北部与中南部区域地表覆被逐渐由原本较为整合的林地景观转变为多种用地类型共存，林区建设用地的持续扩张，社会经济空间与生态空间相互融合。人为因素是引起景观格局变化的主导因素，近年来神农架大力开发的旅游产业造成了人工斑块对自然景观的侵占，其他用地类型向建设用地转化的趋势较明显，将对生态系统安全产生威胁。