毛竹竹秆直径与材积的垂直格局及其模型.docxVIP

毛竹竹秆直径与材积的垂直格局及其模型 竹子有其独特的生物学特性：生长的第一年完成体积的增加。在随后的生长过程中，质量生长继续变化，体积几乎没有变化。在以往竹类植物生物量的相关研究中[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16],构建生物量与胸径关系模型时很少考虑到胸径不变而生物量在变这种特性,直接建立生物量与胸径、竹高的指数模型、线性模型或多项式模型来预测生物量,这种方法所建模型的相关系数偏低,预测效果受到影响,而用竹秆材积建模的研究相对较少。笔者应用竹秆材积这个基本不变的量来反应毛竹的生物量水平,在材积与质量之间转换存在一个中间变量即密度因子,而竹秆的密度是随时间而变的,需要找出毛竹竹秆密度的时间变化函数,材积转换生物量就可完全实现。关于材积模型,在其他物种上也进行过部分研究,针对毛竹的生物学特性,为了消除胸径和高度大小对模型的影响,笔者采用了相对值概念,建立了毛竹秆的相对直径、相对外围累计材积、相对实际累计材积和相对节隔累计材积与相对高度的关系模型,外围材积、实际材积和节隔材积与胸径的关系模型,旨为毛竹秆型结构研究、竹秆材积换算提供理论依据,补充毛竹林生产力研究内容。 1 材料和方法 1.1 样竹与病虫害的关系 福建省永安市竹林面积4.9787万hm2,农民人均拥有竹林面积0.28hm2。研究中选取的毛竹来自福建省永安市天宝岩自然保护区杨梅坪山地竹林地(117°29′08″E,25°53′49″N),海拔705m。按胸径大小选取8~9cm毛竹3株,9~10cm毛竹3株,10cm以上3株共9株。样竹采伐时间为2005年6月5日,选定的样竹要求竹梢保持完好无缺损,顶端不能出现任何断梢现象,竹秆圆而通直,没有任何的病虫害。伐竹时用锯齐地伐倒,并避免伐竹和拖运等作业造成竹秆和竹梢的损伤。 1.2 密度相对值公式 依据微积分理论按长度比把样竹分成若干等份,然后对每一等份进行调查,获取样竹的原数据,测量每一段底部外围圆周长(皮尺,精度为1mm)、壁厚(按4个方向测定,测定工具为游标卡尺,精度为0.2mm)、竹内每个节隔的厚度。考虑到精度要求,每株样竹分成40段,9株样竹共360个竹段。 对于第i竹段而言,底端直径(Di)、底端外半径(r外i)、底端内半径(r内i)、竹段外围材积(V外i)、内腔体积(V内i)、节隔材积(V隔i)、实际材积(V实i)、到第i段的外围累计材积(V′外i)、到第i段的实际累计材积(V′实i)、到第i段的节隔累计材积V′隔i应用以下公式分别求出(其中节隔材积指毛竹内部节隔的体积):Di=Liπ;(1)Di=Liπ;(1) r外i=Di2;(2) r内i=r外i-T壁i; (3) V外i=π3×(r2外i+r外i×r外i+1+r2外i+1)×h;(4) V内i=π3×(r2内i+r内i×r内i+1+r2内i+1)×h;(5) V隔i=π×(r内i+r内i+12)2×Τ隔i;(6) V实i=V外i-V内i+V隔i; (7) V′外i=i∑n=1V外n;(8) V′实i=i∑n=1V实n;(9) V′隔i=i∑n=1V隔n。 (10) 公式(1)~(10)中,Li为第i竹段底端外围周长(i=1,2,3,…,40),T壁i为第i竹段底端壁厚(4个实测值的平均值),T隔i为第i竹段所有节隔累计厚度,h为竹段的高(竹秆全长的1/40)。 用于分析的胸径值取竹秆离地1.3m处直径,直径相对值表示每一竹段底端直径与竹秆胸径的比值,竹秆相对高度表示从竹秆底端至第i竹段上端高度与竹秆全长的比值,相对外围累计材积表示从竹秆第1段到第i段外围累计材积与竹秆外围总材积的比值,相对实际累计材积表示竹秆第1段到i段实际累计材积与竹秆实际总材积的比值,相对节隔累计材积为竹秆第1段到i段节隔累计材积与竹秆节隔总材积的比值,式中相对值均是无量纲值(采用相对值主要目的是消除胸径和高度差异对模型的影响,以取得毛竹共性特征)。 数据处理和作图在Excel和SPSS13.0中进行。 2 结果与分析 2.1 生长参数回归模型 在不同毛竹个体之间,个体大小存在一定的变化,但它们的形态结构具一定规律性,为了消除胸径和高度差异对模型的影响,在分析毛竹秆直径垂直格局时,采用相对值进行数据处理和统计分析,直径的相对值为实测值与胸径的比,相对高度为竹秆基点到竹秆某一高度与竹全高的比。 应用SPSS统计软件对毛竹竹秆直径相对值与竹秆相对高进行回归分析,建立了两者的关系模型: Y1=1.164 — 1.654X+1.645X2— 1.153X3(R2=0.997), (11) 式中:Y1表示竹秆直径相对值,X表示竹秆某高度与竹秆全高的比值。回归模型的估计标准误为0.017,F=46