武秀娟,奥小平,赵育鹏,崔雪晴
(1. 山西省林业和草原科学研究院,山西 太原 030012;2. 国家林业和草原局 调查规划设计院,北京 100714)
林分结构综合反映了林分发展的过程,也是森林未来发展的驱动源泉,已成为分析和管理森林生态系统的重要因子[1]。林分结构包括非空间结构和空间结构,非空间结构用来描述林分平均特征,不受林分相对位置的约束,而空间结构是林木在林地上的分布格局及其属性在空间上的排列方式[2],空间结构在很大程度上决定了林分的稳定性、发展的可能性和经营空间的大小[3],对森林的演替趋势具有决定性作用[4],因此,空间结构已成为森林结构研究的热点[5]。分析森林群落的空间结构特征,有助于科学制定合理的经营措施。芦芽山是汾河和桑干河的发源地,是中国暖温带残存的天然次生林分布区中保存较完整的地区之一,存在大面积完好、典型的以华北落叶松Larixprincipis-rupprechtii和云杉Piceaspp.为建群种组成的属于演替后期相对稳定的群落[6]。前人针对该林分已做了许多研究,其内容主要集中在群落的种间联结性[7]、林内土壤温度及含水量的季节差异[8]、树木径向生长[9−10]及其对气候变化的响应[11−12]等,对其林分空间结构的研究仍鲜见报道,因此,本研究以该地华北落叶松-云杉天然次生林为对象,通过计算其混交度、角尺度、大小比数,以及林分空间结构指数研究不同海拔该林分的结构特征,旨在为该地华北落叶松-云杉天然次生林的科学经营提供依据,同时为相似立地条件下人工林空间结构的优化经营奠定理论基础。
研究地位于山西省忻州市宁武县境内的芦芽山国家级自然保护区内 (38°36′~39°02′N,111°46′~112°54′E),保护区面积21 453 hm2,海拔1 346~2 787 m。该地属暖温带半湿润区,具明显的大陆性气候特点,夏季凉爽多雨,冬季寒冷干燥。据距其最近的五寨气象站1957−2007年资料,其年平均气温为4.9 ℃,1月平均气温为−10.2 ℃,7月平均气温为22.4 ℃,年降水量472.7 mm,年蒸发量1 800 mm,年均相对湿度55%,无霜期130 d。土壤主要是在残积和坡积母质上发育起来的山地褐土和山地淋溶褐土等。海拔1 850~2 650 m是以华北落叶松和云杉属Picea植物白扦Picea meyeri、青扦Picea wilsonii为建群种组成的寒温性针叶林占优势的稳定森林群落。
在芦芽山国家级自然保护区阴坡海拔1 950~2 550 m范围内选取起源相同、林龄相近,未经人为扰动的华北落叶松-云杉天然次生林,按1 950、2 050、2 150、2 350和2 550 m海拔梯度布设样地,每个海拔4个样地,样地面积为20 m×20 m。每个样地内,对胸径5 cm以上的活树进行每木检尺,记录树种、胸径、树高、枝下高、冠幅,并对每株树进行定位。同时在每个海拔的样地内选取主要优势树种的3株优势木,采用生长锥在树干基部钻取髓心,根据3株优势木年龄的均值确定林龄。样地基本概况见表1。
表1 样地基本概况Table 1 Survey of sample plots
芦芽山阴坡不同海拔华北落叶松-云杉天然次生林按相对多度筛选出群落优势度较高的物种仅为华北落叶松和云杉,因此在针对树种进行空间结构计算时,仅针对优势树种进行。群落及各树种平均角尺度及其参数分布见表2,仅海拔2 150 m处群落的平均角尺度为0.51,其余海拔处林分的平均角尺度为0.52~0.53,可见,华北落叶松-云杉天然次生林几乎均为聚集分布型,但聚集程度较低。华北落叶松-云杉天然次生林群落的角尺度整体近似服从正态分布,不同海拔间未表现出明显差异。5个海拔群落中Wi=0的频率很低,甚至不存在,说明该林分中林木很均匀分布的情况极少,而Wi=0.5的频率均为最高,占比均在40%以上,说明该天然林分中林木大多处于随机分布状态,Wi=0.75与Wi=1.00角尺度的频率之和稍大于Wi=0与Wi=0.25角尺度的频率之和,说明不均匀分布的状态稍多于均匀分布的状态。
表2 不同海拔树种和群落角尺度特征Table 2 Species angle index of L. principis-rupprechtii-Picea spp. secondary forests
不同海拔的2种优势树种的角尺度特征表现不尽相同,华北落叶松的角尺度为0.46~0.65,云杉的角尺度为0.44~0.56,华北落叶松在2 350 m处表现为均匀分布,低海拔以聚集分布为主,高海拔以随机分布为主;云杉则在2 050 m处呈均匀分布,低海拔以随机分布为主,高海拔以聚集分布为主。随海拔的升高,华北落叶松分布株数呈先增加后减少的趋势,云杉则呈先减少后增加的趋势,且华北落叶松和云杉的株数分别在海拔2 150和2 550 m达到最大,在其株数较多的海拔均呈聚集分布。2个树种的角尺度均在Wi=0.50时的分布频率最大,优势树种在林分中都属于随机分布。
芦芽山阴坡华北落叶松-云杉天然次生林各海拔群落平均大小比数为0.47~0.52,整体均接近中庸状态,说明林内林木胸径的差异不明显,大小分化不严重,且不同海拔群落各大小比数所占的频率相差不大,约20%(表3)。
表3 不同海拔树种的大小比数特征Table 3 Species neighborhood comparison of L. principis-rupprechtii-Picea spp. secondary forests
2种优势树种大小比数在不同海拔表现不同。随着海拔的升高,华北落叶松先由亚优势状态发展为中庸状态,后由中庸状态发展为优势状态;云杉则先由近亚优势状态发展为劣势状态,后发展为中庸状态,且2种优势树种均在其株数达到最大的海拔处表现为中庸状态。不同海拔华北落叶松在大小比数Ui=0或Ui=0.25时的分布频率均最大,云杉除在1 950 m时大小比数在Ui=0.25时的分布频率最大外,其余海拔均为Ui=0.75或Ui=1.00时分布频率最大。可见,华北落叶松在不同海拔的林分中均处于优势或亚优势状态,云杉处于劣势状态,华北落叶松的竞争力明显大于云杉。
随着海拔高度的上升,华北落叶松-云杉天然次生林群落的平均混交度呈下降趋势,可见随海拔升高混交程度由强变弱,树种的组成趋于单一化。群落的混交度频率分布规律在不同海拔间差异较大(表4),海拔1 950 m时,以中强度混交为主,中强度混交比例之和达到86%以上,海拔 2 050 m以上,以零度和弱度混交为主。零度和弱度混交比例之和达到60%以上,且混交度Mi=0的频率随海拔升高逐渐增加,可见优势树种随海拔增加单种聚集现象加剧,进而导致了高海拔处林分整体处于弱度混交。
表4 不同海拔树种和群落的混交度特征Table 4 Species mingling degree of L. principis-rupprechtii-Picea spp. secondary forests
随着海拔的升高,华北落叶松的混交状态由强变弱,后又变强,云杉则除2 150 m海拔处为强度混交为主外,其余海拔均为中弱度混交为主。可能是随着海拔的升高,华北落叶松的数量先增加后减少,云杉则表现为先减少后增加,进而其数量在增加时的海拔处均表现出了明显的单种聚集现象,因而呈弱度混交状态。
林分空间结构指数能够很好地反映不同林分所处的状态。图1可看出:随着海拔的升高,芦芽山华北落叶松-云杉天然次生林的林分空间结构指数逐渐减小,可见低海拔林分的空间结构较好。
图1 不同海拔的林分空间结构指数Figure 1 Stand spatial structure index of L. principis-rupprechtii-Picea spp. secondary forests
芦芽山阴坡华北落叶松-云杉天然次生林群落的空间分布格局总体上以聚集分布为主,但聚集程度较低,接近随机分布。一般情况下,群落不受严重干扰,经过演替后形成的顶级群落的水平分布格局为随机分布[13−14],可见,芦芽山阴坡华北落叶松-云杉天然次生林正处在演替后期向顶级群落演替的阶段。各海拔的林木均以随机分布的占比最大,比例均大于40%,不均匀分布的比例大于均匀分布的比例,这与柏云龙等[15]及夏富才等[16]研究结果相同。本研究中优势树种华北落叶松和云杉虽在各海拔的空间分布格局不甚相同,角尺度范围为0.44~0.65,但均接近随机分布,进一步证实了该林分正处于向顶级群落类型演替的后期阶段,与赵鸣飞等[6]对芦芽山寒温性针叶林的研究结果一致。
各海拔群落大小分化程度差异不大,且整体均接近中庸状态。林木大小分化程度不大,这与肖欣等[17]和韩金等[18]研究结果相同。2个优势树种的角尺度均在Wi=0.50时分布频率最大,说明优势树种个体大多数呈随机分布的情况。华北落叶松的大小比数在所有海拔均小于云杉,表明该群落中华北落叶松的竞争力大于云杉,但在高海拔华北落叶松的取值仅在0上分布,原因是由于华北落叶松数量较少但胸径都很大,在结构单元中处于明显的优势地位。
华北落叶松-云杉天然次生林群落的平均混交度随海拔的升高呈下降趋势,结果与实际相符。该群落以华北落叶松和云杉为主要优势树种,低海拔时混生有少量的白桦Betula platyphylla、红桦B. albosinensis、辽东栎Quercus wutaishanica等树种,随着海拔的升高,混生树种逐渐减少,甚至消失,进而导致群落的混交度下降。
芦芽山阴坡华北落叶松-云杉天然次生林群落的角尺度、林木大小分化程度在不同海拔间均无显著性差异,而混交度随海拔升高呈明显的下降趋势,可见,该林分不同海拔空间结构的差异主要表现在混交度指数上。陶观护等[19]对不同密度油松Pinus tabulaeformis飞播林的林分空间结构的研究中表明,林分密度对混交度有着显著影响。芦芽山阴坡华北落叶松-云杉天然次生林不同海拔间林分空间结构的差异可能源于其林分密度的差异,2个优势树种在不同海拔间的密度变化导致其空间结构指数的差异也证实了这一观点。但是否存在其他因素(如坡度等)导致不同海拔间林分空间结构的差异仍需进行深入研究。
空间结构指数是反映林分空间结构综合情况的指标,数值越大,说明该林分的空间结构越优化。芦芽山阴坡华北落叶松-云杉天然次生林群落以海拔1 950 m时林分整体结构最好。一般认为,林分空间结构越优,林分的功能越强,稳定性越高[20]。故海拔1 950 m时的华北落叶松-云杉天然次生林的稳定性也最高。因此,在林分经营过程中,针对以混交度为明显差异造成的空间结构指数不同的华北落叶松-云杉天然次生林,可选择以1 950 m处的林分结构为借鉴对高海拔的林分进行结构调整,以提高林分稳定性。高海拔华北落叶松的胸径优势非常明显,占据着群落的最高层,但其株数比例较小,林下更新困难,因此在高海拔,需采取生态抚育措施对混交度为0和0.25的林木进行调整,改善混交比例,降低单种聚集现象,同时促进华北落叶松的更新,以保证物种多样性和群落的持续稳定发展。