森林景观是区域整体景观的重要组成部分，在生态系统和人类游憩活动中发挥着无可替代的作用，是景观生态学研究的重点和热点。研究森林景观格局和空间分布特征对指导森林可持续经营、优化森林景观空间分布结构以及合理发挥生态经济效益具有重要的理论和实践意义^[1]。国外的森林景观格局研究起步较早，其中Forman等最早发起对森林景观组成和格局的研究，全面分析了美国新泽西州濒海平原松枥林景观格局的特点^[2]。Franklin利用棋盘格式模型、集中与分散集合模型对森林景观格局及斑块特征进行深入研究^[3]。Ripple等运用遥感技术分析Oregon州森林景观斑块结构、组成以及动态变化趋势^[4]。Boyden等人采用点格局分析方法对1974—2001年间的美国黄松的历史与当代的林分结构的异质性进行研究^[5]。国内森林景观生态研究引领景观生态研究的前沿，从种群或森林斑块组成结构、年龄结构、空间结构等方面进行了大量的研究，景观格局指数和3S技术等分析法极大地丰富了森林景观定量化的研究^{[6, 7, 8, 9]}。景观分布格局及其生态过程与空间组合规律都具有强烈的尺度依赖性，而传统方法无法全面诠释景观格局的尺度信息，从而使近年来景观格局分析研究陷入困境^［10]。起源于植物生态学的空间点格局分析法，是分析多尺度空间特征的有效手段，与传统景观格局分析法相比显示出更大的优势^{[11, 12, 13]}。

庐山是世界文化景观遗产地和世界地质公园，旅游开发历史悠久，人类活动深刻而持久，其森林植被是世界遗产核心价值的重要体现，也是开展旅游活动的基础和保障。但长期以来受到不同程度的人为破坏，并在保护中破坏、破坏中建设，森林景观空间格局复杂^[14]。本文运用点格局分析方法分别分析庐山森林景观在单一尺度和多尺度下的空间分布和集聚特征，揭示其森林景观空间结构特征，从而为森林抚育、营林造林、景观生态建设和旅游开发等提供理论依据。

庐山位于江西省九江市东南部，北临长江、南傍鄱阳湖，东经115°52′—116°8′，北纬29°26′—29°41′，属于亚热带季风气候区，气候温和湿润，全年平均温度为 11.4℃，年平均降水量1917 mm，年平均相对湿度78%，年平均雾日191 d，山地小气候特征显著。庐山山势雄伟，山体多峭壁悬崖，自东北向西南延伸约25 km，宽约15 km，山体相对高度1200—1400 m，生物资源丰富，森林覆盖率达76.6%。庐山是中国的旅游名山和文化名山，随着旅游业的蓬勃发展，庐山山顶的枯岭镇上聚集着1.5万人的常住人口，旅游活动和其他人类活动历史悠久、频繁深刻，景观格局变化剧烈^[14]。本文以江西省庐山风景名胜区管理局管辖范围为界，研究区域面积121km²。

最邻近距离分析法(Nearest Neighbor Analysis，ANN)是基于空间距离的方法，将其离得最近点的平均距离与随机分布模式下的预期最近邻距离进行比较，用最近邻指数(NNI)来判断“点”的空间聚集性^[15]。

点击浏览公式

Ripley′s K函数是点格局分析的重要方法，最初由Ripley^[16]提出，经Diggle^[17]等人发展，克服了传统景观格局分析方法只能分析单一尺度空间分布格局的缺陷，它可以描述不同空间尺度下的森林景观空间分布情况，最大限度的利用了空间点的信息^{[18, 19]}。目前被广泛应用到多个领域，主要应用于植物种群^[20]、疾病发病点^[21]、地质灾害点^[22]和旅游景点分布^[23]等空间格局分析。Ripey′s K函数公式如下^[24]：

点击浏览公式

点击浏览公式

当L(r)>0，该景观在尺度r服从聚集分布；当L(r)=0服从随机分布；当L(r)<0，服从均匀分布。用 Monte-Carlo 拟合检验计算上下包迹线，即置信区间。假设个体是随机分布的，用随机模型拟合一组点的坐标值，对每一尺度r计算L(r)，这样反复模拟直至达到事先确定的次数，把L(r)的最大值和最小值定义为上、下包迹线的坐标值。再用实际个体分布的点图计算得到真实的L(r)值，若值落在上下包迹线之间，则服从随机分布；若值小于下包迹线，则呈显著均匀分布；若值在上包迹线以上，则显著聚集分布，最大聚集强度处于偏离置信区间最远处，L(r)函数的最大值可以反映典型的特征空间聚集尺度。

庐山森林植被斑块数据来源于全国第二次林相调查资料，并参考2010年庐山土地利用现状，提取林斑中的优势种群，将其分为常绿阔叶林、落叶阔叶林、杉类、松类、竹类、建筑用地、农耕地、水域和裸地9类土地利用景观(考虑到庐山20世纪30年代以来大规模地成片引进种植日本柳杉并迅速遍布山顶，而松类则大多为分散引进，两种森林景观风貌差别较大，因此将两种针叶林景观单独划分)，一共包括1531个斑块数(图 1)。庐山森林景观占据着庐山景观的主体，其中松类、杉类、常绿阔叶林、落叶阔叶林和竹类5种森林景观面积所占百分比为88.62%，5类景观斑块数1389，占总数量的91%，而建筑、农耕地、裸地和水域等斑块镶嵌点缀其中，比例很小。本研究将从中选取常绿阔叶林、落叶阔叶林、松类、杉类和竹类5类森林景观进行深入探讨，不考虑其他非森林景观类型。

图 1 庐山景观斑块图 Fig.1 The landscape patches map of Mount Lushan

图选项

庐山森林植被呈现出古老原生林与次生林共存，天然林与人工林共存，本土物种与外来种共处的生态景象^[13]。森林的年龄结构复杂，大致包括幼龄林、中龄林、近熟林、成熟林、过熟林，本研究将近熟林、成熟林和过熟林统一归为老龄林，即按照森林发育阶段将森林植被划分为幼龄林、中龄林、老龄林3个阶段，分别对应其生命周期的幼年期、壮年期、老年期3个时期。同时又将每一类森林景观划分为天然林和人工林。

在获取庐山森林景观斑块的实地调查矢量数据后，首先运用ARCGIS 10.0软件将1389个森林植被斑块转化为1389个质心点，形成点格局数据；其次，利用平均最近的相邻要素计算最邻近指数并经正态分布检验分别得出各类景观斑块的Z值及P值；最后，利用CrimeStats3.3软件的Ripey′s K函数分析，置信度取99%进行模拟100次，所得结果经软件Grapher 7.0处理出图。

图 2 庐山地名图 Fig.2 The place names landscape of Mount Lushan

图选项

最邻近距离分析通过近邻比和Z值可探测要素的空间分布形态，从而判断区域森林景观是否服从集聚分布、随机分布还是均匀分布。利用ArcGIS10.0软件，将5类森林景观要素转点后，进行平均最近的相邻要素模式分析，得出表 1。由表 1可知，总体森林景观的近邻比接近1，Z-score值落在随机分布区间，呈现随机分布形态，主要是因为庐山属于山地，森林植被广泛分布，集聚程度不明显。而单一景观类型的近邻比均小于1，经正态分布检验，松类、杉类、常绿阔叶林和竹类的集聚特征明显，显著偏离随机分布，这主要受地形、气候、土壤、坡度、海拔等因素影响，其森林斑块在其特定的地域范围内集聚分布；落叶阔叶林分布较广泛，海拔1300m以下都有分布，集聚不显著，集聚程度弱于其它四类森林景观。

在ArcGIS10.0软件的支持下，对5类景观斑块不同发育阶段的幼龄林、中龄林和老龄林进行平均最近的相邻要素模式分析，得出表 2。从表 2中可以看出，常绿阔叶幼龄林显著集聚，中龄林近邻比等于1，服从随机分布，老龄林服从集聚分布但聚集强度弱于幼龄林，主要是因为常绿阔叶林分布在海拔较低的沟谷地带，这里村庄聚落集聚，农业活动干扰强烈，森林植被破坏严重，仅有部分风水林受到较好地保护，以幼龄林为主，残存的中龄林零星分布。落叶阔叶幼龄林和中龄林近邻比小于1，服从随机分布，而老龄林则服从集聚分布但不显著。这主要是因为庐山落叶阔叶林分布海拔偏高，也受农业活动影响，但影响程度较常绿阔叶林低，因此大都演变为次生恢复性植被，幼龄林和中龄林均呈零星分布，部分老龄林则免受破坏呈低集聚分布。松类3个发育阶段均服从聚集分布，其集聚强度以中龄林最显著，老龄林最弱，以台湾松、马尾松为优势种。杉类和松类一样，但集聚强度稍有不同，杉类的不同发育阶段的集聚强度由大到小依次为老龄林、幼龄林、中龄林。松类和杉类的空间分布格局主要受20世纪30年代以来进行大规模的人工造林(至1998年才严格实行封山育林政策)，在空间上主要集中于庐山山顶，目前大部分松类已进入中龄林阶段，杉类大部分进入老龄林阶段。竹类的3个发育阶段各有不同，其中幼龄林服从集聚分布，聚集强度显著，中龄林的近邻比大于1，Z值为正，趋于均匀分布，而成熟林的近邻比为7.35，远大于1，Z值21.05，明显服从均匀分布，竹类主要以毛竹、刚竹林、庐山玉山竹、庐山茶秆竹为主，林龄结构偏年轻化。

总体来看，庐山森林植被在幼龄时期集聚生长，随着树种年龄的增长，分布逐渐趋于随机分布甚至均匀分布。除落叶阔叶林外，其他4类森林景观的幼龄阶段均呈显著集聚分布，中龄阶段为集聚分布的只有松类和杉类，而老龄林阶段只有杉类景观分布最集聚，这与森林植被的自然演替状态基本一致。

在ArcGIS 10.0软件的支持下，运用平均最近的相邻要素分析模式分析5类景观天然林和人工林的分布形态，如表 3。常绿阔叶天然林近邻比小于1，经检验服从集聚分布，人工林则服从随机分布，偏离集聚形态。主要由于庐山常绿阔叶林属于增长型种群，以天然更新为主，呈集聚分布，而部分人工林则为封山育林后形成的，呈随机分布。落叶阔叶天然林服从集聚分布但聚集强度不显著，人工林近邻比接近1，Z值落在随机分布区域范围内，服从随机分布。落叶阔叶林受人为干扰较常绿阔叶林少，天然林和人工林随机性分布比常绿阔叶林明显。松类天然林和人工林均呈现显著集聚分布。杉类天然林服从集聚分布但不显著，其人工林服从集聚分布，集聚强度在各类不同类型景观中最大。松类和杉类受人工栽植影响，人工林集聚现象明显。竹类天然林服从集聚分布，人工林近邻比接近1，服从随机分布，因其天然林面积大且集中，而人工林分布少且分散。

可见，庐山5类森林景观天然林均呈集聚分布形态，偏离随机分布，说明庐山过去森林破坏严重，天然林分布仍然聚集于少数地区。人工林松类和杉类显著集聚分布，而阔叶林和竹类人工林均趋于随机分布，这说明过去人工营林以松类和杉类为主，阔叶林和竹类以自然入侵更新为主。

最邻近分析法只能分析单一尺度下景观空间格局特征，为进一步分析多尺度下的庐山森林景观空间分布格局，运用Ripey′s K函数分析法和Monte Carlo模拟检验，以r(0< r<5200m)为横坐标，模拟的包迹线作为纵坐标，分析结果见图 3，表 4。由此可知，5类森林景观的L(r)曲线在不同的空间尺度上比较光滑，稍有波动，均在小尺度范围内服从随机分布，中、大尺度服从集聚分布，其景观格局以集聚分布为主，随机分布为辅，最大集聚规模范围依次为：松类、竹类、落叶阔叶林、常绿阔叶林、杉类，其中常绿阔叶林和两类针叶林景观在地域分布上太过集中，不利于生态系统的生物多样性发展。

图 3 各类森林景观的Ripley L(d)指数分析 Fig.3 The Ripley L(d) index analysis of forest landscapes in different types

图选项

运用CrimeStats 3.3软件的Ripey′s K函数对森林景观各年龄阶段进行多尺度特征分析，再将分析结果以图表的形式展现，如图 4、表 5。由此可知，除竹类中龄林外，各龄级森林景观空间在小尺度范围内服从随机分布，伴随着空间尺度范围的增大开始服从集聚分布，集聚程度达到最强，随后逐渐减弱趋于随机分布，L(r)曲线随着空间距离的增大也会发生弯曲波动，甚至是一个或多个小集聚峰，多半是受到景观数量以及分布区域影响所致。值得注意的是，本研究是基于植被林斑数据进行空间多尺度特征分析，斑块被视为同一植物群落组成的均质区域，故本研究结果与植物种群不同龄级立木的点格局分析的研究结果有所差异^{[26, 27]}。同一龄级的个体在不同尺度下也有不同的分布格局，主要因为生境异质性所致^[27]。

图 4 不同龄组景观的Ripley L(d)指数分析 Fig.4 The Ripley L(d) index analysis of forest landscapes in different development stages

图选项

运用CrimeStats 3.3软件的Ripey K函数对各类天然林和人工林进行多尺度特征分析，所得结果如图 5，表 6。从图 5，表 6中可知，5种森林景观的天然林以及两种人工针叶林与5类森林景观的空间多尺度集聚特征基本一致，只是空间特征尺度和集聚规模不同，唯有阔叶林的人工林景观分布格局有异，其中常绿阔叶林人工林在460—1860m范围内形成了一个明显的聚集峰，落叶阔叶林人工林随着空间尺度的增大，随机分布与集聚分布形态随机变换，且集聚强度很弱，竹类人工林在所有尺度范围内均服从随机分布。

图 5 庐山不同林分类型景观的Ripley L(d)指数分析 Fig.5 The Ripley L(d) index analysis of forest landscapes in different forest types

图选项

总体上看，庐山森林斑块集聚分布特征明显，大部分森林景观在小尺度内服从随机分布，随着尺度的增大，会出现一个或者多个特征空间尺度，以此为分界线，景观分布格局聚集程度先加强后减弱；从景观类型上说，各类森林景观类型在多尺度下的空间格局多以集聚分布为主，随机分布为辅；而从林龄特征来看，除了竹类中龄林和成熟林两类景观外，其它龄组景观空间格局均在小尺度范围内服从随机分布，而随着尺度范围的扩大开始服从集聚分布特征；从天然林和人工林特征来看，除了落叶阔叶林和竹类人工林外，其它类型景观空间格局小尺度范围内均服从随机分布，随着尺度范围的扩大开始服从集聚分布，以其特征空间尺度为分界，聚集强度先增强后减弱。

庐山森林植被的多尺度特征主要受垂直地带性分异规律和人为活动影响。小尺度范围内庐山森林景观空间格局取决于立地条件和人为活动的影响，森林斑块在区域内普遍分布，表现为随机分布；中等尺度森林景观空间格局主要受到区域地形地貌、土壤、坡度、海拔、水文、地方小气候和人为活动等要素影响，森林植被集聚分布显著；大尺度范围内，气候与地质构造分异决定着整体森林景观空间格局的演变，由于庐山空间面积范围有限，森林植被斑块集聚程度不明显，逐渐减弱并趋向随机分布。

本文运用点格局分析方法分析庐山风景区内各类森林景观以及其不同发育阶段和不同林分结构的空间分布格局和集聚特征。

(1)从单一尺度方面看，庐山常绿阔叶林、落叶阔叶林、松类、杉类和竹类森林景观斑块均服从集聚分布；各类景观的不同发育阶段分布类型以集聚分布为主，随机分布为辅；5类景观天然林均服从集聚分布，人工林只有松类和杉类呈显著集聚分布。在特定的山地气候异质性环境和人为活动的影响下，森林种群在幼年时期集聚合作，共同对抗外界因素影响和侵入，随着树种年龄的增长，对空气、阳光等资源的争夺，促使集聚程度减弱。由于庐山过去植被破坏严重，尚处于植被演替的恢复性阶段，森林植被集聚分布特征明显，这与梁艳艳等人的研究一致^[28]。

(2)从多尺度角度来说，庐山森林景观的空间分布格局与集聚特征随着生态过程和尺度变化而变化。庐山各类森林景观在小尺度内服从随机分布，随着尺度的增大，景观分布格局聚集程度先加强后减弱。这主要受到不同尺度的垂直性地带分布规律和人为活动的影响。但这种多尺度空间特征深受景观斑块个体数量的影响。个体的数量越少，Ripley′s K函数曲线图越不光滑，很容易出现随机分布或者波动；随着个体数量的增多，Ripley′s K函数曲线图越光滑，聚集分布特征显著；当个体数量达到一定的极限，景观很可能会转变为随机分布格局。

(3)无论从单一尺度和多尺度分析都表明，庐山森林景观空间格局集聚特征明显。庐山森林植被长期受到人工栽植、农业开荒、砍伐薪碳林和旅游开发等人为活动的影响，大量原生态森林植被遭到破坏，目前正处于封山育林的恢复阶段，森林植被林龄年轻、以人工林为主，整体生态系统不稳定。纯林现象严重，树种成分单一，森林景观单调，总体质量不高，影响区域景观观赏价值。这一方面说明庐山自1998年实行封山育林政策执行初见成效，另一方面反映了林龄结构失衡，大片同龄林集聚，生态系统脆弱，抗干扰能力差。未来应该继续长久执行封山育林方针，遵循森林植被演替规律朝其地带性植被阔叶林方向演替，保护生态公益林，减少人为干扰，实现森林景观适度集聚。同时减少外来树种成分，增加乡土树种，丰富林分林相结构，增加庐山森林景观生物多样性，从而提高森林景观美景度和观赏价值。

(4)基于庐山2010年各类森林景观的林斑数据进行空间分布格局及多尺度特征研究，尚未进行时间序列的纵向比较，从而无法更加科学地揭示庐山森林景观演替的生态学机制，是本文最大的不足之处，另外将斑块质心转点没有考虑斑块面积、形状等因素，难免会造成一些误差，这些可在今后研究中进行更深入的探讨。