施肥是促进人工林速生丰产的主要技术措施^[1]。长期以来，我国在林地施肥方面存在有机肥利用不足而化肥施用过量等不合理现象^[2]。近年来，随着化肥过量施用导致土壤性状变差、肥力下降以及有机肥尤其是畜便量激增等问题越来越突出，一些学者提出应重新树立有机肥在农、林业生产中的主体地位^[3], 强调采用有机、无机肥相结合的施肥模式^[4]，在有机肥的基础上，配合施用化肥，同时实现作物(林木)增产、品质改善、环境友好等多重目的。

黄梁木Neolamarckia cadamba为茜草科团花属植物，常绿大乔木，生长迅速，材质良好，是热带亚热带地区一种极具发展前景的阔叶树种^[5]，被称为“奇迹树”。目前，黄梁木在广东云浮、新兴、信宜等地以及海南的造林面积很大^[6]，且近年来在广东省造林面积逐年扩大，逐渐成为退耕还林和天然林保护工程的主要造林树种之一^[7-8]，黄梁木作为速生树种，对水分、养分的需求较大，因此，为保证速生丰产，施肥在黄梁木的育苗生产以及造林上逐渐普及。研究发现，合理施肥能明显提高黄梁木的苗高、地径和生物量，且复合肥的施用效果更佳^[9-11]。但盲目、大量的施肥不仅造成资源浪费，增加经济成本，且对林木的生长有一定程度的毒害^[12]。

广东是畜禽生产大省，随着畜禽养殖业集约化、规模化发展，广东省每年的畜禽粪便排放量接近0.74亿t^[13]，而大部分直接排放进入环境^[14-15]，粪污无序排放不仅对环境有害，而且对人类和动物也会产生不利的影响^[16]。因此，找到合适可行的畜禽粪便处置和利用措施是当前亟待解决的问题。充分借鉴国内外一些成功经验可知，畜禽粪便资源化利用应成为畜禽养殖业污染物总量减排的首选，而最主要的方式是作为肥料还田还林。然而，由于饲料添加剂滥用造成许多畜禽粪便重金属含量超标，农用存在很大风险，所以林用逐渐成为首选。基于此，本研究在林地现场条件下，通过正交施肥试验，探究猪粪、鸡粪与无机肥配合施用对黄粱木幼林的生长效应，以期为黄梁木人工幼林营养管理及养殖场有机肥的林地利用提供科学依据。

1 材料与方法 1.1 试验地概况

试验地位于广东省雷州市白沙镇，北纬21°15′~21°20′，东经109°22′~110°27′。属台地地貌类型，海拔35~50 m。年均气温28 ℃，年均日照1 165 h，年均无霜期365 d，年均降水量1 787.5 mm，属热带海洋季风湿润气候区。试验地土壤为浅海沉积物发育的砖红壤，土层厚度在100 cm以上，土壤有机质19.45 g·kg^-1、全氮0.23 g·kg^-1、全磷0.53 g·kg^-1，全钾2.8 g·kg^-1。

1.2 试验材料

试验所用黄梁木苗木为组培苗，种植前长势一致，树高25 cm左右，地径5.6 mm左右。供试化肥过磷酸钙(P₂O₅质量分数为12%)、氯化钾(K₂O质量分数为60%)为市购化肥，有机肥用堆沤发酵腐熟的猪粪、鸡粪。猪粪中养分：w(N)为9.32 g·kg^-1、w(P₂O₅)为7.54 g·kg^-1、w(有机质)为463.1 g·kg^-1，鸡粪中养分：w(N)为18.13 g·kg^-1、w(P₂O₅)为12.41 g·kg^-1、w(有机质)为370.10 g·kg^-1，有机肥含水量600 g·kg^-1。

1.3 试验设计

采用正交试验设计。试验因子为猪粪、鸡粪、过磷酸钙、氯化钾，各因子均设4个水平，具体见表 1。采用L₁₆(4⁵)正交表^[17] (表 2)，其中，以不施肥的试验组1号作为对照(CK)。

试验布置：造林地上除草清杂、堆烧残留物，全面整地，按3 m×3 m株行距挖60 cm×60 cm×60 cm明穴。按照完全随机区组布置试验，25株小区(每个小区5行，每行5株)，造林前按表 2将肥料均匀施入各相应的植树穴，然后覆土，15 d后用相同规格的黄梁木苗木造林。每个小区最外围一圈为保护行，即每个小区有效调查株数为9株(图 1)，每个处理设置4次重复。2012年6月造林，当年12月调查黄梁木成活率，测定树高、地径。

2.4 数据处理方法

试验数据用Excel整理后，用SPSS17.0软件进行数据分析。采用单因素方差分析(One-way ANOVA)和Duncan’s多重比较法分析不同处理间差异，多因素方差分析进行最佳组合的筛选，分别采用Pearson相关系数和Linear Regression进行相关分析和回归分析。

3 结果与分析 3.1 肥料配施对黄梁木幼林生长的影响

造林半年后调查发现，黄粱木成活率为100%，其树高和地径生长情况见表 3。就黄梁木的树高生长而言，与对照(处理1) 相比，除处理2以外，其他不同比例有机、无机肥配施条件下黄梁木树高均显著高于对照，其中，最为显著的是处理13(猪粪14 kg·株^-1+鸡粪0 kg·株^-1+过磷酸钙1.5 kg·株^-1+氯化钾0.25 kg·株^-1)，施肥6个月后，其黄粱木树高为1.78 m，比对照高出109.41%，处理2(猪粪4 kg·株^-1+鸡粪2 kg·株^-1+过磷酸钙0.5 kg·株^-1+氯化钾0.25 kg·株^-1)的树高为1.05 m，虽比对照高出23.53%，但差异未达到统计学显著水平。

同样，不同比例的有机、无机肥配施对黄梁木的地径生长也具有显著的促进作用，施肥条件下的地径为27.17~49.04 mm，比对照(19.69 mm)增加了37.99%~149.06%，促进作用最弱的是处理2，而促进作用最为显著的仍然是处理13。除对照外，其他不同处理间的树高、地径因有机、无机肥比例不同而差异各异，仅从数值来看，处理7、10、13较为突出，这与该处理猪粪、过磷酸钙的施用量较大可能有关。但由于树木生长周期长，肥效反应短时间难以测定，此次结果仅为造林半年后的黄粱木幼树的生长状况，因此施肥是否对此后的黄粱木生长有持续的促进作用，还有待继续观测分析。

3.2 最佳施肥处理组合筛选

正交试验结果的直观分析虽反映了不同处理的差异性，但不能区分该结果的差异是由各因素的水平变化导致还是由试验的随机波动引起。因此，在直观分析的基础上，有必要对试验数据进行方差分析^[18]，从而明确各个因素对试验指标的影响程度，确定各因素的主次顺序，进而找出较好的试验条件或最优参数组合。本研究同时考察了黄梁木树高、地径2个试验指标，将分别以这2个指标作为分析单元，进行最佳施肥处理组合筛选。

表 4反映了黄梁木幼林的树高和地径方差分析情况，F检验结果显示，对黄梁木树高、地径生长而言，各因素施肥效应的主次顺序均为A＞C＞B＞D，即猪粪＞过磷酸钙＞鸡粪＞氯化钾，4个因素中，因素A与因素C对黄梁木的生长具有显著影响(P < 0.05)，而另外2个因素的主效应则不显著(P>0.05)。各因素水平间的多重比较(表 5)发现，因素A、C各水平间具有显著性差异，而因素B、D不同水平间虽有差异，但未达到统计学显著水平。因素A、B、C不同水平对黄梁木树高、地径生长影响的强弱顺序均为4>3>2>1，而因素D对黄梁木树高和地径生长最好的促进水平分别为水平2和3，因此，仅从方差检验结果来看，A₄B₄C₄D₂和A₄B₄C₄D₃分别是对黄梁木树高和地径生长最为有利的组合。但在施肥方案的选优问题上，最佳施肥量和经济效益也是重点考虑因素，即实际生产中一般以符合统计学要求、施肥量少(投资风险低)而经济效应大的作为最优施肥模型。因此，综合考虑经济成本和黄梁木生长状况，最终选择A₄B₁C₄D₁，即猪粪14 kg·株^-1+过磷酸钙1.5 kg·株^-1作为黄梁木幼林施肥方案。

3.3 有机、无机肥施用量与黄梁木生长的关系

对有机、无机肥施用量与黄梁木树高和地径间的关系进行统计分析，从简单相关分析结果看，各类肥料施用量与黄梁木生长呈正相关关系(表 6)，其中，猪粪、过磷酸钙施用量与黄梁木树高、地径相关性均达到极显著，这与前面方差分析结果一致。但从标准回归系数来看，对于黄梁木树高生长而言，猪粪、过磷酸钙贡献要大于其他两者，即回归系数较高，而在地径生长方面，猪粪的贡献略低于氯化钾。采用逐步回归方法对各类型肥料施用量与黄梁木树高和地径的生长关系进行了分析，得到黄梁木地径和树高与各类型肥料施用量间的逐步回归方程:树高(y_H)，y_H=1.016+0.031x₁+0.221x₂(R²=0.835，P=0.000)；地径(y_D)，y_D=26.193+0.991 x₁+6.052x₂(R²=0.832，P=0.000)。式中：x₁和x₂分别代表猪粪和过磷酸钙施用量。

4 讨论与结论 4.1 有机、无机肥配施在林业上的发展前景

我国有机肥资源种类丰富，数量庞大，但利用方式落后且利用率低，有近70%直接排放进入环境中，不仅造成资源浪费，而且会形成新的面源污染^[19-20]。我国农林业生产主要依靠化肥，其增产依赖率达40%^[21]。化肥作为速效肥料，长期大量施用会引起土壤性状变差及肥力下降^[22]。因此，探索利用当地丰富的有机物料资源，结合科学的堆肥技术，生产出优质的有机-无机复合肥产品，成为我国实现农林业可持续发展的重要措施。

有机肥营养均衡，养分全面，可提高土壤生物多样性和土壤生物学活性，改善土壤理化性质、土壤环境及营养循环，提高土壤肥力^[23-24]。近年来，有机肥在农田、菜地、果园等得到广泛应用并取得了一定的成果。然而，在畜牧业生产中，越来越多的金属元素添加剂在规模化禽畜养殖过程中被广泛使用，使得畜禽粪便中重金属大量增加，造成以畜禽粪便为原料的有机-无机复混肥农用后重金属累积，引起农田土壤重金属含量超标，进而通过植物和农产品最终危害动物和人体健康。

相对于有机肥农用而言，有机肥林地利用因具有远离人类食物链、面积广阔、环境容纳量大等多重优势也逐渐被人青睐^[25-26]。当前我国人工林分布广泛，巨大的人工林面积将为畜便利用提供广阔的空间；同时，由于传统粗放的营林方式，已经引起我国林地土壤结构恶化，肥力和生产力下降，亟需采取科学的方法对现有林进行合理经营，维持和提高土壤肥力。因此，从经济、环境、资源角度出发，将有机肥合理利用于林地，可达到畜便总量减排、林地土壤维持和林木增产的目的。

本研究因地制宜，将研究区附近养殖场猪粪、鸡粪等有机肥料与无机肥料配合施用于林地，结果发现，有机、无机肥配施可以显著促进黄粱木幼树的树高和地径生长。基于以往研究发现，有机、无机肥配合施用对于林木生长的促进作用并非偶然，具有广泛性，如有机肥的施用对茶树、橡胶树、杨树^[26-29]等果树或林木生长均具有显著的促进作用。因此，有机肥用于林业生产具有广阔前景。

4.2 黄梁木最佳施肥处理组合筛选

施肥作为提高林地生产力的有效管理手段之一，合理施肥能促进林木生长和提高作物产量，但过量施肥不仅抑制林木生长，还会对环境造成负面影响^[30]。因此，在保证高产的前提下减少生产投入和环境污染是当前农林经营管理所关注的热点问题^[24]。在众多施肥的方法中，设计较为复杂而高效的田间正交施肥试验，通过较少试验得出各因素的主次关系和最优搭配条件，是最常见的方法^[31]。但由于该方法是采取部分试验来代替全面试验的方法，由此可能会产生预想不到的误差^[32]。另外，在多指标正交试验中，不同试验因素及其水平对各试验指标的影响并非同步，某项指标的改善可能使另一项指标恶化^[33]。这可能是本研究中无论是试验还是SPSS理论分析上，未曾出现同时对黄梁木幼林树高、地径生长均有最佳施肥效应组合的一个重要原因。在施肥方案的选优问题上，要同时遵循符合统计学要求、施肥量少而经济效应大等原则。因此，综合考虑黄梁木生长状况(统计结果)和经济成本，最终选用A₄C₄作为黄梁木幼林施肥方案。对于是否还有更好的配方能在林木长势良好的同时继续降低肥料成本，还需进一步的优化试验。

4.3 有机、无机肥施用量与黄梁木生长的关系

肥料效应函数是目前被广泛应用的一类施肥模型^[34]，反映了肥料用量与作物产量之间的统计关系，其方法是通过肥料的单因素、二因素或多因素多水平田间回归试验设计，运用统计学的方法，得出单元肥料效应函数、二元或多元肥料效应函数回归方程，可以在一定程度上预测施肥后作物的生长状况^[35]。本研究分别以黄梁木的树高和地径为目标函数，选取猪粪(x₁)、过磷酸钙(x₂)的施肥水平为调控因子，通过逐步回归模拟黄梁木树高和地径与2个调控因子之间的回归关系，得到树高(y_H)、地径(y_D)与二者的回归方程, 分别为：y_H=1.016+0.031x₁+0.221x₂(R²=0.835)，y_D=26.193+0.991x₁+6.052x₂(R²=0.832)，由此可推测，当最终选用A₄C₄作为黄梁木幼林施肥方案时，造林半年后其树高和地径可达到1.78 m和49.15 mm，同时，由该回归方程中过磷酸钙和猪粪2个因子的系数大小可发现，前者对黄梁木树高和地径生长的贡献远大于后者，由此推断，磷肥可能是黄梁木生长需要重点控制的因子，鉴于本试验中过磷酸钙和猪粪两因素在大试验水平时正效应最大，今后可适当增加磷肥和猪粪的施用。