8年生火炬松施肥效应之数量化分析

2021-01-11张孝谊

绿色科技 2020年23期

张孝谊

(福建省屏南县林业局，福建屏南 352300)

1 引言

火炬松是我国重要的外来引种树种，也是闽东地区重要造林树种之一。屏南县从 20世纪70年代开始小规模引种栽培，目前全县火炬松面积约为0.1019万hm2。施肥是促进林木生长的重要措施之一。开展施肥试验有利于探索、提高施肥技术，最大限度地发挥施肥作用。有关林木施肥试验的报道文献繁多[1～4]，但以往的施肥试验往往忽略树体本身大小这一重要因素，其所得结论值得商榷。为此，笔者尝试运用数量化理论对所开展的火炬松施肥试验进行分析，以期能够得出一个相对客观的试验结论。

2 试验地概括

2016年在屏南县寿山乡太保村选择林相较整齐、林分长势一般、立地条件较一致、造林年度相同的火炬松纯林作为试验林。试验林均为2009年春季造林，初植密度2000株/hm2。现林下植被主要有杂竹、芒萁骨、黄瑞木、檵木等。试验地海拔485～515 m，坡度10°～13°，坡向NE，坡位中上部，土壤类型为山地红壤，土层较厚。属亚热带海洋性季风气候，年平均气温13～19 ℃，年平均降水量1500～2300 mm，降水日数220～220 d，年平均日照时数1500～1700 h,无霜期日数210～280 d。

3 试验设计和方法

试验设计：因素1-施尿素(含有效N46%)，2个水平：A-100g/株、B-200g/株；因素2-施复合肥，3个水平：A-0g/株、B-100g/株、C-200g/株。完全随机区组设计，3次重复，共18个试验小区。每个试验小区面积20×20 m。因素3-树体大小(树高)。

施肥采用沟施，即在每株树的上部，树冠投影外侧挖一弧形沟，宽15 cm、深20 cm、长60 cm，施肥后即覆土，施肥时间2016年5月12～15日。其他经营措施按常规方法进行。2019年1月每木调查各处理小区的树高、胸径、抽高等指标。

4 结果与分析

在数理统计过程中，变量一般分为2类：第一类是定量变量(连续变量)，如温度、树高等，它可以在某个区间内连续取值；第二类是等级变量(定性变量)，如树种、地类等，它的取值只能是几个等级(又叫水平)值中某一个。根据专业原理，有时候定量因子可以变成定性因子，有时候定性因子可以变成定量因子。当因变量与自变量不是线性关系，而使用线性关系的模型来分析他们之间的关系时，使用分级的方法，有利于减少非线性的影响。有时为了简化问题，把定性数据的等级对应一个数值，看成定量数据，数据处理时往往把它当做连续变量。

数量化方法是分析在数据中含有定性因子的一类统计方法。本文采用的数量化方法Ⅰ是分析在自变量中包含定性因子的一元线性模型[5]。其原始数据中的定性因子采用等级计法，如施尿素量、施复合肥量的多少；定量因子采用数字型，如树高、抽高。根据试验调查汇总数据(表1)，经唐守正院士等学者编制ForStat2.0-统计之林软件计算，结果见表2、表3。

表1 8年生火炬松施肥试验调查

表2 8年生火炬松施肥试验数量化分析-方差分析

表3 8年生火炬松施肥试验数量化分析-因子得分

表2中Pr>T表示服从F分布的随机变量§≥F的概率，即§≥F越小，说明该因子影响越显著。从表2可以看出：施尿素、施复合肥以及树体大小对施肥效应均存在极显著的影响，其中影响大小排序为施复合肥>树体大小>施尿素。表3中Pr>T值也是数值越小，该因子影响越显著。同样从表3可以看出Pr>T数值为施尿素>树体大小>施复合肥，故其影响显著性大小排序为施复合肥>树体大小>施尿素。需要指出的是，本文中树体大小的区分以树高做为指标。

根据笔者早年参与的全县林业土壤普查报告，屏南县林业土壤总体有机质含量较高，全氮较为丰富，全磷、全钾含量中等；速效氮含量中等，速效磷和速效钾含量普遍较低；微量元素中硼元素明显缺乏。从本试验结果来看，由于试验地位于山坡中上部，虽然土层较厚，但是碎石较多，土壤蓄水保肥能力较差弱，立地条件较差。复合肥氮磷钾配比较为全面，能针对屏南土壤缺磷、缺钾的状况进行有效补充，对于林木的生长具有明显的促进作用。

林木个体大小与其根系呈正相关，植株越大，根系越发达，吸收土壤水分和养分的能力就越强，生长约旺盛；反之，植株吸收土壤水分和养分的能力就越弱，生长就缓慢。本试验林木个体大大小与施肥效应存在明显的正相关就是对此的一个最好的例证。