研究区域位于江西省中部吉安市(113°48′～115°56′ E，25°58′～27°57′ N)，地貌以丘陵、山地为主，土壤类型主要有红壤、黄壤等，母岩以砂岩、花岗岩为主；属亚热带季风湿润气候，年平均气温17.1～18.6℃，年均降水量1 360～1 577 mm，降水主要集中在春季和初夏，年无霜期281 d。森林资源丰富，主要植被类型有针叶林、常绿阔叶林、针阔混交林和竹林等，森林覆盖率达到67.6%。灌木以木(Loropetalum chinense (R. Br.) Oliv)、杜茎山(Maesa japonica (Thunb.) Moritzi.)为主，草本以铁芒萁(Dicranopteris linearis (Burm.f.) Underw.)、庐山楼梯草(Elatostema stewardii Merr)为主。

根据森林资源二类调查资料及当地林业部门的了解，在研究区的安福县、遂川县、井冈山市的闽楠天然次生林分布地进行踏查，选择人为干扰程度轻且在分布地具有代表性的地块设置标准地，标准地面积依据其分布地形等因素而定，为400 m²(20 m×20 m)或600 m²(20 m×30 m)，共计标准地23块，并采用相邻网格调查方法，将每块标准地划分为4或6个10 m×10 m小样方，共有100个小样方。以每个小样方为调查单元，调查胸径≥2 cm的所有乔木的林木位置(x、y坐标)、树种、胸径、树高等；更新层的调查在10 m×10 m小样方内记录胸径 < 2 cm的幼苗(H < 1 m)、幼树(1 m < H < 3 m)^[20-21]的树种、地径、树高、冠幅等；林下植被在10 m×10 m小样方的上、中、下分别选择具有代表性的样方，灌木样方调查为2 m×2 m，草本样方调查为1 m×1 m，调查记录灌木和草本的种类、数量、高度、盖度等，同时在1 m×1 m小样方内记录凋落物层的盖度和厚度；土壤质地、土层厚度、腐殖质层厚度等均在做土壤剖面时记录；利用GPS测定标准地的地理位置和海拔，同时记录坡位、坡向、坡度和郁闭度等因子。

数量化模型Ⅰ的数学原理是建立某个因子和数量因子或非数量因子之间的线性关系，公式如下：

其中，Y_t为第t点因变量值；B₀为常数项；B(i, j)为第i定性因子第j等级的得分；s_t(i, j)为第t点第i定性因子第j等级的反应值(或0或1)；B(i)为第i定量因子的回归系数；x_ti为第t点第i项定量因子的观测值；n_i为第i定性因子的等级数。为衡量各项目在模型中的贡献，计算各项目的得分范围，得分范围越大，则该项目影响就越大。

重要值计算目前应用较多的是Curtis等修改的重要值指数，是相对密度或相对多度、相对频度和相对显著度3项指数的综合。采用下式计算闽楠幼树幼苗的重要值^[22]：

重要值＝相对密度(%)+相对频度(%)+相对盖度(%)

其中，相对密度=(某种植物的个体数/所有种的个体数)×100%；相对频度=(某种植物的频度/所有种的频度总和)×100%；相对盖度=(某种植物的盖度/所有种盖度总和)×100%。

重要值作为因变量，影响林下植被更新状况的环境因子为自变量，通过多元数量化模型Ⅰ建模，分析各环境因子对闽楠天然次生林闽楠幼树幼苗更新的影响，在此基础上，对影响极显著的定性因子进行控制变量分析，即控制其它影响因子相似的情况下，采用单因素方差分析(one-way ANOVA)方法进一步分析某个因子变化下闽楠幼树幼苗更新状况(平均密度、平均高度)的差异程度。所有建模过程、图表及数据分析过程在SPSS 19.0和Excel 2016中完成。

参考相关文献^[23-24]及结合闽楠天然次生林的特征，选择影响闽楠幼树幼苗自然更新的环境因子并对其进行分解。环境因子包括林分因子(平均胸径、郁闭度和株数密度)、地形因子(坡度、坡位和坡向)、林下植被因子(林下植被盖度、林下植被高度)、土壤因子(土壤厚度、腐殖质层厚度)、凋落物层厚度以及闽楠的状况(闽楠下种母树株数、闽楠密度和闽楠平均胸径)等14个因子。需要说明的是，由于标准地的海拔在200～400 m范围，跨度小，土壤类型均为红壤，其土壤质地等均相差不大，因此没有考虑海拔及其它土壤因子；此外，考虑到能结实的闽楠数量的多少可能对其幼树幼苗更新影响较大，因此，将结实量较为稳定的闽楠视为下种母树，参考文献[18]和闽楠生长特性^[25]，闽楠下种母树定义为胸径≥20 cm以上的树木。各环境因子的分解见表 1。

以闽楠幼树幼苗的重要值为因变量，选取的环境因子(包括定性和定量项目)为自变量，运用多元化数量模型Ⅰ，对14个环境因子进行建模，根据结果进行偏相关检验，把相关系数小的，差异性不显著的项目剔除掉，再对剩下的项目重新建模，反复进行。考虑到减少建模因子的数量，才能增加模型的准确性，最后筛选出闽楠下种母树株数、郁闭度、坡位、腐殖质层厚度、坡向、林下植被盖度、凋落物层厚度、株数密度8个环境因子，其模型如下：

Y=2.284-3.416x₁₁-2.042x₁₂+0.245x₁₃-1.709x₂₁-0.686x₂₂+0.423x₃₁+0.384x₃₂-0.171x₃₃-1.488x₄₁-1.388x₄₂+6.022x₅₁+1.341x₅₂-0.131x₆₁+0.023x₆₂+0.291x₆₃-0.152x₇₁-0.094x₇₂+0.104x₇₃+0.605x₈

在模型检验中，偏相关系数采用t检验，其结果均为极显著或显著；复相关系数用F检验，检验结果为极显著，所以筛选出的8个因子与闽楠天然次生林幼树幼苗自然更新具有显著或极显著的关系。模型运算结果及各类目得分值和项目分值极差见表 2。

从表 2可知，闽楠下种母树株数对闽楠幼树幼苗更新的贡献率最大。自然更新必须有种源，而闽楠下种母树结实比较稳定，可为更新提供种源，对其幼树幼苗的更新起最主要的作用，大树越多，提供的种子越多，越有利于自然更新。因此必需保证林分中有一定数量的闽楠下种母树，这是其自然更新的基础。通过对标准地调查数据的分析，并参考闽楠种子萌发^[26]及其幼树幼苗生长竞争特性^[25]，要保证其天然更新的顺利进行，闽楠下种母树应保留200株·hm^-2 以上。

郁闭度对闽楠幼树幼苗更新影响极显著，其影响程度排名第2位，其中以0.5～0.7的郁闭度最好，其次是0.3～0.5和>0.7， < 0.3最差。进一步通过单因素方差分析表明(表 3)，郁闭度0.5～0.7下闽楠幼树幼苗的平均密度显著高于其他郁闭度(P < 0.05)，而郁闭度为0.3～0.5的幼树幼苗平均高度显著高于郁闭度 < 0.3和>0.7的平均高度(P < 0.05)。林冠层郁闭度的大小直接决定林下层光照的情况。郁闭度过小，会导致光照强度增大，引起土壤湿度变小，不利于种子的萌发；郁闭度过大，会导致光照强度过弱，不利于幼树幼苗生长。

株数密度以2 000～3 000株·hm^-2最好，其次是3 000～4 000株·hm^-2和1 000～2 000株·hm^-2， < 1 000株·hm^-2更新最差。模型运算虽然得出这种结果，但由于株数密度的多少主要是影响光照条件，而在天然林中株数密度对林内光照的影响程度，还与株数按径级的分配比例有很大关系。因此，考虑到株数密度与郁闭度之间存在一定的关系，在经营过程中主要可采用调控郁闭度大小的措施，对郁闭度过大的结合径级结构的调整适当间伐掉一些次要树种，而对郁闭度过小的适当补植树种以提高郁闭度。

根据表 2结果，坡位对幼树幼苗更新特征的影响排名第3位，从得分值来看，下坡位闽楠幼树幼苗更新最好，其次是中坡，最差的是上坡。图 1单因素方差分析结果表明，下坡幼树幼苗的平均密度和平均高度与中坡、上坡之间均存在显著差异(P < 0.05)。这主要是因为下坡土壤的持水性能高，长期的地表径流作用会导致下坡土壤更为肥沃，下坡土壤更有利于乔木种群的繁殖、种子的萌发以及幼苗种群的扩展，且雨水冲刷导致种子冲到下坡，下坡种子库较上坡丰富，因此，下坡更新更好。而上坡位土壤流失水分多，光照比下坡强，同时养分含量具有抑制效应，不利于种子的萌发和幼树幼苗的生长。

图  1  不同坡位对幼树幼苗更新状况的影响

Figure 1.  Effect of different position of slope on natural regeneration

坡向以半阴坡最好，其次是阴坡、半阳坡，阳坡最差。通过单因素方差分析(图 2)可知，阳坡、半阳坡、半阴坡、阴坡闽楠幼树幼苗的平均密度之间的差异均显著(P < 0.05)，其中半阴坡闽楠幼树幼苗的平均密度最大，且其平均高度也显著高于其他坡向(P < 0.05)。不同坡向太阳的辐射强度以及日照时数不同，导致水、热和土壤理化性质有所差异。半阴坡和阴坡林分内太阳辐射及日照少导致温度降低和蒸腾作用减小，土壤水分充足，对种子萌发和幼树幼苗生长起促进作用。而阳坡光照强，导致温度升高和水分蒸腾增大，不利于幼树幼苗生长。闽楠是耐阴性树种，对光热的耐受条件有限。因此，对闽楠幼树幼苗来说，半阴坡最适合闽楠更新。

图  2  不同坡向对幼树幼苗更新状况的影响

Figure 2.  Effect of different aspect of slope on natural regeneration

腐殖质层厚度对闽楠次生林下幼树幼苗更新的影响排在第4位。其中，以>20 cm最好，10～20 cm次之， < 10 cm最差。从单因素方差分析(表 4)可以看出，腐殖质层厚度>20 cm的闽楠幼树幼苗的平均密度和平均高度均显著高于其他腐殖质层厚度的平均密度和平均高度(P < 0.05)。由于土壤中的水分、营养物质等影响着种子的萌发、幼苗的生长和发育，而幼树幼苗的根系较浅，加上腐殖质层的养分、微量元素较为丰富，因此腐殖质层厚度越厚越有利于幼树幼苗的更新。

凋落物层厚度以 < 1 cm最好，其次是1～3 cm，最差是>3 cm。单因素方差分析(表 5)可以看出，凋落物层厚度 < 1 cm闽楠幼树幼苗的平均密度和平均高度均显著高于其他凋落物层厚度的平均密度和平均高度(P < 0.05)。这是由于较厚的凋落物对种子和土壤接触影响较大，不利于种子到达土壤表面，减少了种子萌发的可能性和幼苗定居的机会；同时，凋落物也可能会通过化感作用等抑制幼苗的更新。因此，经营中可通过人为措施对凋落物层厚度进行控制，以促进其更新。

林下植被盖度以 < 30%最好，其次是30%～60%、60%～90%，最差是>90%，从单因素方差分析(图 3)可知，林下植被盖度小于 < 30%的幼树幼苗平均密度和平均高度均显著高于其他林下植被盖度的平均密度和平均高度(P < 0.05)。林下植被盖度的高低直接影响林下幼树幼苗生长空间的大小，其盖度越大林下幼树幼苗与灌草对光照、水分、养分等资源的竞争越强，从而影响幼树幼苗的自然更新。因此在闽楠天然次生林经营中，进行砍杂抚育将有利于其自然更新。

图  3  不同林下植被盖度对幼树幼苗更新状况影响

Figure 3.  Effect of different understory vegetation coverage on natural regeneration of sapling

在更新过程中，土壤种子库是天然更新的基础^[27]，母树的数量和母树的繁殖能力对幼树幼苗的更新有着重要的作用，吴大荣^[25]等研究表明：闽楠更新很大程度上依赖于种子产量，且大径级株数多的种群中其幼树幼苗种子库最丰富，本研究也得出闽楠下种母树株数(结实量较稳定的闽楠)对闽楠幼树幼苗更新的贡献率最大。因此，保证林分中有一定数量的闽楠下种母树对其自然更新至关重要。

林冠层可通过改变土壤和温度条件^[28]以及降低风速和减少蒸腾^[29]等为幼树幼苗提供适宜的环境，因此，郁闭度会影响幼树幼苗的更新和生长。本研究得出郁闭度为0.5～0.7时自然更新状况最好，而刘明国^[24]等对油松的研究认为0.3～0.5的郁闭度下油松更新更好，这可能是由于闽楠和油松有着不同的生物学特性所造成的。同样，林分密度直接导致林下光、热、水分、土壤理化特性等生态因子的异质性，并直接影响到林下幼树幼苗的生长，林分密度过低会抑制大量中生偏阴生性优势植物幼苗的萌发和生长，过大则林内光照强度减弱有助于幼苗的生长却增大了种内和种间的竞争^[30]，本研究表明密度在2 000～3 000株·hm^-2时，闽楠幼树幼苗更新的状况最好，而郑金萍^[5]等对蒙古栎的研究认为林分密度在400～800株·hm^-2更新数量最多，造成这种差异的原因可能是株数按径级的分配比例不一样，且闽楠属于阴性树种，而蒙古栎属于阳性树种，它们对光照的要求不同。此外，坡向也影响着光照条件^[31]，本研究得出半阴坡闽楠幼树幼苗更新状况更好，主要是由闽楠的生物学特性所决定的。

对于坡位及腐殖质层来说，由上坡到下坡，土壤由剥蚀逐渐过渡到堆积，腐殖质层由薄变厚，土壤质地由粗变细，土壤含水量增加，因此，下坡位、厚腐殖质更有助于更新。而对于凋落物层厚度，不少研究均表明其对更新有较大的影响，如张树梓等^[32]研究表明天然更新苗状况与凋落物厚度呈负相关，本研究结论与之相一致，这主要是因为凋落物持续积累增厚会对植被的自然更新形成阻碍，表现在物理的机械阻挡、化学的他感作用、生物的动物侵害和微生物致病作用等方面^[33]。此外，相关研究也表明，林下植被被认为是决定群落结构和动态的重要因素^[34]，但林下植被盖度过大，会加大林下更新树种与灌草对光照、水分、养分等资源的争夺^[35]，从而不利于自然更新，Pham Van Huong^[36]等研究结果及本研究结果均验证了这一点。

影响闽楠天然次生林自然更新的因子主要有密度(闽楠下种母树株数、郁闭度、株数密度)、坡位、腐殖质层厚度、坡向、林下植被盖度、凋落物层厚度。

密度对闽楠天然次生林自然更新影响最大。其中，闽楠下种母树株数起着最主要的作用，其株数多有利于自然更新，而郁闭度与株数密度过高或过低均不利于其自然更新，郁闭度在0.5～0.7之间其幼树幼苗的平均密度均显著高于郁闭度 < 0.3、0.3～0.5和>0.7的(P < 0.05)。凋落物层越厚、林下植被盖度越大均不利于其自然更新，凋落物层厚度 < 1 cm最好，最差是>3 cm，厚度 < 1 cm的幼树幼苗平均密度和平均高度均显著高于厚度1～3 cm及>3 cm的(P < 0.05)；林下植被盖度以 < 30%更新最好，其次是30%～60%、60%～90%，最差是>90%，其 < 30%盖度的幼树幼苗平均密度和平均高度均显著高于其他林下植被盖度(P < 0.05)。腐殖质层越厚越有利于闽楠的自然更新，腐殖质层厚度>20 cm的闽楠幼树幼苗的平均密度和平均高度均显著高于腐殖质层厚度 < 10 cm及10～20 cm的(P < 0.05)。地形对更新的影响主要体现在坡位和坡向，下坡位闽楠幼树幼苗更新最好，最差为上坡，下坡幼树幼苗的平均密度和平均高度均显著高于中坡和上坡(P < 0.05)；坡向则以半阴坡最好，阳坡最差，半阴坡幼树幼苗的平均密度和平均高度显著高于其他坡向(P < 0.05)。

为促进闽楠幼树幼苗的自然更新，林分中闽楠下种母树应保留200株·hm^-2以上，同时，对郁闭度、林下植被盖度及凋落物层厚度等实施相应的调控措施。