Effects of biochar and biochar-based fertilizer on forest growth and soil fertility
-
摘要:目的
探究生物炭及生物炭基肥施用对火力楠Michelia macclurei生长、养分吸收和土壤肥力的影响，为促进生物炭基肥在林地养分管理上的应用提供科学依据。
方法以火力楠幼林为研究对象，设置不施肥(CK)和分别施用复合肥(CF)、稻壳生物炭(RB)、木屑生物炭(WB)、稻壳生物炭基肥(RBF)、木屑生物炭基肥(WBF) 6个处理，测定1年后火力楠的生长指标、Dickson质量指数、养分含量及土壤理化性质。
结果与CK相比，施用生物炭处理的火力楠根生物量、叶面积、叶N含量和根N累积量显著增加，增幅分别为59.1%~62.2%、44.7%~66.0%、29.5%~45.4%和68.6%~92.3%；施用复合肥和生物炭基肥的Dickson质量指数和全株N、P、K累积量显著增加，增幅分别为136.8%~192.2%、167.7%~267.4%、124.4%~224.5%和191.5%~265.5%，对火力楠生长的促进效果优于单施生物炭。与单施复合肥相比，生物炭基肥对植物生长和土壤改良均有促进作用，RBF处理火力楠SPAD值、土壤pH和速效K含量分别显著增加13.4%、6.2%和99.7%；WBF处理火力楠单株总生物量和P、K累积量显著增加31.5%、44.6%和22.1%，土壤容重显著降低17.3%。模糊隶属函数显示，各处理对火力楠的生长及土壤理化性质改善的综合评价排序为WBF > RBF > CF > WB > RB > CK。
结论生物炭和生物炭基肥均能改善土壤理化性质，有效促进火力楠的生长及养分吸收，其中木屑生物炭基肥的综合效果最佳，为林地养分管理提供了新途径。
Abstract:ObjectiveIn order to study the impacts of biochar and biochar-based fertilizer application on the growth, nutrient uptake, and soil fertility of Michelia macclure, and provide a scientific foundation for encouraging the use of biochar in nutrient management on forest land.
MethodTaking the young forest of M. macclure as the research object, six treatments were set up, including no fertilization (CK) and application of compound fertilizer (CF), rice husk biochar (RB), wood chip biochar (WB), rice husk biochar-based fertilizer (RBF), and wood chip biochar-based fertilizer (WBF). After one year, the growth indices, quality indices, nutrient contents, and soil physicochemical properties of M. macclurei were measured.
ResultCompared with CK, the biomass, leaf area, leaf N content, and root N accumulation of M. macclure treated with biochar significantly increased, with the increases of 59.1%−62.2%, 44.7%−66.0%, 29.5%−45.4%, and 68.6%−92.3%, respectively. The Dickson quality index and whole-plant N, P, and K accumulations increased significantly when compound and biochar-based fertilizers were applied, with the increases of 136.8%−192.2%, 167.7%−267.4%, 124.4%−224.5%, and 191.5%−265.5%, respectively. Moreover, the application of compound fertilizer and biochar-based fertilizers had a greater impact on the growth of M. macclurei than the application of biochar alone. Compared to application of compound fertilizer alone, the application of biochar-based fertilizer was found to promote both plant growth and soil improvement. The RBF treatment significantly increased the soil pH, soil quick-acting K content, and SPAD value by 13.4%, 6.2%, and 99.7% respectively, while the WBF treatment significantly increased the whole-plant P and K accumulations and total biomass by 31.5%, 44.6%, and 22.1% respectively, and the soil bulk density significantly decreased by 17.3%. The fuzzy affiliation function revealed that the overall impact of treatments on M. macclurei growth and enhancement of soil physicochemical qualities was evaluated in the following order: WBF > RBF > CF > WB > RB > CK.
ConclusionBoth biochar and biochar-based fertilizer can enhance soil physicochemical qualities and successfully encourage the growth and uptake of nutrients by M. macclurei, and wood chip biochar-based fertilizer yields the greatest overall results, which offers a novel approach for managing nutrients in forest areas.
-
我国是世界第一荔枝生产大国，种植面积和产量分别约占全世界的80%和65%以上，荔枝栽培历史悠久，品种资源和生产经验丰富，其中，广东省是我国荔枝种植大省。荔枝的花量通常都偏多，尤其是大花穗品种，1个花穗其小花量可达数千朵，荔枝常年花量偏大，消耗的养分过多反而影响坐果[1]。荔枝花期的现白点时间、花穗数量、花量、雄花数、雌花数、开花时间、谢花时间等指标要素对荔枝生产精准管理格外重要。根据果树花期长势监测结果，在修剪促萌、灌水、环割作业时间或药剂处理上进行精准控制，以实施精准促花壮花或疏化保花等措施，有利于提高坐果率。近年来，多项研究表明，在花期进行合理的管控有利于果实产量与品质的提升，如疏花加快了无核荔枝果实生长发育，并提高了果实的产量和品质[2]。在花穗生长中期修剪花穗，能明显提高单穗挂果数和产量，综合效果最好[3]。通过“一疏二控三割”控穗疏花技术可以显著提高‘妃子笑’的最终坐果率、株产和果实固酸比[4]。不同程度的疏花方式均能够显著减少花量，缩短开花时期，并增加坐果量[5]。疏花可以提高单个花穗的雌花比例，提前一批雄花和第一批雌花的开放高峰，延长雌花的持续开放时间，进而增加初始坐果量，降低生理落果[6]。此外，随着荔枝种植规模的扩大，人力成本越来越高，基于无人机遥感的大面积果园果树花期长势的智能监测需求日益增加。
近距离拍摄荔枝花的图像，并利用计算机视觉和机器学习的方法进行花期识别和诊断，是目前主流的研究方法。Wang等[7]提出夜间环境下修正的芒果花分割方法，该方法将拍摄效果好的图像作为参照图像，使用归一化放缩法将其他图像灰度放缩到参照图像的灰度范围内，从而起到修正作用。用固定的阈值分割出芒果花、叶片和枝干后，训练识别花的支持向量机，精确分割出芒果花，并计算花密度，其准确率达到85. 8%。Yousefi等[8]得出基于随机森林的油棕榈开花期分类模型在准确率和F度量值(F-measure)方面要高于K最近邻算法 (K-nearest neighbor，K-NN)和支持向量机 (Support vector machine，SVM)算法，证明了利用随机森林模型对油棕榈开花期4个阶段进行分类的可行性，随机森林模型的阶段分类准确率可以达到72.36%。Lu等[9]通过摄像头进行柑橘开花率评估，利用MobileNetv3作为特征提取器，结合深度可分离卷积进一步加速，改进的柑橘开花率统计模型均值平均精度(Mean average percision，mAP)为84.84%，在保证一定精度的前提下，具有内存占用少、检测速度快等优点。朱永宁等[10]提出结合枸杞开花期和果实成熟期的植物学特征，利用Faster R-CNN方法对图像进行训练、分类，构建枸杞开花期和果实成熟期的识别算法，mAP达到74%。杨其晟等[11]提出一种改进 YOLOv5 的苹果花朵生长状态检测方法，对花蕾、半开、全开、凋落这4类苹果树开花期花朵生长状态进行检测。该方法首先对跨阶段局部网络模块进行改进，并调整模块数量，结合协同注意力模块设计主干网络，提高模型检测性能并减少参数，mAP达到92.2%。熊俊涛等[12]提出一种基于深度语义分割网络识别荔枝花、叶片像素并实现分割的方法，构建深度为34层的ResNet 主干网络，引入稠密特征传递方法和注意力模块，提取荔枝花、叶片的特征；最后通过全卷积网络层对荔枝花、叶片进行分割，模型的平均交并比(Mean intersection over union，mIoU)为0. 734，像素识别准确率达到 87%。邓颖等[13]提出基于实例分割的柑橘花朵识别与花量统计方法，以花期的柑橘树冠图像为样本进行花朵实例的识别及分割，通过对 Mask R-CNN 主体卷积部分和掩膜分支部分的优化，实现对复杂结构图像中密集小尺度柑橘花朵目标的高效检测，获取图像中可见花数量识别神经网络的平均精度为 36.3，花量计算误差为 11.9%。
随着无人机行业的发展，低空遥感已经开始应用于农业农事的判断和管理，无人机低空遥感相比于地面摄像头和人工采集图像，可以获得更加丰富的农作物和地理信息，相比于卫星遥感具有获取方式更灵活、时效高、影像分辨率高等特点，在农情监测、精准农业等领域得到越来越多的应用[14-17 ]。杜蒙蒙等[18]在小麦分蘖期内获取可见光波段的无人机遥感图像，建立植被指数与小麦茎蘖密度地面真值之间的预测模型，在田块尺度上定量反演小麦茎蘖密度，小麦茎蘖密度定量反演BP神经网络模型的均方根误差及平均绝对百分比误差分别为19%及3.62%，具有较高的小麦茎蘖密度预测精度。朱圣等[19]通过获取水稻田的低空遥感图像，分析得到杂草分布图，构建CNN网络对图像的分类，精度达到 92.41%，生成的水稻田杂草分布图效果最好。杨蜀秦等[20]提出一种改进的玉米雄蕊识别模型，能够准确地检测无人机遥感图像中尺寸较小的玉米雄蕊，检测精度分别比Faster R CNN 和YOLO v4模型高26. 22%和3. 42%，为玉米抽雄期的农情监测提供参考。
综合研究表明，目前大多通过地面摄像头和人工手持相机采集荔枝花数据，由于角度受限，导致地面拍摄的花期信息不够全面，无法应用于大规模荔枝园。由于荔枝花主要集中在冠层，侧面花较少，且荔枝树均较高，人工观测荔枝冠层耗时长、收效低，因此利用无人机拍摄冠层可以获取更多荔枝花期的信息，同时能够节省部分人力和物力，一定程度上可以提升农业农事的效率。兰玉彬教授研究团队在基于无人机遥感技术进行果树大面积智能监测方面，积累了较好的研究经验，Mo等[21]从无人机获取的图像中对果树冠层进行实例分割，对于果园的精准管理具有重要意义；李继宇等[22]基于文献计量学，得出农用无人机仍处于活跃的发展阶段，应用场景日趋丰富，纵向发展稍有阻滞，未来可能会在畜牧业和渔业上挖掘到新的研究应用方向。现有的研究鲜见通过荔枝冠层进行观测分析，本研究拟采用大疆精灵P4 Multispectral无人机采集荔枝花期冠层数据，重点优选和优化荔枝开花率评估分类模型，并建立更加精准的荔枝花期长势评估模型，以期为后续荔枝花期疏花保果、精准施肥施药提供决策依据。
1. 材料与方法
1.1 研究区域
研究区域位于广东省广州市从化区荔博园(23°58′N,113°62′E)，如图1a所示，该区域属亚热带季风区，气候温和，雨量充沛，全年平均气温22.0 ℃，年降雨量2305.1 mm，年日照1638.0 h，适宜的温度、肥沃的土壤共同组成了荔枝生长的理想环境。
1.2 图像采集与数据集构建
采用大疆精灵P4 Multispectral拍摄园区荔枝树冠层可见光图像，该无人机采用TimeSync时间同步系统，通过将飞控、相机与RTK的时钟系统进行微秒级同步，实现相机成像时刻的毫秒级误差。系统还对每个相机镜头中心点位置与天线中心位置，结合设备姿态信息，进行实时补偿，使影像获得更加丰富的位置信息。本研究于2021年2月底至3月底拍摄园区花期荔枝冠层，飞行高度约45 m，飞行速度约5 m/s，结合地面站软件 DJI GS PRO 进行航线规划。园区占地面积约20 hm2，拍摄园区荔枝共141棵，通过大疆制图软件对获取图像进行拼接。拼接后图像分辨率为10863像素×10974 像素，试验基地荔枝冠层俯拍图如图1b所示。
为获得单棵荔枝冠层图像，将整张园区冠层图像首先进行分块处理，再运用YOLACT实例分割模型对园区内的冠层进行分割，网络训练优化器为随机梯度下降(Stochastic gradient descent，SGD)，学习率为0.001，权重衰减为5×10−4 ，动量因子为0.9，分割最佳精确度(Average precision，AP)达到96.25%，模型推理结束后对图像进行合并处理，去除冗余重复的部分，并结合图片中每棵树的位置信息，拼接后最终获得单棵荔枝树的冠层图像[22]，冠层分割效果如图2所示。综合地面调查，结合园艺专家的综合判断，对每棵冠层图像标记开花率，标注的开花率分别为0、10%~20%、50%~60%、80%及以上，对应荔枝花期的未开花、分化、初花、盛花期：将开花率为0%的冠层图像定为第0类即未开花；开花率10%~20%的冠层图像定为第1类，即分化；开花率50%~60%的冠层图像定为第2类，即初花期；开花率80%及以上的冠层图像定为第3类，即盛花期。得到565张荔枝原始冠层图像，按照6︰2︰2比例随机划分训练集、验证集和测试集。遥感图像开花率示例如图3所示。
1.3 数据集增强
在已有训练集的基础上，在训练过程中利用PyTorch框架进行多项数据增强，包括随机裁剪、垂直翻转、水平翻转、随机旋转一定角度，可以在训练过程中扩充数据集，增加数据的多样性，提升深度学习的鲁棒性，使得模型学习更加全面。数据增强效果如图4所示。
1.4 开花率评估模型构建
1.4.1 深度卷积神经网络
ResNet网络[23]参考了VGG19网络，在其基础上进行了修改，并通过短路机制加入了残差单元解决退化问题，根据研究得出，随着层数的增加，ResNet的精度随之增加，同时计算复杂度也随之增加。后续提出的多维残差网络(Next dimension resNet，ResNeXt)[24]在网络结构上对ResNet 进行了改进，借鉴了Inception分组[25]的思想，区别于Inception需要人工设计每个分支，ResNeXt的每个分支拓扑结构是相同的，通过多分支卷积网络(Multi-branch convolutional networks)、分组卷积(Group convolution)结构在不明显增加参数量级的情况下提升模型的准确率，同时由于拓扑结构相同，超参数减少，便于模型移植。
ShuffleNetv1[26]是一个可用于移动设备的轻量级卷积神经网络，在ResNet的基础上改进，其中Channel shuffle解决了多个分组卷积(Group convolution)叠加出现的边界效应，借鉴了MobileNet的深度可分离卷积、逐点组卷积(Pointwise group convolution)和深度可分离卷积(Depthwise separable convolution)，主要减少了计算量，可以在有限的算力之下实现最优精度，更加适用于高时效、低延迟的场景。ShuffleNetv2[27]中提出了FLOPs不能作为衡量目标检测模型运行速度的标准，ShuffleNetv1大量使用1×1组卷积增加了内存访问成本(Memory access cost)，同时使用过多的组会造成网络碎片化，减少模型的并行度。为了弥补ShuffleNetv1的缺陷引入了通道分割(Channel split)，将1×1组卷积修改为1×1卷积，将2个分支的输出连接在一起，紧接着将2个分支连接的结果进行Channle shuffle，以保证2个分支信息交流，模型基本单元结构如图5a所示。下采样模块不再用Channel split，2个分支都直接复制1份输入，分支均采用步长为2的卷积层进行下采样，接着增加1个 1×1 卷积进行通道数的调整，最后将2个分支进行合并，特征图空间大小减半、通道数翻倍，空间下采样单元结果如图5b所示。
![]()
图 5 ShuffleNetv2的基本单元(a)和下采样单元(b)Figure 5. Basic unit (a) and down sampling unit (b) of ShuffleNetv21.4.2 基于ShuffleNetv2的改进模型
作为轻量级网络，ShuffleNetv2复杂度比ResNet和ResNeXt都低，且适用场景更广，所以选取ShuffleNetv2作为本研究的基础网络模型。ShuffleNetv2的整体结构如图6所示。
本研究在ShuffleNetv2基础上，对尾部结构进行了修改，即将网络最后全连接层1000个类别改为4类，剔除不需要的参数，减少参数量，选取较大的学习率，能够更好地学习拟合本研究的数据，达到更好的效果。训练过程中运用迁移学习的方法，即以在ImageNet数据集训练得到的权重作为本研究的预训练权重，该项操作有助于新模型的训练，避免随机初始化权重导致训练的效果不好、训练时间过长，通过预训练能够使得网络更好地拟合。
此外，CBAM(Convolutional block attention module)是一种轻量型的模块，作为即插即用的模块能够集成到现有的网络中，且能够节约参数和算力，在ResNet和MobileNet中均取得了不错的效果，其中的空间注意力模块(Spatial attention module，SAM)[28]使得网络主要集中学习特征图中的位置信息。因此，本研究在ShuffleNetv2网络结构中加入SAM。SAM的运行流程是将输入的特征图进行基于通道的全局最大池化(Global max pooling)和全局平均池化(Global average pooling)，然后将得到的2个结果进行基于通道的连接(Conact)操作。基于通道的操作可以有效地突出信息区域，经过1个卷积操作，降维后得到1个通道的特征图，再经过Sigmod函数生成空间特征图(Spatial attention feature)，最后将得到的特征图和该模块的输入特征图进行相乘，得到最终生成的特征，如图7所示。
本研究加入的SAM的计算公式如下：
$$\begin{split} {\rm{Ms}}(F)& = \sigma \left\{ {{f^{7 \times 7}}\left[{\rm{AvgPool}}(F);{\rm{MaxPool}}(F)\right]} \right\}\\ & = \sigma \left[{f^{7\times 7}} \left( {{F_{{\rm{avg }}}^{\rm{s}}};F_{{\rm{max }}}^{\rm{s}}} \right)\right]\end{split}$$ (1) 式中：Ms(F)为计算得到的特征；σ为Sigmoid的函数；f 7×7 表示卷积核尺寸是7×7；
${F_{{\rm{avg }}}^{\rm{s}}} $ 、$F_{{\rm{max }}}^{\rm{s}}$ 表示最大池化和平均池化操作。在ShuffleNetv2网络结构中加入SAM，本研究进行多次试验得到加入SAM的最佳位置，即在阶段2后加入SAM。改进后的网络模块如图8所示。
2. 结果与分析
2.1 试验环境
本试验使用的服务器处理器为Intel Core i7-10700 CPU @ 2.90GHZx 16，内存16 GB，显卡为GeForce GTX 3070。软件方面，系统为Ubuntu16.04，CUDA 版本为11.3。基于Python编程语言使用PyTorch框架进行深度学习模型的搭建。
2.2 评价指标
本研究所用的测试指标主要为准确率和损失值。准确率广泛地应用于分类的评估中，准确率是指预测正确的样本数占全部总样本数的比例，针对的是所有样本的统计量，计算较为简单，在分类问题中是最直观的评价指标。
$$ {\rm{Acc}} = \dfrac{{\rm{TP}}+{\rm{TN}}}{{\rm{TP}}+{\rm{FP}}+{\rm{FN}}+{\rm{TN}}} , $$ (2) 式中，Acc表示准确率，TP表示判断正确的正例数，TN 表示判断正确的负例数，FP 表示判断错误的正例数，FN 表示判断错误的负例数。
训练模型时使用交叉熵损失函数来训练优化本研究模型参数，损失函数的公式如下：
$$ \mathrm{L}\mathrm{o}\mathrm{s}\mathrm{s}=-\displaystyle\sum\limits _{i=0}^{{F}_\mathrm{c}}{{y}}_{i}\mathrm{log}{{y}}_{\mathrm{p}\mathrm{r}\mathrm{e}\_i} \text{，} $$ (3) 式中，Loss表示损失值，
${{y}}_{\mathrm{p}\mathrm{r}\mathrm{e}\_i}$ 表示输入数据的类是第i类的概率，i 表示样本标签即真实值，Fc表示类的总数量。2.3 模型训练
针对荔枝冠层花期遥感图像，进行花期评估，本研究进行多项对比试验，将改进后的ShuffleNetv2与深度学习卷积神经网络ResNet18、ResNet34、ResNet50和ResNeXt50进行多方面的对比。均采用迁移学习的方法，加载预训练的权重，模型训练迭代次数为50次，批处理尺寸为32，损失函数为对数交叉熵损失(Cross entropy loss)，优化算法为随机梯度下降(Stochastic gradient descent，SGD) 。动态调整学习率，动量因子为 0. 9，其中ShuffleNetv2属于轻量型网络，相对本研究其他网络较小，ShuffleNetv2和改进后的ShuffleNetv2初始学习率设定为0.1，能够更快地拟合本研究数据，拟合效果更好，其余对比模型较大，学习率过大会造成网络无法收敛，在最优值附近振荡，初始学习率设定为0.001。
2.4 对比分析
图9为本文改进的ShuffleNet v2模型与其他5种主流深度学习卷积神经网络模型的损失值变化曲线的比较。由图9a可以看出，ShuffleNetv2在训练集上损失值降低得较快。损失值可以很好地反映模型的拟合能力，损失越低表示模型拟合能力越强，其中，ResNet18层数过浅导致学习特征不够，模型拟合效果欠佳，而ShuffleNetv2和改进后的ShuffleNetv2最终达到较小的损失值，表明ShuffleNetv2对荔枝花期冠层的特征提取能力强，模型表征充分，和ResNet50、ResNeXt50相比，损失值相差很小。由图9b可以看出，ResNet50损失值最终在0.55左右，ResNeXt50损失值维持在0.51左右，改进后的ShuffleNetv2在验证集上损失值波动很小且达到最低，维持在0.49左右，证明相比其他模型该模型泛化能力强。
由图10a可以看出，在训练集上，ShuffleNetv2和改进后的ShuffleNetv2准确率均高于ResNet，与ResNeXt准确率十分接近。由图10b可以看出，在验证集上，ResNet50的准确率达到0.78左右，ResNeXt50的准确率达到0.80左右，ShuffleNetv2和改进后的ShuffleNetv2准确率分别达到0.83和0.84左右，效果比ResNet、ResNeXt好，表明模型学习效果较好。
表1为不同模型对荔枝开花率评估的性能对比。由表1可知，随着层数的增加，ResNet训练时间延长；ShuffleNetv2的训练时间相比于ResNet18较少。与ResNet、ResNeXt网络相比，ShuffleNetv2的训练速度大幅提高，将有利于后期补充新的数据集，更加迅速地拟合。在单张验证的时间上，ShuffleNetv2和改进的ShuffleNetv2都要比ResNet34短，体现出一个轻量型网络在验证速度上的优势。在测试集上，改进后的ShuffleNetv2模型能够准确地分类出不同花期的冠层，可以准确分析出所处花期为哪一类，尤其是第1类(开花率10%~20%)和第2类(开花率50%~60%)识别效果相比ResNet、ResNeXt有大幅提升。ResNet、ResNeXt模型效果不如改进后的ShuffleNetv2，主要识别错误集中在第2类上，第2类上ResNet34识别准确率只能达到52%，误判率较高，因为第1类和第2类荔枝冠层的花期特征十分接近，所以模型识别困难，本研究关键难点也在于此，改进后的ShuffleNetv2在保持第0类(开花率0%)和第3类(开花率80%及以上)高准确率的基础上，大大提升了第2类的准确率，第2类识别准确率可以达到68%，相比于ResNet、ResNeXt高出约8个百分点，模型整体的识别准确率提升，可以为后期果园精准管控提供参考。
表 1 不同模型的性能对比Table 1. Performance comparison of different models模型 Model t训练/s Training time t验证/ms Validation time 分类准 确率/% Classification accuracy 2类准 确率/% Category 2 accuracy ResNet18 120.4 7.352 83.33 60.0 ResNet34 135.8 7.812 83.33 52.0 ResNet50 175.6 8.802 85.96 60.0 ResNeXt50 199.2 9.116 87.01 60.0 ShuffleNetv2 114.4 7.529 86.84 64.0 改进的 ShuffleNetv2 109.0 7.507 88.60 68.0 Improved ShuffleNetv2 图11为不同网络模型用于评估荔枝开花率的参数量比较。由图11可知，随着层数的增加，ResNet参数量成正比增加，数量级都很大；ResNeXt50相对于ResNet50，参数量有所减少，因为ResNeXt50引入了分组卷积，使得参数减少的情况下效果更好，减少了超参数的数量。ShuffleNetv2的参数量相对ResNet和ResNeXt要少很多，比最少的ResNet18的参数少10866148，而识别效果要好很多。ShuffleNetv2和ResNeXt50相比，参数量相差达到23167106，效果略差一点，改进后的ShuffleNetv2参数量与ShuffleNetv2几乎没有区别，但是识别效果有了较大的提升。荔枝冠层模型分类的最终效果如图12所示。
![]()
图 12 改进的ShuffleNetv2模型荔枝冠层分类效果Figure 12. Effect of improved shufflenetv2 model on litchi canopy classification3. 讨论与结论
3.1 讨论
本研究的荔枝花期遥感图像数据中，单棵荔枝冠层图像分辨率均较低，经过与ResNet、ResNeXt大量的对比试验发现，较深的卷积神经网络应用于本试验，会出现冠层特征消失、过拟合等问题，网络层数过多且没有明显地提升准确率，训练和识别速度比较慢，所以ResNet、ResNeXt网络不适用于本研究。在同样的参数设定下，轻量级网络ShuffleNetv2对本研究数据有着较好的拟合效果，未优化时准确率已达到86.84%，与ResNeXt50相差0.17个百分点，且ShuffleNetv2的训练时间、验证时间有着绝对的优势，参数量最少。试验过程中发现，荔枝冠层花期第2类(开花率50%~60%)误判率最高，识别精度难以有效提升，相比于ResNet、ResNeXt一贯地依靠增加网络层数来提升精度，ShuffleNetv2的基础单元中加入了Channle split，能够巧妙地有效学习第2类花期冠层的特征，且能够保持冠层特征不会消失。
将ShuffleNetv2进行主体结构优化，在阶段2后加入SAM模块，使得网络学习过程中注意面增宽，更加注意输入图像的位置信息，能够对荔枝冠层遥感图像进行更加精准的学习，更加关注荔枝冠层花的位置信息，从而能够更好地判断和识别花期所属的类别，进一步提升网络的泛化性和识别准确率，改进后的ShuffleNetv2分类准确率相比未优化的ShuffleNetv2提升1.76个百分点，相比ResNeXt50提升1.59个百分点，最终能达到88.60%，且参数量和未优化的ShuffleNetv2相差甚微。在增加空间注意力模块和优化主干网络以后，训练时间反而会降低，因为模型更加关注荔枝花的部分，重点学习荔枝冠层的各项特征，有效剔除掉与花期无关的信息，所以训练速度会变快。本研究可以应用于大型荔枝园的花期评估，可以在很大程度上节省人力、物力和财力。方法上充分运用了迁移学习，节省了训练时间，能够充分地学习到冠层花期的特征。
3.2 结论
本文使用无人机采集荔枝冠层的可见光遥感图像，将遥感图像以实例分割的方法进行有效的冠层分割，得到单独每棵荔枝冠层的图像，将图像预处理后送入分类网络中。由于荔枝花期特征较为接近，识别较为困难，通过无人机获取遥感图像，在独特的视角下进行花期的分类。加入空间注意力模块，改进后的ShuffleNetv2能够挖掘学习更为细节的荔枝冠层花期信息，不仅具有较高的识别精度，而且有着很快的训练和识别速度，对荔枝花期的评估有很大的优势。
本研究对荔枝冠层的遥感图像进行处理，实现了大型荔枝园的花期评估，同时降低了计算成本，具有很强的实时性，为荔枝花期的花量提供数据支撑，在荔枝研究领域有着很高的实用价值。同时本研究结果可以更好地被应用与推广，可以在不显著要求算力的情况下进行模型的推理，本研究的方法为其他种类果树的花期评估提供了一个新思路，利用轻量级网络对无人机拍摄的冠层图像进行花期的研究与评估。未来可以加入不同地区的果树冠层数据，提升数据的多样性，并结合地面摄像头数据，进一步提升模型识别效果。将模型精简后部署到边缘智能终端，节约算力资源，可以在无人机拍摄过程中，更快地实现花期的判断。
-
![]()
图 1 生物炭和生物炭基肥对火力楠生长的影响
CK：不施肥，CF：复合肥，RB：稻壳生物炭，WB：木屑生物炭，RBF：稻壳生物炭基肥，WBF：木屑生物炭基肥；各图中，柱子上方的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s 法)。
Figure 1. Effects of biochar and biochar-based fertilizer on the growth of Michelia macclurei
CK: No fertilization, CF: Compound fertilizer, RB: Rice husk biochar, WB: Wood chip biochar, RBF: Rice husk biochar-based fertilizer, WBF: Wood chip biochar-based fertilizer; In each figure, different lowercase letters on columns indicate significant differences among treatments (P<0.05, Duncan’s method).
![]()
图 2 生物炭和生物炭基肥对火力楠生物量的影响
CK：不施肥，CF：复合肥，RB：稻壳生物炭，WB：木屑生物炭，RBF：稻壳生物炭基肥，WBF：木屑生物炭基肥；各图中，柱子上方的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s 法)。
Figure 2. Effects of biochar and biochar-based fertilizer on the biomass of Michelia macclurei
CK: No fertilization, CF: Compound fertilizer, RB: Rice husk biochar, WB: Wood chip biochar, RBF: Rice husk biochar-based fertilizer, WBF: Wood chip biochar-based fertilizer; In each figure, different lowercase letters on columns indicate significant differences among treatments (P<0.05, Duncan’s method).
![]()
图 3 生物炭和生物炭基肥对火力楠Dickson质量指数的影响
CK：不施肥，CF：复合肥，RB：稻壳生物炭，WB：木屑生物炭，RBF：稻壳生物炭基肥，WBF：木屑生物炭基肥；柱子上方的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s 法)。
Figure 3. Effects of biochar and biochar-based fertilizer on the Dickson quality index of Michelia macclurei
CK: No fertilization, CF: Compound fertilizer, RB: Rice husk biochar, WB: Wood chip biochar, RBF: Rice husk biochar-based fertilizer, WBF: Wood chip biochar-based fertilizer; Different lowercase letters on columns indicate significant differences among treatments (P<0.05, Duncan’s method).
![]()
图 4 生物炭和生物炭基肥对火力楠叶面积和SPAD值的影响
CK：不施肥，CF：复合肥，RB：稻壳生物炭，WB：木屑生物炭，RBF：稻壳生物炭基肥，WBF：木屑生物炭基肥；各图中，柱子上方的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s 法)。
Figure 4. Effects of biochar and biochar-based fertilizer on the leaf area and SPAD value of Michelia macclurei
CK: No fertilization, CF: Compound fertilizer, RB: Rice husk biochar, WB: Wood chip biochar, RBF: Rice husk biochar-based fertilizer, WBF: Wood chip biochar-based fertilizer; In each figure, different lowercase letters on columns indicate significant differences among treatments (P<0.05, Duncan’s method).
![]()
图 5 生物炭和生物炭基肥对火力楠各部位养分含量的影响
CK：不施肥，CF：复合肥，RB：稻壳生物炭，WB：木屑生物炭，RBF：稻壳生物炭基肥，WBF：木屑生物炭基肥；各图中，相同植株部位柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s 法)。
Figure 5. Effects of biochar and biochar-based fertilizer on the nutrient contents of various parts of Michelia macclurei
CK: No fertilization, CF: Compound fertilizer, RB: Rice husk biochar, WB: Wood chip biochar, RBF: Rice husk biochar-based fertilizer, WBF: Wood chip biochar-based fertilizer; In each figure, different lowercase letters on the columns of the same plant part indicate significant differences among treatments (P<0.05, Duncan’s method).
![]()
图 6 生物炭和生物炭基肥对火力楠各部位养分累积量的影响
CK：不施肥，CF：复合肥，RB：稻壳生物炭，WB：木屑生物炭，RBF：稻壳生物炭基肥，WBF：木屑生物炭基肥；各图中，相同植株部位柱子上方的不同小写字母表示处理间差异显著，不同大写字母表示全株养分累积量在处理间差异显著 (P<0.05，Duncan’s法)。
Figure 6. Effects of biochar and biochar-based fertilizer on the nutrient accumulations in various parts of Michelia macclurei
CK: No fertilization, CF: Compound fertilizer, RB: Rice husk biochar, WB: Wood chip biochar, RBF: Rice husk biochar-based fertilizer, WBF: Wood chip biochar-based fertilizer; In each figure, different lowercase letters on the columns of the same plant part indicate significant differences among treatments, different capital letters indicate significant differences in nutrient accumulation of the whole plant among treatments (P<0.05, Duncan’s method).
表 1 稻壳生物炭(RB)和木屑生物炭(WB)的基本性质
Table 1 Basic properties of rice husk biochar (RB) and wood chip biochar (WB)
材料
Material容重/(g·cm−3)
Bulk
density毛管持水量/(g·kg−1)
Capillary
capacity总孔隙度/%
Total
porosity毛管孔隙度/%
Capillary
porosity非毛管孔隙度/%
Non-capillary
porosity通气孔隙度/%
Air
porositypH 电导率/(μS·cm−1)
Electrical
conductivityRB 0.50 971.06 81.31 48.08 33.23 72.39 10.17 2320 WB 0.34 1777.55 87.29 59.85 27.44 76.93 9.24 406 
下载: 导出CSV
材料
Materialw/(g·kg−1) w/(mg·kg−1) 有机质
Organic matter全氮
Total N全磷
Total P全钾
Total K碱解氮
Alkali-hydrolyzable N速效磷
Available P速效钾
Available KRB 311.40 12.53 3.28 18.11 56.00 387.36 2761.22 WB 234.19 7.20 1.53 6.42 28.00 13.68 426.12
下载: 导出CSV
表 2 生物炭和生物炭基肥对土壤物理性质的影响1)
Table 2 Effects of biochar and biochar-based fertilizer on soil physical properties
处理
Treatment容重/(g·cm−3)
Bulk density毛管持水量/(g·kg−1)
Capillary capacity总孔隙度/%
Total porosityCK 1.14±0.02a 335.16±21.29a 57.01±0.84b CF 1.10±0.04a 362.78±11.20a 58.31±1.40b RB 1.04±0.06ab 337.04±20.07a 60.58±2.29ab WB 0.93±0.03b 358.75±8.59a 65.09±1.09a RBF 1.01±0.05ab 352.51±16.61a 61.85±1.91ab WBF 0.91±0.06b 356.12±23.39a 65.65±2.44a 1) CK：不施肥，CF：复合肥，RB：稻壳生物炭，WB：木屑生物炭，RBF：稻壳生物炭基肥，WBF：木屑生物炭基肥；同列数据后不同小写字母表示处理间差异显著 (P<0.05，Duncan’s 法)。
1) CK: No fertilization, CF: Compound fertilizer, RB: Rice husk biochar, WB: Wood chip biochar, RBF: Rice husk biochar-based fertilizer, WBF: Wood chip biochar-based fertilizer; Different lowercase letters in the same column indicate significant differences among treatments (P<0.05, Duncan’s method).
下载: 导出CSV
表 3 生物炭和生物炭基肥对土壤化学性质的影响1)
Table 3 Effects of biochar and biochar-based fertilizer on soil chemical properties
处理
TreatmentpH w/(g·kg−1) w/(mg·kg−1) 有机质
Organic matter全氮
Total N全磷
Total P全钾
Total K碱解氮
Available N速效磷
Available P速效钾
Available KCK 4.33±0.04ab 25.35±6.92a 1.99±0.16a 0.17±0.01bc 14.02±0.56b 109.08±9.55a 18.70±5.02a 101.20±15.12b CF 4.21±0.08b 38.19±3.84a 2.09±0.18a 0.21±0.01abc 15.48±1.85ab 118.42±7.65a 29.22±15.93a 101.15±6.43b RB 4.41±0.05a 31.30±6.53a 2.04±0.22a 0.23±0.03abc 16.26±2.05ab 106.75±8.81a 42.27±14.29a 104.52±9.32b WB 4.40±0.03a 38.24±3.41a 2.35±0.23a 0.30±0.04a 14.78±0.80ab 128.92±5.09a 34.64±15.51a 100.76±7.26b RBF 4.47±0.06a 37.42±4.11a 2.17±0.10a 0.25±0.02ab 18.63±1.13ab 110.25±9.26a 54.08±18.19a 201.99±46.49a WBF 4.31±0.02ab 40.16±2.43a 2.17±0.10a 0.14±0.06c 19.43±1.97a 127.75±6.15a 39.54±11.11a 105.20±4.55b 1) CK：不施肥，CF：复合肥，RB：稻壳生物炭，WB：木屑生物炭，RBF：稻壳生物炭基肥，WBF：木屑生物炭基肥；同列数据后不同小写字母表示处理间差异显著 (P<0.05，Duncan’s 法)。
1) CK: No fertilization, CF: Compound fertilizer, RB: Rice husk biochar, WB: Wood chip biochar, RBF: Rice husk biochar-based fertilizer, WBF: Wood chip biochar-based fertilizer; Different lowercase letters in the same column indicate significant differences among treatments (P<0.05, Duncan’s method).
下载: 导出CSV
表 4 生物炭和生物炭基肥对火力楠生长和土壤性质的隶属函数值1)
Table 4 Membership function values of biochar and biochar-based fertilizer on the growth and soil properties of Michelia macclurei
指标 Index CK CF RB WB RBF WBF 株高增长量 Plant height increment 0.05 0.40 0.30 0.19 0.43 0.81 地径增长量 Ground diameter increment 0.13 0.58 0.49 0.41 0.60 0.87 冠幅增长量 Crown width increment 0.14 0.61 0.42 0.38 0.52 0.96 总生物量 Total biomass 0.03 0.49 0.24 0.24 0.53 0.75 Dickson质量指数 Dickson quality index 0.06 0.66 0.32 0.36 0.61 0.84 叶面积 Leaf area 0.07 0.45 0.49 0.69 0.55 0.85 SPAD值 SPAD value 0.09 0.55 0.66 0.62 0.94 0.65 氮含量 N content 0.03 0.66 0.38 0.45 0.38 0.60 磷含量 P content 0.43 0.46 0.38 0.50 0.63 0.63 钾含量 K content 0.28 0.93 0.48 0.33 0.74 0.77 氮累积量 Accumulation of N 0.02 0.45 0.21 0.22 0.43 0.66 磷累积量 Accumulation of P 0.05 0.48 0.23 0.27 0.62 0.83 钾累积量 Accumulation of K 0.06 0.63 0.27 0.23 0.61 0.82 容重 Bulk density 0.11 0.20 0.36 0.67 0.45 0.71 总孔隙度 Total porosity 0.11 0.20 0.36 0.67 0.45 0.71 pH 0.53 0.29 0.68 0.67 0.80 0.50 有机质含量 Organic matter content 0.41 0.80 0.59 0.80 0.77 0.86 碱解氮含量 Available N content 0.36 0.56 0.31 0.79 0.38 0.77 速效磷含量 Available P content 0.15 0.30 0.49 0.38 0.67 0.45 速效钾含量 Available K content 0.13 0.13 0.15 0.13 0.58 0.15 平均值 Average 0.16 0.49 0.39 0.45 0.58 0.71 综合排序 Comprehensive ranking 6 3 5 4 2 1 1) CK：不施肥，CF：复合肥，RB：稻壳生物炭，WB：木屑生物炭，RBF：稻壳生物炭基肥，WBF：木屑生物炭基肥。
1) CK: No fertilization, CF: Compound fertilizer, RB: Rice husk biochar, WB: Wood chip biochar, RBF: Rice husk biochar-based fertilizer, WBF: Wood chip biochar-based fertilizer.
下载: 导出CSV
-
[1] 中国科学院中国植物志编辑委员会. 中国植物志 第30卷[M]. 北京: 科学出版社, 1974: 443. [2] 李清莹, 文珊娜, 姜清彬, 等. 不同营养元素缺乏对火力楠幼苗生长的影响[J]. 生态学杂志, 2017, 36(3): 664-669. [3] MARCIŃCZYK M, OLESZCZUK P. Biochar and engineered biochar as slow- and controlled-release fertilizers[J]. Journal of Cleaner Production, 2022, 339: 130685. doi: 10.1016/j.jclepro.2022.130685
[4] 丁乐, 杨弋, 倪辉, 等. 生物炭配施对芳樟精油产量及品质的影响[J]. 中南林业科技大学学报, 2022, 42(5): 91-100. [5] TAO J C, WAN C X, LENG J J, et al. Effects of biochar coupled with chemical and organic fertilizer application on physicochemical properties and in vitro digestibility of common buckwheat (Fagopyrum esculentum Moench) starch[J]. International Journal of Biological Macromolecules, 2023, 246: 125591. doi: 10.1016/j.ijbiomac.2023.125591
[6] 陈温福, 张伟明, 孟军. 生物炭与农业环境研究回顾与展望[J]. 农业环境科学学报, 2014, 33(5): 821-828. [7] MELO L C A, LEHMANN J, CARNEIRO J S D S, et al. Biochar-based fertilizer effects on crop productivity: A meta-analysis[J]. Plant and Soil, 2022, 472: 45-58. doi: 10.1007/s11104-021-05276-2
[8] YADAV R K, YADAV M R, KUMAR R, et al. Role of biochar in mitigation of climate change through carbon sequestration[J]. International Journal of Current Microbiology and Applied Sciences, 2017, 6(4): 859-866. doi: 10.20546/ijcmas.2017.604.107
[9] 陈志琪, 张海娜, 鲁向晖, 等. 生物炭配施化肥对稀土尾砂基质中樟树幼苗生长及光合作用的影响[J]. 南昌工程学院学报, 2023, 42(3): 34-39. doi: 10.3969/j.issn.1006-4869.2023.03.006 [10] EL-SYED N M M, HELMY A M, FOUDA S E E, et al. Biochar with organic and inorganic fertilizers improves defenses, nitrogen use efficiency, and yield of maize plants subjected to water deficit in an alkaline soil[J]. Sustainability, 2023, 15(16): 12223. doi: 10.3390/su151612223
[11] YE L, CAMPS ARBESTAIN M, SHEN Q, et al. Biochar effects on crop yields with and without fertilizer: A meta-analysis of field studies using separate controls[J]. Soil Use and Management, 2020, 36(1): 2-18. doi: 10.1111/sum.12546
[12] 黄庆, 刘忠珍, 朱根发, 等. 生物质炭基肥料及作物施用技术研究进展[J]. 广东农业科学, 2021, 48(1): 26-34. [13] 吴志庄, 王道金, 厉月桥, 等. 施用生物炭肥对黄连木生长及光合特性的影响[J]. 生态环境学报, 2015, 24(6): 992-997. [14] LIU M, LINNA C, MA S, et al. Biochar combined with organic and inorganic fertilizers promoted the rapeseed nutrient uptake and improved the purple soil quality[J]. Frontiers in Nutrition, 2022, 9: 997151.
[15] ADEKIYA A O, ADEBIYI O V, IBABA A L, et al. Effects of wood biochar and potassium fertilizer on soil properties, growth and yield of sweet potato (Ipomea batata)[J]. Heliyon, 2022, 8(11): e11728. doi: 10.1016/j.heliyon.2022.e11728
[16] SCHULZ H, GLASER B. Effects of biochar compared to organic and inorganic fertilizers on soil quality and plant growth in a greenhouse experiment[J]. Journal of Plant Nutrition and Soil Science, 2012, 175(3): 410-422. doi: 10.1002/jpln.201100143
[17] 刘士玲, 贾宏炎, 陈琳, 等. 容器规格和添加生物炭的基质配方对西南桦幼苗生长的影响[J]. 生态学杂志, 2019, 38(9): 2875-2882. [18] 国家林业局. 火力楠培育技术规程: LY/T 2461—2015[S]. 北京: 中国标准出版社, 2015. [19] 鲍士旦. 土壤农化分析[M]. 3版. 北京: 中国农业出版社, 2000: 39-110. [20] 陆洁, 冯嘉仪, 盛晗, 等. 生物炭对污泥混合基质性质和植物生长的影响[J]. 农业环境科学学报, 2023, 42(5): 1082-1090. doi: 10.11654/jaes.2022-0936 [21] 彭钟通, 林佳娜, 刘丽, 等. 水氮耦合对辣木生长特性的影响[J]. 生态学杂志, 2021, 40(2): 363-372. [22] FALOYE O T, AJAYI A E, ALATISE M O, et al. Maize growth and yield modelling using AquaCrop under deficit irrigation with sole and combined application of biochar and inorganic fertiliser[J]. Journal of Soil Science and Plant Nutrition, 2020, 20(4): 2440-2453. doi: 10.1007/s42729-020-00310-1
[23] METE F Z, MIA S, DIJKSTRA F A, et al. Synergistic effects of biochar and NPK fertilizer on soybean yield in an alkaline soil[J]. Pedosphere, 2015, 25(5): 713-719. doi: 10.1016/S1002-0160(15)30052-7
[24] FALOYE O T, ALATISE M O, AJAYI A E, et al. Effects of biochar and inorganic fertiliser applications on growth, yield and water use efficiency of maize under deficit irrigation[J]. Agricultural Water Management, 2019, 217: 165-178. doi: 10.1016/j.agwat.2019.02.044
[25] 文慧宝, 杜双江, 程贵庆, 等. 生物炭与有机肥配施对白菜生长品质、土壤养分及酶活性的影响[J]. 江苏农业科学, 2023, 51(20): 224-230. [26] 刘大翔, 高贤, 许亚坤, 等. 施加生物炭对植被混凝土理化性质及狗牙根生物量的影响[J]. 应用基础与工程科学学报, 2021, 29(1): 1-14. [27] 王高飞, 刘鸿雁, 邢丹, 等. 育苗基质添加辣椒秆生物炭对辣椒幼苗生长与养分的影响[J]. 西南农业学报, 2022, 35(3): 543-549. [28] YE Z X, ZHANG L M, HUANG Q Y, et al. Development of a carbon-based slow release fertilizer treated by bio-oil coating and study on its feedback effect on farmland application[J]. Journal of Cleaner Production, 2019, 239: 118085. doi: 10.1016/j.jclepro.2019.118085
[29] SAHA A, BASAK B B, GAJBHIYE N A, et al. Sustainable fertilization through co-application of biochar and chemical fertilizers improves yield, quality of Andrographis paniculata and soil health[J]. Industrial Crops and Products, 2019, 140: 111607. doi: 10.1016/j.indcrop.2019.111607
[30] BAI S H, OMIDVAR N, GALLART M, et al. Combined effects of biochar and fertilizer applications on yield: A review and meta-analysis[J]. Science of the Total Environment, 2022, 808: 152073. doi: 10.1016/j.scitotenv.2021.152073
[31] HE K, XU Y, HE G, et al. Combined application of acidic biochar and fertilizer synergistically enhances Miscanthus productivity in coastal saline-alkaline soil[J]. Science of the Total Environment, 2023, 893: 164811. doi: 10.1016/j.scitotenv.2023.164811
[32] 邓洪涛, 肖琳, 冯为迅, 等. 施用生物炭和复合肥对桢楠生长及养分吸收的影响[J]. 广西科学, 2023, 30(5): 869-875. [33] 吴志庄, 高贵宾, 欧建德, 等. 生物炭肥对毛竹林下三叶青叶绿素含量、光合与荧光特性的影响[J]. 西北林学院学报, 2017, 32(5): 59-63. [34] XIANG Y Z, DENG Q, DUAN H L, et al. Effects of biochar application on root traits: A meta-analysis[J]. GCB Bioenergy, 2017, 9(10): 1563-1572. doi: 10.1111/gcbb.12449
[35] 王耀锋, 刘玉学, 吕豪豪, 等. 水洗生物炭配施化肥对水稻产量及养分吸收的影响[J]. 植物营养与肥料学报, 2015, 21(4): 1049-1055. [36] 张令. 生物炭与氮肥配施对巨尾桉幼苗生长生理及养分吸收利用的影响[D]. 长沙: 中南林业科技大学, 2017. [37] 夏浩, 张梦阳, 刘波, 等. 生物炭对作物氮肥利用率影响的整合分析[J]. 华中农业大学学报, 2021, 40(3): 177-186. [38] 凌宏伟, 李小英, 周月惠, 等. 不同生物炭混施对土壤养分和蚕豆生长的影响[J]. 江苏农业科学, 2019, 47(17): 110-112. [39] 李青山, 王德权, 杜传印, 等. 有机无机肥与生物炭配施对烤烟生长发育和烟叶质量的影响[J]. 土壤通报, 2021, 52(6): 1393-1401. [40] 袁金华, 徐仁扣. 生物质炭的性质及其对土壤环境功能影响的研究进展[J]. 生态环境学报, 2011, 20(4): 779-785. [41] 胡京钰, 杨红军, 刘大军, 等. 酒糟生物炭与化肥配施对土壤理化特性及作物产量的影响[J]. 植物营养与肥料学报, 2022, 28(9): 1664-1672. [42] 丛铭, 张梦阳, 夏浩, 等. 施用生物炭对红壤中不同形态钾含量及小白菜生长的影响[J]. 华中农业大学学报, 2020, 39(4): 22-28. [43] 胡华英, 殷丹阳, 曹升, 等. 生物炭对杉木人工林土壤养分、酶活性及细菌性质的影响[J]. 生态学报, 2019, 39(11): 4138-4148. [44] 尧芳, 任天宝, 徐敏, 等. 生物质炭改善土壤矿质营养吸收的研究进展及作用机制分析[J]. 江苏农业科学, 2020, 48(10): 46-51. [45] 严陶韬, 丁子菊, 朱倩, 等. 生物质炭对黄棕壤理化性质及龙脑樟幼苗生长的影响[J]. 土壤, 2018, 50(4): 681-686.
计量
- 文章访问数: 241
- HTML全文浏览量: 70
- PDF下载量: 24
