Effects of nitrogen and phosphorus additions on carbon storage and allocation of Cinnamomum camphora seedlings under different planting densities
-
摘要:目的
对氮(N)、磷(P)添加条件下4种种植密度的樟树Cinnamomum camphora幼苗各器官碳(C)含量、储量和分配比例进行研究，以期为氮沉降和磷添加背景下森林碳储量分配格局的变化提供参考。
方法以1年生樟树幼苗为试验材料，选择氯化铵(NH4Cl)作为氮肥模拟氮沉降，以二水合磷酸二氢钠(NaH2PO4·2H2O)作为磷添加，设置4个水平：不加N和P(对照，CK)，加N，加P，加N和P(N+P)。N、P及N+P每年的添加量分别为NH4Cl 40 g·m−2、NaH2PO4·2H2O 20 g·m−2和NH4Cl 40 g·m−2+ NaH2PO4·2H2O 20 g·m−2；种植密度设置4个水平，即10、20、40和80株·m−2。
结果各氮、磷添加和密度处理下幼苗的根、茎和枝的C含量基本上差异不显著，而添加N和N+P能够促使樟树幼苗叶的C含量上升。随着种植密度的增大，樟树幼苗叶片C含量表现出下降的趋势；N、P添加处理基本上能够促进幼苗单株C储量和单位面积C储量的增加；随着种植密度的增大，单株幼苗C储量呈现下降的趋势。
结论樟树幼苗叶的单株C储量和单位面积C储量分配比例随着种植密度的增大逐渐减小。高密度种植有利于茎的分配比例增加。N+P添加处理对幼苗C储量的促进效果大于单一N或P添加处理。
Abstract:ObjectiveTo study carbon (C) content, storage and allocation in seedling organs of Cinnamomum camphora cultivated in four different densities under nitrogen (N) and phosphorus (P) additions, and provide information for forest C storage and allocation under the background of N deposition and P addition.
MethodThe 1-year-old C. camphora seedlings were used as test materials. NH4Cl and NaH2PO4·2H2O were selected to simulate atmospheric N deposition and P addition, respectively. N and P additions were performed with four different levels (control, N, P, and N+P). The N and P addition amounts per year in N, P, and N+P treatments were 40 g·m−2 NH4Cl, 20 g·m−2 NaH2PO4·2H2O and 40 g·m−2 NH4Cl + 20 g·m−2 NaH2PO4·2H2O, respectively. Seedlings were planted in four different densities (10, 20, 40 and 80 seedlings·m−2).
ResultC contents in roots, stems and branches of seedlings in all treatments had no significant difference. N and N+P treatments increased C content of leaves. With the increase of planting density, C content of leaves tended to decrease. The N and P additions increased C storage per seedling and C storage in unit area. C storage per seedling decreased with the increase of planting density.
ConclusionC storage of leaves per seedling and C storage in unit area decreases with the increase of planting density. C storage percentage of stems increases in high cultivation density treatment. The effect of N+P treatment on C storage per seedling and C storage in unit area is greater than that of single N or P addition.
-
Keywords:
- N and P additions /
- planting density /
- carbon storage /
- Cinnamomum camphora
-
奶牛养殖业在畜牧业发展中具有极其重要的地位[1]。集成化、智能化和标准化养殖已成为保证奶牛进食量、提高牛奶产量及保障乳制品品质的必然选择[2]。对于奶牛饲养，自由采食能提高饲料摄入量和牛奶产量，但奶牛进食时易将饲料拱出饲喂区，造成饲料的分散，不利于奶牛的进食。这样不仅浪费饲料，最重要的是影响奶牛的营养摄入量，最终影响产奶量和产奶质量[3-4]。目前，国外很多牛场采用以Lely、Valmetal公司的PRO-FEED为代表的推料机器人，Lely推料机器人[5]通过旋转的圆柱形外壳将饲料推回到栏杆下，PRO-FEED推料机器人[6]是螺旋式结构，通过滚动的螺旋式叶片将饲料推回到栏杆下。这两种机器人在推进过程中能对饲料进行一定的搅拌，推料效果较好，但结构复杂、加工难度大。现如今国内牛场养殖的推料方式有三种，一是人工扫料，二是人工驾驶的刮板式推料车，三是自走式推料机器人[7]。3种推料工作中，传统的人工扫料方式劳动强度大，需要投入大量的人力[8]。张勤等[9]研发的基于信息融合的智能推料机器人是采用刮板的推料方式，虽然刮板式推料车刮板小、结构简单、便于操作，但推料效果一般，且在纵向方向上有推料现象，影响奶牛采食。现有的自旋式推料机器人[5]推料效果较好，没有漏料情况，但推料效率较低；螺旋式推料机器人推料效率较高，但推料效果不好，有漏料现象。
对于皮带式推料机，具有代表性的是德国RMH-futterschieber皮带式推料机，针对含水率高、撒料量和撒料面积大的青饲料的推料效果较好。因此，本文对皮带式推料机结构进行设计，并对其推料过程进行仿真试验，为自动推料机实际推料作业参数提供参考。该研究利用离散元软件EDEM创建颗粒实体的参数化模型，此外，将材料力学特性和其他物理特性添加到EDEM中，形成颗粒模型和力学分析过程；基于Design-Expert软件进行正交试验设计，运用EDEM在不同试验条件下对皮带式推料机进行虚拟试验仿真和分析，探究推料机推料过程中物料颗粒群的运动规律，分析推料机不同推料角度、皮带速度、前进速度3个试验因素对推送率和推送效率的影响，以期为寻求推料机推料最优参数组合提供一定的指导。
1. 推料机器人结构及工作原理
1.1 整机结构
可变角度皮带式推料机器人主要由大容量锂电池、机器人主体部分、控制系统、精料抛撒装置、TMR搅拌撒料装置、激光雷达、RGB相机和IMU以及皮带式推料机组成，整机结构如图1所示。
![]()
图 1 推料机器人结构简图1：机器人主体部分，2：控制系统，3：精料抛撒装置，4：TMR饲料搅拌装置，5：RGB相机和IMU，6：激光雷达，7：皮带式推料机Figure 1. Structure sketch of pushing robot1: The main part of the robot, 2: Control system, 3: Concentrate throwing device, 4: TMR feed stirring device, 5: RGB camera and IMU, 6: LiDAR, 7: Belt type pusher工作时，首先，推料机器人融合RGB和雷达信息智能识别牛舍饲料带，控制机器人沿着饲料带自主行走；其次，根据饲料带是否存在饲料，机器人控制推料电机带动推料皮带转动，皮带上的推板将草料推回饲料带，从而完成推料作业。
1.2 皮带式推料机结构
皮带式推料机采用的是立式结构，外形尺寸长×宽×高为1 000 mm × 460 mm × 550 mm。如图2所示，皮带式推料机主要由1个伺服电机、2个滚筒、软质PVC带、推料机支撑架和可变角度连接架组成。
![]()
图 2 皮带式推料机结构图1：滚筒，2：皮带，3：伺服电机，4：张紧装置，5：皮带机支撑架Figure 2. Structure sketch of belt pusher1: Roller, 2: Belt, 3: Servo motor, 4: Tensioning device, 5: Belt conveyor support frame其中，滚筒的尺寸直接影响其与皮带的接触面积与摩擦力，这将影响皮带推料效率。滚筒的直径越大，皮带与滚筒的接触面积越大，两者之间的摩擦力越大，皮带运动越稳定。然而，滚筒直径过大会导致设备整机尺寸增大，增大能量消耗。因此，本设计选用的滚筒直径为155 mm、高度为450 mm、材质为碳钢、壁厚为2.5 mm，滚筒表面采用车床粗加工，以增加表面的摩擦力，满足推料作业的实际需求，三维模型如图3所示。
2. TMR饲料特性参数测量与接触参数仿真标定
2.1 试验材料
TMR饲料于2024年4月取自山东省泰安市岱岳区满庄镇泥沟村的金兰奶牛养殖牧场，分别在撒料1、2、3 h后进行取料，含水率(w)为55%~70%。
2.2 饲料特性参数测量
2.2.1 TMR饲料休止角测量
利用TMR饲料休止角对饲料间的恢复系数进行仿真标定试验，间接得出饲料间的恢复系数。采用排出法对TMR饲料的休止角进行测定，如图4所示。将一定量的饲料装入无底圆筒中，向上提升圆筒使饲料缓慢排出，待料堆稳定后，测量料堆竖直高度(h)与底面直径(d)的尺寸大小，通过公式(1)计算其休止角[10]，试验重复3次取平均值，测得休止角为45°~47°，平均46.75°。为降低摩擦因数对测量结果的影响，圆筒材质选用钢制。
![]()
图 4 休止角测定h：料堆竖直高度，d：料堆底面直径Figure 4. Determination of rest angleh: Vertical height of feed pile, d: Basal diameter of feed pile休止角计算公式：
$$ \varphi = {\tan ^{ - 1}}\dfrac{{2h}}{d} 。 $$ (1) 2.2.2 滑动摩擦角测量
通过饲料与皮带、饲料与地面间的滑动摩擦角，可以仿真标定饲料与皮带、饲料与瓷砖地面间的恢复系数。饲料滑动摩擦角测量所用仪器为斜面仪，摩擦面材料为PVC皮带和瓷砖。将所取饲料样本放置在水平摩擦面上，缓慢匀速转动手柄，使摩擦面一端缓慢上升，避免振动。如图5所示，当饲料试样开始向下滑动时，记录此时的摩擦面倾角，即为饲料滑动摩擦角(θ)。测得饲料与皮带间的滑动摩擦角平均值为37°，饲料与瓷砖间的滑动摩擦角平均值为42°。
![]()
图 5 滑动摩擦角测量θ：饲料滑动摩擦角Figure 5. Measurement of sliding friction angleθ: Sliding friction angle of feed2.2.3 摩擦系数测定
使用MXD-2型摩擦因数试验仪，测定饲料之间、饲料与皮带之间以及饲料与瓷砖地面之间的动、静摩擦系数，如图6所示。测定饲料之间的动、静摩擦系数时，将饲料用胶水平整均匀地固定在水平试验台上，并在试验仪滑块底面也均匀粘固一层饲料，使2个面上的饲料平整接触，启动试验仪后2个接触面上的饲料产生相对移动，滑块产生一定距离的相对移动后，仪器根据摩擦力的大小自动计算摩擦系数。试验重复5次取平均值，结果如表1所示。
表 1 试验测定摩擦系数Table 1. Determination of friction coefficient by tests参数
Parameter饲料与饲料
Feed and
feed饲料与皮带
Feed and
belt饲料与地面
Feed and
ground静摩擦系数
Static friction
coefficient0.54 0.75 0.80 滚动摩擦系数
Rolling friction
coefficient0.85 0.40 0.48 2.3 TMR饲料接触参数仿真标定
采用试验测试与离散元仿真相结合的方法对饲料进行仿真参数标定及优化。试验时，以实测的休止角、滑动摩擦角为响应值设计中心组合试验方案，利用Design-Expert软件分析各因素的编码值和数据，对各项进行方差分析和回归系数显著性检验，分析仿真参数对仿真结果的影响，确定饲料EDEM的最优参数。TMR饲料仿真离散元模型建立方法如下。
1)TMR饲料模型。由作物秸秆、球形米粒和青贮块组成奶牛饲料，并进行仿真建模。采用理想球形颗粒组合法对实际物料进行替代。对混合物料几何尺寸进行实测并取平均值[11]，3种球形填充物料颗粒的直径分别设为12、10和8 mm。物料离散元模型如图7a~7c所示。
2)接触模型。利用EDEM软件对推料作业进行仿真。由于不考虑传热及磨损等问题，选用Hertz Mindlin(no slip)模型作为离散元仿真中饲料颗粒与皮带的接触模型[12]。
3)碰撞恢复系数仿真标定试验。利用EDEM离散元参数标定对饲料堆积过程进行仿真，如图8所示；对饲料滑动摩擦角测定过程进行仿真，如图9所示。通过优化碰撞恢复系数，使饲料仿真堆积过程的休止角与实际休止角一致，以及饲料仿真过程的滑动摩擦角与实际滑动摩擦角一致，从而预测碰撞恢复系数，为推料仿真分析提供基础数据[13]。
![]()
图 9 TMR饲料滑动摩擦角测定过程模拟θ：饲料滑动摩擦角Figure 9. Simulation of TMR feed sliding friction angle measurement processθ: Sliding friction angle of feed4)仿真参数。TMR饲料模型参数包括材料参数和接触参数。其中，材料参数包括TMR饲料、皮带、地面的密度、泊松比、剪切模量等[5, 10, 14-15]，具体数值如表2所示。
表 2 物料特性参数Table 2. Property parameters of materials材料
Material密度/(kg·m−3)
Density泊松比
Poisson’s ratio剪切模量/Pa
Shear modulus饲料 Feed 700 0.30 2.30×107 皮带 Belt 1 380 0.45 2.40×108 地面 Ground 1 500 0.20 1.25×1010 通过以上测量试验、仿真测定及查阅相关文献，得到饲料与饲料之间、饲料与推料皮带装置之间、饲料与瓷砖地面之间的碰撞恢复系数均为0.30。
3. EDEM推料仿真试验
3.1 EDEM推料仿真模型及参数设置
为便于推料仿真的计算，本研究去除了与物料运动无接触关系的部件。推料机的皮带与饲料直接接触，其材质为软质PVC材料，将其导入EDEM中，如图10所示。
根据表1、2和标定的碰撞恢复系数设置仿真参数，通过动态生成的方式在颗粒工厂中生成物料颗粒，设定颗粒总质量为180 kg，作物秸秆、球形米粒、青贮块的质量比为35∶10∶55。
根据相关文献[14]可知，推送率随着饲料含水率的增加而增大。本研究推料机皮带速度、推料角度以及前进速度的定义图如图11所示。
在推料机构作用下，饲料颗粒间相互作用，饲料颗粒沿着皮带表面速度方向运动。如图12所示，当推料角度为0°时，饲料运动轨迹垂直运动方向，若运动速度增大会将未及时推送到饲喂区的饲料向斜前方推送，出现饲料分布不均匀现象；随着推送角度增大，饲料运动轨迹开始向斜后方运动，在一定前进速度下，饲料可以更好地被推到原位置；但当推料角度超过30°时，会出现将已推过去的饲料带到皮带后面的现象。
![]()
图 12 物料颗粒在各个角度的运动方向各小图中的箭头指饲料运动方向Figure 12. Movement direction of the material particles at various anglesThe arrows in each figure indicate the movement direction of feed选取推送率(Y1)和推送效率(Y2)作为推料效果的评价指标，计算方式如公式(2)和(3)。设置图13中阴影区域的颗粒质量为150 kg，仿真完成后对未被推到撒料带的饲料进行统计。
$$ {Y}_{1}=\dfrac{m-{m}_{_{{0}}}}{m}\times 100{\text{%}} \text{，} $$ (2) $$ {Y}_{2}=\dfrac{m-{m}_{_{{0}}}}{{t}} \text{，} $$ (3) 式中，m：推送前阴影区域饲料的质量，kg；
$m_{_{{0}}} $ ：推送后未推到撒料带饲料的质量，kg；t：推料时间，min。3.2 推料试验仿真分析
以推料机推料角度(x1)、皮带速度(x2)和前进速度(x3)为研究因素，设计三因素三水平试验表，如表3所示。
表 3 试验因素水平表1)Table 3. Table of test factor and level水平
Levelx1/(°) x2/(m·s−1) x3/(m·s−1) −1 0 1 0.5 0 15 2 1.0 1 30 3 1.5 1)x1：推料角度，x2：皮带速度，x3：前进速度
1) x1: Pushing angle, x2: Belt speed, x3: Forward speed以各因素编码值x1、x2、x3为自变量，Y1和Y2为目标值，试验设计与结果如表4所示。
表 4 推料仿真试验设计与结果1)Table 4. Design and results of pushing simulation experiment序号
Serial numberx1 x2 x3 Y1/% Y2/(kg·min−1) 1 0 1 1 92.6 833.4 2 1 1 0 96.5 579.0 3 −1 0 1 91.8 826.2 4 1 −1 0 93.3 559.8 5 1 0 1 92.0 828.0 6 −1 1 0 94.8 568.8 7 0 0 0 93.5 561.0 8 0 −1 −1 93.9 281.7 9 −1 −1 0 91.1 546.6 10 1 0 −1 95.8 287.4 11 0 0 0 94.0 564.0 12 −1 0 −1 95.1 285.3 13 0 −1 1 89.7 807.3 14 0 0 0 93.9 563.4 15 0 0 0 94.3 565.8 16 0 1 −1 98.1 294.3 17 0 0 0 95.0 570.0 1) x1：推料角度，x2：皮带速度，x3：前进速度，Y1：推送率，Y2：推送效率
1) x1: Pushing angle, x2: Belt speed, x3: Forward speed, Y1: Pushing rate, Y2: Pushing efficiency对各因素与推送率和推送效率的关系进行方差分析和回归系数显著性检验，结果如表5、6所示。依据表5，利用响应面法建立推送率与试验因素编码值的关系模型为
表 5 推料仿真试验推送率响应面方差分析Table 5. Variance analysis of response surface of pushing rate in material pushing simulation test来源
Source平方和
Square
sum自由度
Freedom
degree均方
Mean
squareF P1) 模型
Model64.08 9 7.12 15.44 0.000 8** x1 2.88 1 2.88 6.25 0.042 0* x2 24.50 1 24.50 53.15 0.000 2** x3 35.28 1 35.28 76.53 <0.000 1** x1x2 0.06 1 0.06 0.14 0.723 6 x1x3 0.06 1 0.06 0.14 0.723 6 x2x3 0.42 1 0.42 0.92 0.370 3 x12 0.01 1 0.01 0.03 0.867 0 x22 0.10 1 0.10 0.23 0.648 6 x32 0.70 1 0.70 1.52 0.257 9 残差
Residual error3.23 7 0.46 失拟项
Misfit term1.98 3 0.66 2.10 0.240 0 纯误差
Pure error1.30 4 0.31 总和
Sum67.30 16 1)“*”“**”分别表示在P<0.05和P<0.01水平影响显著(方差分析)
1) “*” and “**” indicate significant effects at P<0.05 and P<0.01 levels, respectively (Analysis of variance)Y1=94.140 0+0.600 0x1+1.750.0 x2−2.100 0x3−0.125 0x1x2−0.125 0x1x3+0.325 0x2x3−0.057 5x12−0.157 5x22−0.408 0x32。
由表5可知，推送率与试验因素编码值之间的关系是极显著的(P<0.01)，失拟项影响不显著(P>0.05)，表示模型合理。其中x1对推料率影响显著，x2和x3对推送率影响极显著，其余各项都影响不显著。从单因素分析来看，各因素对推送率的影响由大到小为x3>x2>x1。从各因素交互作用来看，交互作用项对推送率的影响由大到小为x2x3>x1x3>x1x2。
各因素的交互作用对推送率的响应面如图14。当推料角度在中间水平时，推送率随皮带速度增加而增大；当皮带速度在中间水平时，推送率随推料角度增加而增大。当推料角度在中间水平时，推送率随前进速度增加而明显降低；当前进速度在中间水平时，推送率随推送角度增加而增加。当前进速度在中间水平时，推送率随皮带速度增加而明显增高。
![]()
图 14 各因素对推送率的响应面图Figure 14. Response surface graph of various factors to pushing rate依据表6，利用响应面法建立推送效率与试验因素编码值之间的关系模型为
表 6 推料仿真试验推送效率响应面方差分析Table 6. Variance analysis of response surface of pushing efficiency in material pushing simulation test来源
Source平方和
Square
sum自由度
Freedom
degree均方
Mean
squareF P1) 模型
Model5.77×105 9 64 117 3 484 <0.000 1** x1 93.16 1 93.16 5.06 0.060 0 x2 802.00 1 802.00 43.60 0.000 3** x3 5.76×105 1 5.76×105 31 293 <0.000 1** x1x2 2.25 1 2.25 0.12 0.736 9 x1x3 0.02 1 0.02 0.00 0.973 1 x2x3 45.56 1 45.56 2.48 0.159 6 x12 1.67 1 1.67 0.09 0.771 9 x22 15.52 1 15.52 0.84 0.388 9 x32 322.00 1 322.00 17.50 0.004 1** 残差
Residual error128.79 7 18.40 失拟项
Misfit term83.72 3 27.91 2.48 0.200 8 纯误差
Pure error45.07 4 11.27 总和
Sum5.77×105 16 1)“**”表示在P<0.01水平影响显著(方差分析)
1) “**” indicates significant effect at P<0.01 level (Analysis of variance)Y2=564.840+3.410x1+10.010x2+268.280x3−0.750x1x2−0.075x1x3+3.380x2x3+0.630x12−1.920x22−8.740x32。
由表6可知，推送效率与试验因素编码值之间的关系是极显著的(P<0.01)，失拟项影响不显著(P>0.05)；这表明模型合理，可以用来预测推料机推送效率。其中，x2、x3和x32对推料机推送效率影响极显著，其余各项均不显著。从单因素分析来看，各因素对推送效率的影响由大到小为x3>x2>x1。从各因素交互作用来看，交互作用项对推送效率的影响由大到小为x2x3>x1x2>x1x3。
各因素的交互作用对推送效率的响应面如图15。推料角度和皮带速度对推料效率的影响不明显。当推料角度和皮带速度分别在中间水平时，推料效率都会随着前进速度增加而明显增高。
![]()
图 15 各因素对推送效率的响应面图Figure 15. Response surface graph of various factors to pushing efficiency以推料机的推送率和推料效率最大为目标值，通过响应面法进行优化求解，得到推料机理论最优作业参数水平组合：推料角度30°、皮带速度2.99 m/s、前进速度0.83 m/s。在此作业参数下，推送率和推送效率分别为94.6%、724.5 kg/min。
4. 样机的推料试验
4.1 试验方案
如图16所示，在泰安市金兰奶牛场进行试验，选取推料机理论最优作业参数水平组合和仿真结果较好的2组试验参数进行试验。每次推料前对图13阴影区域TMR饲料的质量进行记录，推料作业完成后收集未推到撒料区的TMR饲料并称质量，推料过程中TMR饲料水分损失忽略不计，每组试验重复3次。选取推送率和推送效率作为推料效果的评价指标，将结果与EDEM仿真试验结果做对比，测试皮带式推料机器人在真实工作情况下的推料效果。
4.2 试验结果与分析
利用TMR饲料进行试验，推料机器人3次推料试验推送率为91%~96%，推送效率为720~820 kg/min。当推料前进速度为1.5 m/s时，机器人推料效率最高；但是当皮带速度低于2 m/s时，皮带无法及时完成推料作业，导致大量饲料堆积在皮带前方。当推料机器人推送到位的饲料较少，剩余饲料会导致饲料流反向运动，造成饲料分布不均匀。前进速度和皮带速度一定时，当推料角度为0°时，饲料运动轨迹是斜前方，当推送角度逐渐增大时，饲料运动轨迹逐渐向斜后方运动，但是推送角度过大会导致出现将已推过去的饲料带到皮带后面的现象，即漏料现象。当皮带速度为3 m/s时，无论推料机器人的前进速度是高速还是低速，在3个工作角度下的推料效果都是良好的；在皮带高速转动的情况下，皮带通过摩擦力和皮带表面凸起将饲料从一侧带动到另一侧，高速转动的皮带能快速地将饲料推送到指定位置，从而验证了本文的推料机推料效果可以满足牧场需求。实际试验结果与仿真结果基本一致，验证了推料机结构和作业参数的合理性。
5. 结论
针对目前国内采用螺旋式推料机器人对含水率高、撒料量多的TMR饲料推料效果不好等问题，研制了一种皮带式推料机，采用旋转的皮带将饲料推回到奶牛采食区域。皮带式推料机试验表明，推送效率达到牧场的要求，仿真模拟结果为皮带式推料机的动力学分析和参数优化提供了参考。
-
表 1 幼苗处理前各器官单株生物量和碳含量(平均值±标准误差)
Table 1 Biomass and carbon content per seedling before treatment (mean ± SE)
器官
Organ单株生物量/g
Biomass per seedlingw(C)/(g·kg−1) 根 Root 2.18±1.12 469.57±0.34 茎 Stem 1.11±0.25 487.02±1.92 枝 Branch 0.11±0.08 460.53±1.49 叶 Leaf 0.88±0.31 515.13±2.56 
下载: 导出CSV
表 2 不同种植密度及氮、磷添加处理幼苗各器官的碳含量1)(平均值±标准误差)
Table 2 Carbon contents of organs per seedling under different planting densities and nitrogen and phosphorus additions (mean±SE)
种植密度/(株·m−2)
Planting densityN、P添加
N and P additionsw(C)/(g·kg−1) 根 Root 茎 Stem 枝 Branch 叶 Leaf 10 CK 452.75±10.52Aa 492.11±21.67Aa 478.82±7.42Aa 477.81±3.18Bc N 447.34±4.82Ba 494.13±2.16Aa 477.82±3.31Ba 499.36±4.91Ab P 446.27±3.87Aa 500.56±1.57Aa 481.53±4.88ABa 501.17±1.56Ab N+P 462.95±2.60Aa 489.97±1.04Ba 481.98±5.13Aa 512.34±3.61Aa 20 CK 460.81±5.39Aa 500.11±1.83Aa 481.11±3.66Aa 494.30±2.87Aab N 457.91±8.47ABa 503.78±7.88Aa 480.64±3.53ABa 502.38±4.51Aa P 466.52±3.48Aa 496.91±2.59Aa 474.63±4.16Ba 484.32±4.44ABb N+P 460.04±4.44Aa 500.79±1.46Aa 482.60±5.00Aa 503.61±4.71Aa 40 CK 459.78±1.86Aab 497.80±2.11Aa 474.88±5.83Aab 474.64±5.10Bab N 469.60±5.10Aab 499.21±8.71Aa 488.58±3.32Aa 480.58±6.92Bab P 447.74±12.65Ab 485.44±5.03Bab 471.16±4.35Bb 465.51±4.06BCb N+P 471.79±4.12Aa 472.21±1.47Cb 475.74±4.70Aab 484.62±4.63Ba 80 CK 457.01±5.11Aa 490.15±7.01Aa 492.16±6.15Aa 480.11±2.77Ba N 457.79±6.02ABa 492.82±2.25Aa 482.77±2.22ABab 483.05±1.17Ba P 449.37±3.65Aa 499.93±3.63Aa 490.58±3.17Aa 448.61±14.61Cb N+P 452.35±14.07Aa 496.53±3.24Aa 475.81±5.28Ab 486.18±4.52Ba 1) 同列数据后，不同大写字母表示相同N、P添加处理下不同种植密度间差异显著，不同小写字母表示相同种植密度下不同N、P添加处理间差异显著(P<0.05，Duncan’ s法)
1)Different capital letters in the same column indicate significant differences among different planting densities under the same N/P addition treatment, and different lowercase letters in the same column indicate significant differences among different N/P addition treatments under the same planting density (P<0.05, Duncan’ s test)
下载: 导出CSV
表 3 不同种植密度及氮、磷添加处理下单株幼苗各器官的碳储量和分配比例1)
Table 3 Carbon storage of organs and its allocation per seedling under different planting densities and nitrogen and phosphorus additions
种植密度/
(株·m−2)
Planting densityN、P添加
N and P additions根 Root 茎 Stem 枝 Branch 叶 Leaf 全株 Total C储量/g
C storage分配比例/%
RatioC储量/g
C storage分配比例/%
RatioC储量/g
C storage分配比例/%
RatioC储量/g
C storage分配比例/%
RatioC储量/g
C storage分配比例/%
Ratio10 CK 1.19Aa 41.66 0.97Aa 33.95 0.16Bb 5.47 0.54Ac 18.92 2.86Ab 100 N 2.04Aa 35.06 1.30Aa 22.32 0.78Aa 13.37 1.70Ab 29.25 5.82Aa 100 P 1.90Aa 32.78 1.04ABa 17.94 0.93Aa 16.00 1.93Ab 33.28 5.79Aa 100 N+P 2.13Aa 29.27 1.30Aa 17.86 1.05Aa 14.50 2.79Aa 38.37 7.27Aa 100 20 CK 1.00Aa 39.16 0.91ABa 35.33 0.15Bb 5.78 0.51Ac 19.73 2.57ABb 100 N 1.37Ba 32.94 1.07Aa 25.74 0.52Ba 12.38 1.21Ab 28.93 4.17Ba 100 P 1.42Ba 39.72 1.16Aa 32.53 0.26Bb 7.20 0.73Bc 20.55 3.58Ba 100 N+P 1.19Ba 28.11 0.98Aa 23.18 0.51Ba 11.94 1.56Ba 36.76 4.24Ba 100 40 CK 0.79Ab 37.97 0.75Ba 35.96 0.14Bb 6.83 0.40Ab 19.23 2.09BCb 100 N 1.36Ba 46.43 0.96Aa 32.86 0.23Cab 7.74 0.38Bb 12.97 2.93Ca 100 P 1.34Ba 47.16 0.80Ca 28.24 0.27Ba 9.34 0.43BCb 15.26 2.84BCa 100 N+P 1.27Ba 39.60 1.00Aa 31.36 0.31BCa 9.72 0.62Ca 19.33 3.20BCa 100 80 CK 0.79Aa 42.11 0.75Ba 39.84 0.22Aa 11.85 0.12Ba 6.20 1.88Ca 100 N 1.00Ba 43.05 0.90Aa 38.54 0.23Ca 9.74 0.20Ba 8.66 2.33Ca 100 P 0.71Ca 36.03 0.85BCa 42.94 0.21Ba 10.71 0.20Ca 10.33 1.98Ca 100 N+P 1.07Ba 42.72 0.99Aa 39.66 0.18Ca 7.01 0.27Ca 10.61 2.50Ca 100 1) 同列数据后，不同大写字母表示相同N、P添加处理下不同种植密度间差异显著，不同小写字母表示相同种植密度下不同N、P添加处理间差异显著(P<0.05，Duncan’ s法)
1) Different capital letters in the same column indicate significant differences among different planting densities under the same N/P addition, and different lowercase letters in the same column indicate significant differences among different N/P addition under the same planting density (P<0.05, Duncan’ s test)
下载: 导出CSV
表 4 不同种植密度和氮、磷添加对单株幼苗各器官碳储量交互作用的方差分析
Table 4 Variance analyses of interaction effects of planting density and nitrogen and phosphorus additions on carbon storage of organs per seedling
器官碳储量
Carbon storage of organs种植密度
Planting densityN、P添加
N and P additions种植密度与N、P添加交互作用
Interaction effect of planting density and N and P additionsF P F P F P 根 Root 17.219 0.000 6.261 0.001 1.003 0.441 茎 Stem 6.800 0.000 4.157 0.007 0.745 0.667 枝 Branch 33.459 0.000 13.035 0.000 5.281 0.000 叶 Leaf 93.367 0.000 28.212 0.000 8.978 0.000 全株 Total 54.661 0.000 18.687 0.000 3.565 0.001
下载: 导出CSV
表 5 不同种植密度及氮、磷添加处理下幼苗各器官的单位面积碳储量和分配比例1)
Table 5 Carbon storage of organs and its allocation in unit area under different planting densities and nitrogen and phosphorus additions
种植密度/
(株·m−2)
Planting densityN、P添加
N and P additions根 Root 茎 Stem 枝 Branch 叶 Leaf 全株 Total C储量/
(g·m−2)
C storage分配比
例/%
RatioC储量/
(g·m−2)
C storage分配比
例/%
RatioC储量/
(g·m−2)
C storage分配比
例 /%
RatioC储量/
(g·m−2)
C storage分配比
例/%
RatioC储量/
(g·m−2)
C storage分配比
例/%
Ratio10 CK 11.89Ca 41.66 9.70Da 33.95 1.56Cb 5.47 5.40Cc 18.92 28.55Db 100 N 20.40Ca 35.00 12.98Ca 22.32 7.78Ba 13.37 17.02Bb 29.55 58.18Ca 100 P 18.99Ba 32.78 10.39Ca 17.94 9.27Ca 16.00 19.28Ab 33.28 57.93Ca 100 N+P 21.28Ca 29.27 12.99Ca 17.86 10.54Aa 14.50 27.90Aa 38.37 72.71Ca 100 20 CK 20.07BCa 39.16 18.11Ca 35.33 2.96Cb 5.78 10.11Bc 19.73 51.25Cb 100 N 27.47Ca 32.94 21.46BCa 25.74 10.33Ba 12.38 24.13Bb 28.93 83.39BCa 100 P 28.41Ba 39.72 23.27Ba 32.53 5.15BCb 7.20 14.70Ac 20.55 71.53Ba 100 N+P 23.84Ca 28.11 19.66Ca 23.18 10.13Aa 11.94 31.18Aa 36.76 84.80BCa 100 40 CK 31.80Bb 37.97 30.11Ba 35.96 5.72Bb 6.83 16.11Ab 19.23 83.74Bb 100 N 54.37Ba 46.43 38.47Ba 32.86 9.07Bb 7.74 15.18Bb 12.97 117.10Ba 100 P 53.62Aa 47.16 32.10Ba 28.24 10.62Ba 9.34 17.34Ab 15.26 113.68Ba 100 N+P 50.62Ba 39.60 40.08Ba 31.36 12.42Aa 9.72 24.70Aa 19.33 127.82Ba 100 80 CK 63.24Aa 42.11 59.83Aa 39.84 17.80Aa 11.85 9.31Ba 6.20 150.19Aa 100 N 80.25Aa 43.05 71.86Aa 38.54 18.17Aa 9.74 16.15Aa 8.66 186.43Aa 100 P 57.14Aa 36.02 68.12Aa 42.94 16.99Aa 10.71 16.39Aa 10.33 158.64Aa 100 N+P 85.55Aa 42.72 79.42Aa 39.66 14.03Aa 7.01 21.25Aa 10.61 200.24Aa 100 1) 同列数据后，不同大写字母表示相同N、P添加处理下不同种植密度间差异显著，不同小写字母表示相同种植密度下不同N、P添加处理间差异显著(P<0.05，Duncan’ s法)
1)Different capital letters in the same column indicate significant differences among different planting densities under the same N/P addition, and different lowercase letters in the same column indicate significant differences among different N/P additions under the same planting density (P<0.05, Duncan’ s test)
下载: 导出CSV
表 6 不同种植密度及氮、磷添加对幼苗各器官的单位面积碳储量交互作用的方差分析
Table 6 Variance analyses of interaction effects of planting density and nitrogen and phosphorus additions on carbon storage of organs in unit area
器官单位面积碳储量
Carbon storage of organs in
unit area种植密度
Planting densityN、P添加
N and P additions种植密度与N、P添加交互作用
Interaction effect of planting density and N and P additionsF P F P F P 根 Root 51.780 0.000 3.743 0.013 1.003 0.313 茎 Stem 127.420 0.000 2.741 0.047 0.676 0.730 枝 Branch 29.407 0.000 6.744 0.000 2.670 0.007 叶 Leaf 1.959 0.123 26.787 0.000 2.340 0.017 全株 Total 201.698 0.000 3.062 0.030 0.874 0.550
下载: 导出CSV
-
[1] 辛小娟, 王刚, 杨莹博, 等. 氮、磷添加对亚高山草甸地上/地下生物量分配的影响[J]. 生态科学, 2014, 33(3): 452-458. [2] 陈瑶, 王法明, 莫其锋, 等. 氮磷添加对华南热带森林尾叶桉木质残体分解和养分动态的影响[J]. 应用与环境生物学报, 2015, 21(4): 747-753. [3] CHEN H, ZHANG W, GILLIAM F, et al. Changes in soil carbon sequestration in Pinus massoniana forests along an urban-to-rural gradient of southern China[J]. Biogeosciences, 2013, 10(10): 6609-6616. doi: 10.5194/bg-10-6609-2013
[4] 刘顺, 罗达, 刘千里, 等. 川西亚高山不同森林生态系统碳氮储量及其分配格局[J]. 生态学报, 2017, 37(4): 1074-1083. [5] LUO Y, ZHANG X, WANG X, et al. Dissecting variation in biomass conversion factors across China’s forests: Implications for biomass and carbon accounting[J]. PLoS One, 2014, 9(4): e94777. doi: 10.1371/journal.pone.0094777
[6] LIE Z Y, XUE L, JACOBS D F. Allocation of forest biomass across broad precipitation gradients in China’ s forests[J]. Sci Rep, 2018, 8: 10536. doi: 10.1038/s41598-018-28899-5
[7] YAN G Y, XING Y J, WANG J Y, et al. Sequestration of atmospheric CO2 in boreal forest carbon pools in Northeastern China: Effects of nitrogen deposition[J]. Agric Forest Meteorol, 2018, 248: 70-81. doi: 10.1016/j.agrformet.2017.09.015
[8] SPARRIUS L B, SEVINK J, KOOIJMAN A M. Effects of nitrogen deposition on soil and vegetation in primary succession stages in inland drift sands[J]. Plant Soil, 2012, 353(1/2): 261-272.
[9] GAO Y, HAO Z, YANG T T, et al. Effects of atmospheric reactive phosphorus deposition on phosphorus transport in a subtropical watershed: A Chinese case study[J]. Environ Pollut, 2017, 226: 69-78. doi: 10.1016/j.envpol.2017.03.067
[10] 王鑫瑶, 王旭, 朱美玲, 等. 海南主要森林类型植被碳贮量与固碳价值评价[J]. 中南林业科技大学学报, 2017, 37(7): 92-98. [11] 明安刚, 刘世荣, 莫慧华, 等. 南亚热带红锥、杉木纯林与混交林碳贮量比较[J]. 生态学报, 2016, 36(1): 244-251. [12] 郭剑芬, 杨玉盛, 陈光水, 等. 火烧对森林土壤有机碳的影响研究进展[J]. 生态学报, 2015, 35(9): 2800-2809. [13] DEZI S, MEDLYN B E, TONON G, et al. The effect of nitrogen deposition on forest carbon sequestration: a model-based analysis[J]. Global Change Biol, 2010, 16(5): 1470-1486. doi: 10.1111/j.1365-2486.2009.02102.x
[14] ZAEHLE S, CIAIS P, FRIEND A D, et al. Carbon benefits of anthropogenic reactive nitrogen offset by nitrous oxide emissions[J]. Nat Geosci, 2011(9): 601-605.
[15] 朱仕明, 肖玲玲, 薛立, 等. 密度对乐昌含笑幼苗的生长和生物量的影响[J]. 中南林业科技大学学报, 2015, 35(8): 77-80. [16] 肖玲玲, 朱仕明, 胡继文, 等. 不同密度条件下樟树幼苗生长和幼苗重量分配格局[J]. 安徽农业大学学报, 2015, 42(3): 353-356. [17] 董喜光, 张越, 薛立, 等. 不同密度的山杜英幼苗生长分析[J]. 生态科学, 2016, 35(4): 86-90. [18] 伍艳芳, 肖复明, 徐海宁, 等. 樟树全基因组调查[J]. 植物遗传资源学报, 2014, 15(1): 149-152. [19] 张赐成, 韩广, 关华德, 等. 樟树和桂花树光合最适温度对环境温度改变的响应[J]. 生态学杂志, 2014, 33(11): 2980-2987. [20] 王卓敏, 郑欣颖, 薛立. 樟树幼苗对干旱胁迫和种植密度的生理响应[J]. 生态学杂志, 2017, 36(6): 1495-1502. [21] 文丽, 雷丕锋, 戴凌. 不同林龄樟树林土壤碳氮贮量及分布特征[J]. 中南林业科技大学学报, 2014, 34(6): 106-111. doi: 10.3969/j.issn.1673-923X.2014.06.021 [22] 文汲, 闫文德, 刘益君, 等. 施氮对亚热带樟树人工林土壤氮矿化的影响[J]. 中南林业科技大学学报, 2015, 35(5): 103-108. [23] 赵晶, 闫文德, 郑威, 等. 樟树人工林凋落物养分含量及归还量对氮沉降的响应[J]. 生态学报, 2016, 36(2): 350-359. [24] AERTS R, VANLOGTESTIJN R, KARLSSON P S. Nitrogen supply differentially affects litter decomposition rates and nitrogen dynamics of sub-arctic bog species[J]. Oecologia, 2006, 146(4): 652-658. doi: 10.1007/s00442-005-0247-5
[25] 佘汉基, 郑欣颖, 薛立, 等. 外源性氮和磷对尾叶桉凋落叶分解的影响[J]. 安徽农业大学学报, 2017, 44(3): 409-414. [26] 鲍士旦. 土壤农化分析[M]. 北京: 中国农业出版社, 2000: 30-106. [27] 陈莉, 宋敏, 宋同清, 等. 广西不同林龄软阔林碳储量及其分配格局[J]. 生态学杂志, 2017, 36(3): 592-600. [28] CURREY P M, JOHNSON D, DAWSON L A, et al. Five years of simulated atmospheric nitrogen deposition have only subtle effects on the fate of newly synthesized carbon in Calluna vulgaris and Eriophorum vaginatum[J]. Soil Biol Biochem, 2011, 43(3): 495-502. doi: 10.1016/j.soilbio.2010.11.003
[29] 张慧, 郭卫红, 杨秀清, 等. 麻栎种源林叶片碳、氮、磷化学计量特征的变异[J]. 应用生态学报, 2016, 27(7): 2225-2230. [30] 翁俊, 顾鸿昊, 王志坤, 等. 氮沉降对毛竹叶片生态化学计量特征的影响[J]. 生态科学, 2015, 34(2): 63-70. [31] 蒋思思, 魏丽萍, 杨松, 等. 不同种源油松幼苗的光合色素和非结构性碳水化合物对模拟氮沉降的短期响应[J]. 生态学报, 2015, 35(21): 7061-7070. [32] 吴家兵, 井艳丽, 关德新, 等. 氮沉降对森林碳汇功能影响的研究进展[J]. 世界林业研究, 2012, 25(2): 12-16. [33] HOFMOCKEL K S. Increases in the flux of carbon belowground stimulate nitrogen uptake and sustain the long-term enhancement of forest productivity under elevated CO2[J]. Ecol Lett, 2011, 14(4): 349-57. doi: 10.1111/j.1461-0248.2011.01593.x
[34] 黎磊, 周道玮, 盛连喜. 密度制约决定的植物生物量分配格局[J]. 生态学杂志, 2011, 30(8): 1579-1589. [35] 沈杰, 蔡艳, 何玉亭, 等. 种植密度对烤烟养分吸收及品质形成的影响[J]. 西北农林科技大学学报(自然科学版), 2016, 44(10): 51-58. [36] LU X K, MO J M, GILLIAM F S, et al. Nitrogen addition shapes soil phosphorus availability in two reforested tropical forests in southern China[J]. Biotropica, 2012, 44(3): 302-311. doi: 10.1111/j.1744-7429.2011.00831.x
[37] 苗惠田, 吕家珑, 张文菊, 等. 潮土小麦碳氮含量对长期不同施肥模式的响应[J]. 植物营养与肥料学报, 2015, 21(1): 72-80. doi: 10.11674/zwyf.2015.0108 [38] ELSER J J, BRACKEN M E S, CLELAND E E, et al. Global analysis of nitrogen and phosphorus limitation of primary producers in freshwater, marine and terrestrial ecosystems[J]. Ecol Lett, 2007, 10(12): 1135. doi: 10.1111/j.1461-0248.2007.01113.x
[39] 何利元, 胡中民, 郭群, 等. 氮磷添加对内蒙古温带草地地上生物量的影响[J]. 应用生态学报, 2015, 26(8): 2291-2297. [40] 黄钰辉, 张卫强, 甘先华, 等. 南亚热带杉木林改造不同树种配置模式的土壤质量评价[J]. 中国水土保持科学, 2017, 15(3): 123-130. [41] 樊后保, 刘文飞, 李燕燕, 等. 杉木(Cunninghamia lanceolata)人工林生长与土壤养分对氮沉降的响应[J]. 生态学报, 2007, 27(11): 4630-4642. doi: 10.3321/j.issn:1000-0933.2007.11.031 [42] 詹书侠, 陈伏生, 胡小飞, 等. 中亚热带丘陵红壤区森林演替典型阶段土壤氮磷有效性[J]. 生态学报, 2009, 29(9): 4673-4680. doi: 10.3321/j.issn:1000-0933.2009.09.010 [43] 苏波, 韩兴国, 渠春梅, 等. 森林土壤氮素可利用性的影响因素研究综述[J]. 生态学杂志, 2002, 21(2): 40-46. doi: 10.3321/j.issn:1000-4890.2002.02.011 [44] XIA J Y, WAN S Q. Global response patterns of terrestrial plant species to nitrogen addition[J]. New Phytolog, 2008, 179(2): 428-439. doi: 10.1111/j.1469-8137.2008.02488.x
[45] REICH P B, LUO Y, BRADFORD J B, et al. Temperature drives global patterns in forest biomass distribution in leaves, stems, and roots[J]. Proc Nat Acad Sci USA, 2014, 111(38): 13721-13726. doi: 10.1073/pnas.1216053111
[46] XUE L, HAGIHARA A. Density effects on organs in self-thinning Pinus densiflora Sieb. et Zucc. stands[J]. Ecol Res, 2008, 23(4): 689-695. doi: 10.1007/s11284-007-0427-3
[47] XUE L, PAN L, ZHANG R, et al. Density effects on the growth of self-thinning Eucalyptus urophylla stands[J]. Trees, 2011, 25(6): 1021-1031. doi: 10.1007/s00468-011-0576-4
[48] POORTER H, NIKLAS K J, REICH P B, et al. Biomass allocation to leaves, stems and roots: Meta-analyses of interspecific variation and environmental control[J]. New Phytolog, 2012, 193(1): 30-50. doi: 10.1111/j.1469-8137.2011.03952.x
[49] GOWER S T, GHOLZ H L, NAKANE K, et al. Production and carbon allocation patterns of pine forests[J]. Ecol Bull, 1994, 43(43): 115-135.
[50] FANG Y R, ZOU X J, LIE Z Y, et al. Variation in organ biomass with changing climate and forest characteristics across Chinese forests[J]. Forests, 2018, 9: 521. doi: 10.3390/f9090521
计量
- 文章访问数: 1151
- HTML全文浏览量: 3
- PDF下载量: 1388
