天目山常绿阔叶林优势种群胸径的空间连续性分析

方国景，汤孟平

（1.浙江省森林资源监测中心，浙江杭州 310020；2.浙江农林大学环境与资源学院，浙江临安 311300）

天目山常绿阔叶林优势种群胸径的空间连续性分析

方国景1，汤孟平2

（1.浙江省森林资源监测中心，浙江杭州 310020；2.浙江农林大学环境与资源学院，浙江临安 311300）

研究常绿阔叶林优势种群胸径的空间连续性以揭示常绿阔叶林的空间结构特征。在天目山国家级自然保护区内选择典型常绿阔叶林设置100 m×100 m样地，用全站仪测定每株树木坐标，选择胸径在5 cm以上的乔木，用优势度分析法确定群落优势种群，并用GS+软件进行空间统计分析。结果显示：样地整体结构比为0.7，平均空间连续性范围为5.7m，表明天目山国家级自然保护区常绿阔叶林的优势树种林分胸径具有较高的空间相关性，但林分整体的空间连续性范围较小，这也体现了常绿阔叶林的结构复杂性。图2表3参15

森林生态学；常绿阔叶林；优势种群；地统计学；空间连续性；天目山

常绿阔叶林是中国亚热带地区最复杂、生产力最高、生物多样性最丰富的地带性植被类型之一，对保护环境、维持全球性碳循环的平衡和人类的持续发展等都具有极重要的作用［1］。对常绿阔叶林的研究是维护亚热带地区生物多样性和构建稳定森林结构的基础工作，也是近年来生态学研究的热点之一［2-4］。地统计学（GS，geostatistics）于20世纪50年代初开始形成，60年代在法国著名统计学家Matheron［5］的大量理论研究工作基础上形成一门新的统计学分支。地统计学如今已经被广泛用于地理学、生态学、环境科学、土壤学等诸多领域的研究中［6］。在林业上，通常把林分中的树木视为离散点，分析林木空间分布格局［7］。实际上，树木的属性特征（胸径、冠幅、树高等）的空间分布被认为是具有空间连续性的［8］。在树木的各属性特征中，胸径是最常用的指标，也是最容易准确测定的指标。本研究通过对浙江省天目山国家级自然保护区内的常绿阔叶林优势种群的胸径空间分布进行地统计分析，以便揭示常绿阔叶林的空间结构特征，为常绿阔叶林的保护和恢复重建提供理论依据。

1 研究区概况

浙江省天目山国家级自然保护区位于浙江省西北部临安市境内的西天目山，距省会杭州94 km，30°18′30″～30°24′55″N，119°23′47″～119°28′27″E。区内年平均气温为8.8～14.8℃［9］；≥10℃积温2 500～5 100℃；年降水量1 390～1 870 mm，相对湿度76%～81%。自然保护区受海洋暖湿气候影响较深，具有中亚热带向北亚热带过渡特征，森林植被十分茂盛。由于区内地势较为陡峭，海拔上升快，气候差异大，植被的分布有着明显的垂直界限，自山麓到山顶垂直带谱为：海拔870m以下为常绿阔叶林；870～1 100 m为常绿、落叶阔叶混交林；1 100～1 380 m为落叶阔叶林；1 380～1 506 m为落叶矮林。其中，常绿阔叶林是自然保护区内重点保护植被类型。

2 研究方法

2.1 调查方法

选择保存较完好的常绿阔叶林设置样地，样地大小为100 m×100 m（图1）。用相邻格子调查方法，把样地划分为100个10 m×10 m的调查单元。在每个调查单元内，对胸径大于或等于5 cm的木本植物进行每木调查，记录树木种类，测定每株树木的胸径、树高、活枝下高、冠幅等因子，采用激光对中全站仪（徕卡TCR702Xrange）测定每株树木坐标（x，y，z）。

2.2 数据分析方法

2.2.1 优势种确定优势树种按优势度分析法确定［10］。其方法是首先计算群落乔木层每个树种的相对胸高断面积，并作为优势度，按优势度从大到小排序。然后，通过下式确定优势树种数：

式（1）中：xi为排序在前的上位种（T）的相对胸高断面积，x为优势树种所占的理想百分比，xj为上位种以外的剩余种（U）的相对胸高断面积，N为总种数。如果群落只有1个优势树种，则优势树种的理想百分比为100.0%。如果有2个优势树种，则它们的理想百分比为50.0%，如果有3个优势树种，则理想百分比为33.3%，依次类推，分别计算d值。当d为最小值时的上位种数为群落乔木层优势树种数。

2.2.2 空间统计分析地统计学的核心是半变异函数。半变异函数分析是通过测定区域化变量分隔等距离的样点间的差异来研究变量的空间相关性和空间结构。分隔距离h的2点x0和x的区域化变量Z（xi）和Z（xi+h）之间的变异，可用增量值Z（xi）-Z（x+h）平方的数学期望（即区域化变量增量的方差）表示：

图1 样地三维地形图Figure 1 3D map of the plot

式（2）中：γ（h）称为半变异函数。半变异函数既是距离h函数，又是方向a的函数。半变异函数曲线图（semivariogram）是半变异函数值γ（h）对距离h的函数的图形，它有3个特征参数：基台值（sill），变程（range）和块金值（nugget），可作为各个方向的平均值，也可作为某一特定方向的值。基台值表示变量在研究系统中最大的变异程度，包括空间结构方差（c）和块金方差（c0）。空间结构方差表示非随机的结构原因形成的变异，块金方差则反映的是由实验误差和小于最小取样尺度所引起的随机变异。变程表示研究变量在空间上自相关的范围。对于观测的数据系列Z（xi）（其中i=1，2，3，…，n），样本半变异函数值可由下式计算：

式（3）中：γ（h）是相隔距离为h的半变异函数曲线图的估计值，N（h）为被h分割的数据对（xi，xi+h）的总数即是相隔距离为h的所有点的配对数，Z（xi）和Z（xi+h）分别为在点x和（xi+h）处样本的测量值，h为样点间距离（图2）。当用变异函数定量描述整个区域时，需给变异曲线配制相应的数学模型。本研究拟合了有基台值的3个模型：球状模型、指数模型和高斯模型（表1）。Journel等［11］推荐遵循间距和最远配对样品之间的距离。间距距离4m保证有1个足够的间距配对数。最大搜索距离使用相应样地宽的一半，即50 m［12］。地统计学计算分析使用的是GS+软件。

图2 1个典型的半变异函数Figure 2 A typical Simivariogram

表1 地统计模型Table 1 Models of geostatistics

3 结果与分析

3.1 优势种群

根据调查结果进行统计，样地内胸径5 cm以上的乔木共计1 603株，72个种。计算群落乔木层各树种的相对胸高断面积，并按降序排列，根据式（1）确定优势树种。结果显示：群落有13个优势树种（表2）。13个优势树种的株数和胸高断面积分别占样地的78.66%和81.13%，表明这些树种在森林群落乔木层中占明显的优势。其中细叶青冈Cyclobalanopsismyrsinaefolia（A），杉木Cunninghamia lanceolata（B），短尾柯Lithocarpus brevicaudatus（C），青冈Cyclobalanopsis glauca（E）和豹皮樟Litsea coreana（M）是常绿树种，占样地总株树的61.01%。细叶青冈（A），短尾柯（C）和青冈（E）三者的株树占了总株树的52.71%，胸高断面积占了总样地胸高断面积的34.04%。

3.2 优势种群胸径空间连续性

不同方向上的半变异函数往往会表现出不同的性质，这种不同方向上的差异称为各向异性［13］，反之称为各向同性。各向同性是各向异性的特例。对本次调查的数据进行方向差异性检验时发现树木胸径的空间相关性仅仅与2点间的距离有关，即具有各向同性，因此不考虑方向效应，用各向同性变异函数进行分析。在各优势树种中，蓝果树（H）和金钱松（L）的株数太少，达不到最低配对数，不符合地统计分析要求，不予计算。黄檀（K）函数变程大于50m，超出了最大搜索距离，在空间分布的函数假设中，认为超出最大距离的树是没有意义的，也不予计算。地统计模型的拟合，首要考虑决定系数r2的大小，然后考虑变程和块金值的大小，从而确定最佳模型。结果如表3。

从表3可见：豹皮樟（M）块金值为0，半变异函数值具有纯金块效应，即在该样地空间尺度上不存在空间相关性。杉木（B），枫香（D），白栎（F），黄连木（G），天目木姜子（I）和榉树（J）的结构比大于0.75，空间相关性较强；细叶青冈（A），短尾柯（C），青冈（E）结构比都在0.55以上，具有中度的空间相关性。短尾柯（C），青冈（E）和细叶青冈（A）的变程较大，空间连续性范围分别为35.70 m，45.55 m和11.10 m，其余优势树种的变程都在7.50m以内。样地整体结构比为0.70，具有中度的空间相关性，平均空间连续性范围为5.70m。

4 结论与展望

表2 优势种群及其比例Table 2 Dominant tree species and mingling

天目山国家级自然保护区常绿阔叶林的优势树种中，以细叶青冈、杉木、短尾柯和枫香为主，形成多优势种结构特征，林分胸径特征具有较高的空间相关性，但林分整体的空间连续性范围较小，仅在6 m以内，这也体现了常绿阔叶林的结构复杂性。汤孟平等［14］采用Ripley's K（d）函数分析优势种群空间分布格局和种间关联关系，认为天目山常绿阔叶林群落的优势种均呈显著聚集分布，多数优势种间有较强的种间关联性。这与本研究的结果是一致的。本研究仅对树木胸径这一属性特征，分析其空间连续性，林分其他因子，包括土壤性质、竞争等的空间连续性有待进一步研究。地统计分析是研究林分特征的有力工具，有必要在林业领域开发出具有林业特色的实用空间统计软件，使空间统计学理论、方法、程序系统一体化［15］，以促进空间统计分析在林业领域的广泛应用。

表3 优势种群胸径变异函数计算结果Table 3 Dominant tree speciesbasal area and height variograms

参考文献：

［1］包维楷，刘照光，刘朝禄，等.中亚热带湿性常绿阔叶次生林自然恢复15年来群落乔木层的动态变化［J］.植物生态学报，2000，24（6）：702-709.

BAOWeikai，LIU Zhaoguang，LIU Chaolu，etal.Fifteen-year changes of tree layer in secondaryCastanopsis-Schimahumid evergreen broad-leaved forest in central subtropics of western China［J］.ActaPhytoecolSin，2000，24（6）：702-709.

［2］丁圣彦，宋永昌.常绿阔叶林植被动态研究进展［J］.生态学报，2004，24（8）：1769-1774.

DING Shengyan，SONG Yongchang.Research advances in vegetation dynamic of evergreen broad-leaved forest［J］.ActaEcolSin，2004，24（8）：1769-1774.

［3］蔡锡安，彭少麟，曹洪麟，等.广州罗岗村边次生常绿阔叶林群落分析［J］.应用与环境生物学报，1998，4（2）：107-114.

CAIXi'an，PENG Shaolin，CAO Honglin，etal.Community analysis of a secondary evergreen broadleaved forest at Luogang，Guangzhou［J］.ChinJApplEnvironBiol，1998，4（2）：107-114.

［4］石胜友，尚进，田海燕，等.缙云山风灾迹地常绿阔叶林生态恢复过程中优势种群分布格局和动态［J］.武汉植物学研究，2003，21（4）：321-326.

SHIShengyou，SHANG Jin，TIAN Haiyan，etal.Distribution pattern and dynamics of dominant population in the progression of ecological restoration of evergreen broadleaved forest after wind-damage in Jinyun Mountain［J］.JWuhanBotRes，2003，21（4）：321-326.

［5］MATHERONG.Principles of geostatistics［J］.EconGeol，1963，58：1246-1266.

［6］肖斌，赵鹏大，侯景儒.地统计学新进展［J］.地球科学进展，2000，15（3）：293-296.

XIAO Bin，ZHAO Pengda，HOU Jingru.New development of geostatistics［J］.AdvEarthSci，2000，15（3）：293 -296.

［7］MARK R TDALE.SpatialPatternAnalysisinPlantEcology［M］.Cambridge：Cambridge University Press，1999.

［8］CHEN Jiquan，SONG Bo，RUDNICKIM，etal.Spatial relationship of biomass and species distribution in an oldgrowthPseudotsuga-Tsugaforest［J］.ForSci，2004，50（3）：364-375.

［9］张志华，韦新良，汤孟平，等.天目山针阔混交林中枫香的结构特征［J］.浙江农林大学学报，2012，29（6）：867-874.

ZHANG Zhihua，WEIXinliang，TANG Mengping，etal.Structural characteristics ofLiquidambarformosanafor a mixed coniferous-broadleaf forest in Mount Tianmu［J］.JZhejiangA&FUniv，2012，29（6）：867-874.

［10］OHSAWA M.Difference of vegetation zones and species strategies in the subalpine region of Mt.Fuji［J］.Vegetation，1984，57（1）：15-52.

［11］JOURNEL A G，HUIJBREGTSCh J.MiningGeostatistics［M］.Salt Lake：Academic Press，1978.

［12］LIEBHOLD A M，ROSSIR E，KEMPW P.Geostatistics and geographic information systems in applied insect ecology［J］.AnnuRevEntomol，1993，38（1）：303-327.

［13］王政权.地统计学及在生态学中的应用［M］.北京：科学出版社，1999.

［14］汤孟平，周国模，施拥军，等.天目山常绿阔叶林优势种群及其空间分布格局研究［J］.植物生态学报，2006，30（5）：743-752.

TANGMengping，ZHOU Guomo，SHIYongjun，etal.Study of dominant plant populations and their spatial patterns in evergreen broadleaved forest in Tianmu Mountain，China［J］.ActaPhytoecolSin，2006，30（5）：743-752.

［15］冯益明，唐守正，李增元.空间统计分析在林业中的应用［J］.林业科学，2004，40（3）：149-155.

FENG Yiming，TANG Shouzheng，LIZengyuan.Application of spatial statistic analysis in forestry［J］.SciSilv Sin，2004，40（3）：149-155.

Spatial continuity for DBH in dominant populations of an evergreen broadleaved forest in NationalNature ReserveofMount Tianmu，China

FANGGuojing1，TANGMengping2
（1.Monitoring Center for Forest Resources in Zhejiang Province，Hangzhou 310020，Zhejiang，China；2.School of Environmental and Resource Sciences，Zhejiang A&FUniversity，Lin'an 311300，Zhejiang，China）

To reveal structural characteristics，a study on spatial continuity for DBH of dominant plant populations in an evergreen broadleaved forest was conducted.A typical plot，100 m×100 m，in National Nature Reserve of Mount Tianmu，Zhejiang was surveyed for tree diameter at breast height（DBH）above 5 cm with Total Station（Leica TCR702Xrange）that is used to measure tree coordinates.Also，dominant plant populations were first determined using the Dominance Index and then completed using a spatial statistical analysis with GS+software.Results showed that DBH of the dominant plant populations of evergreen broadleaved trees in the reserve had a spatial dependence of 0.7 and a range of spatial continuity of 5.7 m. Thus，the structural complexity for this evergreen broadleaved forest showed a high spatial dependence and a small range of spatial continuity.［Ch，2 fig.3 tab.15 ref.］

forest ecology；evergreen broadleaved forest；dominant species，geostatistics；spatial continuity；Mount Tianmu

S718.54

A

2095-0756（2014）05-0663-05

2013-12-09；

2014-03-13

“十二五”国家科技支撑计划项目（2012BAD22B0503）

方国景，工程师，从事森林经理及林业工程规划等研究。E-mail：fangguojing2008@163.com。通信作者：汤孟平，教授，博士，博士生导师，从事森林经理学等研究。E-mail：goodtmp@yahoo.com.cn