森林是陆地生态系统中重要的养分储存库之一，有关森林碳(C)、氮(N)、磷(P)储量和潜力的研究在全球和区域尺度上都十分广泛^[1-4]。地处中国亚热带地区的森林生态系统具有极高的碳密度和养分含量，准确评估其碳、氮、磷储量具有重要现实和理论意义^[5-6]。植被养分储量评估的准确程度与生物量和养分含量密切相关。以往研究所涉及的平均生物量密度法、平均生物量转换因子法、改善模型连续生物量转换因子法和异速生长方程模型等^[7-9]，都是通过提高生物量估算精度以提升养分储量评估准确性；针对树木养分含量计算误差和不确定性的研究相对较少^[10-11]。虽有部分研究关注了不同树种、部位、树龄等因素间养分含量的差异^[12-14]，但受取样方法所限，对树干养分含量的研究仍较匮乏。此外，林学和生态学领域通常采用的小样本取样法(取样量n=3~5)难以客观代表总体养分含量，增加样本数可在何种程度上提高估算结果的精度值得开展进一步研究。因此，本研究以浙江省五洩国家森林公园中5种具有代表性的亚热带常绿阔叶树种为研究对象，采用有放回抽样方法(Bootstrap)模拟并比较小样本取样(n=3~5)和完整样本取样(n=18~32)的估算结果，评估小样本取样产生的养分质量分数估算误差。权衡取样量和估算精度，分别提出估算碳、氮、磷储量的最优取样量，为相关野外研究工作提供取样方法的建议。

1.   研究方法

1.1   样地设置与采样

研究样地位于浙江省绍兴市五洩国家森林公园(29°43′19.08″N，120°03′06.09″E)，在森林公园中部无人为干扰的典型亚热带常绿阔叶林区域，设3个50 m × 50 m的样地作为长期固定监测样地。样地内优势树种为马尾松Pinus massoniana、杉木Cunninghamia lanceolata、木荷Schima superba、青冈Quercus glauca和毛竹Phyllostachys edulis。2018年11月于样地内随机选取一定数量的5种优势树种，分别截取1 m长度的树干圆盘，具体取样数为马尾松32、杉木21、木荷30、青冈30和毛竹18。将取样圆盘于45 ℃烘干至恒量，选取部分组织粉碎待用。采用元素分析仪测定碳和氮质量分数(Elementar Analysen system, Bran & Luebb GmbH，德国)，采用Auto Analyzer测定磷质量分数(AA3, Bran+Luebb GmbH，德国)。

1.2   数据分析

采用单因素方差分析检测树种间和树种内部不同圆盘样本之间的碳、氮、磷质量分数差异。在不同树种的碳、氮、磷样本中(n=18~32)，采用Bootstrap估算随取样量增加各树种养分估算值总平方误差变化情况。预计估算值的总平方误差将随取样量增加而减小，取样量较小时(n=3~5)估算误差最大，可能检测到一个处于中间区间建议取样量(n_i)。利用随机抽样方式，比较常用的小样本(n=3~5)、n_i与全样本导致的碳、氮、磷质量分数变异情况。所有统计和绘图工作利用R软件完成，其中利用“segmented”软件包进行阈值点的估算。

2.   结果与分析

2.1   不同树种树干碳、氮、磷质量分数变异

由图1和表1可知：不同树种树干碳、氮、磷质量分数之间具有极显著差异(P＜0.001)；杉木碳质量分数极显著高于其他4个树种(P＜0.001)，木荷和青冈最低；马尾松氮质量分数最高(P＜0.001)，青冈和毛竹次之，杉木和木荷最低，但青冈和毛竹、杉木和木荷之间无显著差异；马尾松和毛竹的磷质量分数极显著高于其他树种(P＜0.001)，杉木和木荷较低，两者之间差异不显著。

图 1  亚热带常绿阔叶林5个优势树种的树干木材碳、氮、磷质量分数差异

Figure 1  Carbon (C), Nitrogen (N), and Phosphorus (P) concentration variation in the tree boles of five tree species in subtropical evergreen broadleaf forest

下载: 全尺寸图片 幻灯片

由表1可知：马尾松树干碳变异系数为7%，其他树种为2%~3%；马尾松氮变异系数最低，为12%，杉木最高，为31%，其他树种为20%~25%；青冈的磷变异系数最低，仅为14%，木荷最高，为46%，其他树种为20%~33%。总体上，树种之间的碳、氮、磷变异系数差异较大。马尾松碳变异幅度高于其他树种，氮变异系数低于其他树种；青冈磷变异系数较小，木荷较高，其他树种处于中间水平。

表 1  5个优势树种的树干碳、氮、磷质量分数和变异系数

Table 1  Mean (standard deviation) and CV of C, N and P concentrations for 5 dominant tree species

树种	样本数/个	碳			氮			磷
		质量分数/(g·kg−1)	变异系数/%		质量分数/(g·kg−1)	变异系数/%		质量分数/(g·kg−1)	变异系数/%
马尾松	32	449.1±29.4 b	7		3.4±0.4 a	12		1.10±0.29 a	26
杉木	21	487.6±10.7 a	2		1.3±0.4 c	31		5.80±0.19 c	33
木荷	30	431.1±12.5 c	3		1.2±0.3 c	25		0.52±0.24 c	46
青冈	30	435.9±11.5 c	3		2.0±0.4 b	20		0.98±0.14 b	14
毛竹	18	438.6±14.1 bc	3		2.2±0.5 b	23		1.28±0.26 a	20
说明：数值为平均值±标准差。同列不同小写字母表示树种间差异极显著(P＜0.001)

下载: 导出CSV 

| 显示表格

2.2   不同取样量下的树干碳、氮、磷质量分数估算

由图2可知：5个亚热带常绿阔叶林优势树种碳的建议取样量为5~10。其中马尾松建议取样量最大，为10，木荷为8，杉木和青冈为6，毛竹为5。亚热带常绿阔叶林优势树种氮的建议取样量为4~7。其中杉木建议取样数最大，为7，马尾松为6，木荷和青冈为5，毛竹最少，为4。亚热带常绿阔叶林优势树种磷的建议取样量为5~9。其中马尾松建议取样量最大，为9，青冈为8，杉木和木荷为7，毛竹最小，为5。

图 2  不同取样量下树干碳、氮、磷质量分数估算的总平方误差变化趋势

Figure 2  Trends of total squared error in the estimation of C, N, P concentration of tree boles with different sample sizes

下载: 全尺寸图片 幻灯片

2.3   比较小样本和建议取样量抽样引起的估算误差

将该领域常用的小取样量(n=3~5)与上述估算的取样量进行比较(表2)发现：对于碳估算变化，马尾松由于树干碳质量分数变异较大，n=5小样本条件下质量分数估算误差为6%~7%，采样样本增加至n=10后估算误差为4%~5%；而对于其他树种，n=3即可保证5%的估算误差。对于氮质量分数变化，马尾松n=5小样本条件下氮质量分数为14%~15%，采样样本增加至n=10不能显著提高估算精度；其他树种也表现出一致规律，增加取样量不能有效降低估算误差。磷质量分数变异较大，增加样本没有显著降低取样误差，n=5小样本可以保证马尾松、杉木、木荷估算误差为30%~50%，而青冈和毛竹估算误差为15%~20%。

表 2  5个优势树种建议取样量和小样本取样量(n=3~5)对总体的估算误差比较

Table 2  Comparison of the estimated extent of the population for different tree species suggested for n_i and small sample size (n =3−5)

养分	树种	总样本/建议取样量	建议取样量估算误差/%	估算误差/%
				n=3	n=4	n=5
碳	马尾松	32/10	−3.99~5.12	−8.41~11.75	−7.47~9.46	−6.14~7.23
	杉木	21/6	−1.90~2.00	−3.36~2.92	−2.81~2.58	−2.09~1.87
	木荷	30/8	−2.95~2.22	−4.70~5.01	−3.98~4.40	−3.13~3.62
	青冈	30/6	−2.89~2.78	−4.32~3.97	−3.20~3.02	−3.14~2.88
	毛竹	18/5	−3.37~3.30	−5.15~4.48	−3.37~3.76	−3.37~3.30
氮	马尾松	32/6	−12.42~17.11	−19.30~23.01	−16.60~18.83	−14.97~13.96
	杉木	21/7	−26.46~27.03	−41.27~51.44	−31.67~44.95	−29.82~37.00
	木荷	30/5	−22.39~33.72	−31.16~43.34	−20.94~40.85	−22.39~33.72
	青冈	30/5	−21.63~20.36	−31.89~20.93	−23.86~20.70	−21.63~20.36
	毛竹	18/4	−26.04 ~25.65	−30.48~30.12	−26.04~25.65	−19.74~19.71
磷	马尾松	32/9	−21.94~22.49	−41.23~44.80	−31.65~30.50	−28.66~37.46
	杉木	21/7	−30.70~26.10	−42.85~46.81	−40.42~49.80	−38.35~43.98
	木荷	30/7	−41.40~42.89	−58.95~73.93	−54.89~66.40	−42.95~51.08
	青冈	30/8	−11.79~13.29	−18.32~20.15	−18.76~17.18	−15.07~17.55
	毛竹	18/5	−18.06~22.52	−24.53~26.72	−22.01~24.37	−18.06~22.52

下载: 导出CSV 

| 显示表格

3.   讨论

相比于其他树种，马尾松的碳质量分数变异系数最大，建议采样样本数为10。马尾松是广泛分布于亚热带地区材用和造林树种，有关其碳含量的报道较多^[15-18]。常绿阔叶林中，马尾松碳质量分数为41.78%~51.49%，平均为46.38%，变异系数约10%^[18]；而全国范围内，马尾松树干中碳平均质量分数为44.83%，变异系数约12%^[11]；全球尺度范围内，马尾松及同纬度地区陆生植物碳质量分数为46.4%，变异幅度约14%^[19]。而在本研究中，样地尺度上马尾松树干碳质量分数为425.4~475.8 g·kg⁻¹，平均值为450.0 g·kg⁻¹，整体变异系数为7%。总体而言，碳质量分数变异幅度随尺度变大而升高。这些变异幅度指标可以为不同尺度水平下总碳质量分数估算的不确定性提供借鉴，即在全国范围内，关于马尾松的碳储量估算可能存在7%~12%的不确定性，而其他树种n=3~5小样本的估算误差在5%以下。李永进等^[20]研究杉木等8种乡土树种碳质量分数时发现：树干、树枝、树叶和树根碳质量分数分别为454.7~500.0、426.0~474.0、448.4~495.8、436.0~490.4 g·kg⁻¹，变异系数均在8%以上。唐学君等^[21]研究内蒙古地区油松林碳储量变异情况时发现：油松人工林平均碳质量分数为502.49 g·kg⁻¹，最大值为枝(514.20 g·kg⁻¹)，最小值为树干(485.58 g·kg⁻¹)，变异系数高达20%。赵瑞等^[22]研究了东北常见的10个阔叶树种的碳质量分数，每个树种采样3个，各树种碳变异系数均约7%。以上研究结果与本研究结果相似。

蔡国俊等^[23]分析了黔南喀斯特峰丛洼地建群树种不同器官氮和磷化学计量特征，发现取样量为3时，氮和磷平均质量分数分别为(8.90±7.68)和(0.52±0.32) mg·g⁻¹，变异系数分别为86.29%和61.54%，氮和磷质量分数变异较大。周永姣等^[24]研究了亚热带59个常绿与落叶树种细根养分及化学计量特征，发现取样量为3时，根氮质量分数变异系数为42.14%~48.87%；而将采样样本数提高到4或5时，亚热带常绿阔叶林11个树种细根氮质量分数的总体变异系数可降至41.29%^[25]。本研究也发现将采样强度由3提升至4或5时，树干氮质量分数估算误差显著降低。

傅梓奇^[26]测定了亚热带地区18个树种的根系磷质量分数(取样量为4或5)，总体磷质量分数为0.3~0.8 g·kg⁻¹ ，与本研究树干木质部的研究结果一致。刘璐等^[27]对神农架常绿落叶阔叶混交林植物不同器官磷质量分数的研究发现：加权不同器官生物量(取样量为3)可知磷在不同整株植物中的总体变异系数为47.67%，这与本研究14%~46%的变异系数相吻合。另外，本研究还发现增加采样样本不能显著降低5个树种磷的估计误差。

4.   结论

总体上，马尾松树干碳质量分数变异较大，在估算和评估马尾松林碳储量时，建议适当增加取样量(n=5~10)以降低估算误差，或在仍采用n=3~5小样本取样量的同时，将±10%估算误差范围考虑在内。对于其他树种的碳质量分数及所有树种的氮和磷质量分数，增加采样样本对估算误差的降低效果不明显，建议仍使用取样量n=4或n=5。