人类足迹对东北虎豹国家公园景观破碎化的影响

张晓宇,宁晓刚,王 浩,刘娅菲 ,刘若文

1 辽宁工程技术大学测绘与地理科学学院,阜新 123000

2 中国测绘科学研究院,北京 100036

建立国家公园体制是新时代中国特色社会主义生态文明建设的重要内容,对于有效保护国家重要自然生态系统的原真性、完整性[1],促进人与自然和谐共生,具有极其重要的意义。尽管目前我国国家公园体制建设取得了相当大的进展[2],但依然普遍存在能源建设、矿产开采、居民生活、农业开垦、道路建设等人类活动,致使自然生境破碎、生物多样性降低,严重影响着国家公园的生态系统服务功能[3—4]。因此,量化评价国家公园范围内人类活动对自然生态系统的影响程度、识别空间分异规律显得尤为迫切且意义重大。

人类足迹(Human Footprint)是从人类活动压力变化角度来量化评价人类活动影响程度的一种方法[5],相较于人类活动强度指数[6]、人为干扰度[7]、陆地表层人类活动强度[8]等从土地利用状态变化角度出发的评价方法,人类足迹通过选取与人类活动直接相关的空间因子对人类活动进行表达,能够更加全面客观的反映人类活动影响程度。2002年,Sanderson等[9]首次提出了人类足迹指数(Human Footprint Index,HFI)模型,利用人口密度、土地利用变化、通达性及电力基础设施4个方面的9个评价指标进行影响力赋值和叠加计算,最终得到全球陆地范围内1km格网的人类足迹数据。该方法能够直观体现人类活动在区域上的空间分布差异,因此在生态保护区域中应用广泛。如Woolmer等[10]利用90m×90m分辨率的空间数据,在Sanderson的研究基础上引入工矿、大坝和其他基础设施等指标,制备了适用于北阿巴拉契亚/阿卡迪亚生态保护区的人类足迹数据,更加精细化的体现出生态保护区域中人类活动的复杂性。随后人类足迹方法在苏格兰的凯恩戈姆斯国家公园和洛蒙德湖与特罗萨克斯山国家公园[11]、厄瓜多尔东南部的罗汉松国家公园[12]等地均有应用。我国的研究学者分别将人类足迹指数模型应用于横断山区[13]、青藏高原[14]、台湾日月潭国家公园[15]等生态保护区,并通过分析保护区内外人类足迹的空间分布差异,评估建立保护区对于减少人类活动影响的有效性。综上,目前关于生态保护区人类足迹的研究均从评价指标的选择和细化赋值入手,然后将各指标层进行空间叠加,最后根据叠加得到的人类足迹指数值直接判断人类活动的影响程度,但并未探讨各人类足迹指标对自然生态系统的影响。

作为景观生态学的重要研究内容,景观破碎化(Landscape Fragmentation)是指由于受到外界干扰所导致的景观由单一、均质和连续的整体趋向复杂、异质和不连续的斑块镶嵌体[16],在一定程度上能够直观反映自然生态系统受到人类活动影响的程度。鉴于此,本文基于人类足迹理论,以东北虎豹国家公园为研究区域,以2019年为研究时间点,构建适地化的人类足迹指标体系,并借助景观格局分析、回归分析、GIS空间分析等方法,从全局和局部两个角度上深入探讨各人类足迹指标对东北虎豹国家公园及不同管控分区景观破碎化的影响范围及其空间分异特征,以期为东北虎豹国家公园生态修复与人类活动的合理管控提供参考依据。

1 研究区概况

东北虎豹国家公园(简称“虎豹公园”)体制试点正式启动于2017年1月,是我国首个由中央直管的国家公园,其园区位于吉林省和黑龙江省交界处(129°05′—131°19′E,42°39′—44°15′N),东部、东南部与俄罗斯豹地国家公园接壤,西南部隔图们江与朝鲜相邻,行政区划涉及吉林省延边朝鲜族自治州珲春、汪清、图们和黑龙江省牡丹江市东宁、穆棱、宁安6个县(市),总面积14612km2(图1)。

图1 东北虎豹国家公园地理位置Fig.1 Geographical location of NCTL-NPNCTL-NP：东北虎豹国家公园,Northeastern China Tiger and Leopard National Park

按照国土空间和自然资源用途管制要求,遵循生态系统完整性、原真性和保护发展协调性三大原则,虎豹公园试点区划分为核心保护区与一般控制区[17],分别占虎豹公园试点区总面积的52%和48%。其中,核心保护区是维护现有虎豹种群正常繁衍、迁移的关键区域,主要包括近年来监测确认的东北虎豹繁殖家域、定居区及活动频繁的区域。一般控制区是恢复虎豹适宜栖息地及探讨人与自然和谐发展的关键区域,该区域将东北虎豹种群扩散栖息地、迁移廊道等划为虎豹潜在栖息地；将部分乡镇驻地,以及人口数量较大的林场场部、村屯及周边农业耕种区划为人口聚集区。长久以来,虎豹公园内传统的生产生活方式使得人类活动空间与东北虎豹栖息地和潜在栖息地空间高度重叠,农业用地、居民点、工矿企业、道路铁路穿插,严重压缩和分割了东北虎豹的生存空间,虎豹栖息地破碎化问题突出。

2 研究数据与方法

2.1 数据来源及预处理

本文采用的数据包括虎豹公园试点区及管控分区矢量边界数据、地表覆盖数据、人口密度数据、夜间灯光数据。其中矢量边界数据来源于吉林省延边朝鲜族自治州汪清县林草部门；地表覆盖数据来源于2019年基础性地理国情监测成果数据,数据生产所使用的影像空间分辨率优于2m,并通过外业核实和质量检查保证数据质量；人口密度数据采用美国国家橡树岭实验室(https://landscan.ornl.gov/landscan-datasets)开发的LandScan全球人口动态统计分析数据库,空间分辨率为30″；夜间灯光数据来源于美国科罗拉多矿业大学EOG小组官网(https://eogdata.mines.edu/products/vnl/)提供的NPP/VIIRS年度产品,空间分辨率为15″。

为充分体现研究区内的景观异质性,以2019年《基础性地理国情监测内容与指标》(CH/T 9029—2019)为基础,并参考前人对于自然保护区土地利用分类方法[18—19],将种植土地中的水田、旱地和构筑物中的温室、大棚合并为耕地；将种植土地中的果园、其他经济苗木合并为园地；将林草覆盖拆分为林地和草地；将构筑物中除温室、大棚外的其他地类以及房屋建筑(区)、铁路与道路、人工堆掘地合并为建设用地,最终得到耕地、园地、林地、草地、建设用地、水域、荒漠与裸露地7种土地利用类型(图2),用于景观格局指数的计算。

图2 东北虎豹国家公园2019年土地利用覆盖分布Fig.2 Land use coverage distribution of NCTL-NP in 2019

人口密度、夜间灯光、重新分类的土地利用覆盖,以及基础性地理国情监测成果中的道路要素数据用于构建人类足迹指标体系。由于各数据源在分辨率、投影和空间范围等方面存在不一致,因此在开展研究之前,需要根据研究区范围统一进行投影转换和裁剪。此外,人口密度和夜间灯光数据还需要进行降噪等校正处理,并综合考虑各数据的分辨率,最终采用最邻近法生成100m×100m分辨率数据。

2.2 研究方法

2.2.1人类足迹指标体系的构建

本文在Sanderson等人[9]构建的人类足迹指数方法的基础上,结合虎豹公园内的实际人类活动影响情况[20],考虑数据的可获取性,最终选用人口密度、土地利用状况、工矿活动、交通可达性、夜间灯光5个与人类活动直接相关的指标来定量评价公园内的人类足迹,具体选取原则如表1所示。

表1 人类足迹指标选取原则Table 1 Principles for the selection of human footprint index

2.2.2人类足迹指标的离散化赋值

为便于统一比较各人类足迹指标之间的大小关系,本文对各指标的栅格数据进行重新赋值,赋值范围为0—10,值越大表示人类足迹影响程度越大,最终得到各人类足迹指标空间分布图(图3)。具体赋值方法如下：

图3 各人类足迹指标空间分布Fig.3 Spatial distribution of each human footprint index

(1)人口密度

本文参考Venter等[24]对人口密度的赋值方法,将人口密度大于1000人/km2的栅格赋值为10,对于小于1000人/km2的栅格按照对数方程进行计算赋值,具体计算方程如下：

popscore=3.333×log(popdensity+1)

(1)

式中,popscore表示该栅格数据重新赋值的分数,popdensity表示各栅格像元的人口密度值。

(2)土地利用状况

本文结合虎豹公园土地利用分类标准,参考段群滔等[27]的研究成果,为公园内7种地类匹配相应的影响值。其中,建设用地受人类活动影响最大,赋值为10；其次为耕地,影响力赋值为8；园地影响力赋值为6；草地影响力赋值为2；林地、水域及荒漠与裸露地的影响力赋值为0。

(3)工矿活动

提取基础性地理国情监测数据中的露天采掘场、尾矿堆放物及工业设施来表征工矿活动。本文设置工矿活动对生态系统的最大影响距离为10000m[10],将距离大于10000m的栅格赋值为0,对于距离小于10000m的栅格参考InVEST模型生境质量模块中的距离衰减函数[28—29]来表示工矿活动得分,具体计算公式如下：

(2)

式中,gkscore表示该栅格数据重新赋值的分数,ddensity表示经过欧氏距离计算的栅格数据值,dmax表示工矿活动的最大影响距离,即10000m。

(4)交通可达性

本文参考Woolmer等[10]对道路的赋值方法,结合公园内各等级道路及铁路的分布情况,将交通可达性分为5个数据层,同时根据两侧距离进行重分类赋值(表2)。最后将重分类后的5个数据层进行等权叠加,并采用Min-Max标准化法将叠加后的栅格结果进行标准化处理,使其分数范围与其他人类足迹指标统一,具体公式如下：

表2 交通可达性赋值Table 2 Traffic accessibility assignment

(3)

式中,dlscore表示叠加后栅格数据经过标准化处理的分数,dldensity表示叠加后的栅格数据值,dlmax-density表示叠加后栅格数据中的最大值,dlmin-density表示叠加后栅格数据中的最小值。

(5)夜间灯光

将经过预处理后的NPP/VIIRS夜间灯光数据中DN值等于0的栅格赋值为0,再利用分位数分级法将DN值大于0的栅格赋值为1—10。

2.2.3景观破碎化分析方法

综合考虑研究区域的景观类型及分布特点,本文在景观尺度上,从面积、形状、聚散性、多样性4个方面分别选取斑块密度(PD)、边缘密度(ED)、景观分离度(DIVISION)、香农多样性指数(SHDI)来表征研究区景观破碎化程度,各指数计算方法及生态学意义参见《景观生态学—格局、过程、尺度与等级》[30]。根据研究区范围大小,本文分别设定边长为300m、600m、900m、1200m、1500m、1800m、2100m、2400m、2700m的正方形移动窗口计算并生成景观指数栅格图。从图4中可以看出,各景观指数曲线在1200—1500m尺度趋于缓和。由于窗口设置过小将导致生成的栅格图像缺乏连续性；窗口设置过大则会导致生成的栅格图像模糊,综合考虑最终选取1500m为最佳研究尺度。

图4 景观指数随窗口尺度变化曲线Fig.4 Variation curve of landscape index with window scale

2.2.4人类足迹对景观破碎化影响分析方法

(1)全局水平影响分析方法

普通最小二乘(Ordinary Least Squares,OLS)模型是一种常用于探索解释变量与因变量之间全局关系的线性回归模型。本文以景观破碎化作为因变量,以5种人类足迹指标作为解释变量,建立回归方程并得到每个解释变量的回归系数,通过对比该系数大小确定各指标影响程度。其方程形式如下：

(4)

式中,yi是空间i位置的景观破碎化指数；β为二乘法空间截距；βk为第k项自变量(人类足迹指标)的回归系数；xik为第k项人类足迹指标在空间i位置的取值；εi为误差修正项。

(2)局部水平影响分析方法

1)空间自相关检验

在构建地理加权回归模型之前,需要对研究区域内各变量进行空间自相关检验,确定各变量值之间是否存在相互依赖性。本文采用应用较为广泛的全局Moran′sI指数来检验区域内各变量是否具有空间相关性。其计算公式如下：

(5)

式中,n为研究区域内空间单元的总数；xi和xj分别代表第i个和第j个空间单元上的变量值,wij是单元i和单元j之间的空间权重矩阵；S0是所有空间权重矩阵之和。

2)地理加权回归模型

地理加权回归(Geographically Weighted Regression,GWR)模型将数据的空间位置嵌入到回归方程中,通过建立空间范围内每个位置点的局部回归方程,来探索研究对象在某一尺度下的空间变化及相关驱动因素。本文通过构建GWR模型以获得不同人类足迹指标随空间地理位置变化而变化的回归系数,该系数代表了人类足迹指标对景观破碎化的影响程度。其方程形式如下：

(6)

式中,yi是空间i位置的景观破碎化指数；(ui,vi)为i位置点的空间坐标；β0和βk分别为空间i位置的截距和局部回归系数；xik为第k项人类足迹指标在空间i位置的取值；εi为误差修正项。

3 结果与分析

3.1 景观破碎化空间分布格局

为避免表征景观破碎化的4个指数存在信息重叠,故通过主成分分析法将4个指数综合成一组新的综合指标来表征景观破碎化,具体结果如表3所示。

表3 景观指数主成分分析结果Table 3 Results of principal component analysis of landscape index

其中KMO检验指数达到0.795,满足因子分析的条件,第一成分的方差贡献率达到了93.54%,有且只有一个主成分。进一步计算各指数的综合模型得分和归一化权重,最终得到景观破碎化综合指数表达式为：

LFCI=ZPD×0.2428+ZED×0.2547+ZDIVISION×0.2485+ZSHDI×0.2539

(6)

式中,LFCI为景观破碎化综合指数值,ZPD、ZED、ZDIVISION、ZSHDI为PD、ED、DIVISION、SHDI标准化后的数值。利用自然间断点法将LFCI分为五个等级并按照管控分区进行统计,以定量化表示各管控分区内的景观破碎化程度及空间分布差异(表4、图5)。

表4 各分区中景观破碎化程度面积占比/%Table 4 Percentage of area of landscape fragmentation in each zoning district

从表4、图5中可以看出,虎豹公园景观破碎化整体上呈现“北高南低”的空间分布格局,且以低、较低破碎化为主,但仍有8.96%的区域为较高、高度破碎化。其中,核心保护区内较高、高度破碎化地区面积占该区域总面积的4.08%,主要分布于公园东南部,呈条带状分布；虎豹潜在栖息地内较高、高度破碎化地区面积占该区域总面积的5.09%,主要分布于该区域内的人类居住点周围,呈点状分布；人口聚集区内较高、高度景观破碎化地区则占到该区域面积的60.54%,且以罗子沟镇、老黑山镇、复兴镇、春化镇、春阳镇等乡镇驻地为中心向外围扩散。

3.3 人类足迹对景观破碎化影响分析

结合景观破碎化分析确定的1500m正方形格网最佳研究尺度成果,将研究区域划分为7006个网格单元,以分区统计5种人类足迹指标和景观破碎化综合指数。

3.3.1全局影响因素分析

OLS模型计算结果如表5所示,其中VIF为方差膨胀系数,一般认为当VIF大于10时,各变量间存在明显的多重共线性[31]。本文中5种人类足迹指标的VIF值均小于10,因此各指标不存在多重共线性问题。从显著性检验水平来看,5种人类足迹指标中只有夜间灯光未通过0.05水平显著性检验,这表明在全局水平上,LFCI与夜间灯光的相关关系不具有统计学意义。由于虎豹公园内以第一产业为主,社会经济水平较低,电力基础设施建设相对薄弱,夜间灯光集中分布在各乡镇驻地及周围地区,不足以对当地景观格局产生影响,导致夜间灯光对景观破碎化的解释能力不显著。

从回归系数来看,LFCI与其余4项人类足迹指标均呈正向相关,且土地利用状况之间的回归系数最大,达0.332。表明在全局水平上,土地利用状况影响力每增加1%,景观破碎化综合指数增加33.2%,其他3项人类足迹指标影响效应从高到低依次为交通可达性(0.289)、人口密度(0.224)、工矿活动(0.148)。

3.2.2局部影响因素分析

(1)空间自相关检验结果分析

通过OLS模型的计算结果可知,夜间灯光与景观破碎化之间的回归关系不显著,因此在进行GWR模型计算之前,首先将夜间灯光这一指标剔除,以保证GWR模型拟合结果的准确性。从图6可以看出,因变量景观破碎化综合指数LFCI和剩余解释变量(人口密度、土地利用状况、工矿活动、交通可达性)的Moran′sI指数散点图均沿一、三象限角平分线分布,聚集特征(高-高值或低-低值)较为明显,其结果在0.01极显著水平上通过检验(Z>2.58),说明上述各变量的空间分布存在较强的空间自相关性,因此可以构建地理加权回归模型进一步分析各人类足迹指标对景观破碎化的影响。

图6 各变量Moran′s I指数散点图Fig.6 Scatter plot of Moran′s I index for each variable

(2)各人类足迹指标对景观破碎化的定量影响

采用自然间断点法对GWR模型回归系数进行分级,得到LFCI与各人类足迹指标相关性空间格局(图7)。总体来看,LFCI与各人类足迹指标回归系数的空间格局呈现集聚状态,且以正向为主,表明各项人类足迹指标与LFCI均呈现出较强的正向相关性。其中,土地利用状况的回归系数高值区明显高于其他指标,为1.17—2.41,占公园全域面积的6.80%,且大部分位于核心保护区和虎豹潜在栖息地,仅有0.34%位于人口聚集区。空间上主要分布于公园南部、东部；人口密度的回归系数高值区(0.52—0.81)占全域面积的8.96%,其中有53.45%位于核心保护区,37.41%位于虎豹潜在栖息地,9.14%位于人口聚集区。空间上主要分布于公园中部罗子沟镇及复兴镇周边地区,东南部、北部和西北部的居民点附近有也少量分布；工矿活动的回归系数高值区(0.46-0.73)占全域面积的2.01%,且有44.22%位于人口聚集区。空间上主要分布于公园中部以罗子沟镇为中心的人口聚集区及其周围；交通可达性的回归系数高值区(0.47-0.86)占全域面积的13.03%,且有55.81%和37.40%位于核心保护区和虎豹潜在栖息地。空间上主要分布于公园东部春化镇及周边地区,东北部老黑山镇、大肚川镇及周边地区,在北部、西部也有一定分布。

图7 LFCI与各人类足迹指标相关性空间格局Fig.7 Spatial pattern of regression coefficients of LFCI and each human footprint indexLFCI：景观破碎化综合指数,Landscape Fragmentation Composite Index

(3)景观破碎化的主导人类足迹指标区域面积占比

从各回归系数的显著性水平来看,在7006个网格单元中,仅有1.78%未通过0.05水平显著性检验,主要分布于数据质量较差的边界地区。为确定虎豹公园景观破碎化的最大影响因素,分析通过显著性检验的剩余网格单元中回归系数正向最大值,并以管控分区进行统计分析。从图8中可以看出,土地利用状况影响为主的区域占公园全域面积的61.04%,其次为交通可达性(18.70%)、人口密度(14.72%)、工矿活动(4.26%),这也表明GWR模型与全局水平上OLS模型回归结果一致。从各管控分区来看,核心保护区与虎豹潜在栖息地中景观破碎化的各影响因素面积占比与公园全域基本一致,而在人口聚集区中,除土地利用状况外,工矿活动的影响范围较大,占人口聚集区面积的19.37%,其次为交通可达性(15.21%)、人口密度(14.69%)。

图8 LFCI最大影响因素在各管控分区中的面积占比/%Fig.8 Area share of LFCI maximum impact factors in each control zone

(4)景观破碎化的主导人类足迹指标空间格局分析

从图9中可以看出,景观破碎化以土地利用状况影响为主的区域最为广泛,主要分为三大集中连片区,其中东部连片区范围最大,涵盖了核心保护区大部,虎豹潜在栖息地的部分地区以及以复兴镇、老黑山镇、大肚川镇为中心的人口聚集区,而西北部和西部连片区则主要位于虎豹潜在栖息地中。受交通可达性影响为主的景观破碎化区域主要集中分布于公园东部春化镇及周边地区、东北部大肚川镇部分地区、西部春阳镇及周边地区以及公园北部地区,其中东部春化镇附近及公园北部地区的交通可达性对核心保护区影响较大,西部春阳镇附近的交通可达性对虎豹潜在栖息地影响较大。以人口密度影响为主导的区域呈现大集中、小分散的特征,主要分布在公园中部罗子沟镇、复兴镇附近地区,此外在东部边境地区、核心保护区南部及其他地区一些小型居民点周围均有分布。工矿活动的影响区域最小,主要分布在公园中部工矿业较为集中的罗子沟镇及周边地区,且对核心保护区东北部地区也存在一定影响。

图9 LFCI最大影响因素空间分布Fig.9 Spatial distribution of the most influential factors of LFCI

4 讨论与结论

4.1 讨论

我国东北虎豹目前主要分布于老爷岭、张广才岭和完达山3个区域,其中虎豹公园覆盖了老爷岭分布区的核心区域,对俄罗斯西南滨海边疆区域和我国老爷岭区域分布的东北虎豹种群的保护有着重大意义[32]。 总体来看,虎豹公园内景观破碎化程度较轻,尤其是核心保护区和虎豹潜在栖息地,主要是由于园区内以森林覆盖为主,斑块较为完整,但人类的过度开发活动导致部分地区景观破碎化严重,且存在从人口聚集区向虎豹潜在栖息地和核心保护区扩散的趋势,严重阻碍东北虎豹的生存繁衍和自由迁徙。本文通过分析各人类足迹指标与虎豹公园及其管控分区景观破碎化的空间关系,在一定程度上揭示了各人类足迹指标对虎豹公园景观破碎化影响程度的空间差异性,同时充分证明了人类的过度开发活动是虎豹公园景观破碎化加剧的重要驱动因素之一,对生态廊道建设工程存在巨大威胁。其一,随着近代以来人口的不断流入,各种规模的林场、牧场、耕地、园地以及大量用于人类生产生活的建设用地散布其中,导致森林破碎化严重,园区内落后的生产生活方式至今仍未发生改变,极大地威胁着东北虎豹栖息地的稳定以及种群数量的发展；其二,由于经济发展和国防建设需要,虎豹公园内由各等级道路组成的交通道路网交错纵横,总长度达5000余公里,特别是边境公路国道G331贯穿公园东部,过往车辆和潜在较高的交通可达性在一定程度上阻隔了中俄边境地区的虎豹种群向内陆迁徙的通道；其三,虎豹公园内人口相对较多且分散,其中500人以上的村屯有38个,其他规模较小的林场和农场人口则分散于其他地区,此外还有31个边境村屯分布于中俄、中朝边境线附近,较为分散的人口分布对虎豹栖息地破碎化影响巨大；其四,工矿活动作为虎豹公园内传统的生产方式,是当地居民主要的经济来源,然而工矿用地的过度开发对地表破坏严重,是虎豹栖息地景观破碎化加剧的重要影响因素。

综上所述,为保护东北虎豹种群资源及其所在的自然生态系统,亟待对虎豹公园内的人类活动进行合理管控,公园管理部门应当遵循东北虎豹种群可持续生存与发展的保护生态学规律,优先清理园区内不符合保护和规划要求的各类人工设施,对已有交通道路实行车辆通行管控,严格控制人口增量,逐渐减少农业生产、矿产开发等活动,最终实现人与虎豹的和谐共存。

此外,本文仍存在不足之处。首先,尽管本文在局部1500m格网尺度上分析了人类活动对虎豹公园内景观破碎化的影响,但仍需进一步把人类活动落到更加精细的空间位置上,因此,在未来的研究中将通过大数据等手段来细化人类活动的表达和空间定位；其次,本文仅分析了各人类足迹指标与景观破碎化之间的空间关系,并未深入探讨人类足迹对景观破碎化影响的阈值效应,即人类活动对自然生态系统产生风险的阈值。因此,这也将是未来重要的研究方向之一,旨在为国家公园的规划和发展提供更加直观的数据支持。

4.2 结论

本文通过构建适地化的人类足迹指标体系,并借助景观格局分析、回归分析、GIS空间分析等方法,从全局和局部两个尺度上探讨各人类足迹指标对虎豹公园及其管控分区景观破碎化的影响及其空间分异特征。研究结果表明：

(1)从景观破碎化空间分布来看,虎豹公园内景观破碎化程度以低、较低为主,但高度和较高破碎化区域仍占据公园全域面积的8.96%。其中在核心保护区占比为4.08%,主要分布于东南部,呈条带状分布；在虎豹潜在栖息地占比为5.09%,主要分布于该区域内人类居住点周围,呈点状分布；而在人口聚集区中的占比则高达60.54%,且以乡镇驻地为中心向外围扩散。

(2)从公园全域范围来看,对景观破碎化的影响从高到低依次为土地利用状况(0.332)、交通可达性(0.289)、人口密度(0.224)、工矿活动(0.148),而夜间灯光与景观破碎化的相关性不显著。

(3)从各管控分区来看,核心保护区与虎豹潜在栖息地内各人类足迹指标对景观破碎化的影响规律与全局一致,而在人口聚集区中,工矿活动成为该区域中第二影响因素。

(4)从各人类足迹指标主导的空间格局来看,土地利用状况对景观破碎化的影响区域最为广泛,涵盖了园区的大部分地区；交通可达性的影响区域集中分布于春化镇、大肚川镇、春阳镇周边地区及公园北部地区；人口密度的影响区域主要以罗子沟镇、复兴镇附近地区为主,此外在东部边境地区、核心保护区南部及其他小型居民点周围也有分布；工矿活动的影响区域最小,主要分布在工矿业较为集中的罗子沟镇及周边地区,且对核心保护区东北部地区也存在一定影响。