Effects of biochar amendment with phosphorus fertilizer on organic phosphorus transformation and soybean phosphorus uptake in red soil
-
摘要:目的
研究不同磷水平下添加生物炭对红壤磷周转与作物磷吸收的影响及潜在机制，为生物炭在减磷增效和促进大豆磷吸收的应用提供科学依据。
方法以酸性红壤为供试土壤，研究4种磷肥水平(0、30、60、90 kg·hm−2)下配施4%(w)稻秆生物炭对土壤磷组分含量、磷酸酶活性、微生物量磷含量和大豆磷吸收的影响及相关关系。
结果磷水平、生物炭及二者的交互作用对大豆磷吸收、土壤磷组分及周转有显著影响。添加生物炭促进大豆磷吸收、提升磷利用效率，且在30 kg·hm−2磷施用条件下作用更显著。不同磷水平下配施生物炭均提高土壤无机磷含量，尤其NH4F−Pi和HCl−Pi含量增加最显著；不施磷肥条件下添加生物炭使NaOH提取态有机磷降低48.5%。此外，不施磷肥和施30 kg·hm−2磷条件下添加生物炭分别使土壤碱性磷酸酶活性增加71.7%和46.0% (P<0.05)，土壤微生物量磷增加458.3%和84.0% (P<0.05)。相关分析表明，植株磷吸收量与土壤NH4F−Pi、HCl−Pi、碱性磷酸酶活性和微生物量磷含量呈显著正相关。
结论不施磷肥或施30 kg·hm−2磷条件下配施生物炭对植物磷吸收和磷利用效率的促进更显著，同时促进土壤活性有机磷的生物学转化。
Abstract:ObjectiveTo investigate the effect of biochar addition at different phosphorus (P) levels on P transformation and crop P uptake in red soil, as well as the underlying processes and mechanisms, and provide a scientific basis for the application of biochar in reducing P input and increasing efficiency, as well as promoting soybean P uptake.
MethodAn acidic red soil was used as the test soil to explore the effects of 4%(w) rice straw biochar under four phosphate fertilizer levels (0, 30, 60 and 90 kg·hm−2) on contents of soil components, phosphatase activity, microbial biomass P content and soybean P uptake and their correlations.
ResultP level, biochar and their interaction had significantl effects on soybean P absorption, soil P composition and turnover. The biochar addition promoted soybean P uptake and P use efficiency, with a more significant effect at 30 kg·hm−2 P treatment. Biochar application at different P levels increased soil inorganic P content, especially the NH4F-Pi, HCl-Pi. Without phosphate fertilizer, adding biochar reduced NaOH-extractable organic P by 48.5%. Furthermore, adding biochar with 0 and 30 kg·hm−2 P treatment increased soil alkaline phosphatase activity by 71.7% and 46.0% respectively, and soil microbial biomass P content by 458.3% and 84.0% respectively.
ConclusionBiochar addition under low P conditions significantly promotes plant P uptake and P use efficiency, while also enhancing the biological transformation of soil labile organic P.
-
Keywords:
- Soybean /
- Red soil /
- Biochar /
- P fraction /
- Phosphatase /
- P uptake
-
随着农业机械化水平的不断提高，越来越多的新技术应用于农业。针对农作物病虫害防治，利用农业机器人实现自主作业，可以解决从事农业人口比例逐年降低、劳动者工作强度大等问题，提高作业效率和质量，降低生产成本[1-2]。自动导航是农业机器人实现自主作业重要的一步，目前应用于农业的自动导航技术主要有RTK-GPS导航、机器视觉导航、地磁导航和无线电导航等，但GPS导航容易受到农田作物遮挡影响，信号不稳定，而多传感器融合导航具有实时性好、价格低廉、适用范围广等优点，因此相机与激光雷达融合导航已经逐步成为国内外研究热点。
Søgaard等[3]通过计算图像中玉米的重心代替分割步骤，从而减少了图像处理的耗时；杨洋等[24]利用动态感兴趣区域去除图像中干扰部分，采用最小二乘法直线拟合获取了玉米导航线。常昕等[4]以激光雷达点云数据修正跟踪，解决图像中受光照、阴影和背景干扰的问题。胡丹丹等[5]以玉米收割机器人为载体，采用机器视觉获取玉米植株离散点特征，基于Hough变换进行直线拟合得到导航基准线。Marchant等[6]通过植物的颜色特征和相机标定对作物行进行Hough变换处理，提高了检测效率。梁习卉子等[7]对区域内各色域分量权值进行调整，强调绿色像素而淡化红色和蓝色像素，去除枯草、阴影对检测结果造成的干扰。
目前农业机器人在田间行走主要依靠单一传感器识别，利用相机获取图像数据，将植株绿色像素投影曲线峰值特性，通过算法拟合峰值点可以获得导航基准线。俞毓锋等[8]结合激光雷达的深度信息和图像的颜色纹理信息，构建在时序帧间的特征点匹配关系，终将定位问题转化为特征点对的加权重投影误差优化问题。但实际上存在光照和地面起伏等因素，在信息融合过程中容易出现信息不匹配问题，降低识别准确率。宋宇等[9]提出特征点处植株绿色特征列像素累加值大于左右相邻像素累加值，通过设定峰值点距离阈值来获取准确根茎定位点，最终利用拟合定位点来获取导航线。但实际上玉米田中普遍存在缺苗、一穴多株现象，并且相机拍摄角度或焦距成像不同，在设定距离阈值进行筛选根茎定位点时会出现错误，降低拟合导航线的准确率。冯娟等[10]利用特征点附近像素值梯度变化程度，通过对形态学处理后的图像进行扫描获取土壤与根茎交界线。该方法在根茎附近有杂草干扰的情况时，会在相邻根茎定位点之间产生伪特征点，对导航线拟合产生严重影响。
本文针对单一传感器无法有效识别导航线问题，提出一种基于离散因子的多传感器特征检测和数据融合方法，离散因子是由激光雷达和相机采集到的孤立的信息点集，因此构建多传感器目标识别的相对一致性和加权一致性函数，根据多传感器特征检测的权重以及相关的一致性函数，建立多传感器特征识别的数据融合支持度数学模型以获取特征点。最后利用正确的玉米特征中心点进行拟合，可以获得导航基准线，为农业机械在玉米田间行走提供导航基准。
1. 数据采集和预处理
1.1 图像采集与特征提取
视觉系统采集到的3~5叶期玉米植株图片为RGB彩色图像(图1)，图像分辨率为640像素×480像素。测试环境为阴天，环境中有石子、树枝及行人等干扰因素，用
$ 2G - B - R $ 可以提取绿色特征分量[11]，采用超绿化算法处理后，灰度图像仍为三维向量数据组，因此需要对图像进行二值化处理，可以有效减少计算量，提高效率。使用最大类间方差法[12-14](OTSU)以及形态学滤波中的开运算[15-17]进行二值化处理，处理效果如图2所示，该算法有效屏蔽了干扰因素并提取植物特征，同时对每个玉米植株连通域计算其重心位置并保存备用。1.2 点云数据采集与聚类分析
采用Velodyne Lidar VLP-16PUCK的激光雷达对视觉采集的同一目标进行数据扫描。激光雷达扫描的数据以点云形式存储，每个扫描点包含该点的三维坐标和反射率等信息。但由于复杂环境存在干扰以及设备存在精度的缺陷，原始点云数据存在伪特征点。针对原始数据中的伪特征点，采用基于半径滤波去噪的方法进行剔除，即设置规定半径范围内点云的最小数量，当给定点周围点云数量大于给定值时保留该点，反之则剔除该点，依序迭代可获得最密集的点，从而去除原始数据中的伪特征点，处理后的点云俯视图如图3所示。
预处理后的玉米植株点云数据依然存在模糊、分散等问题，因此进行部分聚类操作，导出最大密度相连的样本集合，即为最终聚类的一个类别，从而加强玉米植株的特征点云数据。采用DBSCAN[18-19]聚类算法进行点云数据处理，通过将紧密相连的样本划为一类，这样就得到一个聚类类别，通过将所有各组紧密相连的样本划为各个不同的类别，则可以获得最终所有聚类类别的结果。经过DBSCAN聚类操作处理的点云数据如图4所示，玉米植株的特征被划分得更加明显，解决了点云数据杂乱、噪音多等问题。
2. 多传感器数据融合
基于相机和激光雷达的融合导航识别算法是将预处理后点云数据的反射率信息加入到预处理后图像数据中，降低光照不足、遮挡等干扰对可通行区域提取的影响，提高识别精度，算法流程如图5所示。
2.1 相机与激光雷达联合标定
为达到数据融合的目的，首先要将相机和激光雷达数据进行空间匹配，因此需要对相机和激光雷达进行联合标定[20-21]，即得到两者之间的数据转换关系，找到同一时刻激光雷达点云数据对应图像中的像素点。
建立激光雷达坐标系
${O_{\rm{L}}} - {X_{\rm{L}}},{Y_{\rm{L}}},{Z_{\rm{L}}}$ ；相机坐标系${O_{\text{C}}} - {X_{\rm{C}}},{Y_{\rm{C}}},{Z_{\rm{C}}}$ ；图像坐标系$ O - x,y $ ；像素坐标系$ O - u,v $ ，如图6所示。假设$ {\theta _X},{\theta _Y},{\theta _Z} $ 为激光雷达坐标系相对于相机坐标系在$ x,y,z $ 坐标轴方向的转角，转动矩阵为R；$ {\boldsymbol{T}}({t_1},{t_2},{t_3}) $ 为相机坐标系到激光雷达坐标系的平移向量；$ ({u_0},{v_0}) $ 为图像坐标系原点在像素坐标系下的位置。传感器融合的关键就是将物体在激光雷达下的坐标转化为像素坐标，从而达到激光点云与图像融合的目的。因此，假设空间有一点P在激光雷达坐标系下的坐标为
$({X_{\rm{L}}},{Y_{\rm{L}}},{Z_{\rm{L}}})$ ，在相机坐标系下的坐标为$({X_{\rm{C}}},{Y_{\rm{C}}},{Z_{\rm{C}}})$ ，在图像坐标系下的坐标为$ (x,y) $ ，在像素坐标系下的坐标为$ (u,v) $ ，t1,t2,t3分别为相机到激光雷达$ x,y,z $ 方向的距离。具体过程为：将P点在激光雷达坐标系下的坐标转为相机坐标，再从相机坐标转为图像坐标，最后从图像坐标转为像素坐标。1) 从激光雷达坐标系到相机的坐标系转化为：
$$ \left[ {\begin{array}{*{20}{c}} {{X_{\rm{C}}}} \\ {{Y_{\rm{C}}}} \\ {{Z_{\rm{C}}}} \end{array}} \right] = {\boldsymbol{R}}\left[ {\begin{array}{*{20}{c}} {{X_{\rm{L}}}} \\ {{Y_{\rm{L}}}} \\ {{Z_{\rm{L}}}} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{t_1}} \\ {{t_2}} \\ {{t_3}} \end{array}} \right] ,$$ (1) 式中，
$$ \begin{split} &{\boldsymbol{R}} = {{\boldsymbol{R}}_x}{{\boldsymbol{R}}_y}{{\boldsymbol{R}}_z} = \left[ {\begin{array}{*{20}{c}} 1&0&0 \\ 0&{\cos {\theta _x}}&{ - \sin {\theta _x}} \\ 0&{\sin {\theta _x}}&{\cos {\theta _x}} \end{array}} \right]\\ &\quad\left[ {\begin{array}{*{20}{c}} {\cos {\theta _y}}&0&{\sin {\theta _y}} \\ 0&1&0 \\ {\sin {\theta _y}}&0&{\cos {\theta _y}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\cos {\theta _z}}&{ - \sin {\theta _z}}&0 \\ {\sin {\theta _z}}&{\cos {\theta _z}}&0 \\ 0&0&1 \end{array}} \right]。 \end{split} $$ (2) 2) 从相机坐标系到图像坐标系的转换为：
$$ Z_{\rm{C}}\left[ {\begin{array}{*{20}{c}} x \\ y \\ 1 \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} f&0&0 \\ 0&f&0 \\ 0&0&1 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{X_{\rm{C}}}} \\ {{Y_{\rm{C}}}} \\ {{Z_{\rm{C}}}} \end{array}} \right] ，$$ (3) 式中，
$ f $ 是相机的焦距。3) 从图像坐标到像素坐标的转换为：
$$ \left[ {\begin{array}{*{20}{c}} u \\ v \\ 1 \end{array}} \right]{\text{ = }}\left[ {\begin{array}{*{20}{c}} {\dfrac{1}{{d_x}}}&0&{{u_0}} \\ 0&{\dfrac{1}{{d_y}}}&{{v_0}} \\ 0&0&1 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} x \\ y \\ 1 \end{array}} \right]， $$ (4) 式中：
$d_x$ 为单个像素点在像素坐标x方向的长度，单位mm；$d_y$ 为单个像素点在像素坐标y方向的长度，单位mm；u0为像素坐标系和图像坐标在x方向的平移偏量，v0为像素坐标系和图像坐标在y方向的平移偏量，单位mm。2.2 数据融合模型
在相机和激光雷达对地面作物进行识别的过程中，2种传感器会探测到地面作物的不同特征，通过信号处理将2个传感器收集到的目标相应特征输出，其输出数据对目标识别结果的支持程度以及输出数据的一致性函数[22]是系统准确识别作物的重要依据。
假设
$ \alpha$ 是需要识别的作物特征组成的集合，则将$ \alpha$ 中一个模糊集定义为一个隶属函数：$$ {m}_{j}(w)=1-\dfrac{|E-{\mu }_{j}|}{2{\sigma }_{j}}\text{；}{m}_{j}(w): \alpha \to [0,1],w\in \alpha ，$$ (5) 式中，
$ {m_j}(w) $ 表示第$ j $ 个传感器对模糊命题的隶属度，E表示地面目标特征，$ {\mu _j} $ 、$ {\sigma _j} $ 分别表示地面目标特征的均值和标准差，w表示识别的作物目标特征，不同的传感器获取作物的不同特征。因为采用多传感器进行作物识别，所以采用关联性结果来表示2个传感器对同一个支持模糊命题的支持程度。具体操作是先假设2个传感器
$ i $ 和$ j $ 对模糊命题的共同支持程度，即$ {m_i}(w) $ 和$ {m_j}(w) $ 之间的关联性，称之为一致性函数：$$ {S_{ji}} = \dfrac{{{m_j}(w) \wedge {m_i}(w)}}{{{m_j}(w) \vee {m_i}(w)}} = \dfrac{{\min \left\{ {{m_j}(w),{m_i}(w)} \right\}}}{{\max \left\{ {{m_j}(w),{m_i}(w)} \right\}}}， $$ (6) 式中，
$ \wedge $ 表示取两者中数值较小者，$ \vee $ 表示取两者中数值较大者。然后进行判别，当2个传感器对模糊命题的支持度相同时，取
$ {S_{ji}} = 1 $ ；当2个传感器对模糊命题的支持度将近时，表明2个传感器观测值相互支持度高，则$ 0 < {S_{ji}} < 1 $ ；当2个传感器对模糊命题的支持度相差较大，表明2个传感器观测值相互支持度低，甚至相互背离，因此对其信息不采用，则取$ {S_{ji}} = 0 $ 。如果出现2个传感器支持度持续不同的结论，就需要进一步判别这2个传感器是否出现故障。例如采用一组随机变量
$ {E_1},{E_2}, \cdots ,{E_h} $ 来表示通过多传感器获得地面目标的h个特征，通过随机变量的标准差与均值计算隶属度$ {m_j}(w) $ ，从而计算传感器对模糊命题的支持度$ {S_{ji}} $ 。假设随机变量的标准差和均值如表1所示。
表 1 玉米植株特征随机变量及其标准差和均值Table 1. Random variables, standard deviations and means of maize plant characteristics传感器编号
Sensor number目标1 Target 1 目标2 Target 2 目标3 Target 3 随机变量
Random variable标准差
Standard deviation均值
Mean标准差
Standard deviation均值
Mean标准差
Standard deviation均值
Mean1 50.0 246.0 68.0 372.0 23.0 289.0 318.0 2 1.1 1.8 0.6 2.6 0.7 2.9 3.4 将数值带入公式(5)计算的隶属度
$ {m_j}(w) $ 结果如表2所示。表 2 各传感器对不同目标的隶属度Table 2. Membership of each sensor to different targets传感器编号
Sensor number目标1
Target 1目标2
Target 2目标3
Target 31 0.28 0.59 0.37 2 0.27 0.33 0.64 根据公式(6)可以求得
$ {S}_{12}(目标1)\text{=}0.964 $ ，$ {S}_{12} (目标2)\text{=}0.559 $ ，$ {S}_{12}(目标3)\text{=}0.578 $ 。根据相机和激光雷达目标识别输出的一致性函数，构建传感器采集的所有特征对模糊命题的支持度一致性矩阵：
$$ {\boldsymbol{M}} = \left[ {\begin{array}{*{20}{c}} 1&{{S_{12}}}&{...}&{{S_{1i}}}&{...}&{{S_{1l}}} \\ \vdots & \vdots &{}& \vdots &{}& \vdots \\ {{S_{j1}}}&{{S_{j2}}}&{...}&{{S_{ji}}}&{...}&{{S_{jl}}} \\ \vdots & \vdots &{}& \vdots &{}& \vdots \\ {{S_{l1}}}&{{S_{l2}}}&{...}&{{S_{li}}}&{...}&1 \end{array}} \right] ，$$ (7) 式中，
$ l $ 表示采集到的目标的特征数量。利用上式中传感器采集的所有特征对模糊命题的支持度一致性计算出平均一致性，通过归一化处理，计算出第
$ j $ 个传感器作物识别输出的相对一致性函数：$$ M'{'_j} = \dfrac{{M{'_j}}}{{\displaystyle\sum\limits_{i = 1}^l {\dfrac{1}{{l - 1}}\displaystyle\sum\limits_{i = 1,i \ne j}^l {{S_{ji}}} } }} = \dfrac{{\dfrac{1}{{l - 1}}\displaystyle\sum\limits_{i = 1,i \ne j}^l {{S_{ji}}} }}{{\displaystyle\sum\limits_{i = 1}^l {\dfrac{1}{{l - 1}}\displaystyle\sum\limits_{i = 1,i \ne j}^l {{S_{ji}}} } }}， $$ (8) 式中，
$ M{'_j} $ 为第$ j $ 个传感器作物识别输出的一致性函数。通过一致性函数，即可将不同物理量的传感器数据进行匹配，保留图像特征重心和激光点云共同的作物特征点，剔除伪特征点。由于测试环境的目标干扰，行人、测试区域动态变化性等不确定因素，利用2种传感器采集目标的多类型特征数据，并依据获得的目标离散信息点集之间的离散状态确定多传感器的当前权重，离散因子根据不同时间段、多测试区域传感器测量的目标信息进行融合，有效避免了传感器权重选取的经验性和局限性；同时，将传感器当前权重引入到多传感器目标识别的平均加权一致性函数数学模型，通过计算处理给出多传感器目标识别的相对加权一致性计算函数；结合多传感器目标识别的相对一致性、多传感器目标识别的当前权重，构建基于离散因子的多传感器目标识别的数据融合支持度计算模型。
3. 导航线提取
将预处理后的图像与激光雷达特征点进行融合，剔除了伪特征点，对其获得的特征点进行直线拟合即可得到导航基准线。为减少异常点对作物特征点进行直线拟合的影响，引入随机采样一致(Random sample consensus, RANSAC)算法[23]，通过反复选择数据集估计出模型，一直迭代到估计出较好的模型。具体步骤为：
1) 首先筛选可以估计出模型的最小数据，使用这个数据计算出数据模；
2) 将所有数据带入这个模型，计算出组成模型参数数据的数量即内点数量；
3) 比较当前模型和之前推出最好模型的内点数量，记录最大内点数的模型参数；
4) 重复上述步骤，直到迭代结束或者当前模型内点数量达到要求。
此时需要对迭代次数进行选择，假设内点在数据中的占比为t。则：
$$ {t} = \dfrac{{{n_{\rm{inliers}}}}}{{{n_{\rm{inliers}}} + {n_{\rm{outliers}}}}}, $$ (9) 式中，
${n_{{\rm{inliers}}}}$ 为内点数量，${n_{{\rm{outliers}}}}$ 为外点数量。在每次计算模型使用N个点的情况下，选取的点至少有1个不是内点的概率是
$ 1 - {t^N} $ ，在迭代k次的情况下，$ {(1 - {t^N})^k} $ 为模型至少1次没有采样到内点时计算出模型的概率，即计算出错误模型的概率。那么能采样到正确的N个点并计算出正确模型的概率为：$$ P = 1 - {(1 - {t^N})^k} ,$$ (10) 通过上式，可以求得：
$$ k = \dfrac{{\lg (1 - P)}}{{\lg (1 - {t^n})}} ，$$ (11) 内点占比t通常为先验值，P是我们希望RANSAC算法得到正确模型的概率，此时可以求出最优的迭代次数k。通过RANSAC算法分别得到左、右导航基准线，2条导航基准线的角平分线即玉米田导航线，如图7所示。设左导航基准线斜率为
$ \;{\beta _{\text{1}}} $ ，右导航基准线斜率为$ \;{\beta _{\text{2}}} $ ，则导航线斜率$\; \beta $ 满足：$$ \dfrac{{\beta - {\beta _1}}}{{1 + \beta {\beta _1}}} = \dfrac{{{\beta _{\text{2}}} - \beta }}{{1 + \beta {\beta _2}}} 。$$ (12) 4. 试验及可靠性分析
试验采用英特尔Pentium(R) CPU G3250 @ 3.20 GHZ、4 GB内存、Windows 7(64位)操作系统的计算机，在Python3.7(Anaconda)的集成开发环境下编程完成。试验中的农业机器人宽45 cm、长60 cm，装有森云智能SG2-AR0233C-5200-GMSL2相机、16线激光雷达、4 G模块等多种传感器，采集前方环境和玉米植株信息，如图8所示。
为了更精确地分析该方法在玉米田中的植保机器人行走区域识别精度，本试验通过布置标定好行距和株距的玉米植株，模拟真实场景开展精度分析试验。结果表明，本文算法可以准确提取导航线。由于实时性和正确性是导航线提取的正确指标，表3统计了本文算法不同传感器对随机900帧图像的平均处理时间及正确率，通过人工标定图像最合理的导航线，定义人工与本文算法提取的导航线之间的夹角为误差角，误差角超过5°时，视为导航线错误。
表 3 不同传感器性能评价指标Table 3. Performance evaluation index of different sensors传感器
Sensor平均处理时间/ms
Average processing time处理帧数
No. of processing frames正确帧数
No. of correct frames正确率/%
Correct rate平均误差角度/(°)
Mean error angle相机 Camera 49.76 300 272 90.67 3.457 三维激光雷达
3D LiDAR58.15 300 278 92.67 2.670 传感器融合
Sensor fusion95.62 300 286 95.33 1.595 由表3可知，分别用相机、三维激光雷达和传感器融合的识别方法，平均处理时间逐渐增加，正确率也随之增加，当将相机和激光雷达融合时，平均处理时间和误差角度均满足要求。将该算法与候补定位点二次判别算法、动态迁移算法进行对比，本算法单帧平均处理时间为95.62 ms，优于候补定位点二次判别算法的单帧平均处理时间200 ms[9]和动态迁移算法的单帧平均处理时间97.56 ms[24]，具有更高的实时性；本算法正确率为95.33%，优于候补定位点二次判别算法的正确率90%[9]和动态迁移算法的正确率95%[24]，为农业机械的视觉导航提供了可靠保障。
5. 结论
本文基于半径的滤波方法依序迭代剔除点云数据噪音，利用DBSCAN算法进行点云数据聚类，导出最大密度相连的样本集合，采用激光雷达和相机信息融合的方式进行可通行区域提取，在图像数据中加入点云数据的反射率信息，降低了阴影、遮挡对导航线提取的影响。
根据多时间段、多空间区域传感器采集的目标特征量进行系统建模分析，并依据获得的目标特征值之间的离散因子确定多传感器权重选择，根据权重引入多传感器目标识别的平均加权一致性函数，获得具有冗余性的数据融合支持度模型。最终试验平均处理时间仅为95.62 ms，正确率高达95.33%，提高了检测结果的稳定性和准确性，为农业机器人在农田环境中作业提供了可靠的导航路径。
-
![]()
图 1 不同磷水平添加生物炭对大豆生物量和磷吸收的影响
P0：不施磷，P30：施用30 kg·hm−2磷，P60：施用60 kg·hm−2磷；P90：施用90 kg·hm−2磷；不同大写字母表示相同生物炭处理不同磷水平间差异显著(P<0.05，Duncan’s法)，不同小写字母表示相同磷水平不同量生物炭处理差异显著(P<0.05，t检验)。
Figure 1. Effect of biochar addition at different P levels on soybean biomass and P uptake
P0: No phosphorus, P30: 30 kg·hm−2 phosphorus applied, P60: 60 kg·hm−2 phosphorus applied, P90: 90 kg·hm−2 phosphorus applied; Different uppercase letters indicate significant differences among different phosphorus levels under the same biochar treatment (P<0.05, Duncan’s method), different lowercase letters indicate significant differences among different biochar treatments under the same phosphorus level (P<0.05, t test).
![]()
图 2 不同磷水平添加生物炭对土壤磷组分含量的影响
P0：不施磷；P30：施用30 kg·hm−2磷；P60：施用60 kg·hm−2磷；P90：施用90 kg·hm−2磷；BP0：施4%生物炭；BP30：施用4%的生物炭和30 kg·hm−2磷；BP60：施用4%生物炭和60 kg·hm−2磷；BP90：施用4%的生物炭和90 kg·hm−2磷。NH4F-Pi：活性无机磷；NaOHI-Pi、HCl-Pi：中等活性无机磷；NaOHII-Pi：稳定态无机磷；TOP：总有机磷；Residual-P：残余态磷。不同大写字母表示相同生物炭处理不同磷水平间差异显著(P<0.05，Duncan’s法)，不同小写字母表示相同磷水平不同生物炭处理间差异显著(P<0.05，t检验)。
Figure 2. Effect of biochar addition at different P levels on soil P fraction
P0: No phosphorus; P30: 30 kg·hm−2 phosphorus applied; P60: 60 kg·hm−2 phosphorus applied; P90: 90 kg·hm−2 phosphorus applied; BP0: 4% biochar; BP30: 4% biochar and 30 kg·hm−2 phosphorus are applied; BP60: 4% biochar and 60 kg·hm−2 phosphorus are applied; BP90: 4% biochar and 90 kg·hm−2 phosphorus are applied. NH4F-Pi: Labile inorganic phosphorus; NaOHI-Pi, HCl-Pi: Moderately labile inorganic phosphorus; NaOHII-Pi：Stable inorganic phosphorus; TOP: Total organic phosphorus; Residual-P: Residual phosphorus. Different uppercase letters indicate significant differences among different phosphorus levels under the same biochar treatment (P<0.05, Duncan’s method), different lowercase letters indicate significant differences among different biochar treatments under the same phosphorus level (P<0.05, t test).
![]()
图 3 不同磷水平添加生物炭对土壤磷酸酶活性和微生物量磷含量的影响
P0：不施磷；P30：施用30 kg·hm−2磷；P60：施用60 kg·hm−2磷；P90：施用90 kg·hm−2磷。不同大写字母表示相同生物炭处理不同磷水平间差异显著(P<0.05，Duncan’s法)，不同小写字母表示相同磷水平不同生物炭处理间差异显著(P<0.05，t检验)。
Figure 3. Effect of biochar addition at different P levels on soil phosphatase activity and microbial biomass P content
P0: No biochar and phosphorus; P30: 30 kg·hm−2 phosphorus applied; P60: 60 kg·hm−2 phosphorus applied; P90: 90 kg·hm−2 phosphorus applied. Different uppercase letters indicate significant differences among different phosphorus levels under the same biochar treatment (P<0.05, Duncan’s method), different lowercase letters indicate significant differences among different biochar treatments under the same phosphorus level (P<0.05, t test).
![]()
图 4 植株磷吸收量与不同磷组分含量、相关酶活性的相关性
P0：不施生物炭和磷；P30：施用30 kg·hm−2磷；P60：施用60 kg·hm−2磷；P90：施用90 kg·hm−2磷；BP0：施4%生物炭；BP30：施用4%的生物炭和30 kg·hm−2磷；BP60：施用4%生物炭和60 kg·hm−2磷；BP90：施用4%的生物炭和90 kg·hm−2磷。
Figure 4. Correlation between plant P uptake and different P components contents/related enzyme activities
P0: No biochar and phosphorus; P30: 30 kg·hm−2 phosphorus applied; P60: 60 kg·hm−2 phosphorus applied; P90: 90 kg·hm−2 phosphorus applied; BP0: 4% biochar; BP30: 4% biochar and 30 kg·hm−2 phosphorus are applied; BP60: 4% biochar and 60 kg·hm−2 phosphorus are applied; BP90: 4% biochar and 90 kg·hm−2 phosphorus are applied.
表 1 不同磷水平添加生物炭对土壤化学性质的影响1)
Table 1 Effects of biochar addition at different P levels on soil chemical properties
生物炭(w/%)
Biochar磷水平
P levelpH w(交换性阳离子)/(mg·kg−1)
Exchangeable cation content盐基饱和度/%
Salt base
saturation交换性酸
Exchangeable acidK+ Na+ Ca2+ Mg2+ b(H+)/
(cmol·kg−1)w(Al3+)/
(mg·kg−1)0 P0 5.06±0.04Ab 31.2±4.7Ab 89.2±4.2Ab 918.6±9.4Ab 47.2±4.4Ab 76.7±0.6Ab 0.43±0.02Ca 127.6±11.9Ba P30 5.04±0.01ABb 32.9±2.3Ab 92.0±4.1Ab 922.0±32.0Ab 44.4±2.6Ab 72.3±1.7BCb 0.53±0.03ABa 165.4±17.0ABa P60 4.94±0.05Bb 34.9±6.6Ab 95.0±2.6Ab 942.7±27.4Ab 44.6±4.3Ab 74.7±1.0ABb 0.46±0.03BCa 155.2±2.3ABb P90 4.78±0.03Cb 28.2±2.4Ab 98.7±2.6Ab 912.9±14.5Ab 37.6±2.6Ab 69.6±1.4Cb 0.58±0.02Aa 193.0±28.1Aa 4 P0 5.28±0.03Aa 230.5±11.6Aa 121.0±5.1Aa 1163.0 ±19.1ABa151.9±3.5Aa 94.8±0.7Aa 0.16±0.02BCb 26.3±7.4Bb P30 5.31±0.04Aa 231.0±13.6Aa 121.9±2.4Aa 1216.5 ±33.4Aa150.2±3.1ABa 94.8±0.5Aa 0.14±0.01Cb 32.4±7.1Bb P60 5.27±0.01Aa 219.8±1.8Aa 127.2±2.6Aa 989.5±26.6Ca 131.6±2.2Ca 89.4±0.9Ba 0.23±0.01Ab 70.9±8.3Ab P90 5.31±0.02Aa 221.9±15.0Aa 129.4±5.8Aa 1100.6 ±17.4Ba141.2±3.4Ba 93.5±0.2Aa 0.21±0.03ABb 27.0±7.6Bb 生物炭 Biochar *** *** *** *** *** *** *** *** 磷 Phosphorus ** ** ** ** ** 生物炭×磷
Biochar×Phosphorus*** *** * *** ** * 1)P0：不施磷，P30：施用30 kg·hm−2磷，P60：施用60 kg·hm−2磷，P90：施用90 kg·hm−2磷；同列数据后的不同大写字母表示相同生物炭处理不同磷水平间差异显著(P<0.05，Duncan’s法)，不同小写字母表示相同磷水平不同量生物炭处理间差异显著(P<0.05，t检验)；*、**和***分别表示在0.05、0.01和0.001水平差异显著(双因素方差分析)。
1) P0: No phosphorus, P30: 30 kg·hm−2 phosphorus applied, P60: 60 kg·hm−2 phosphorus applied, P90: 90 kg·hm−2 phosphorus applied; Different uppercase letters in the same column indicate significant differences among different phosphorus levels under the same biochar treatment (P<0.05, Duncan’s method), different lowercase letters indicate significant differences among different biochar treatments under the same phosphorus level (P<0.05, t test); *, ** and *** indicate significant differences at 0.05, 0.01 and 0.001 respectively (Two-way ANOVA).
下载: 导出CSV
表 2 添加生物炭对土壤磷组分的影响1)
Table 2 Effect of biochar addition on soil P fraction
处理
TreatmentNH4F-Pi NaOHI-Pi HCl-Pi NaOHII-Pi TOP Residual-P 生物炭 Biochar *** *** *** ** 磷 Phosphorus *** 生物炭×磷 Biochar × Phosphorus * 1)NH4F-Pi：活性无机磷；NaOHI-Pi、HCl-Pi：中等活性无机磷；NaOHII-Pi：稳定态无机磷；TOP：总有机磷；Residual-P：残余态磷。*、**和***分别表示在0.05、0.01和0.001水平差异显著(双因素方差分析)。
1) NH4F-Pi: Labile inorganic phosphorus; NaOHI-Pi, HCl-Pi: Moderately labile inorganic phosphorus; NaOHII-Pi：Stable inorganic phosphorus; TOP: Total organic phosphorus; Residual-P: Residual phosphorus. *, ** and *** indicate significant differences at 0.05, 0.01, and 0.001 respectively (Two-way ANOVA).
下载: 导出CSV
-
[1] 刘佳, 陈晓芬, 刘明, 等. 长期施肥对旱地红壤细菌群落的影响[J]. 土壤学报, 2020, 57(2): 468-478. doi: 10.11766/trxb201810090507 [2] 王蕾, 王艳玲, 李欢, 等. 长期施肥下红壤旱地磷素有效性影响因子的冗余分析[J]. 中国土壤与肥料, 2021(1): 17-25. doi: 10.11838/sfsc.1673-6257.19549 [3] 张宝贵, 李贵桐. 土壤生物在土壤磷有效化中的作用[J]. 土壤学报, 1998(1): 104-111. doi: 10.11766/trxb199508270115 [4] 李力, 陆宇超, 刘娅, 等. 玉米秸秆生物炭对Cd(Ⅱ)的吸附机理研究[J]. 农业环境科学学报, 2012, 31(11): 2277-2283. [5] ZHANG L, WU Z, ZHOU J, et al. Meta-analysis of the response of the productivity of different crops to parameters and processes in soil nitrogen cycle under biochar addition[J]. Agronomy, 2022, 12(8): 1857. doi: 10.3390/agronomy12081857
[6] SIGUA G C, NOVAK J M, WATTS D W, et al. Efficacies of designer biochars in improving biomass and nutrient uptake of winter wheat grown in a hard setting subsoil layer[J]. Chemosphere, 2016, 142: 176-183. doi: 10.1016/j.chemosphere.2015.06.015
[7] 刘玉学, 唐旭, 杨生茂, 等. 生物炭对土壤磷素转化的影响及其机理研究进展[J]. 植物营养与肥料学报, 2016, 22(6): 1690-1695. doi: 10.11674/zwyf.14585 [8] 才吉卓玛. 生物炭对不同类型土壤中磷有效性的影响研究[D]. 北京: 中国农业科学院, 2013. [9] GUL S, WHALEN J K. Biochemical cycling of nitrogen and phosphorus in biochar-amended soils[J]. Soil Biology and Biochemistry, 2016, 103: 1-15. doi: 10.1016/j.soilbio.2016.08.001
[10] XU G, SUN J, SHAO H, et al. Biochar had effects on phosphorus sorption and desorption in three soils with differing acidity[J]. Ecological Engineering, 2014, 62: 54-60. doi: 10.1016/j.ecoleng.2013.10.027
[11] 张晓娜, 陈平, 庞婷, 等. 玉米/豆科间作种植模式对作物干物质积累、分配及产量的影响[J]. 四川农业大学学报, 2017, 35(4): 484-490. [12] 任嘉宇. 大豆玉米带状复合种植技术[J]. 现代化农业, 2024(3): 14-17. doi: 10.3969/j.issn.1001-0254.2024.03.007 [13] SATTARI S Z, VAN ITTERSUN M K, GILLER K E, et al. Key role of China and its agriculture in global sustainable phosphorus management[J]. Environmental Research Letters, 2014, 9(5): 054003. doi: 10.1088/1748-9326/9/5/054003
[14] 邓金环. 生物炭改良酸性土壤及提高大豆硅磷吸收转化的机理研究[D]. 广州: 华南农业大学, 2019. [15] 鲍士旦. 土壤农化分析[M]. 3版. 北京: 中国农业出版社, 2000. [16] 鲁如坤. 土壤农业化学分析方法[M]. 北京: 中国农业科技出版社, 2000. [17] KUO S. Phosphorus[M]//SPARKS D L. Methods of Soil Analysis. Madison, WI: Soil Science Society of America. 1996: 869-919.
[18] WEAVER R W, ANGLE J S, BOTTOMLEY P S. Methods of Soil Analysis[M]//Microbiological and Biochemical Properties. Madison, WI: Soil Science Society of America. 1994.
[19] HEDLEY M J, STEWART J W B, CHAUHAN B S. Changes in inorganic and organic soil phosphorus fractions induced by cultivation practices and by laboratory incubations[J]. Soil Science Society of America Journal, 1982, 46(5): 970-976. doi: 10.2136/sssaj1982.03615995004600050017x
[20] TESFAYE F, LIU X, ZHENG J, et al. Could biochar amendment be a tool to improve soil availability and plant uptake of phosphorus? A meta-analysis of published experiments[J]. Environmental Science and Pollution Research, 2021, 28(26): 34108-34120. doi: 10.1007/s11356-021-14119-7
[21] 张登晓, 高雅, 介红彬, 等. 生物质炭对农田土壤磷有效性的影响研究进展[J]. 河南农业大学学报, 2021, 55(2): 199-205. [22] XU G, ZHANG Y, SUN J, et al. Negative interactive effects between biochar and phosphorus fertilization on phosphorus availability and plant yield in saline sodic soil[J]. Science of the Total Environment, 2016, 568: 910-915. doi: 10.1016/j.scitotenv.2016.06.079
[23] 武毅昶. 土壤调理剂对酸性土壤改良效果初探[D]. 呼和浩特: 内蒙古农业大学, 2017. [24] LATATI M, BLAVET D, ALKAMA N, et al. The intercropping cowpea-maize improves soil phosphorus availability and maize yields in an alkaline soil[J]. Plant and Soil, 2014, 385(1): 181-191.
[25] ZHU J, LI M, WHELAN M. Phosphorus activators contribute to legacy phosphorus availability in agricultural soils: A review[J]. Science of the Total Environment, 2018, 612: 522-537. doi: 10.1016/j.scitotenv.2017.08.095
[26] YANG C, LU S. Straw and straw biochar differently affect phosphorus availability, enzyme activity and microbial functional genes in an Ultisol[J]. Science of the Total Environment, 2022, 805: 150325. doi: 10.1016/j.scitotenv.2021.150325
[27] ZHANG H, CHEN C, GRAY E M, et al. Roles of biochar in improving phosphorus availability in soils: A phosphate adsorbent and a source of available phosphorus[J]. Geoderma, 2016, 276: 1-6. doi: 10.1016/j.geoderma.2016.04.020
[28] JIANG J, YUAN M, XU R, et al. Mobilization of phosphate in variable-charge soils amended with biochars derived from crop straws[J]. Soil and Tillage Research, 2015, 146: 139-147. doi: 10.1016/j.still.2014.10.009
[29] DELUCA T H, GLANVILLE H C, HARRIS M, et al. A novel biologically-based approach to evaluating soil phosphorus availability across complex landscapes[J]. Soil Biology and Biochemistry, 2015, 88: 110-119. doi: 10.1016/j.soilbio.2015.05.016
[30] RICHARDSON A E, SIMPSON R J. Soil microorganisms mediating phosphorus availability update on microbial phosphorus[J]. Plant Physiology, 2011, 156(3): 989-996. doi: 10.1104/pp.111.175448
[31] NANNIPIERI P, GIAGNONI L, LANDI L, et al. Role of phosphatase enzymes in soil[M]//BÜNEMANN E, OBERSON A, FROSSARD E. Soil Biology. Berlin, Heidelberg: Springer Berlin Heidelberg, 2010: 215-243.
[32] LEHMANN J, RILLIG M C, THIES J, et al. Biochar effects on soil biota: A review[J]. Soil Biology and Biochemistry, 2011, 43(9): 1812-1836. doi: 10.1016/j.soilbio.2011.04.022
[33] 和文祥, 朱铭莪, 张一平. 土壤酶与重金属关系的研究现状[J]. 土壤与环境, 2000, 9(2): 139-142. [34] 刘海旭, 吴俊江, 王金生, 等. 大豆耐低磷研究进展[J]. 大豆科学, 2017, 36(4): 639-644. [35] WANG W, CHEN Y, ZHANG F, et al. Cotton-maize intercropping increases rhizosphere soil phosphorus bioavailability by regulating key phosphorus cycling genes in Northwest China[J]. Applied Soil Ecology, 2023, 182: 104734. doi: 10.1016/j.apsoil.2022.104734
[36] LUO G, SUN B, LI L, et al. Understanding how long-term organic amendments increase soil phosphatase activities: Insight into phoD-and phoC-harboring functional microbial populations[J]. Soil Biology and Biochemistry, 2019, 139: 107632. doi: 10.1016/j.soilbio.2019.107632
[37] 易观文. 生物炭与磷肥配施对盐碱土磷素转化的影响及微生物学机制[D]. 上海: 华东师范大学, 2023. [38] LIANG J L, LIU J, JIA P, et al. Novel phosphate-solubilizing bacteria enhance soil phosphorus cycling following ecological restoration of land degraded by mining[J]. The ISME Journal, 2020, 14(6): 1600-1613. doi: 10.1038/s41396-020-0632-4
[39] HAN B, LI J, LIU K, et al. Variations in soil properties rather than functional gene abundances dominate soil phosphorus dynamics under short-term nitrogen input[J]. Plant and Soil, 2021, 469(1): 227-241.
[40] SILES J A, STARKE R, MARTINOVIC T, et al. Distribution of phosphorus cycling genes across land uses and microbial taxonomic groups based on metagenome and genome mining[J]. Soil Biology and Biochemistry, 2022, 174: 108826. doi: 10.1016/j.soilbio.2022.108826
[41] ZHOU L L, CAO J, ZHANG F S, et al. Rhizosphere acidification of faba bean, soybean and maize[J]. Science of the Total Environment, 2009, 407(14): 4356-4362. doi: 10.1016/j.scitotenv.2009.02.006
-
期刊类型引用(7)
1. 李志刚，彭兆伟，宁琳如，李亚美，金海勇，田鹏. 多传感器数据融合下气体密度继电器自动温控技术研究. 粘接. 2025(05): 139-142 . 
百度学术
2. 王俊鹏，黎良坤，冯杨. 卷烟条包端部贴标技术设计. 自动化仪表. 2025(04): 17-20+25 . 百度学术
3. 刘洋，冀杰，赵立军，冯伟，贺庆，王小康. 基于激光雷达与相机融合的树干检测方法. 西南大学学报(自然科学版). 2024(02): 183-196 . 百度学术
4. 江雪峰. 传感器信息融合下新能源汽车动力电池信号故障检测方法. 东莞理工学院学报. 2024(03): 94-99 . 百度学术
5. 王张夫，汤显峰. 基于协方差交叉融合的多传感器数据融合研究. 电子测量技术. 2024(08): 78-85 . 百度学术
6. 朱哲人，张伟. 基于传感器信息融合的运动员上肢运动姿态识别方法. 河北北方学院学报(自然科学版). 2024(07): 4-10 . 百度学术
7. 焦萍萍，周显春，高华玲，杨真. 复杂环境下多模态交通警示信息推荐算法仿真. 计算机仿真. 2023(07): 121-125 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 95
- HTML全文浏览量: 22
- PDF下载量: 44
- 被引次数: 11
