随着由全球生态环境剧变带给人类生存与发展的压力越来越大, 生态安全问题已成为目前人类实现可持续发展的重要障碍。森林生态系统作为陆地上最大、最复杂的生态系统, 在生态安全的保护过程中具有根本性和不可替代的作用^{[ 1]}。然而, 我国的森林资源面临诸多威胁：盗砍滥伐、毁林开荒、非法占用林地、火灾、病虫害、环境污染等, 这些有害干扰导致森林正常结构的破坏、生态平衡的失调和生态功能的退化, 严重威胁着森林生态健康和森林生态安全^{[ 2]}。因此, 通过建立评估指标体系对森林生态安全进行评价与监测具有十分重要的意义。

国外对生态安全问题研究起步较早, 1996年《地球公约》中首次提出了生态安全的定义, 随后该定义被国际应用系统分析研究所(IIASA)进行了广义和狭义上的区分。针对生态安全的研究分为宏观和微观两个领域, 宏观方面主要集中于概念以及与国家安全、军事战略、可持续发展和全球化的关系^{[ 3- 5]}, 微观方面体现在化学品的使用对农业生态系统健康及生态安全的影响^{[ 6]}。在生态安全评价方面, 研究较少, 主要集中在生态风险评价^{[ 7- 8]}和生态系统健康评价两方面^{[ 9]}, 评价模型主要有PSR、DPSEEA和ECCO等^{[ 9- 11]}。国内研究中, 米峰等在基于广义森林生态安全的基础上, 提出了森林生态安全既指森林生态系统的自身健康和安全, 又指在受到人类行为干扰和影响时保持安全的状态。强调了森林所具有的经济、社会、生态三大效益^{[ 12]}。张智光等在前人研究基础上首次提出了林业生态安全的概念, 创新性地将森林生态安全与林业产业结合论证, 并对两者之间的关系进行了探讨。同时, 他认为对森林生态安全进行监测预警非常重要^{[ 13]}。此后, 国内学者开始对生态安全进行方法上的探索, 此阶段研究也呈现出研究领域较为广泛的特征。如赵宏波采用TOPSIS模型研究了吉林省耕地生态安全时空变化^{[ 14]}。类似的还有对用地^{[ 15]}、海洋^{[ 16]}以及草原^{[ 17]}生态系统的研究。而针对森林生态系统进行评价的研究还比较少, 刘心竹等人在前人研究的基础上创新性地从森林有害干扰的视角构建评价指标体系, 并运用该体系对全国31个省进行了森林生态安全评价^{[ 18]}。张智光等采用压力-状态-影响-响应结构模型及结构方程模型对林业生态安全的共生耦合情况进行了测度^{[ 19]}。但这两位学者并未将研究区域缩小至县, 并进行时空演变分析。汤旭等虽然从空间相关性的角度考察了湖南省各县的森林生态安全状况, 得出了各地存在显著相关性的结论^{[ 20]}, 但没有考虑生态区位因素对森林生态安全造成的影响, 也没有建立模型深入阐释森林生态安全空间相关性的内在效应机制。

基于此, 本文以江苏省80个县为研究样本, 尝试引入生态区位模型结合GIS技术对2000—2015年江苏省县域森林生态安全指数(ecological security index, ESI)及其时空分布进行研究( 图 1), 并在此基础上建立空间滞后模型(SLM)与空间误差模型(SEM)探讨森林ESI的空间扩散效应与回流效应, 旨在定量化客观评价江苏省县域森林生态安全状况并深入探讨县域间森林生态发展的时空演化规律及内在效应机制, 为生态保护政策制定提供依据。

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

江苏省地处中国大陆东部沿海地区, 是长江三角洲重要组成部分, 地跨116°18′—121°57′E, 30°45′—35°20′N, 内陆面积10.72万km², 占中国土地总面积的1.12%。江苏省生态区位条件较好, 属于温带向亚热带的过度性气候, 气候温和, 雨量适中, 四季分明。但江苏省森林资源较为贫乏, 截至2016年, 江苏省森林覆盖率仅为15.8%, 居全国28位。江苏省人口非常稠密, 2016年人口密度为每平方千米744人, 居全国之首；实现地区生产总值76086.2亿元, 居全国第二位。由此可见, 江苏是经济大省, 但同时又是森林资源小省。人口密度大, 经济发展水平高, 但森林资源短缺, 这就给经济、社会、资源、环境可持续发展带来了巨大压力。因此研究该省森林生态安全状况具有十分重要的现实意义。

1.2 数据处理

本研究的样本为江苏省80个区县, 涉及2000年、2005年、2010年、2015年4个年份的数据, 所有数据来自生态安全(森林)指数研究项目。为方便数据的收集与统计, 课题组设计了数据投放器, 以方便各区县在线输入数据, 其设计思路如 图 2所示。

数据收集工作分几个步骤完成：首先组织江苏省林业厅、规划院进行了数据投放器的使用培训会议；江苏省林业厅将数据收集工作细化到各县林业局、规划院, 各单位完成数据填报, 江苏省林业厅对数据进行审核；课题组专家对数据进行修正。所有数据处理工作于2017年4月完成。

2 研究方法

文章采用PSR模型建立森林生态安全指标体系, 熵权法与专家法相结合的方法来确定权重, 选用模糊物元法与欧式贴近度来评价森林ESI。然后借助生态区位模型对森林ESI进行修正。最后运用GIS技术和空间计量模型对森林ESI的空间分布进行深入分析。

2.1 指标体系的建立 2.1.1 森林生态安全指标体系

文章以PSR模型为基础对江苏省县域森林ESI评价指标体系进行研究。同时对PSR模型进行改进, 将模型中三类指标合并为压力指标和状态指标, 其中压力指标是把原模型中的压力指标与响应指标结合起来, 因为这两类指标反映了人类对森林施加的压力以及人类保护森林的努力, 都属于人类行为指标。课题组经过文献综述及专家咨询首先选定了对森林生态安全影响较大的多个指标。然后通过STATA软件的主成分分析功能, 经方差最大化旋转后生成10个主因子, 累计方差贡献率达到90%, 在旋转后的主成分因子负荷矩阵中初步选取负荷系数和≥1.00的30余指标。2016年4月课题组又进一步对江苏省盱眙县、宜兴市进行了实地调研, 发现天然林蓄积量、森林生物量等指标存在年份连续性差、县域之间数据差别较大或统计口径不一致等问题, 因此剔除了数据之间相关性强以及数据可获得性差的指标, 随后通过电话咨询以及查询统计年鉴进行数据填补, 最终形成了包括森林覆盖率、森林单位面积蓄积量在内的16个具体指标, 如 表 1所示。

2.1.2 森林生态安全指标权重

文章采用熵权法与专家打分法相结合的方法确立综合权重, 计算结果如 表 1所示。这样既克服了熵权法极值点对计算结果影响过大的缺陷, 也克服了专家法主观性过强的缺陷, 实现了良好的评估效果。综合权重中, 人口密度、单位面积能源消耗量、退耕还林面积占比、森林自然保护区面积占比、森林生态建设与保护投资强度等指标权重较大(>0.08), 因此提高江苏省森林生态安全状况要降低人类活动对森林生态系统的负面影响, 增加县域森林覆盖面积, 加大投资力度与保护力度。

2.2 模糊物元法

为了解决各类指标的不相容的问题, 本文决定采用模糊物元法。模糊物元法用M、C、x分别表示模糊物元、特征、模糊量值。若有m个样本, n项指标(本文有80个样本, 16项指标), 则可用M_i表示第i个样本(i=1, 2…m), C_j表示第i个事物的第j项指标(j=1, 2…n), x_ij表示第i个样本第j个指标的模糊量值(i=1, 2…m；j=1, 2…n)^{[ 21- 22]}。则江苏省县域森林生态安全的复合模糊物元R_mm如下所示：

(1)

2.2.1 计算从优隶属度

对于森林生态安全的每项指标C_j, 可以计算其相对于标准事物的从优隶属度。采用以下两个公式：越大越优型指标：

(2)

越小越优型指标：

(3)

式中, u_xij表示森林生态安全各指标的从优隶属度(i=1, 2…m；j=1, 2…n), 各指标中的最大值和最小值分别用maxx_ij和minx_ij表示。基于以上计算结果得出新的模糊物元R_mm：

(4)

2.2.2 标准模糊物元与差平方模糊物元

因为森林生态安全指标体系各特征值中u_xij的最大值均为1, 所以得出标准模糊物元R_0n：

(5)

依据Δ_ij=(u_x0j-u_xij)²和R_0n来计算差平方, 得出差平方模糊物元R_Δ：

(6)

2.3 欧式贴近度与ESI综合评价

经过标准化之后, 文章采用欧式贴近度公式来计算森林ESI值, 该方法计算时首先确定最优方案, 然后计算其他方案与该方案的贴近程度, 越接近则该方案越优^{[ 23]}。本文中复合模糊物元R_PH的计算采用M(*, +)方法：

(7)

式中, 是第j个指标的熵权权重, Δ_ij为第i个样本县第j个指标的平方差。R_PH为从第一个到第m个样本县森林ESI的集合, PH_i即为第i个样本县的森林ESI值。

2.4 生态区位模型

生态区位模型是一个综合反映气象、地形等区域自然条件因素的系统, 它能够通过影响物种分布和植被生产力来对森林生态安全发生作用, 因此可以用它来对森林生态安全指标体系进行修正和补充。本文运用ArcGIS软件对8个生态区位指标做栅格化处理, 再对各评价指标进行分级赋值, 最后运用专家打分法确定了气象类系数与地形类系数的权重分别为0.625和0.375。在确定指标的权重后, 再对江苏省各区县的生态区位系数进行评估。根据指标的具体特征, 将每个指标分为5个级别, 等级越高对于林业的发展越有利。目前尚缺少对每一指标分级区间划定的研究, 因此对于每个指标分属不同级别的区间划定, 主要通过查阅参考文献^{[ 24]}, 再结合区位系数的定义和应用范围最终确定( 表 2)。在计算生态区位系数时, 首先应依据 表 2的评级标准和参评单元评价因子的现状值, 确定现状区位系数评分空间分布图, 然后计算指标层区位系数分值:

(8)

式中, F_ij、S_ijk、W_ijk分别为j指标的区位系数分值、评价指标现状值、权重值, n为指标个数^{[ 25]}。同理可以计算得到准则层与目标层的区位系数分值。结合前面已经计算出的森林生态安全指数, PH_i可以通过公式：计算修正以后的森林生态安全指数。ESI为经生态区位系数修正后的森林生态安全值, PH_i为最初森林生态安全值, F为区位系数综合分值, F为全省区位系数平均值。

2.5 探索性空间数据分析 2.5.1 全局空间自相关

全局空间相关性考察属性值在研究区域范围内空间相关性的整体趋势。常用统计量Moran′s I。当Moran′s I为正, 表明森林ESI在空间上正相关, 在空间上趋于聚集；当Moran′s I为负, 表明森林ESI在空间上负相关, 在空间上趋于分散。当值等于0, 表示森林ESI在空间单元之间空间不相关^{[ 26]}。莫兰指数的计算公式为：

(9)

式中, n是研究区域内地域单元总数；w_ij是空间权重矩阵的元素值；x_i是地域单元i的x变量值；y_i是地域单元i的y变量值。

2.5.2 局部空间自相关

局部空间自相关统计量Local Moran′s I是Moran′s I的分解。局部莫兰指数高值表明有相似变量值的面积单元在空间集聚(高值或低值), 低值表明不相似变量值的面积单元在空间集聚。地域单元i的局部莫兰指数用来度量它和周围地域单元之间的相关性, 定义为：

(10)

式中, 正的I_i表示一个高值被高值所包围(高-高), 或者是一个低值被低值所围(低-低)。负的I_i表示一个低值被高值所包围(低-高), 或者一个高值被低值所包围(高-低)。

2.5.3 空间滞后模型与空间误差模型

传统计量分析在研究过程中没有将空间因素纳入到分析模型, 从而忽视了空间效应, 使模型设定与现实情况存在一定偏差, 进而可能导致研究得出的各种结果和推论不够完整、科学, 缺乏应有的解释力^{[ 27]}。因此本文引入空间滞后模型与空间误差模型对森林ESI的扩散效应与回流效应进行探讨与分析。

空间滞后模型(SLM)测度局部地区之间是否存在扩散效应, 其模型表达式为：

(11)

式中, ρ是向量空间自回归参数；ε_{i, t}是随机误差项。

空间误差模型(SEM)测度空间依赖作用是否由于误差扰动项引起, 其模型表达式为：

(12)

式中, X_{i, t}为第i个地区在t年的森林ESI, X_{i, t+k}为第i个地区在t+k年的森林ESI；α、β为待估的参数；I是单位矩阵；ξ是误差项空间自回归参数；W为空间权重矩阵；μ(i, t)为随机误差项。

3 结果与分析 3.1 生态区位系数计算结果分析

文章首先基于公式(8)及 表 2的8个指标和江苏省80个区县的数据, 计算出江苏省各县域的区位系数。本文运用ArcGIS软件将江苏省县域生态区位系数采用五分位图展示出来, 如 图 3所示。

本文基于自然断裂法将区位系数分为五个等级。从图中可以直观地看出, 区位系数最高(30.5310—32.8123)区域主要分布在江苏省南部少数地区, 如南京市、宜兴市、无锡市、东台市、昆山市、太仓市等。区位系数次高(29.3620—30.5309)区域主要分布在江苏省中南部, 如苏州市、丹江市、盱眙县、宝应县等。区位系数最低(23.0436—28.2045)的区域主要分布在江苏省北部, 如响水县、涟水县、丰县、沛县等地。

从总体来看, 江苏省的生态区位系数平均值为29.27, 生态区位较好地区主要分布在江苏省南部, 而最高值也出现在苏南地区的高淳县(32.8123)。这些地区自然基础条件较好, 适合森林树木生长。而生态区位较差区域则集中在苏北, 尤其是江苏省东北部地区, 如响水县、滨海县、涟水县、射阳县等, 说明这些地区自然条件较不适合森林植被生长。与此同时, 江苏省内生态区位较好地区大多集中在苏南经济发达地区, 而生态区位较差的地区则一般聚集在苏北地区, 这些地方经济相对落后, 可见生态区位系数与当地经济发展水平具有一定的耦合性。

3.2 森林ESI修正指数计算与ArcGIS分析

计算出生态区位系数之后, 结合模糊物元法与欧式贴近度方法, 并利用公式(1)到(7)可计算出2000—2015年江苏省森林ESI指数的修正值, 同时用ArcGIS绘制五分位分布图进行展示, 如 图 4所示。

从2000年来看, 江苏省森林ESI指数最高区域集中在苏南地区以及苏北少部分地区, 如东海县、沭阳县等。ESI指数次高区域则分布相对均匀。而ESI指数最低区域较为集中在江苏省中东部及北部地区个别县市, 如滨海县、射阳县、邳州市、新沂市等。

从2005年来看, 江苏省森林ESI指数最高区域较多分布在苏南地区, 江苏中部及北部地区也有少部分分布, 如沛县、睢宁县、宝应县、江阴市、无锡市等。相较2000年来看, ESI指数最高区域分布更为均衡。ESI指数次高区域分布较上年更为广泛。而ESI指数最低区域则仍集中在苏北及中东部地区。

从2010年来看, 江苏省南部地区森林ESI指数变化不大, 只有江阴和无锡退出最高层级, 同时新增溧水县和高淳县, 苏北及中部地区退步较为明显, 沛县、沭阳县、宝应县和兴化市退出该档。森林ESI指数最低区域新增泗洪县、宝应县、兴化县、海安县、启东县、如皋市等地, 而苏北的宿迁、新沂、邳州则有所改善。

从2015年来看, 江苏省森林ESI指数最高区域仍集中在苏南地区, 全省新增金坛市、常州市、如东县、启东县, 而宜兴市、常熟市、太仓市、昆山市、张家港市则退出该档, 苏北地区则仅剩灌云县；森林ESI指数最低区域则集中在中东部地区, 苏南也有部分地区出现ESI指数过低的情况, 如句容市、宜兴市、常熟市。

从2000—2015年江苏省县域森林ESI等级变化来看, 67.5%的区县都呈现出较明显的下降趋势, 仅徐州市泉山区、铜山区、新沂市一直保持上升势头。而苏北很多处于ESI最低档县市的状况却得到了改善。究其原因, 这与2003年“绿色江苏”的建设实施关系密切。自2003年以来, 全省林地面积、森林面积、森林覆盖率和林地利用率明显增加, 活立木总蓄积和森林蓄积大幅度增长, 江苏省森林绿化建设在“量”上有了突破性进展。但2010年以后, 苏北地区大部分县市森林ESI变化不大, 同时苏南地区部分县市森林ESI出现急剧下降, 如宜兴市、苏州市、常熟市、句容市等。原因是2010年以后, 虽新增造林不断增加, 但受诸多因素影响, 全省森林面积净增量有所下降, 且由于树种单一造成的森林生态功能降低问题也越来越突出。与此同时, 苏南与苏中等经济发达地区不断加大开发建设力度, 进一步加大了对当地森林生态安全的压力, 导致苏南地区许多县市森林ESI有所下降。可见, 江苏省在森林生态安全建设过程中, 不仅要注重大规模植树造林, 更要注重“质”的提升, 优化树种结构, 协调区域生态发展。

3.3 空间相关分析 3.3.1 全局自相关分析

在得出2000—2015年江苏省县域森林ESI值后, 可根据全局自相关公式(9), 并运用GeoDA软件计算得出森林ESI的莫兰指数及检验结果, 如 表 3所示。

从检验结果来看, 都通过了P≤0.01的检验, 表明全局自相关实验在99.9%置信度下显著, 说明江苏省县域森林ESI存在较为明显的空间相关性。从时间变化看, Moran′s I指数在不同阶段呈现出不同的特征。2000—2005年, 空间相关性有较小幅度上升；2005—2010年, 空间相关性则出现大幅上升；2010—2015年, 空间相关性则出现急剧下降。但从总体来看, Moran′s I指数还是有一定程度的降低, 说明江苏省县域间的森林ESI空间聚集状况有所下降。

如 图 5所示, 2000—2015年江苏省全局Moran′s I散点大部分位于一、三象限内, 即高-高和低-低类型聚集区占主导地位, 表现出较明显的空间依赖性, 也就是森林ESI较高的区县相互聚集, 而较低值的区县也相互聚集。高值区县互相聚集表明各区县能够形成相互促进作用, 而低值区县聚集说明各区县会形成一个相互的负面影响, 进而陷入一个恒差的恶性循环。2000、2005年, 江苏省位于第一象限与第三象限的县域个数基本相当, 而2010、2015年, 第三象限县域个数明显多于第一象限, 说明江苏省县域存在明显的森林生态安全较好和较差的县域, 而森林生态状况较差区县的聚集数量明显多于较好区县的数量, 并且有逐渐加重的趋势。究其原因, 在于2010年之后, 苏北平原早期发展的杨树林大面积进入采伐期, 且杨树木材价格相比农作物明显下降, 许多林农对杨树进行采伐之后就不再栽植。由此可见, 苏北杨树林面积大幅度减少是造成低-低聚集区增多的最主要原因。

3.3.2 局部自相关分析

为了更加深入地研究森林ESI指数在江苏省各县域之间的空间分布情况, 利用公式(10)及GeoDA软件可绘制出2000、2005、2010、2015年江苏省县域森林生态安全指数LISA聚类图, 且通过P≤0.05的显著性检验, 如 图 6所示。

图 6反映了各区县之间的局部自相关关系, High-High表示森林ESI高值区县被周围高值区县包围, Low-Low表示低值区县被低值区县包围, Low-High表示低值区县被高值区县包围, High-Low表示高值区县被低值区县包围。通过分析LISA聚集图可以发现：(1)4个年份中, 森林生态安全指数High-High集聚区县一般分布在苏南地区, 以苏州、无锡一带为中心。(2)Low-Low集聚区县则集中在江苏东北部一带, 以阜宁县、建湖县为中心；(3)从各年的显著聚类结果来看, 2000、2005、2010和2015年High-High集聚区县分别有9个、8个、8个、3个, Low-Low集聚区县则变化不大, 分别有8个、5个、6个、5个, 2005年之后基本稳定在5个左右, 但分布呈现分散化的趋势。为了更加直观地表示2015年江苏省县域森林ESI各类聚集情况, 本文将各类聚集区域制成 表 4。

3.3.3 江苏省县域森林ESI空间计量分析

上述探索性空间数据分析(ESDA)仅证明了江苏省各县市之间存在空间相关性, 且空间分异特征明显。然而ESDA在本质上是“数据驱动”的, 单纯依靠ESDA在对森林ESI进行时空分析时是缺乏理论模型解释力的。因此, 文章在ESDA的基础上, 运用CSDA(证实性空间数据分析)对江苏省各县市森林ESI指数的时空演进模式进行深入研究, 探讨县域之间森林ESI指数的空间结构模式。文章将2005年、2010年设为趋势分界点, 运用公式(11)、(12)及Geoda软件, 得出空间滞后模型(SLM)与空间误差模型(SEM)估计结果, 如 表 5所示。

在 表 5中, 反映模型拟合效果的指标有LIK、AIC及SC指标, 其中LIK越大, AIC和SC指标越小, 则说明模型的拟合效果越好。通过对比2000—2005、2005—2010和2010—2015 3个研究时间段LIK、AIC及SC指标的大小, 可以发现空间滞后模型的解释力比空间误差模型更强, 即森林ESI在江苏省县域之间存在扩散效应。同时, 空间滞后模型结果显示, 估计值β在3个研究时期内均为正, 表明森林ESI在江苏省县域之间也存在回流效应。因此, 我们可以得到, 森林ESI在江苏省县域之间并不只存在单一的某种效应, 而是扩散效应与回流效应并存。这也证明了江苏省内各个县市之间的森林生态安全状况并不只与自身的森林生态状况有关, 也与临近地区的森林生态安全状况存在密切的关系。

通过ESDA和CSDA分析发现, 江苏省县域森林ESI确实具有空间依赖性, 因为县域森林生态安全往往受周边地区森林生态状况的影响, 但不同区域之间表现差异较大, 如 图 7所示。江苏省县域森林ESI扩散效应主要存在于苏南与苏北之间, 苏南以苏州、无锡一带为中心点对临近地区产生正向扩散效应, 但随着时间推移, 影响范围逐渐缩小, 到2015年甚至出现了中心点向西南转移的现象；而苏北则以阜宁县、建湖县为中心点对临近地区动态地产生负向扩散效应, 如临近的响水县从2005年的高低孤立区变为了2010年的低值萧索区。而江苏中南部个别地区受到苏南的正向扩散效应并不明显, 反而由于受到回流效应的影响与苏南的差距有所扩大。

4 结论与讨论 4.1 结论

本文首先运用熵权法、专家法及模糊物元法计算出江苏省县域森林ESI指数, 然后计算出江苏各区县的生态区位系数, 并用来修正森林ESI指数。本文还独创性地运用了ArcGIS技术、ESDA和CSDA相结合的方法来分析江苏省各区县森林ESI指数的空间分布规律, 结论如下：

(1) 从16个评价指标的权重来看, 人口密度的权重最高, 其次是单位面积能源消耗量与退耕还林面积占比, 而土地沙化强度所占的权重最小。

(2) 从生态区位模型来看, 江苏省生态区位系数最高的地区主要分布在江苏省南部少数地区, 如南京市、宜兴市等；区位系数最低的区域主要分布在江苏省东北部, 如响水县、涟水县、丰县等地。

(3) 从ESI指数的空间分布来看, 2000—2015年森林ESI等级最高的区域集中在江苏省南部, 而ESI等级最差的区域则较为集中在江苏省中部及东北部地区。这些地区一方面区位系数较低, 反映了苏北地区的光、水、热等条件较差, 不利于森林植被生长, 另一方面这些地区大面积种植杨树林, 树种单一, 森林的生态功能也比较脆弱, 故ESI指数较低。

(4) 从ESI指数的时间变化来看, 与2000年相比, 2015年江苏省67.5%的区县呈现出比较明显的下降趋势, 反映出江苏省较为严峻的森林生态安全发展状况。其中下降幅度最大的是宜兴市, 其次是江阴市和苏州市吴中区, 这些地区虽然生态区位优势明显, 但近年来由于经济发展, 人口增长较快, 对森林生态系统造成了较大压力。而增长幅度最为迅猛的是新沂市, 其次是邳州市, 这些地区虽然生态区位不佳, 但由于2003年绿色江苏造林工程实施以后, 其森林生态安全状况得到了明显改善。

(5) 从空间相关性来看, 江苏省各区县森林ESI呈现出比较明显的聚集性, 但2000—2015年空间聚集程度有所下降, 且Low-Low集聚区多于High-High集聚区。森林ESI的High-High集聚区县一般分布在江苏省南部地区, 但2015年区域范围大幅度缩小并发生了向西南方向转移的现象；Low-Low集聚区县则大多集中在江苏东北部一带。

(6) 从SLM与SEM模型来看, 森林ESI在江苏省县域之间是扩散效应与回流效应并存的。扩散效应主要存在于苏南与苏北之间, 苏南以苏州、无锡一带为中心点对临近地区产生正向扩散效应, 但随着时间推移, 影响范围逐渐缩小；苏北以阜宁县、建湖县为中心点对临近地区产生负向扩散效应。江苏中南部个别地区回流效应发挥主导作用, 且不具有收敛性, 与苏南的差距有所扩大。

4.2 讨论

从森林ESI计算结果来看, 苏南地区森林生态安全状况整体上好于苏北地区。但处于生态区位系数最高的苏南地区在2015年却出现了ESI值过低的区域, 说明经济发展、人口压力、能源消耗、废水废气排放等压力因素制约了该区域森林生态系统的发展。而苏北地区虽然生态区位系数不高, 但随着造林工程的实施, 森林生态安全状况也得到了一定程度的改善, 但也有诸如树种单一, 森林生态稳定性差等问题亟待解决。因此本文建议该地区相关政府、企业在控制人口、能耗及排污强度的同时, 不仅要提高森林植被的“量”, 更要加大造林树种的多样化选择, 彻底扭转森林结构不合理, 珍贵化和彩叶化水平低的不利局面。从空间分析结果来看, 江苏省县域森林ESI呈现出较明显的南北分异特征, 且扩散效应与回流效应并存。因此该地区政府在进行森林生态治理时应注意省内地域平衡, 防止南北差异进一步扩大, 同时加大对阜宁县、建湖县的治理力度, 防止其对周边地区森林ESI的负面扩散效应。而对苏州、无锡一带应加大其对周边地区的正向扩散效应, 且防止其产生回流效应从而对周边地区产生不利影响。

从森林生态安全的研究方法来看, 本文评价指标忽略了一些管理因素的影响, 比如自然保护区级别和各地林场等, 保护区的级别不同保护效果也会有所差别^{[ 28]}, 此外各类林场对森林的保护在客观上比个体农户的效果要好^{[ 29]}, 因此后续研究将逐渐加入这些指标。其次本文生态区位系数模型仅考虑了气象类指标和地形类指标, 但在影响森林植被的影响因素中, 还有水文、土壤、环境敏感性^{[ 30]}等因素没有考虑进来。水文因素中, 汛期、径流量等对流域内的地表水和地下水有较大影响, 而这些水源是植物生长的重要因素。土壤能够给森林植物提供大量的养分, 其中的有机质、微生物和水分等是植物生长的必要条件。环境敏感性因素中, 土壤盐碱化敏感性、水污染敏感性等对森林生态环境都有较大影响。因此, 本文今后将在水文指标中引入汛期、结冰期、径流量等指标, 在土壤类指标中引入土壤有机质、土壤类型等指标, 在环境敏感性指标中引入土壤盐碱化敏感性、水污染敏感性等指标。后续研究应引入更多理论方法^{[ 31- 33]}, 对不同模型^{[ 34- 37]}与技术^{[ 38- 40]}之间进行横向比较，更大程度地丰富森林生态安全评价理论。