Prediction of Carbon, Nitrogen, and Phosphorus Contents of Chinese Fir Based on Ecological Factors and Artificial Neural Networks

研究区位于我国亚热带地区12个杉木主产省份。地理坐标为103.30°～119.32° E，22.25°～33.55° N，海拔高度60～1150 m。该区为亚热带季风气候，年平均气温为10.7～ 22.4℃，年平均降水量为773～782 mm。研究区土壤类型变化很大，区域分布可大致划分为南部红壤和砖红壤，中部黄壤和黄棕壤，以及北部褐土和黄褐土。林下植被较为稀疏，主要包括蕨类、草本类和菝葜类（Smilax spp.）等。

2018年7月底、8月初，在杉木主产区调查19个平均林龄为16 a的杉木中龄林。为减轻南北方物候差异造成的影响，野外采样从南到北进行。所选林分的经营管理方式相似，且较少受外界干扰。在每个人工林内建立3个20 m × 20 m的样地作为重复。相邻样地间垂直距离超过20 m。

在每个样地内，随机、均匀地选取5～6株长势良好的植株进行采样。在向阳面中部树冠随机摘取1年生和2年生成熟叶片，随后将所摘取的叶片混合成1个样本。去除下层植物和表面凋落物后，运用五点采样法采集0～20 cm深度的表层土壤。

此外，记录每个采样点地理坐标和海拔高度数据。从世界气候网站（www.worldclim.org）获得分辨率为1 km × 1 km的年平均气温（mean annual temperature, MAT）和年平均降水量（mean annual precipitation, MAP）数据。

植物样品经120℃杀青，70℃烘干至恒质量后，用球磨仪（Retsch MM400，德国）对植物样品进行粉碎，保存备测。叶片有机碳（C）含量采用重铬酸钾法测定；叶片氮（N）、磷（P）含量采用硫酸-高氯酸消煮法测定。土壤样品经风干后过2 mm筛，备测。土壤有机碳（SOC）含量采用重铬酸钾氧化-外加热法测定；土壤全氮（SN）含量采用凯氏定氮法测定；土壤全磷（SP）含量碱熔法测定；土壤pH采用电导法测定。

实地野外采样获得57组杉木叶片C、N、P含量数据，57组土壤SOC、SN、SP含量和pH值数据，以及19组经度、纬度（LAT）、海拔高度（ALT）、MAT和MAP数据（表1）。此外，运用文献收集法，从已发表文献中得到21组生态因子（土壤和地理气候）数据（数据表略）。

在Matlab环境下构建我国亚热带地区杉木叶片C、N、P含量预测的RBF神经网络模型。RBF神经网络模型构建按以下步骤进行：

（1）设置RBF输入和输出向量：选择LAT、ALT、MAT、MAP、土壤pH、SOC、SN和SP作为输入层。以杉木叶片C、N、P含量为输出层。

（2）数据归一化：将所有测量数据进行归一化处理，以消除量纲、指标数量级差异过大带来的影响。利用Matlab2012软件中自带函数mapminmax ( )对数输入值进行归一化处理，进而对RBF神经网络输出的模拟值进行反归一化处理。公式如下：

其中，I_ij是样本x_ij的归一值；x代表训练样本；i是输入指标数目，n=8；j是训练样本数目，n=50；x_min和x_max分别是输入因子的第i个指标的样本最小值和最大值；本研究中，α和β为参数，分别被设置为0.90和0.05。所有的归一化值都在[0,1]之间。

利用Matlab2012中的newrb函数建立3个3层RBF神经网络。通过迭代的方法设计隐含层，每迭代1次就增加1个隐含层神经元，迭代一致持续到误差低于目标误差或神经元数目达到最大值为止。newrb函数格式为：

其中，net为网络对象；P为输入向量；T为输出向量；GOAL是均方误差；SPREAD为径向基函数的分布函数，需要反复尝试以确定最佳值；MN为神经元的最大数目；DF为两次迭代之间添加的神经元数目。

训练样本包括来自我国亚热带杉木人工林的50组生态因子数据和叶片相关数据，使用newrb函数进行网络训练，经反复训练以获得以下参数：

（1）预测杉木叶片C含量的RBF神经网络各参数的值分别为：GOAL设为10⁻⁸；SPREAD取0.027；MN为100；DF为1。训练样本为1～57号中除了1、9、25、28、32、47、56号的50组样本。

（2）预测杉木叶片N含量的RBF神经网络各参数的值分别为：GOAL设为10⁻⁸；SPREAD取0.0152；MN为100；DF为1。训练样本为1～57号中除了6、7、10、30、32、41、45号的50组样本。

（3）预测杉木叶片P含量的RBF神经网络各参数的值分别为：GOAL设为10⁻⁸；SPREAD取0.0088；MN为100；DF为1。训练样本为1～57号中除了11、14、28、41、44、51、54号的50组样本。

由此构成了3个3层RBF神经网络模型用于分别预测杉木叶片C、N、P含量，其中输入神经元8个，隐藏神经元49个，输出神经元1个（图1）。

Figure 1.  RBF artificial neural networks structure for predicting leaf carbon, nitrogen, and phosphorus contents of Chinese fir

杉木叶片C、N、P含量分布用7组检验样本对训练后的RBF神经网络模型进行检验。7个检验样本在RBF神经网络中的预测结果和实测值的对比见表2。

可以看出，应用RBF人工神经网络模拟杉木叶片C、N、P含量结果与实测值拟合较好。其中杉木叶片C含量的预测值最小误差为0.40%，最大误差为4.32%，平均相对误差为1.82%，相关系数较低，为0.66；杉木叶片N含量的预测值误差范围为3.89%～15.08%，平均相对误差为9.88%，相关系数高，为0.85；杉木叶片P含量的预测值误差范围为2.43%～12.38%，平均相对误差为7.02%，相关系数较高，达到0.99。以上结果表明，RBF人工神经网络预测模型较合理，精度较高，预测值具有较高的可信度。

应用上述经训练的RBF神经网络模型，通过文献收集法获得的21组（编号58～78）生态因子数据对杉木叶片C、N、P含量进行预测。21组叶片C含量预测结果在356.40～532.17 mg·g⁻¹，平均值为453.05 mg·g⁻¹；叶片N含量预测结果在6.00～12.53 mg·g⁻¹，平均值为10.05 mg·g⁻¹；叶片P含量预测结果在0.94～1.48 mg·g⁻¹，平均值为1.27 mg·g⁻¹（表3）。

杉木叶片C、N、P含量特征见表4。杉木叶片C含量的变化范围为356.40～532.17 mg·g⁻¹，算术平均值为476.68 mg·g⁻¹，高于我国森林植物叶片C平均含量（455.10 mg·g⁻¹）^[16]，标准误为3.70 mg·g⁻¹，变异系数为6.86%，变异水平极低。杉木叶片N含量的变化范围为6.00～19.80 mg·g⁻¹，算术平均值为12.27 mg·g⁻¹（几何平均值为12.02 mg·g⁻¹），远低于全球1280种、我国753种陆地植物叶片N含量几何平均值（18.34 mg·g⁻¹、18.60 mg·g⁻¹）^{[4, 17]}，标准误为0.27 mg·g⁻¹，变异系数为19.52%，属于弱变异。杉木叶片P含量的变化范围为0.73～2.96 mg·g⁻¹，算术平均值为1.24 mg·g⁻¹（几何平均值为1.19 mg·g⁻¹），与我国753种陆地植物叶片P含量几何平均值（1.21 mg·g⁻¹）^[4]相当，略低于全球1280种陆地植物叶片P含量几何平均值（1.42 mg·g⁻¹）^[17]，标准误为0.04 mg·g⁻¹，变异系数为31.48%，属于中等变异。杉木叶片N含量空间变异性低于P含量，与Han等^[4]对我国753种陆地植物研究结果相一致，且两者又明显高于叶片C含量的空间变异性。

杉木叶片生态化学计量比结果见表5。杉木叶片C/N均值为40.28，明显高于全球水平^[18]，略高于三峡库区库首常绿乔木^[19]，但低于浙江天童山常绿针叶树C/N值^[20]，分布范围为24.75～49.22，标准误为0.94，变异系数为20.62%，属中等变异；叶片C/P均值为412.01，远高于全球平均水平^[18]，但明显低于三峡库区库首常绿乔木和浙江天童山常绿针叶树C/P值^[19-20]，分布范围为169.59～706.85，标准误为11.54，变异系数为24.74，属中等变异；叶片N/P均值为10.50，低于中国、全球及三峡库区库首常绿乔木和浙江天童山常绿针叶树N/P值^{[4, 18-20]}，分布范围为4.46～15.62，标准误为0.32，变异系数为26.64，属中等变异。

植物C、N、P生态化学计量特征能够反映植物养分利用效率以及养分限制类型^[21-22]。叶片C/N和C/P常被用于反映植物吸收营养元素时所能同化碳的能力，是体现植物营养元素利用效率的重要指标。杉木叶片C/N和C/P值均偏高，说明杉木对N、P的利用效率均处于较高水平，这很可能是杉木对瘠薄生境的一种生态适应策略。植物N/P值可以用作植物生长所受养分限制类型的判断，有研究表明，植物叶片N/P<14.00时，植物生长受N限制，叶片14<N/P<16.00时，同时受N、P限制或不受限制，叶片N/P>16.00时，受P限制^[23]。本研究中杉木叶片N/P值仅为10.50，显著小于14.00（One-sample T test，p<0.001），表明我国亚热带地区杉木生长主要受N元素的限制。

荟萃分析（Meta-analysis）是目前在全球水平或大尺度上对植物与生态因子间的生态化学计量学研究中最常用的方法^{[4, 18, 24]}。近年来，实地野外取样的方法逐渐被广泛运用，克服了荟萃分析存在的样点分布不均，采样时间不统一等诸多问题^[25]。然而，上述两种方法都是基于传统的破坏性取样手段，包括叶片、枝条、茎干、根系甚至全株取样等，对植物生长、群落结构稳定以及生态系统生产力等都会造成极大的负面影响。因此，从生态化学计量学角度出发，选取典型数据建立合适的植物化学性质预测模型，对于揭示植物性状对外界环境变化的响应与适应性具有重要意义。

在以往的研究中，估算植物营养元素含量通常是通过运用传统的回归或支持向量机回归（SVMR）方法对各类光谱技术所获得数据进行分析实现的^{[15, 26]}。因此，在掌握土壤基本理化性质本底值的前提下，在大的区域尺度上对植物化学性状进行预测成为了可能。Li等^[27]利用植物-土壤系统所具有的养分循环特征，构建了基于土壤性质预测养分含量水平的数学模型。本研究使用的RBF神经网络具有操作简单、训练速度快、精度高等优点，且能够准确地分析复杂系统的非线性关系。这就为预测杉木叶片C、N、P含量提供了一种经济、快捷的手段。本研究除将土壤性质运用到RBF神经网络模型的构建外，还将与植物营养元素含量密切相关的地理、气候等生态因子一并纳入，研究结果具有较强普适性，且准确性较高，为杉木人工林可持续经营、适宜种植区划分以及大径材培育等科学决策提供了一种新的方法和手段。

通常情况下，N素被认为是温带地区森林生态系统的限制因子，P素则是热带、亚热带地区森林生态系统的限制因子^[28-30]。在全球尺度的研究中，Reich等^[17]认为随着纬度增加（或温带降低），叶片N、P含量增加，而叶片N/P减小。Han等^[4]对我国753种植物叶片N、P含量的地理分布格局进行研究，同样发现叶片N、P含量随纬度的增加而增加，而叶片N/P并未表现出显著的变化趋势。一般认为，叶片和土壤P素在生态系统各尺度上均具有较强的耦合关系^{[24, 30-31]}。土壤P含量相对匮乏可能是导致我国植物叶片N/P高于全球水平的重要原因^{[4, 17]}。因此，相较于全球植物而言，我国植物生长更大程度上受到P的限制，尤其是热带和亚热带地区。

杉木是我国亚热带地区广泛栽培的人工林树种，目前对其生长所受的限制元素研究已较多^{[10, 32-33]}。然而，由于研究区地理位置、生长阶段、取样时间和经营管理模式等存在差异，并未有一致的结论。此外，尽管判断林木生长限制元素的方法众多^[34]，但由于野外操作难度大等诸多原因限制，导致这些方法的应用并不普遍。因此，将更多的养分循环利用机理和树种特性等与叶片N/P值结合起来对杉木生长限制元素进行综合判断是很有必要的。

目前，针对陆地植物而言，判断限制元素的N/P阈值有很多^[23]。研究者普遍认为N/P值越低，则植物生长越可能受N素限制。本研究中，叶片N/P值处于较低水平，因此杉木生长受到N素限制的可能性大。元素变异水平常被用于植物限制元素的判断，Han等^[35]认为元素变异水平越低，则植物生长受其限制的可能性越大。本研究中，叶片N含量变异系数高于P含量，同样表明杉木生长受N素限制的可能性大。植物养分重吸收是指植物将营养元素从凋落叶片中转移到其他它组织的过程，被认为是植物的一种重要的营养保存机制^[35-37]。大多数研究发现，杉木P重吸收效率均不同程度的高于N重吸收效率^{[7, 38]}。相对重吸收率假说认为植物倾向于吸收更受限制的元素^[39]，但Mariana等^[40]则认为重吸收效率高是对相应土壤养分匮乏的一种适应策略。相比之下，作者更倾向于后者，高的P重吸收效率使得植物-土壤系统P循环处于相对健康的状态，而低的N重吸收效率并未改善杉木受土壤N素匮乏的现状。此外，郑璐嘉等^[41-42]通过对不同林龄杉木叶片、细根氮稳定同位素（δ¹⁵N）组成进行测定，同样得出N素是处于速生期杉木人工林生长限制因素的结论。

RBF神经网络是预测杉木叶片C、N、P含量的一种新方法，具有合理的模型结构和较高的可信度。本研究可为利用机器学习方法预测除C、N、P元素外的其它植物养分含量提供有价值的参考。杉木叶片N含量和N/P均偏低，表明在我国亚热带地区土壤普遍缺磷的现状下，杉木生长仍然受到N素的限制。然而，N素限制还需要从杉木人工林养分循环利用以及生理生态学特性等方面进一步确认。