文章信息
- 仇建习, 汤孟平, 娄明华, 沈利芬, 庞春梅, 邹奕巧
- QIU Jianxi, TANG Mengping, LOU Minghua, SHEN Lifen, PANG Chunmei, ZOU Yiqiao
- 基于Hegyi改进模型的毛竹林空间结构和竞争分析
- Analysis of the spatial structure and competition with a Phyllostachys edulis standbased on an improved Hegyi model
- 生态学报, 2016, 36(4): 1058-1065
- Acta Ecologica Sinica, 2016, 36(4): 1058-1065
- http://dx.doi.org/10.5846/stxb201406131227
-
文章历史
- 收稿日期: 2014-06-13
2. 浙江省森林生态系统碳循环与固碳减排重点实验室, 临安 311300;
3. 莲都区林业局, 丽水 323000;
4. 天目山国家级自然保护区管理局, 临安 311300;
5. 天台县林业特产局, 台州 318000
2. Zhejiang Provincial Key Laboratory of Carbon Cycling in Forest Ecosystems and Carbon Sequestration, Linan 311300, China;
3. Liandu District Bureau of forestry, Lishui 323000, China;
4. Authority in National Nature Reserve of Tianmu Mountain, Linan 311300, China;
5. Tiantai County Bureau of forestry specialty, Taizhou 318000, China
在有限资源空间里,竞争无处不在。树木竞争意味着资源不足以供给2株及以上树木的最佳生长能力[1]。学者们已提出多种衡量树木竞争能力大小的单木竞争指数[2,3,4,5,6,7,8,9]。根据是否考虑距离因素,可分为2类:距离无关竞争指数和距离相关竞争指数[1,10]。距离无关竞争指数未考虑空间位置因素[11,12]。通常用与距离有关的竞争指数定量描述树木竞争状况[13]。竞争指数的选取应兼顾两个方面即理论上的合理性和实际中的可操作性[14]。Hegyi竞争指数[4](以下简称Hegyi指数,用CH表示)容易测定、计算简单,是最常用的竞争指数[15,16]。
汤孟平等[16]已将Hegyi指数引入毛竹林的空间结构研究,但该指数的应用有待改进。采用Hegyi指数计算乔木林分竞争时,林木胸径作为反映林木竞争能力的主要参数,可以直接反映了林木竞争能力随年龄的动态变化。但是,毛竹(Phyllostachys edulis)属于禾本科竹亚科竹种,无次生形成层组织,胸径在毛竹秆形生长结束后停止增长。因此,毛竹的胸径是一个相对静态的结构因子,不能反映毛竹竞争能力的动态变化。为了说明运用Hegyi指数分析毛竹林竞争的不足,建立2个典型的对象竹与其4株最近相邻竹的具有相同胸径和相等距离,但年龄不同的空间结构单元,使用Hegyi指数分别计算2个不同年龄结构下空间结构单元中对象竹的竞争强度。图1中,每个空间结构单元中间1株表示对象竹,周围4株表示竞争竹,不同的形状表示不同年龄毛竹,圆表示9年生毛竹,五边形表示6年生毛竹。假定毛竹的胸径均为10cm,对象竹与竞争竹之间的距离均为1m。采用Hegyi指数计算图1中2个不同年龄结构下空间结构单元的对象竹竞争强度为CHA=CHB=4,即对象竹所受的竞争压力相同。说明Hegyi指数不能区分在对象竹与竞争竹的距离相等和胸径相同的条件下,不同年龄结构下空间结构单元之间的竞争差异。因此,直接引用Hegyi指数分析毛竹竞争关系是不全面的,应对其进行改进。
事实上,毛竹年龄作为毛竹林的主要结构因子,是反映毛竹自身生理活动强度动态变化的重要指标,是影响毛竹竞争能力变化的重要因素。毛竹从出笋成竹、成长、壮大、衰落和死亡,各个时期的毛竹生理活动强度存在较大差异[17],即吸收矿质元素和水分、制造和积累有机物质等生理活动的能力有强弱之分。依据不同年龄毛竹生理活动强度不同的特点,可将毛竹的材质生长周期分为3个阶段,即增进期、稳定期和下降期。3个阶段的年龄范围分别为:2—5a为增进期,6—8a为稳定期,9a以上为下降期[17]。由于毛竹秆形生长结束即进入材质生长期,因此1年生也属于增进期。本文采用以上结论。
本研究在Hegyi指数基础上,根据毛竹的生物学特性,从毛竹的生长过程和机理角度来分析毛竹间的竞争关系,改进Hegyi指数,提出增加毛竹年龄的竞争分析指数——竞争势(CP),为毛竹林空间结构竞争分析提供了新的分析指标。
1 研究区概况浙江天目山国家级自然保护区位于东经119°23′47″—119°28′27″,北纬30°18′30″—30°24′55″。地处东海之滨,浙江省西北部杭州市临安市境内的西天目山(图2),距杭州市区94km。保护区位于中亚热带北缘至北亚热带过渡地带,受海洋暖湿气候影响较深,雨量充沛,年平均相对湿度71%—76%,年平均降水量1390—1870mm,年平均气温8.8—14.8 ℃,≥10 ℃年积温2500—5100 ℃。植被垂直分布明显,自山麓到山顶垂直带谱为:海拔870m以下为常绿阔叶林;870—1100 m为常绿、落叶阔叶混交林;1100—1380m 为落叶阔叶林;1380—1506m为落叶矮林。近自然毛竹林主要分布于海拔350—900m,其林下植被稀少,主要有茶(Camellia sinensis)、微毛柃(Eurya hebeclados)、豹皮樟(Litsea coreana var. sinensis)、毛花连蕊茶(Camellia fraterna)、细叶青冈(Cyclobalanopsis gracilis)、小叶青冈(Cyclobalanopsis myrsinifolia)、马银花(Rhododendron ovatum)和短尾柯(Lithocarpus brevicaudatus)等。
2 研究方法 2.1 样地调查2 009年7月,在浙江省天目山国家级自然保护区中的禁伐区内选择少受干扰的毛竹纯林(无采伐、垦翻、施肥等经营活动),设置1个100m×100m的固定样地,样地中心海拔840m,主坡向为南偏东30°。利用南方全站仪NTS355,将样地划分100个10m×10m的网格单元,在每个网格单元内,对所有毛竹用钢围尺测胸径和竹秆皮色法识别年龄,利用南方全站仪NTS355测定每竹的3维坐标(X,Y,Z)。以后每年进行1次复查,复查老竹的状态,并对新竹进行每竹调查及测定3维坐标。至2013年,已连续调查5次。
2.2 边缘校正Stage and Ledermann[11]和Ledermann[12]认为距离有关竞争指数不包括样地边缘木,因此必须进行边缘校正。采用距离缓冲区法进行校正,就是在原样地四条边界线的内侧设置带状缓冲区,中间核心区的林木均作为对象木参与林分信息的计算[18]。本文设定缓冲区宽度为5m,90m×90m为中间核心区。
2.3 Hegyi指数的改进竞争势(CP)是在Hegyi指数(CH)基础上,增加考虑毛竹材质生长周期中各个阶段不同的生理活动强度。实际上,竞争势反映了对象竹所承受来自竞争竹的竞争压力。对某一对象竹而言,对象竹的材质生长期越晚,其取值越小,生理活动能力越弱,或对象竹的胸径越小,均导致竞争势越大,来自周围竞争竹的竞争压力越大,在竞争中越处于不利地位。计算公式如下:
竞争势可定义为立竹竞争强度与对象竹所对应材质生长期及胸径成反比,与竞争竹所对应材质生长期以及胸径成正比,且随对象竹与竞争竹之间距离增大而减小。利用公式(1)分析图1中两个不同的空间结构单元竞争强度,得CPA=8,CPB=2,而CHA=4,CHB=4,说明与Hegyi指数相比,竞争势可以有效区分不同年龄结构空间结构单元的竞争差异。
合理的空间结构单元有助于深入分析毛竹林空间结构动态。惠刚盈等[19]指出由对象木与其4株最近相邻木可以构成最佳空间结构单元(简称“1+4”结构),汤孟平等[20]认为对象竹周围有4株最近邻竹时,最有可能获得高产。因此,本文采用“1+4”结构即空间结构单元中最近邻竹株数即竞争竹株数n=4。
3 结果与分析 3.1 毛竹林分生长因子动态变化2 009—2013年毛竹主要生长因子动态变化表见1。可见,毛竹林大小年现象较明显,毛竹大年(2010、 2 012)新竹株数远高于小年(2009、2011、2013)。2009年,毛竹的总株数、断面积、平均胸径均最低,这是因为除2009年是毛竹小年的内在因素以外,2008年雪灾造成大量毛竹死亡的环境因素[21],综合导致2009年毛竹林分主要生长因子普遍较低。2010年的新竹株数、总株数和断面积均为最高。毛竹年龄划分为龄级,即1周年毛竹为1龄级,2—3a毛竹为2龄级,4—5a毛竹为3龄级,以此类推。由表2可知,毛竹各龄级所占比例随年份发生动态变化。2009—2010年毛竹林包含8个龄级,2011—2012年为9个龄级,2013年则为10个龄级,其中2009—2011年有4个龄级毛竹所占比例大于10%,2012—2013年则有5个龄级毛竹所占比例大于10%,说明毛竹林龄级趋于多样化,龄级比例趋于均匀化。2009—2012年中存在9龄级或10龄级毛竹比例为0的现象,原因可能在于老龄竹生活力较差、力学强度较低,雪灾进一步加剧了老龄竹的死亡。
年份 Year | 新竹株数 New growth number | 总株数 All number | 断面积 Basal area/(m2/hm2) | 胸径Diameter at breast height/cm | |
均值 Mean value | 标准差 Standard Deviation | ||||
2009 | 44 | 5488 | 51.484 | 10.8 | 1.58 |
2010 | 1304 | 6561 | 63.325 | 11.0 | 1.61 |
2011 | 4 | 5981 | 58.604 | 11.1 | 1.57 |
2012 | 763 | 6094 | 60.926 | 11.2 | 1.57 |
2013 | 4 | 5757 | 57.984 | 11.2 | 1.55 |
年份Year | 龄级分布Age class distribution/% | |||||||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | |
2009 | 0.8 | 11.1 | 14.4 | 30.5 | 3.6 | 6.9 | 32.5 | 0.2 | 0 | 0 |
2010 | 19.9 | 9.6 | 11.9 | 24.2 | 3.0 | 5.6 | 25.6 | 0.2 | 0 | 0 |
2011 | 0.1 | 21.4 | 9.5 | 12.4 | 24.3 | 3.0 | 5.6 | 23.5 | 0.2 | 0 |
2012 | 12.5 | 19.7 | 9.5 | 11.6 | 20.6 | 2.7 | 4.9 | 18.4 | 0.1 | 0 |
2013 | 0.1 | 13.1 | 20.6 | 9.4 | 11.8 | 20.2 | 2.6 | 4.6 | 17.5 | 0.1 |
利用竞争势(CP)分析2009—2013年5个调查年度的竹林竞争动态,并用C#程序进行编程计算,绘制得到的毛竹林竞争强度随年份的动态变化图(图3)。可见毛竹大年(2010,2012)的林分竞争强度明显大于小年(2009,2011,2013)。采用单因素方差分析方法检验(One-way ANOVA)各年份林分竞争强度的差异,结果表明,毛竹大年与小年之间的林分竞争强度存在极显著性差异(P<0.01),毛竹大年之间林分竞争强度相差不大,方差检验结果表现为无显著性差异(P>0.05),毛竹小年中2011与2013年的林分竞争强度同样无显著性差异(P>0.05),但是,2009年与其它年份的林分竞争强度均存在极显著性差异(P<0.01)。该结果的原因在于,2008年雪灾造成大量毛竹死亡,直接导致林分的竞争强度降低,且2009年是毛竹小年,对林分竞争强度的增强影响不大。因此,新竹和自然灾害(如雪灾)是导致毛竹林竞争强度变化的主要原因。
3.3 林分不同径阶、龄级竞争动态对象竹个体的大小和年龄能反映其竞争能力[17]。对象竹胸径按2cm径阶整化,年龄按龄级划分,分别绘制2009—2013年5个调查年度对象竹竞争强度与其对应径阶和对应龄级的关系图(图4,图5)。各年份毛竹林分竞争强度均呈现随对象竹径阶的增大而呈现逐渐减小的趋势,且各年份对象竹的竞争强度与径阶均服从幂函数关系,这个结果与采用Hegyi指数计算得到的乔木林竞争强度与胸径的关系基本一致[22,23,24]。不同年份对象竹竞争强度与径阶之间的拟合决定系数R2的值不同。4径阶毛竹常年处于竞争强度较大的状态,主要原因在于其杆形弱小,受到周围竞争竹的竞争压力较大,在营养空间的竞争中处于绝对劣势。因此可根据这一特性,在竹林抚育过程中,对4径阶毛竹进行采伐,调节和改善立竹之间的竞争关系。
各年份毛竹林分竞争强度与对象竹龄级的关系呈现随龄级的增大而逐渐增大的趋势,且各年份对象竹的竞争强度与龄级均服从线性函数关系,各年份竞争强度与龄级之间的拟合决定系数R2的值不同。5龄级对象竹竞争强度随年份大致呈现逐渐下降的趋势,到2013年基本降低至各龄级对象竹竞争强度的平均水平,而5龄级以上对象竹的竞争强度常年处于较高状态,说明5龄级竹以上对象竹始终受周围相邻竹的竞争压力较大,在竞争中处于极度不利地位。这一结果可为毛竹林空间结构优化调控中将5龄级以上毛竹全部采伐的择伐方式提供参考依据[25]。
4 结论与讨论(1)毛竹大小年现象明显,毛竹大年新竹所占比例明显高于毛竹小年。2012—2013年毛竹龄级数多于2009—2011年,且各龄级株数比例更加均匀。
(2)毛竹林竞争强度随年份发生动态变化,呈现出毛竹小年向大年过渡年份林分竞争强度增加,大年向小年过渡年份林分竞争强度减小的现象。毛竹小年的林分竞争强度明显弱于大年,大小年之间林分竞争强度存在极显著性差异(P<0.01)。毛竹大年之间的林分竞争强度无显著性差异(P>0.05),小年2011和2013之间的林分竞争强度无显著性差异(P>0.05),但2009年与其它年份的林分竞争强度均存在极显著性差异(P<0.01)。新竹株数和自然灾害是导致毛竹林分竞争强度变化的主要原因。
(3)采用竞争势计算得到的毛竹林竞争强度,可以有效全面地反映毛竹林竞争强度与立竹胸径及立竹年龄的密切关系。各年份毛竹林竞争强度均随着对象竹径阶的增大而减小,两者之间服从幂函数关系,这与学者们在研究乔木竞争时采用Hegyi指数得到的结论相同[22,23,24]。各年份毛竹林竞争强度均随着年龄的增加而增加,两者之间服从线性关系。上述结论均符合竞争势所描述的对象竹竞争强度与对象竹径阶成反比,与对象竹年龄呈正比的特性,进而说明立竹胸径越大,并且年龄越小,在林分中所处的竞争位置越有利。
毛竹不同于树木,它是克隆植物,地上分离的竹子在地下是通过竹鞭相连的,而且相连的竹子间具有资源共享作用。因此,相邻竹子之间不仅有竞争关系,还有相互的生理整合影响。这两方面作用的结果可以通过毛竹个体的大小差异和生长状况反映。Hegyi指数是一种简便易测的指标,它仅考虑了个体大小和立竹之间的距离,,同时还从毛竹生长机理角度出发,根据不同年龄阶段毛竹生理活动强度的不同,增加了材质生长期因子,更深层次地反映立竹之间的竞争关系,是分析毛竹林竞争关系的较理想指标,可为毛竹林经营过程中优化调控竹林竞争关系提供一定的理论指导。
毛竹的生长包括竹笋地下生长、笋-幼竹生长和成竹生长3个阶段[17,26],本文仅分析了毛竹林成竹生长阶段的竞争及其动态特征,关于地下生长阶段和笋-幼竹生长阶段的竞争及其动态有待深入研究。
[1] | Rivas J J C, González J G A, Aguirre O, Hernández F J. The effect of competition on individual tree basal area growth in mature stands of Pinus cooperi Blanco in Durango (Mexico). European Journal of Forest Research, 2005, 124(2):133-142. |
[2] | Bella I E. A new competition model for individual trees. Forest Science, 1971, 17(3):364-372. |
[3] | Hegyi F. A simulation model for managing jack pine stands//Fries J. Growth Models for Tree and Stand Simulation. Stockholm:Royal College of Forestry, 1974:74-90. |
[4] | Alemdag I S. Evaluation of Some Competition Indexes for the Prediction of Diameter Increment in Planted White Spruce. Ottawa:Forest Management Institute, 1978. |
[5] | Wykoff W R, Crookston N L, Stage A R. User's Guide to the Stand Prognosis Model. Ogden:Intermountain Forest and Range Experiment Station, 1982. |
[6] | Holmes M J, Reed D D. Competition indices for mixed species northern hardwoods. Forest Science, 1991, 37(5):1338-1349. |
[7] | Reineke L H. Perfecting a stand density index for even-aged forest. Journal Agricultural Research, 1993, 46(7):627-638. |
[8] | Rouvinen S, Kuuluvainen T. Structure and asymmetry of tree crowns in relation to local competition in a natural mature Scots pine forest. Canadian Journal of Forest Research, 1997, 27(6):890-902. |
[9] | Schröder J, von Gadow K. Testing a new competition index for Maritime pine in northwestern Spain. Canadian Journal of Forest Research, 1999, 29(2):280-283. |
[10] | Contreras M A, Affleck D, Chung W. Evaluating tree competition indices as predictors of basal area increment in western Montana forests. Forest Ecology and Management, 2011, 262(11):1939-1949. |
[11] | Stage A R, Ledermann T. Effects of competitor spacing in a new class of individual-tree indices of competition:semi-distance-independent indices computed for Bitterlich versus fixed-area plots. Canadian Journal of Forest Research, 2008, 38(4):890-898. |
[12] | Ledermann T. Evaluating the performance of semi-distance-independent competition indices in predicting the basal area growth of individual trees. Canadian Journal of Forest Research, 2010, 40(4):796-805. |
[13] | 汤孟平. 森林空间结构研究现状与发展趋势. 林业科学, 2010, 46(1):117-122. |
[14] | 关玉秀, 张守攻. 竞争指标的分类及评价. 北京林业大学学报, 1992, 14(4):1-8. |
[15] | 刘彤, 李云灵, 周志强, 胡海清. 天然东北红豆杉(Taxus cuspidata)种内和种间竞争. 生态学报, 2007, 27(3):924-929. |
[16] | 汤孟平. 森林空间经营理论与实践. 北京:中国林业出版社, 2007. |
[17] | 江泽慧. 世界竹藤. 沈阳:辽宁科学技术出版社, 2002. |
[18] | 周红敏, 惠刚盈, 赵中华, 胡艳波. 林分空间结构分析中样地边界木的处理方法. 林业科学, 2009, 45(2):1-5. |
[19] | 惠刚盈, Gadow K V. 森林空间结构量化分析方法. 北京:中国科学技术出版社, 2003. |
[20] | 汤孟平, 徐文兵, 陈永刚, 邓英英, 赵明水. 天目山近自然毛竹林空间结构与生物量的关系. 林业科学, 2011, 47(8):1-6. |
[21] | 尹新华, 翁益明, 董云富, 傅红梅, 吴礼栋. 毛竹受雨雪冰冻危害的受损特点. 浙江林学院学报, 2008, 25(6):823-827. |
[22] | 邹春静, 徐文铎. 沙地云杉种内、种间竞争的研究. 植物生态学报, 1998, 22(3):269-274. |
[23] | 段仁燕, 王孝安. 太白红杉种内和种间竞争研究. 植物生态学报, 2005, 29(2):242-250. |
[24] | 史宇, 余新晓, 岳永杰, 王小平, 秦永胜, 陈俊崎. 北京山区天然侧柏林种内竞争研究. 北京林业大学学报, 2008, 30(2):36-40. |
[25] | 汤孟平, 徐文兵, 陈永刚, 娄明华, 仇建习, 庞春梅, 赵明水. 毛竹林空间结构优化调控模型. 林业科学, 2013, 49(1):120-125. |
[26] | 何奇江. 毛竹三笋栽培新技术. 杭州:浙江科学技术出版社, 2010. |